poj1753Flip Game<高斯消元>

链接 http://poj.org/problem?id=1753

 

View Code 
  1 #include <stdio.h>
  2 #include <stdlib.h>
  3 #include <string.h> 
  4 #include <algorithm>
  5 using namespace std;
  6 const int N=16;
  7 int d[20][20];
  8 int slove( int n )   //枚举自由解 求最小值 
  9 {
 10     int x[20], temp, ans=100, m[4], t, p, k; 
 11     for(  t=1; t<=16; ++t ){
 12         memset( x, 0, sizeof x ); 
 13         temp=0; 
 14         for(  p=0; p<4; ++p ){
 15             m[p]=(t&( 1<<p ));
 16             if( m[p] ){
 17                 x[13+p]=1;    
 18             }    
 19         }
 20         for( p=N; p>0; --p ){
 21             for(  k=N; k>0; --k ){
 22                 if( p!=k )
 23                     x[p] ^= (x[k]*d[p][k]);    
 24             }
 25             x[p] ^= d[p][n];
 26             if (x[p]) {
 27                 ++temp;
 28             }    
 29         }    
 30         ans=min( ans, temp ); 
 31     } 
 32     return ans;    
 33 } 
 34 void Gauss( )
 35 {
 36     int i=1, j, k, p, t;
 37      
 38     for( j=1; j<=N; ++i,++j ){
 39         for( p=i; p<=N; p++ ){
 40                 if(d[p][j])break;
 41             }
 42         if( p==N+1 ){--i;continue;} 
 43         if( p!=i ){
 44             for( k=j; k<=N+2; ++k ){
 45                 t=d[i][k];
 46                 d[i][k]=d[p][k];
 47                 d[p][k]=t;    
 48             }    
 49         } 
 50         for( k=i+1; k<=N; ++k ){
 51             if( d[k][j] ){ 
 52                 for( t=j; t<=N+2; ++t ){
 53                     d[k][t]^=d[i][t]; 
 54                 }    
 55             } 
 56         } 
 57     
 58     }  
 59     int flag1=0, flag2=0, ans, ans1, ans2; 
 60     for( int p=i; p<=N; ++p ){
 61         if( d[p][N+1]&& (!d[p][N+2]) )flag1=1;
 62         if( (!d[p][N+2])&& d[p][N+2] )flag2=1;
 63         if( d[p][N+1]&& d[p][N+2] )    {
 64             puts( "Impossible" );
 65             return;    
 66         } 
 67     } 
 68     if( ! flag1 ){
 69         ans1=slove( N+1 ); 
 70     } 
 71     if( !flag2 ){
 72         ans2=slove( N+2 );    
 73     } 
 74     if( !flag1 && ! flag2 ){
 75         ans=min( ans1, ans2 );
 76         printf( "%d\n", ans );
 77         return;    
 78     } 
 79     if( (!flag1) && flag2 ){
 80         printf( "%d\n", ans1 );
 81         return;    
 82     }
 83     if( flag1 && (!flag2) ){
 84         printf( "%d\n", ans2 );
 85         return;    
 86     } 
 87 }
 88 int main( )
 89 {
 90     char s[4][4];
 91     while( scanf( "%s", s[0])!= EOF ){
 92         for(int i=1; i<4; ++ i )
 93             scanf( "%s", s[i] );
 94         memset( d, 0, sizeof d );
 95         for( int i=1; i<=16; ++ i )
 96             d[i][i]=1;
 97         for( int i=0; i<4;++ i ){
 98             for(int j=0; j<4; ++ j){
 99                 if( s[i][j]=='b')d[i*4+j+1][N+1]=1;
100                 if( s[i][j]=='w')d[i*4+j+1][N+2]=1;
101                 if( j>0 )d[4*i+j][4*i+1+j]=1; 
102                 if( i>0 )d[4*(i-1)+j+1][4*i+1+j]=1;
103                 if( j<3 )d[4*i+j+2][4*i+1+j]=1;
104                 if( i<3 )d[4*(i+1)+j+1][4*i+1+j]=1; 
105             }
106         }
107         
108         Gauss( );
109     }
110     return 0;
111 }

 

 

 

posted @ 2012-07-28 10:42  淡墨æ末央  阅读(151)  评论(0)  编辑  收藏  举报