wyiu
Follow.
C++博客
::
首頁
::
新隨筆
::
聯系
::
聚合
::
管理
posts - 100, comments - 15, trackbacks - 0
<
2009年5月
>
日
一
二
三
四
五
六
26
27
28
29
30
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
1
2
3
4
5
6
常用鏈接
我的隨筆
我的評論
我參與的隨筆
留言簿
(1)
給我留言
查看公開留言
查看私人留言
隨筆分類
(84)
Design Pattern(1)
POJ(67)
常用模板和函數(3)
數據結構(2)
數值分析(3)
算法(8)
隨筆檔案
(100)
2010年10月 (8)
2010年3月 (11)
2009年10月 (16)
2009年9月 (1)
2009年8月 (1)
2009年7月 (20)
2009年5月 (16)
2009年4月 (27)
向高手學習
cai0715
RyanWang
wingyiu
搜索
積分與排名
積分 - 28269
排名 - 675
最新評論
1.?re: poj 3368 rmq(st)
有錯誤~~
index[] 可能會以負數為下標~
--tongjiantao
2.?re: pku poj 3009
@BOAT
行列搞錯了吧?
--yaoyaozii
3.?re: pku poj 3009
@ww
郁悶啊。。我怎么也是3 啊。。。郁悶!!!
--BOAT
4.?re: pku poj 3009 [未登錄]
excit就是這個游戲
--X
5.?re: pku2777
這個只能用c++交嗎?用G++交的結果很神奇喔,居然CE。。
--share4
閱讀排行榜
1.?POJ 線段樹題(2000)
2.?pku poj 3009 (1015)
3.?矩陣轉置_十字鏈表(960)
4.?全主元高斯消元法(805)
5.?關于“逆序數”[轉](776)
評論排行榜
1.?pku poj 3009 (10)
2.?pku 1019 poj(2)
3.?pku2777(1)
4.?pku 2151(1)
5.?poj 3368 rmq(st)(1)
列主元高斯消元算法
#include
<
iostream
>
#include
<
math.h
>
using
namespace
std;
#define
MAX 100
double
A[MAX
+
1
][MAX
+
1
];
double
B[MAX
+
1
];
double
X[MAX
+
1
];
double
e;
int
n;
void
ColGaussianElimination()
{
int
i,j,k,kmi;
double
T;
for
(k
=
1
;k
<=
n
-
1
;k
++
)
{
//
選主元
T
=
0
;
for
(i
=
k;i
<=
n;i
++
)
if
( fabs(A[i][k])
>
T )
{ T
=
A[i][k];kmi
=
i;}
if
( T
<=
e)
{ printf(
"
Error!\n
"
);
return
;}
if
(kmi
!=
k)
{
T
=
B[k];B[k]
=
B[kmi];B[kmi]
=
T;
//
swap B[k] and B[kmi]
for
(j
=
k;j
<=
n;j
++
)
//
swap row kmi and k of A
{
T
=
A[k][j];
A[k][j]
=
A[kmi][j];
A[kmi][j]
=
T;
}
}
//
消元
for
(i
=
k
+
1
;i
<=
n;i
++
)
{
T
=
A[i][k]
/
A[k][k];
B[i]
-=
T
*
B[k];
for
(j
=
k;j
<=
n;j
++
)
A[i][j]
-=
T
*
A[k][j];
}
}
//
回代
if
( fabs(A[n][n])
<=
e )
{ printf(
"
Error!\n
"
);
return
;}
X[n]
=
B[n]
/
A[n][n];
double
S_Aij_Xj;
for
(i
=
n
-
1
;i
>=
1
;i
--
)
{
S_Aij_Xj
=
0
;
for
(j
=
i
+
1
;j
<=
n;j
++
)
S_Aij_Xj
+=
A[i][j]
*
X[j];
X[i]
=
(B[i]
-
S_Aij_Xj)
/
A[i][i];
}
}
void
print(
double
X[])
{
int
i;
printf(
"
X[]:\n
"
);
for
(i
=
1
;i
<=
n;i
++
)
printf(
"
%f\n
"
,X[i]);
}
int
main()
{
int
i,j;
printf(
"
n:
"
);
scanf(
"
%d
"
,
&
n);
printf(
"
A[][]:\n
"
);
for
(i
=
1
;i
<=
n;i
++
)
for
(j
=
1
;j
<=
n;j
++
)
scanf(
"
%lf
"
,
&
A[i][j]);
printf(
"
B[]:\n
"
);
for
(i
=
1
;i
<=
n;i
++
)
scanf(
"
%lf
"
,
&
B[i]);
printf(
"
e:
"
);
scanf(
"
%lf
"
,
&
e);
ColGaussianElimination();
print(X);
system(
"
pause
"
);
return
0
;
}
posted on 2009-05-16 16:34
wyiu
閱讀(616)
評論(0)
編輯
收藏
引用
所屬分類:
數值分析
只有注冊用戶
登錄
后才能發表評論。
【推薦】100%開源!大型工業跨平臺軟件C++源碼提供,建模,組態!
相關文章:
全主元高斯消元法
列主元高斯消元算法
順序高斯消元算法
網站導航:
博客園
IT新聞
BlogJava
博問
Chat2DB
管理
Powered by:
C++博客
Copyright ©2025 wyiu
国产精品免费久久
|
天天爽天天狠久久久综合麻豆
|
精品蜜臀久久久久99网站
|
97久久精品国产精品青草
|
国产精品亚洲美女久久久
|
奇米影视7777久久精品
|
国产精品成人久久久久久久
|
久久人做人爽一区二区三区
|
www久久久天天com
|
亚洲国产成人精品无码久久久久久综合
|
亚洲欧洲中文日韩久久AV乱码
|
久久99国内精品自在现线
|
7777精品伊人久久久大香线蕉
|
国产69精品久久久久777
|
久久婷婷五月综合成人D啪
|
久久精品这里热有精品
|
久久午夜羞羞影院免费观看
|
久久久这里有精品中文字幕
|
伊人色综合久久天天
|
久久久免费精品re6
|
久久人做人爽一区二区三区
|
久久AⅤ人妻少妇嫩草影院
|
国产精品久久久久9999高清
|
少妇久久久久久被弄高潮
|
久久久久久久91精品免费观看
|
精品无码久久久久久久久久
|
99久久综合狠狠综合久久止
|
热re99久久精品国99热
|
无码人妻久久一区二区三区免费丨
|
国产麻豆精品久久一二三
|
一级a性色生活片久久无少妇一级婬片免费放
|
久久99久久99精品免视看动漫
|
亚洲国产精品人久久
|
日韩人妻无码精品久久免费一
|
日韩人妻无码一区二区三区久久99
|
狠狠色丁香婷婷综合久久来来去
|
97久久久精品综合88久久
|
国产精品欧美久久久天天影视
|
97久久超碰国产精品2021
|
久久se精品一区二区
|
色综合久久天天综合
|