wyiu
Follow.
C++博客
::
首頁(yè)
::
新隨筆
::
聯(lián)系
::
聚合
::
管理
posts - 100, comments - 15, trackbacks - 0
<
2025年7月
>
日
一
二
三
四
五
六
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
7
8
9
常用鏈接
我的隨筆
我的評(píng)論
我參與的隨筆
留言簿
(1)
給我留言
查看公開留言
查看私人留言
隨筆分類
(84)
Design Pattern(1)
POJ(67)
常用模板和函數(shù)(3)
數(shù)據(jù)結(jié)構(gòu)(2)
數(shù)值分析(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)
向高手學(xué)習(xí)
cai0715
RyanWang
wingyiu
搜索
積分與排名
積分 - 28270
排名 - 675
最新評(píng)論
1.?re: poj 3368 rmq(st)
有錯(cuò)誤~~
index[] 可能會(huì)以負(fù)數(shù)為下標(biāo)~
--tongjiantao
2.?re: pku poj 3009
@BOAT
行列搞錯(cuò)了吧?
--yaoyaozii
3.?re: pku poj 3009
@ww
郁悶啊。。我怎么也是3 啊。。。郁悶!!!
--BOAT
4.?re: pku poj 3009 [未登錄](méi)
excit就是這個(gè)游戲
--X
5.?re: pku2777
這個(gè)只能用c++交嗎?用G++交的結(jié)果很神奇喔,居然CE。。
--share4
閱讀排行榜
1.?POJ 線段樹題(2000)
2.?pku poj 3009 (1015)
3.?矩陣轉(zhuǎn)置_十字鏈表(960)
4.?全主元高斯消元法(805)
5.?關(guān)于“逆序數(shù)”[轉(zhuǎn)](776)
評(píng)論排行榜
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)
評(píng)論(0)
編輯
收藏
引用
所屬分類:
數(shù)值分析
只有注冊(cè)用戶
登錄
后才能發(fā)表評(píng)論。
【推薦】100%開源!大型工業(yè)跨平臺(tái)軟件C++源碼提供,建模,組態(tài)!
相關(guān)文章:
全主元高斯消元法
列主元高斯消元算法
順序高斯消元算法
網(wǎng)站導(dǎo)航:
博客園
IT新聞
BlogJava
博問(wèn)
Chat2DB
管理
Powered by:
C++博客
Copyright ©2025 wyiu
怡红院日本一道日本久久
|
精品国产一区二区三区久久久狼
|
97精品久久天干天天天按摩
|
久久精品国产亚洲欧美
|
久久93精品国产91久久综合
|
东方aⅴ免费观看久久av
|
色婷婷综合久久久久中文
|
99久久国产综合精品成人影院
|
精品久久久久久无码国产
|
伊人久久大香线焦AV综合影院
|
丰满少妇高潮惨叫久久久
|
亚洲国产成人乱码精品女人久久久不卡
|
久久久久久免费视频
|
国内精品久久久久
|
久久婷婷五月综合97色直播
|
久久精品无码一区二区三区
|
久久精品无码一区二区WWW
|
国产午夜精品理论片久久
|
人妻无码αv中文字幕久久琪琪布
|
国产亚洲精午夜久久久久久
|
久久久免费精品re6
|
久久久久亚洲AV无码观看
|
欧美久久久久久精选9999
|
精品久久久久久综合日本
|
亚洲AV成人无码久久精品老人
|
久久人人爽人人澡人人高潮AV
|
国内精品伊人久久久久777
|
久久91精品国产91久久户
|
亚洲日韩中文无码久久
|
国产精品成人久久久
|
亚洲精品无码久久久久AV麻豆
|
精品久久久久一区二区三区
|
亚洲国产综合久久天堂
|
7国产欧美日韩综合天堂中文久久久久
|
久久99热这里只频精品6
|
久久国产乱子伦精品免费午夜
|
久久九九亚洲精品
|
国产国产成人久久精品
|
91精品婷婷国产综合久久
|
国产精品日韩深夜福利久久
|
国产精品内射久久久久欢欢
|