var
n:longint;
a:array[0..101,0..101]of extended;
ans:array[0..101]of extended;
procedure init;
var
i,j:longint;
begin
readln(n);
for i:=1 to n do
begin
for j:=1 to n+1 do
read(a[i,j]);
readln;
end;
end;
procedure swap(i,j:longint);
var
c:array[0..101]of extended;
begin
c:=a[i];
a[i]:=a[j];
a[j]:=c;
end;
procedure work;
var
i,j,k,maxi:longint;
p:extended;
begin
for i:=1 to n-1 do
begin
maxi:=i;
for j:=i+1 to n do
if abs(a[j,i])>abs(a[i,maxi]) then maxi:=j;
swap(i,maxi);
for j:=i+1 to n do
begin
p:=-(a[j,i]/a[i,i]);
for k:=i to n+1 do
a[j,k]:=a[j,k]+a[i,k]*p;
end;
end;
for i:=n downto 1 do
begin
p:=a[i,n+1];
for j:=n downto i+1 do
p:=p-a[i,j]*ans[j];
ans[i]:=p/a[i,i];
end;
for i:=1 to n do
write(round(ans[i]),' ');
writeln;
end;
begin
init;
work;
end.
3 条评论
-
twd2 LV 9 MOD @ 2013-10-06 00:13:42
希望以后题解能发到各题目的题解区下，谢谢合作！
I wish that the problem solutions can be sent to each topic area under the problem solution, thank you!
Решение Надежда проблема может быть отправлен в каждой тематической области под решение проблемы, спасибо!
希望以後題解能發到各題目的題解區下，謝謝合作！
ホープ問題解決、問題の解決策の下で、各トピックエリアに送信することができ、ありがとうございました！
La solution du problème de l'espoir peut être envoyé à chaque sujet dans la solution du problème, merci!
Rozwiązanie problemu Nadzieja mogą być wysyłane do każdego obszaru tematu w ramach rozwiązania problemu, dziękuję! -
@ 2013-10-05 22:54:22
希望以后题解能发到各题目的题解区下，谢谢合作！
I wish that the problem solutions can be sent to each topic area under the problem solution, thank you!
Решение Надежда проблема может быть отправлен в каждой тематической области под решение проблемы, спасибо!
希望以後題解能發到各題目的題解區下，謝謝合作！
ホープ問題解決、問題の解決策の下で、各トピックエリアに送信することができ、ありがとうございました！
La solution du problème de l'espoir peut être envoyé à chaque sujet dans la solution du problème, merci!
Rozwiązanie problemu Nadzieja mogą być wysyłane do każdego obszaru tematu w ramach rozwiązania problemu, dziękuję! -
@ 2013-10-05 18:55:36
为什么
- 1
