bzoj 1013

又做了一道水题

实在是太水啦！

直接列方程高斯消元即可

#include<bits/stdc++.h>using namespace std;const double eps = 1e-8 ;int n;double k[20][20],a[20][20],x[20][20],r,sum[20];int main(){    scanf("%d",&n);    for(int i=1;i<=n+1;i++)     for(int j=1;j<=n;j++)       scanf("%lf",&k[i][j]),sum[i]+=k[i][j]*k[i][j];    for(int i=1;i<=n;i++)    {     for(int j=1;j<=n;j++)       a[i][j]=(k[i][j]-k[i+1][j])*2;     a[i][n+1]=sum[i]-sum[i+1];    }    for(int i=1;i<=n;i++)    {       int temp=i;       for(temp=i;temp<=n;temp++)         if(fabs(a[temp][i])>0) break;       if(temp!=i)        for(int j=1;j<=n+1;j++) swap(a[i][j],a[temp][j]);       double key=a[i][i];       for(int j=1;j<=n+1;j++) a[i][j]/=key;       for(int j=1;j<=n;j++) if(j!=i)       {          key=a[j][i];          for(int k=1;k<=n+1;k++)           a[j][k]-=key*a[i][k];       }    }       for(int i=1;i<n;i++) printf("%.3lf ",a[i][n+1]);    printf("%.3lf",a[n][n+1]);}

心得：
1.高斯消元是很重要的，是必须要会的，这个模板是非常好的