Null__Cat假假
by muyang_233 @ 2019-08-22 20:56:22
@[muyang_233](/space/show?uid=113521) 您先别在这fAKe，为什么连样例都过不了qwq
by Null_Cat @ 2019-08-22 20:57:59
tql
by End_donkey @ 2019-08-22 21:02:24
就写球形空间产生器的时候写过高斯消元
by End_donkey @ 2019-08-22 21:03:06
@[Null_Cat](/space/show?uid=213346) 那好吧 ~~我来现场学一下高斯消元~~
by muyang_233 @ 2019-08-22 21:04:06
@[muyang_233](/space/show?uid=113521) ~~现场切orz~~
by Null_Cat @ 2019-08-22 21:04:38
@[Null_Cat](/space/show?uid=213346)
我截了一段来，~~太久没写了，忘记咋敲来着了~~
b是系数
c是常数
```cpp
for(int i=1;i<=n;++i){
for(int j=i;j<=n;++j){
if(b[j][i]!=0){
for(int k=1;k<=n;++k){
swap(b[j][k],b[i][k]);
}
swap(c[i],c[j]);
}
}
}
for(int i=1;i<=n;++i){
for(int j=i;j<=n;++j){
if(i==j)continue;
double tmp=b[j][i]/b[i][i];
for(int k=1;k<=n;++k){
b[j][k]=b[j][k]-tmp*b[i][k];
}
c[j]=c[j]-tmp*c[i];
}
}
for(int i=n;i>=1;--i){
for(int j=i+1;j<=n;++j){
c[i]=c[i]-c[j]*b[i][j];
}
c[i]=c[i]/b[i][i];
}
```
by End_donkey @ 2019-08-22 21:07:19
@[End_donkey](/space/show?uid=141637) ~~为什么还要在开数组存常数~~求解释一下代码qwq小蒟蒻看不太懂qwq
by Null_Cat @ 2019-08-22 21:11:25
@[Null_Cat](/space/show?uid=213346) 您跟我讲下原理吧 我写下qaq
by muyang_233 @ 2019-08-22 21:11:46
@[muyang_233](/space/show?uid=113521) ~~出门右转百度“线性方程组”~~
by Null_Cat @ 2019-08-22 21:12:57