高斯消元 poj-1222、1681、1830、3185,sicily-1048

最新推荐文章于 2020-12-31 17:20:39 发布
最新推荐文章于 2020-12-31 17:20:39 发布
阅读量471 收藏
点赞数
分类专栏： 高斯消元 文章标签： 数学
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u013113958/article/details/41450123
版权

这几个题目都是翻转之后关联其他点的翻转，并且都是采用的01翻转，实属良心。

另外建立增广矩阵真的很费劲

  poj 1222 
 
 
  1. #include<iostream> 
  2. #include<cstring> 
  3. #include <algorithm> 
  4. #include<cstdlib> 
  5. #include<vector> 
  6. #include<cmath> 
  7. #include<stdlib.h> 
  8. #include<iomanip> 
  9. #include<list> 
  10. #include<deque> 
  11. #include<map> 
  12. #include <stdio.h> 
  13. #include <queue> 
  14.  
  15. const int maxn=1000+5; 
  16.  
  17. #define inf 0x3f3f3f3f 
  18.   #define INF 0x3FFFFFFFFFFFFFFFLL 
  19. #define rep(i,n) for(i=0;i<n;i++) 
  20.  #define reP(i,n) for(i=1;i<=n;i++) 
  21.  
  22. #define ull unsigned long long 
  23.  #define ll long long 
  24.  
  25. #define cle(a) memset(a,0,sizeof(a)) 
  26.  
  27. using namespace std; 
  28. //对2取模的01方程组 
  29.  
  30. const int MAXN = 40; 
  31. //有equ个方程，var个变元。增广矩阵行数为equ,列数为var+1,分别为0到var 
  32. int equ,var;//方程数和变量数,注意对其赋值 
  33. int a[MAXN][MAXN]; //增广矩阵 
  34. //增广矩阵:在系数矩阵的右边添上一列，这一列是线性方程组的等号右边的值。 
  35. int x[MAXN]; //解集 
  36. int free_x[MAXN];//用来存储自由变元（多解枚举自由变元可以使用） 
  37. int free_num;//自由变元的个数 
  38. void init(){//初始化 
  39.     cle(a),cle(x),cle(free_x); 
  40. } 
  41. //返回值为-1表示无解，为0是唯一解，否则返回自由变元个数 
  42. void out(){ 
  43.     int i,j; 
  44.     rep(i,20){ 
  45.         rep(j,21){ 
  46.             cout<<a[i][j]<<" "; 
  47.         } 
  48.         cout<<endl; 
  49.     } 
  50.     cout<<"_____bug______"<<endl; 
  51. } 
  52. int ans[MAXN]; 
  53. int now; 
  54.  
  55. int dfs(int v) 
  56. { 
  57.     int i,j; 
  58.     int cop_a[MAXN][MAXN]; 
  59.     if(v==0){ 
  60.         int tt=0; 
  61.         for(i=0;i<20;i++){ 
  62.             x[i]=ans[i]; 
  63.         } 
  64.         memcpy(cop_a,a,sizeof(a)); 
  65.         for(i=now-1;i>=0;i--){ 
  66.             for(j=i;j<var;j++){ 
  67.                 cop_a[i][20]^=(x[j]&cop_a[i][j]); 
  68.             } 
  69.             x[i]^=cop_a[i][20]; 
  70.         } 
  71.         rep(i,20)tt+=x[i]; 
  72.         return tt; 
  73.     } 
  74.     else{ 
  75.         ans[free_x[v-1]]=0; 
  76.         int tt=dfs(v-1); 
  77.         ans[free_x[v-1]]=1; 
  78.         tt=min(tt,dfs(v-1)); 
  79.         return tt; 
  80.     } 
  81. } 
  82. int Gauss() { 
  83.     int max_r,col,k; 
  84.     free_num = 0; 
  85.     for(k = 0, col = 0 ; k <equ&& col <var ; k++, col++) { 
  86.         max_r = k; 
  87.         for(int i = k+1; i<equ; i++) { 
  88.             if(abs(a[i][col]) > abs(a[max_r][col])) 
  89.                 max_r = i; 
  90.         } 
  91.         if(a[max_r][col] == 0) { 
  92.             k--; 
  93.             free_x[free_num++] = col;//这个是自由变元 
  94.             continue; 
  95.         } 
  96.         if(max_r != k) { 
  97.             for(int j = col; j < var+1; j++) 
  98.                 swap(a[k][j],a[max_r][j]); 
  99.         } 
  100.         for(int i = k+1; i<equ; i++) { 
  101.             if(a[i][col] != 0) { 
  102.                 for(int j = col; j < var+1; j++) 
  103.                     a[i][j] ^= a[k][j]; 
  104.             } 
  105.         } 
  106.     } 
  107.     now=k; 
  108. //    cout<<k<<":"<<var<<endl; 
  109.     for(int i = k; i<equ; i++) 
  110.         if(a[i][col] != 0) 
  111.             return -1;//无解 
  112.     if(k<var)  
  113.         //自由变元个数 
  114.     { 
  115.         int tt; 
  116.         tt=dfs(var-k); 
  117.         now=k; 
  118.         cout<<tt<<endl; 
  119.         return -1; 
  120.     } 
  121.     //唯一解，回代 
  122.    else{ 
  123.        for(int i = var-1; i>= 0; i--) { 
  124.         x[i] = a[i][var]; 
  125.         for(int j = i+1; j <var; j++) 
  126.             x[i] ^= (a[i][j] && x[j]); 
  127.        } 
  128.    } 
  129.     return 0; 
  130. } 
  131. int main() 
  132. { 
  133. #ifndef ONLINE_JUDGE 
  134.      freopen("in.txt","r",stdin); 
  135.      //freopen("out.txt","w",stdout); 
  136. #endif 
  137.     int n;init(); 
  138.     while(cin>>a[0][20]){ 
  139.         int i; 
  140.         for(i=1;i<=19;i++){ 
  141.             scanf("%d",&a[i][20]); 
  142.         } 
  143.         equ=var=20; 
  144.         rep(i,20){ 
  145.             a[i][i]=1; 
  146.             if(i-1>=0)a[i][i-1]=1; 
  147.             if(i+1<20)a[i][i+1]=1; 
  148.         } 
  149.         int ans=Gauss(); 
  150.         if(ans==0){ 
  151.             rep(i,20){ 
  152.             ans+=x[i]; 
  153.          } 
  154.             cout<<ans<<endl; 
  155.  
  156.         }init(); 
  157.     } 
  158.     return 0; 
  159. }
  poj 1681
 
 
  1.    #include<iostream> 
  2. #include<cstring> 
  3. #include <algorithm> 
  4. #include<cstdlib> 
  5. #include<vector> 
  6. #include<cmath> 
  7. #include<stdlib.h> 
  8. #include<iomanip> 
  9. #include<list> 
  10. #include<deque> 
  11. #include<map> 
  12. #include <stdio.h> 
  13. #include <queue> 
  14.  
  15. const int maxn=1000+5; 
  16.  
  17. #define inf 0x3f3f3f3f 
  18.   #define INF 0x3FFFFFFFFFFFFFFFLL 
  19. #define rep(i,n) for(i=0;i<n;i++) 
  20.  #define reP(i,n) for(i=1;i<=n;i++) 
  21.  
  22. #define ull unsigned long long 
  23.  #define ll long long 
  24.  
  25. #define cle(a) memset(a,0,sizeof(a)) 
  26.  
  27. using namespace std; 
  28. const int MAXN = 250; 
  29. //有equ个方程，var个变元。增广矩阵行数为equ,列数为var+1,分别为0到var 
  30. int equ,var;//注意对其赋值 
  31. int a[MAXN][MAXN]; //增广矩阵 
  32. //增广矩阵:在系数矩阵的右边添上一列，这一列是线性方程组的等号右边的值。 
  33. int x[MAXN]; //解集 
  34. int free_x[MAXN];//用来存储自由变元（多解枚举自由变元可以使用） 
  35. int free_num;//自由变元的个数 
  36. void init(){//初始化 
  37. cle(a),cle(x),cle(free_x); 
  38. } 
  39. int n; 
  40. int dir[4][2]={1,0,-1,0,0,1,0,-1}; 
  41. int getbool(int x,int y){ 
  42.     if(x<0||y<0||x>=n||y>=n)return 0; 
  43.     return 1; 
  44. } 
  45. void getcon(int x,int y){ 
  46.     int newx,newy; 
  47.     int i; 
  48.     rep(i,4){ 
  49.         newx=x+dir[i][0],newy=y+dir[i][1]; 
  50.         if(getbool(newx,newy)){ 
  51.             a[x*n+y][newx*n+newy]=1; 
  52.         } 
  53.     } 
  54. } 
  55. void out(){ 
  56.     int i,j; 
  57.     rep(i,n*n){ 
  58.         rep(j,n*n+1){ 
  59.             cout<<a[i][j]<<" "; 
  60.         } 
  61.         cout<<endl; 
  62.     } 
  63. } 
  64. //返回值为-1表示无解，为0是唯一解，否则返回自由变元个数 
  65. int Gauss() { 
  66.     int max_r,col,k; 
  67.     free_num = 0; 
  68.     for(k = 0, col = 0 ; k <equ&& col <var ; k++, col++) { 
  69.         max_r = k; 
  70.         for(int i = k+1; i<equ; i++) { 
  71.             if(abs(a[i][col]) > abs(a[max_r][col])) 
  72.                 max_r = i; 
  73.         } 
  74.         if(a[max_r][col] == 0) { 
  75.             k--; 
  76.             free_x[free_num++] = col;//这个是自由变元 
  77.             continue; 
  78.         } 
  79.         if(max_r != k) { 
  80.             for(int j = col; j < var+1; j++) 
  81.                 swap(a[k][j],a[max_r][j]); 
  82.         } 
  83.         for(int i = k+1; i<equ; i++) { 
  84.             if(a[i][col] != 0) { 
  85.                 for(int j = col; j < var+1; j++) 
  86.                     a[i][j] ^= a[k][j]; 
  87.             } 
  88.         } 
  89.     } 
  90.     for(int i = k; i<equ; i++) 
  91.         if(a[i][col] != 0) 
  92.             return -1;//无解 
  93.    // if(k <var) return var-k;//自由变元个数 
  94.     //唯一解，回代 
  95.     for(int i = var-1; i>= 0; i--) { 
  96.         x[i] = a[i][var]; 
  97.         for(int j = i+1; j <var; j++) 
  98.             x[i] ^= (a[i][j] && x[j]); 
  99.     } 
  100.     return 0; 
  101. } 
  102. int main() 
  103. { 
  104. #ifndef ONLINE_JUDGE 
  105.      freopen("in.txt","r",stdin); 
  106.      //freopen("out.txt","w",stdout); 
  107. #endif 
  108.     int T; 
  109.     cin>>T; 
  110.     while(T--){ 
  111.         cin>>n; 
  112.         init(); 
  113.         equ=var=n*n; 
  114.         int i,j; 
  115.         for(i=0;i<n*n;i++){ 
  116.             char c; 
  117.             cin>>c; 
  118.             if(c=='w'){ 
  119.                 a[i][n*n]=1; 
  120.             } 
  121.         } 
  122.         for(i=0;i<n;i++){ 
  123.             for(j=0;j<n;j++){ 
  124.                 a[i*n+j][i*n+j]=1; 
  125.                 getcon(i,j); 
  126.             } 
  127.         } 
  128.         //out(); 
  129.         int ans=Gauss(); 
  130.         if(ans==0){ 
  131.             int sum=0; 
  132.             rep(i,n*n){ 
  133.                 sum+=x[i]; 
  134.             } 
  135.             cout<<sum<<endl; 
  136.         } 
  137.         else if(ans==-1){ 
  138.             cout<<"inf"<<endl; 
  139.         } 
  140.     } 
  141.     return 0; 
  142. }
   poj 1830
 
 
  1. #include<iostream> 
  2. #include<cstring> 
  3. #include <algorithm> 
  4. #include<cstdlib> 
  5. #include<vector> 
  6. #include<cmath> 
  7. #include<stdlib.h> 
  8. #include<iomanip> 
  9. #include<list> 
  10. #include<deque> 
  11. #include<map> 
  12. #include <stdio.h> 
  13. #include <queue> 
  14.  
  15. const int maxn=1000+5; 
  16.  
  17. #define inf 0x3f3f3f3f 
  18.   #define INF 0x3FFFFFFFFFFFFFFFLL 
  19. #define rep(i,n) for(i=0;i<n;i++) 
  20.  #define reP(i,n) for(i=1;i<=n;i++) 
  21.  
  22. #define ull unsigned long long 
  23.  #define ll long long 
  24.  
  25. #define cle(a) memset(a,0,sizeof(a)) 
  26.  
  27. using namespace std; 
  28. const int MAXN = 910; 
  29. //有equ个方程，var个变元。增广矩阵行数为equ,列数为var+1,分别为0到var 
  30. int equ,var;//注意对其赋值 
  31. int a[MAXN][MAXN]; //增广矩阵 
  32. //增广矩阵:在系数矩阵的右边添上一列，这一列是线性方程组的等号右边的值。 
  33. int x[MAXN]; //解集 
  34. int free_x[MAXN];//用来存储自由变元（多解枚举自由变元可以使用） 
  35. int free_num;//自由变元的个数 
  36. void init(){//初始化 
  37.     cle(a),cle(x),cle(free_x); 
  38. } 
  39. int n; 
  40. void out(){ 
  41.     int i,j; 
  42.     rep(i,n){ 
  43.         rep(j,n+1){ 
  44.             cout<<a[i][j]<<" "; 
  45.         } 
  46.         cout<<endl; 
  47.     } 
  48. } 
  49. //返回值为-1表示无解，为0是唯一解，否则返回自由变元个数 
  50. int Gauss() { 
  51.     int max_r,col,k; 
  52.     free_num = 0; 
  53.     for(k = 0, col = 0 ; k <equ&& col <var ; k++, col++) { 
  54.         max_r = k; 
  55.         for(int i = k+1; i<equ; i++) { 
  56.             if(abs(a[i][col]) > abs(a[max_r][col])) 
  57.                 max_r = i; 
  58.         } 
  59.         if(a[max_r][col] == 0) { 
  60.             k--; 
  61.             free_x[free_num++] = col;//这个是自由变元 
  62.             continue; 
  63.         } 
  64.         if(max_r != k) { 
  65.             for(int j = col; j < var+1; j++) 
  66.                 swap(a[k][j],a[max_r][j]); 
  67.         } 
  68.         for(int i = k+1; i<equ; i++) { 
  69.             if(a[i][col] != 0) { 
  70.                 for(int j = col; j < var+1; j++) 
  71.                     a[i][j] ^= a[k][j]; 
  72.             } 
  73.         } 
  74.     } 
  75.     for(int i = k; i<equ; i++) 
  76.         if(a[i][col] != 0) 
  77.             return -1;//无解 
  78.     if(k <var) 
  79.         return var-k;//自由变元个数 
  80.     //唯一解，回代 
  81.     for(int i = var-1; i>= 0; i--) { 
  82.         x[i] = a[i][var]; 
  83.         for(int j = i+1; j <var; j++) 
  84.             x[i] ^= (a[i][j] && x[j]); 
  85.     } 
  86.     return 0; 
  87. } 
  88. int main() 
  89. { 
  90. #ifndef ONLINE_JUDGE 
  91.      freopen("in.txt","r",stdin); 
  92.      //freopen("out.txt","w",stdout); 
  93. #endif 
  94.     int T; 
  95.     cin>>T; 
  96.     while(T--){ 
  97.         cin>>n; 
  98.         equ=var=n; 
  99.         int i,j; 
  100.         int be[maxn]; 
  101.         int en[maxn]; 
  102.         rep(i,n){ 
  103.             cin>>be[i]; 
  104.         } 
  105.         init(); 
  106.         rep(i,n){ 
  107.             cin>>en[i]; 
  108.             if(be[i]!=en[i]){ 
  109.                 a[i][n]=1; 
  110.             } 
  111.         } 
  112.         int x,y; 
  113.         while(cin>>x>>y){ 
  114.             if(x==0&&y==0){ 
  115.                 break; 
  116.             } 
  117.             a[y-1][x-1]=1; 
  118.         } 
  119.         rep(i,n)a[i][i]=1; 
  120.         int ans=Gauss(); 
  121.        // cout<<ans<<endl; 
  122.         if(ans==-1){ 
  123.             cout<<"Oh,it's impossible~!!"<<endl; 
  124.         } 
  125.         else{ 
  126.             cout<<(1<<ans)<<endl; 
  127.         } 
  128.     } 
  129.     return 0; 
  130. }
       poj 3185


 
 
  1. #include<iostream> 
  2. #include<cstring> 
  3. #include <algorithm> 
  4. #include<cstdlib> 
  5. #include<vector> 
  6. #include<cmath> 
  7. #include<stdlib.h> 
  8. #include<iomanip> 
  9. #include<list> 
  10. #include<deque> 
  11. #include<map> 
  12. #include <stdio.h> 
  13. #include <queue> 
  14.  
  15. const int maxn=1000+5; 
  16.  
  17. #define inf 0x3f3f3f3f 
  18.   #define INF 0x3FFFFFFFFFFFFFFFLL 
  19. #define rep(i,n) for(i=0;i<n;i++) 
  20.  #define reP(i,n) for(i=1;i<=n;i++) 
  21.  
  22. #define ull unsigned long long 
  23.  #define ll long long 
  24.  
  25. #define cle(a) memset(a,0,sizeof(a)) 
  26.  
  27. using namespace std; 
  28. //对2取模的01方程组 
  29.  
  30. const int MAXN = 40; 
  31. //有equ个方程，var个变元。增广矩阵行数为equ,列数为var+1,分别为0到var 
  32. int equ,var;//方程数和变量数,注意对其赋值 
  33. int a[MAXN][MAXN]; //增广矩阵 
  34. //增广矩阵:在系数矩阵的右边添上一列，这一列是线性方程组的等号右边的值。 
  35. int x[MAXN]; //解集 
  36. int free_x[MAXN];//用来存储自由变元（多解枚举自由变元可以使用） 
  37. int free_num;//自由变元的个数 
  38. void init(){//初始化 
  39.     cle(a),cle(x),cle(free_x); 
  40. } 
  41. //返回值为-1表示无解，为0是唯一解，否则返回自由变元个数 
  42. void out(){ 
  43.     int i,j; 
  44.     rep(i,20){ 
  45.         rep(j,21){ 
  46.             cout<<a[i][j]<<" "; 
  47.         } 
  48.         cout<<endl; 
  49.     } 
  50.     cout<<"_____bug______"<<endl; 
  51. } 
  52. int ans[MAXN]; 
  53. int now; 
  54.  
  55. int dfs(int v) 
  56. { 
  57.     int i,j; 
  58.     int cop_a[MAXN][MAXN]; 
  59.     if(v==0){ 
  60.         int tt=0; 
  61.         for(i=0;i<20;i++){ 
  62.             x[i]=ans[i]; 
  63.         } 
  64.         memcpy(cop_a,a,sizeof(a)); 
  65.         for(i=now-1;i>=0;i--){ 
  66.             for(j=i;j<var;j++){ 
  67.                 cop_a[i][20]^=(x[j]&cop_a[i][j]); 
  68.             } 
  69.             x[i]^=cop_a[i][20]; 
  70.         } 
  71.         rep(i,20)tt+=x[i]; 
  72.         return tt; 
  73.     } 
  74.     else{ 
  75.         ans[free_x[v-1]]=0; 
  76.         int tt=dfs(v-1); 
  77.         ans[free_x[v-1]]=1; 
  78.         tt=min(tt,dfs(v-1)); 
  79.         return tt; 
  80.     } 
  81. } 
  82. int Gauss() { 
  83.     int max_r,col,k; 
  84.     free_num = 0; 
  85.     for(k = 0, col = 0 ; k <equ&& col <var ; k++, col++) { 
  86.         max_r = k; 
  87.         for(int i = k+1; i<equ; i++) { 
  88.             if(abs(a[i][col]) > abs(a[max_r][col])) 
  89.                 max_r = i; 
  90.         } 
  91.         if(a[max_r][col] == 0) { 
  92.             k--; 
  93.             free_x[free_num++] = col;//这个是自由变元 
  94.             continue; 
  95.         } 
  96.         if(max_r != k) { 
  97.             for(int j = col; j < var+1; j++) 
  98.                 swap(a[k][j],a[max_r][j]); 
  99.         } 
  100.         for(int i = k+1; i<equ; i++) { 
  101.             if(a[i][col] != 0) { 
  102.                 for(int j = col; j < var+1; j++) 
  103.                     a[i][j] ^= a[k][j]; 
  104.             } 
  105.         } 
  106.     } 
  107.     now=k; 
  108. //    cout<<k<<":"<<var<<endl; 
  109.     for(int i = k; i<equ; i++) 
  110.         if(a[i][col] != 0) 
  111.             return -1;//无解 
  112.     if(k<var)  
  113.         //自由变元个数 
  114.     { 
  115.         int tt; 
  116.         tt=dfs(var-k); 
  117.         now=k; 
  118.         cout<<tt<<endl; 
  119.         return -1; 
  120.     } 
  121.     //唯一解，回代 
  122.    else{ 
  123.        for(int i = var-1; i>= 0; i--) { 
  124.         x[i] = a[i][var]; 
  125.         for(int j = i+1; j <var; j++) 
  126.             x[i] ^= (a[i][j] && x[j]); 
  127.        } 
  128.    } 
  129.     return 0; 
  130. } 
  131. int main() 
  132. { 
  133. #ifndef ONLINE_JUDGE 
  134.      freopen("in.txt","r",stdin); 
  135.      //freopen("out.txt","w",stdout); 
  136. #endif 
  137.     int n;init(); 
  138.     while(cin>>a[0][20]){ 
  139.         int i; 
  140.         for(i=1;i<=19;i++){ 
  141.             scanf("%d",&a[i][20]); 
  142.         } 
  143.         equ=var=20; 
  144.         rep(i,20){ 
  145.             a[i][i]=1; 
  146.             if(i-1>=0)a[i][i-1]=1; 
  147.             if(i+1<20)a[i][i+1]=1; 
  148.         } 
  149.         int ans=Gauss(); 
  150.         if(ans==0){ 
  151.             rep(i,20){ 
  152.             ans+=x[i]; 
  153.          } 
  154.             cout<<ans<<endl; 
  155.  
  156.         }init(); 
  157.     } 
  158.     return 0; 
  159. }
     sicily 1048 是找的题目给大二小朋友做搜索的，结果意外发现可以用高斯定理，然后就给顺带A了
 
 
  1. // Problem#: 1048 
  2. // Submission#: 3258134 
  3. // The source code is licensed under Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License 
  4. // URI: http://creativecommons.org/licenses/by-nc-sa/3.0/ 
  5. // All Copyright reserved by Informatic Lab of Sun Yat-sen University 
  6. /************************************************************************* 
  7.     > File Name: 3.cpp 
  8.     > Author: cy 
  9.     > Mail: 1002@qq.com  
  10.     > Created Time: 14/11/23 14:23:15 
  11.  ************************************************************************/ 
  12.  
  13. #include<iostream> 
  14. #include<cstring> 
  15. #include <algorithm> 
  16. #include<cstdlib> 
  17. #include<vector> 
  18. #include<cmath> 
  19. #include<stdlib.h> 
  20. #include<iomanip> 
  21. #include<list> 
  22. #include<deque> 
  23. #include<map> 
  24. #include <queue> 
  25.  
  26. const int maxn=1000+5; 
  27.  
  28. #define inf 0x3f3f3f3f 
  29.   #define INF 0x3FFFFFFFFFFFFFFFLL 
  30. #define rep(i,n) for(i=0;i<n;i++) 
  31.  #define reP(i,n) for(i=1;i<=n;i++) 
  32.  
  33. #define ull unsigned long long 
  34.  #define ll long long 
  35.  
  36. #define cle(a) memset(a,0,sizeof(a)) 
  37.  
  38. using namespace std; 
  39. const int MAXN = 40; 
  40. //有equ个方程，var个变元。增广矩阵行数为equ,列数为var+1,分别为0到var 
  41. int equ,var;//方程数和变量数,注意对其赋值 
  42. int a[MAXN][MAXN]; //增广矩阵 
  43. //增广矩阵:在系数矩阵的右边添上一列，这一列是线性方程组的等号右边的值。 
  44. int x[MAXN]; //解集 
  45. int free_x[MAXN];//用来存储自由变元（多解枚举自由变元可以使用） 
  46. int free_num;//自由变元的个数 
  47. void init(){//初始化 
  48.     cle(a),cle(x),cle(free_x); 
  49.         equ=var=9; 
  50. } 
  51. //返回值为-1表示无解，为0是唯一解，否则返回自由变元个数 
  52. void out(){ 
  53.     int i,j; 
  54.     rep(i,9){ 
  55.         rep(j,10){ 
  56.             cout<<a[i][j]<<" "; 
  57.         } 
  58.         cout<<endl; 
  59.     } 
  60.     cout<<"_____bug______"<<endl; 
  61. } 
  62. int ans[MAXN]; 
  63. int now; 
  64.  
  65. int Gauss() { 
  66.     int max_r,col,k; 
  67.     free_num = 0; 
  68.     for(k = 0, col = 0 ; k <equ&& col <var ; k++, col++) { 
  69.         max_r = k; 
  70.         for(int i = k+1; i<equ; i++) { 
  71.             if(abs(a[i][col]) > abs(a[max_r][col])) 
  72.                 max_r = i; 
  73.         } 
  74.         if(a[max_r][col] == 0) { 
  75.             k--; 
  76.             free_x[free_num++] = col;//这个是自由变元 
  77.             continue; 
  78.         } 
  79.         if(max_r != k) { 
  80.             for(int j = col; j < var+1; j++) 
  81.                 swap(a[k][j],a[max_r][j]); 
  82.         } 
  83.         for(int i = k+1; i<equ; i++) { 
  84.             if(a[i][col] != 0) { 
  85.                 for(int j = col; j < var+1; j++) 
  86.                     a[i][j] ^= a[k][j]; 
  87.             } 
  88.         } 
  89.     } 
  90.     now=k; 
  91. //  cout<<k<<":"<<var<<endl; 
  92.     for(int i = k; i<equ; i++) 
  93.         if(a[i][col] != 0) 
  94.             return -1;//无解 
  95.     if(k<var)  
  96.         //自由变元个数 
  97.     { 
  98.         int tt; 
  99.         tt=dfs(var-k); 
  100.         now=k; 
  101.         cout<<tt<<endl; 
  102.         return -1; 
  103.     } 
  104.     //唯一解，回代 
  105.    else{ 
  106.        for(int i = var-1; i>= 0; i--) { 
  107.         x[i] = a[i][var]; 
  108.         for(int j = i+1; j <var; j++) 
  109.             x[i] ^= (a[i][j] && x[j]); 
  110.        } 
  111.    } 
  112.     return 0; 
  113. } 
  114. int dir[8][2]={-1,-1,-1,0,-1,1,0,-1,0,1,1,-1,1,0,1,1}; 
  115. bool getbool(int x,int y){ 
  116.     if(x<0||y<0||x>=3||y>=3)return false; 
  117.     return true; 
  118. } 
  119. int main() 
  120. { 
  121.  
  122.     int n; 
  123.     int T; 
  124.     cin>>T; 
  125.  
  126.     while(T--){ 
  127.         init(); 
  128.         int i; 
  129.         int num=0; 
  130.         for(i=0;i<9;i++){ 
  131.             char c; 
  132.             cin>>c; 
  133.             if(c=='b')num++,a[i][9]=1; 
  134.         } 
  135.         if(!num){ 
  136.             cout<<11<<endl; 
  137.             continue; 
  138.         } 
  139.         int j,k; 
  140.         for(i=0;i<3;i++){ 
  141.             rep(j,3){ 
  142.                 a[i*3+j][i*3+j]=1; 
  143.                 rep(k,8){ 
  144.                     int newx=i+dir[k][0]; 
  145.                     int newy=j+dir[k][1]; 
  146.                     if(getbool(newx,newy)){ 
  147.                         a[i*3+j][newx*3+newy]=1; 
  148.                     } 
  149.                 } 
  150.             } 
  151.         } 
  152.     //  out(); 
  153.         int ans=Gauss(); 
  154.             rep(i,9){ 
  155.                 if(x[i]){ 
  156.                     cout<<i+1; 
  157.                 } 
  158.             } 
  159.             cout<<endl; 
  160.     } 
  161.     return 0; 
  162. }
