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 ;

 }