POJ 3185 The Water Bowls （高斯消元 求最小步数）

题目链接

题意：有20个数字，0或1。如果改变一个数的状态，它左右两边的两个数的状态也会变反。问从目标状态到全0，至少需要多少次操作。

分析：

和上一题差不多，但是比上一题还简单，不多说了，但是在做题的时候犯了一个非常二的错误。。看图吧。

先输入了a[0]又，初始了a[]【】数组

 #include <iostream>
 #include <cstdio>
 #include <cstring>
 #include <cstdlib>
 #include <cmath>
 #include <algorithm>
 #define LL __int64
 const int maxn = +;
 const int INF = <<;
 using namespace std;
 int equ, var, fn;
 int a[maxn][maxn], x[maxn];
 int free_x[maxn];
 int gcd(int a, int b)
 {
     return b==?a:gcd(b, a%b);
 }
 int lcm(int a, int b)
 {
     return a*b/gcd(a, b);
 }
 int Gauss()
 {
     int x_mo;
     x_mo = ;
     int i, j, k, max_r, col;
     int ta, tb, LCM, fx_num = ;
     col = ;
 
     for(k = ; k<equ && col<var; k++, col++)
     {
         max_r = k;
         for(i = k+; i < equ; i++)
             if(abs(a[i][col])>abs(a[max_r][col]))
                 max_r = i;
 
         if(max_r != k)
             for(j = k; j < var+; j++)
                 swap(a[k][j], a[max_r][j]);
 
         if(a[k][col]==)
         {
             free_x[fx_num++] = col; //求自由变元所在的列
             k--;
             continue;
         }
         for(i = k+; i < equ; i++)
         {
             if(a[i][col] != )
             {
                 LCM = lcm(abs(a[i][col]), abs(a[k][col]));
                 ta = LCM/abs(a[i][col]);
                 tb= LCM/abs(a[k][col]);
                 if(a[i][col]*a[k][col] < ) tb = -tb;
 
                 for(j = col; j < var+; j++)
                     a[i][j] = ((a[i][j]*ta - a[k][j]*tb)%x_mo+x_mo)%x_mo;
             }
         }
     }
     for(i = k; i < equ; i++)
         if(a[i][col] != )
             return INF;
 
     int stat=<<(var-k);
     int res=INF;
     for(i=; i<stat; i++)
     {
         int cnt=;
         int index=i;
         for(j=; j<var-k; j++)
         {
             x[free_x[j]]=(index&);
             if(x[free_x[j]]) cnt++;
             index>>=;
         }
         for(j=k-; j>=; j--)
         {
             int tmp=a[j][var];
             for(int l=j+; l<var; l++)
                 if(a[j][l]) tmp^=x[l];
             x[j]=tmp;
             if(x[j])cnt++;
         }
         if(cnt<res)res=cnt;
     }
     return res;
 }
 
 void init()
 {
     int i;
     memset(a, , sizeof(a));
     memset(x, , sizeof(x));
     for(i = ; i < ; i++)
     {
         a[i][i] = ;
         if(i+<)
         a[i+][i] = ;
         if(i->=)
         a[i-][i] = ;
     }
 }
 int main()
 {
     int i;
         equ = ;
         var = ;
         init();
         for(i = ; i < ; i++)
         scanf("%d", &a[i][]);
         fn = Gauss();
         printf("%d\n", fn);
     return ;
 }