uvalive 6185

高斯消元，以前从来没写过，今天的模拟比赛里面，添琦给了我一个模板！

虽然是个裸的，但是因为从来没写过，一个小细节竟然搞了我几个小时；

终于最后在小珺同志的帮助下成功a掉了，太开心了！

存一下，作为模板！

代码：

 #include<cstdio>
 #define maxn 10
 #include<cstring>
 #include<cmath>
 #include<algorithm>
 #define eps 0.00001
 using namespace std;
 double matrix[maxn][maxn+];
 double ans[maxn];
 double a[maxn];

 void exchange_col(int p1,int p2,int n)
 {
     double t;
     int i;
     for(int i=; i<=n; i++)
         t=matrix[p1][i],matrix[p1][i]=matrix[p2][i],matrix[p2][i]=t;
 }

 bool gauss(int n)
 {
     int i,j,k;
     int p;
     double r;
     for(i=; i<n-; i++)
     {
         p=i;
         for(j=i+; j<n; j++)
             if(abs(matrix[j][i])>abs(matrix[p][i]))
                 p=j;
         if(p!=i)
             exchange_col(i,p,n);
         if(matrix[i][i]==) return false;
         for(j=i+; j<n; j++)
         {
             r=matrix[j][i]/matrix[i][i];
             for(k=i; k<=n; k++)
                 matrix[j][k]-=r*matrix[i][k];
         }
     }
     for(i=n-; i>=; i--)
     {
         ans[i]=matrix[i][n];
         for(j=n-; j>i; j--)
             ans[i]-=matrix[i][j]*ans[j];
         if(abs(matrix[i][i])<eps)return false;
         ans[i]/=matrix[i][i];
     }
     return true;
 }

 int main()
 {
     int d;
     while(scanf("%d",&d)&&d)
     {
         memset(matrix,,sizeof matrix);
         for(int i=; i<=d+; i++)
             scanf("%lf",&a[i]);
         for(int k=; k<=d+; k++)
         {
             int cnt=;
             for(int j=; j<=d+; j++)
             {
                 if(k==j) continue;
                 matrix[cnt][d+]=a[j];
                 for(int i=; i<=d; i++)
                 {
                     if(i==) matrix[cnt][i]=;
                     else matrix[cnt][i]=matrix[cnt][i-]*j;
                 }
                 cnt++;
             }
             if(!gauss(d+))
             {
                 printf("%d\n",k);
                 break;
             }
         }
     }
     return ;
 }