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国产精品久久
|
久久精品桃花综合
|
久久精品国产一区二区三区日韩
|
人人狠狠综合久久亚洲88
|
久久人人爽人人爽人人片AV东京热
|
欧美精品一区二区久久
|
久久婷婷五月综合97色
|
久久99亚洲综合精品首页
|
日本WV一本一道久久香蕉
|
久久男人Av资源网站无码软件
|
国产2021久久精品
|
久久亚洲私人国产精品vA
|
色播久久人人爽人人爽人人片aV
|
日韩人妻无码精品久久免费一
|
精品久久人人爽天天玩人人妻
|
天堂久久天堂AV色综合
|
久久亚洲国产成人影院网站
|
精品久久久久久久久中文字幕
|
亚洲综合久久久
|
久久精品二区
|
国内精品久久久久久久亚洲
|
久久久无码人妻精品无码
|
一级a性色生活片久久无
|
国产成人无码精品久久久免费
|
久久国产高清字幕中文
|
亚洲精品无码专区久久同性男
|
久久99免费视频
|
99久久国语露脸精品国产
|
日韩久久久久久中文人妻
|
伊人久久综合成人网
|
久久久国产视频
|
久久久久亚洲av综合波多野结衣
|
久久夜色精品国产亚洲av
|
亚洲精品高清久久
|
国产福利电影一区二区三区,免费久久久久久久精
|
久久精品国产99久久丝袜
|
大蕉久久伊人中文字幕
|
AA级片免费看视频久久
|
国产精品99久久久久久www
|
94久久国产乱子伦精品免费
|
久久中文娱乐网
|