caoyu01
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
开关问题 POJ - 1830 高斯消元
九大于七的博客
04-11 3084
开关问题 POJ - 1830 有N个相同的开关，每个开关都与某些开关有着联系，每当你打开或者关闭某个开关的时候，其他的与此开关相关联的开关也会相应地发生变化，即这些相联系的开关的状态如果原来为开就变为关，如果为关就变为开。你的目标是经过若干次开关操作后使得最后N个开关达到一个特定的状态。对于任意一个开关，最多只能进行一次开关操作。你的任务是，计算有多少种可以...
POJ---3185：The Water Bowls【高斯消元】
KobeDuu
03-06 462
题意： 一排有20个碗，给出每个碗的状态，改变一个碗的状态，其左右两边的碗的状态也随之改变(0 -> 1，1 -> 0)，求全部碗的状态为0时的最少操作数 分析: （1）每个碗的状态只有0或1（所以mod 2），设第i个碗旋转了 xi 次，得到以下方程： 说明：如果转动第j个碗影响第i个碗，aij = 1，否则aij = 0；vi 为第 i 碗的初末状态的异...
KMP算法解释
蔡军帅
02-12 280
给定两个字符串A,B，判断T是否为S的子串（变式：寻找子串B在串A中的位置）。  要求一个O(|A|+|B|)的做法。  通常称A为目标串（或主串），B为模式串。  算法过程：  我们假设串A的长度为n，串B的长度为m，每个字符串的开头下标默认为1。     定义两个变量i和j，这两个变量共同表示：A[i-j+1~i]与B[1~j]均匹配，即：A中以第i个字符结尾的、长度为j的字符...
hdu 3416 Marriage Match IV
绝望的乐园
08-10 814
#include #include #include #include #include #include #include #include #include #include #define inf 2139062143 #define Inf 0x3FFFFFFFFFFFFFFFLL #define eps 1e-9 #define pi acos(-1.0) using namespac
SRM 622 div1
绝望的乐园
05-29 1016
好久没写
HDU 2865 Birthday Toy(Polya+矩阵乘法+dp)
绝望的乐园
04-16 857
题意：要构成一个长度为n
记录
acmer的随记
09-27 347
long long INF:const ll INF=0x3FFFFFFFFFFFFFFFLL;
高斯消元 POJ 1830 1222 1753
Baileys0530的专栏
07-05 772
做了高斯消元题目若干，总结： AX = B A是已知的状态矩阵，a[i][j]表示物j对物i有关联作用。如果作用是相互的，则是对称矩阵。 B的设置根据初始状态和末状态的差值状态。B0+AX = B1，B = B1 - B0. A是n*n矩阵，秩为r；增广矩阵（A|B）是n*(n+1)矩阵，秩为s Ⅰ.  s != r 方程无解，即没有方案满足； Ⅱ.  s == r
高斯消元入门②-异或线性方程组-POJ1222
qq_35577488的博客
12-31 511
题目大意: 给你一个5∗65 * 65∗6大小的二维01矩阵代表灯的开关状态。改变任意一个灯的状态，会导致其上下左右的灯的状态的改变。问一种灯的开关方案使得最终所有灯都关闭. 题目思路: 首先，一个灯只会被按0次或1次。两次等于没按。不同灯按下的效果是叠加起来的. 解释:想象按下某个位置,等效于原矩阵异或上一个新矩阵。 例如按下(2,2)(2,2)(2,2)产生的效果. A2,2=[010111010]A_{2,2} = \left[\begin{array}{lll} 0 & 1 & 0
高斯消元入门③-异或线性方程,无解,自由元,计数-POJ1830
qq_35577488的博客
12-31 951
题目大意: 给你n(n≤29)n(n \leq 29)n(n≤29)个开关。每个开关改变状态会跟着改变其他若干个开关的状态。给你一个开始状态，一个结束状态。问你有多少种操作方案.(每个开关至多改变一次状态，且无视操作顺序). 题目思路: 套路与POJ1222几乎一样。将每个开关按或不按设成未知数,列出异或线性方程求解即可。 这里关键对解的存在性进行讨论: 答案是2free2^{free}2free. freefreefree代表自由元的个数. 1.如何求自由元：（自由元的含义是：可以取定义域中的任何值的未知
算法-开关问题（POJ-1830）（包含源程序）.rar
09-16
标题中的“算法-开关问题（POJ-1830）”是指一个编程竞赛题目，源自编程平台POJ（Programming Online Judge）。这类题目通常要求参赛者编写程序来解决特定的数学或逻辑问题。POJ-1830是这个特定问题的编号，方便用户...
POJ 2947 高斯消元--判断一解多解无解
Jackyguo1992的专栏
04-23 3654
又是偷来的代码，但是第一道高斯消元，纪念下 知道每条式子的结果取模后的值，求方程组，这不同于一般的高斯消元，而且答案必须为整数并在一定区间内 1.由于是整数要用到gcd； 2.解的判断： 有解的情况都是必须满足row一下全为0 一解：化简后是严格的三角矩阵 多解：自由变量个数较少 无解：row以下有非0的 #include #include #include #inclu
codeforces　round 261
caoyu01
08-17 422
A :给出正方形的两个点求出另外两个点，没有的话输出-1 B
UVA 11995 I Can Guess the Data Structure!
绝望的乐园
06-06 753
题意：根据所给的操作，猜测数据结构。 思路：直接模拟即可。   #include #include #include #include #include #include #include #include #include #include #define inf 0x3f3f3f3f #define Inf 0x3FFFFFFFFFFFFFFFLL #define eps 1e-9 #
1446-基于51单片机的气压检测（MPX4115,上位机，LCD）proteus，原理图、流程图、物料清单、仿真图、源代码
09-06
1446-基于51单片机的气压检测（MPX4115,上位机，LCD）proteus，原理图、流程图、物料清单、仿真图、源代码 资料介绍： 1、使用气压传感器MPX4115检测气压 2、利用LCD1602显示气压 3、把气压值上传到上位机 有哪些资料： 1、仿真工程文件 2、源代码工程文件 3、原理图工程文件 4、流程图 5、功能介绍 6、元件清单 7、上位机软件 8、演示视频
8051Proteus仿真c源码用定时器T0的中断实现长时间定时
09-06
8051Proteus仿真c源码用定时器T0的中断实现长时间定时提取方式是百度网盘分享地址
威纶通HMI登陆方式（二）：索引密码登陆
09-06
通过用户参数设置、高级安全模式控制地址的使用以及用户登陆界面的设计实现威纶通HMI索引密码登陆。
基于yolov8的抽烟检测系统python源码+onnx模型+评估指标曲线+精美GUI界面.zip
09-06
【测试环境】 windows10 anaconda3+python3.8 torch==1.9.0+cu111 ultralytics==8.2.70 【模型可以检测出类别】 cig-pack smoke 更多实现细节参考博文：https://blog.csdn.net/FL1623863129/article/details/141964642
epoll+线程池实现的高并发服务器
最新发布
09-06
epoll+线程池实现的高并发服务器
POJ - 1661
09-18
根据提供的引用内容，这个问题讨论的是POJ-1661题，这是一个dp问题。题目给定了一系列平台的高度和人的起始位置，要求从起始位置跳到地面的最短时间。通过设置状态dp[i]表示从第i个平台的左边跳到地面所需的最小时间...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值