高斯消元裸题 记录模板 BZOJ 1013

两遍过

第一遍是PE了 最后一个数不能输出空格 气死啦

                                                        By——WXH

//我的变量是需要就加一个 所以就有很多定义 惨啊

#include<bits/stdc++.h>using namespace std;int n;double st[20];double s,tmp;double gs[20][20];double ans[20];void gauss(){for(int i=1;i<n;i++){double t=gs[i][i];for(int j=i;j<=n+1;j++)gs[i][j]/=t;for(int j=i+1;j<=n;j++){double tt=-gs[j][i];for(int k=i;k<=n+1;k++){gs[j][k]+=tt*gs[i][k];}}}for(int i=n;i;i--){double t=gs[i][n+1];for(int j=n;j>i;j--){t-=gs[i][j]*ans[j];}ans[i]=t/gs[i][i];}}int main(){scanf("%d",&n);for(int i=1;i<=n;i++){scanf("%lf",&st[i]);s+=st[i]*st[i];}for(int i=1;i<=n;i++){double sum=0;for(int j=1;j<=n;j++){scanf("%lf",&tmp);sum+=tmp*tmp;gs[i][j]=2*(tmp-st[j]);}gs[i][n+1]=sum-s;}gauss();for(int i=1;i<n;i++){printf("%.3lf ",ans[i]);}printf("%.3lf",ans[n]);return 0;}