poj1753Flip Game<高斯消元>

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

#include <stdio.h>
#include <stdlib.h>
#include <string.h> 
#include <algorithm>
using namespace std;
const int N=16;
int d[20][20];
int slove( int n )   //枚举自由解 求最小值 
{
    int x[20], temp, ans=100, m[4], t, p, k; 
    for(  t=1; t<=16; ++t ){
        memset( x, 0, sizeof x ); 
        temp=0; 
        for(  p=0; p<4; ++p ){
            m[p]=(t&( 1<<p ));
            if( m[p] ){
                x[13+p]=1;    
            }    
        }
        for( p=N; p>0; --p ){
            for(  k=N; k>0; --k ){
                if( p!=k )
                    x[p] ^= (x[k]*d[p][k]);    
            }
            x[p] ^= d[p][n];
            if (x[p]) {
                ++temp;
            }    
        }    
        ans=min( ans, temp ); 
    } 
    return ans;    
} 
void Gauss( )
{
    int i=1, j, k, p, t;
     
    for( j=1; j<=N; ++i,++j ){
        for( p=i; p<=N; p++ ){
                if(d[p][j])break;
            }
        if( p==N+1 ){--i;continue;} 
        if( p!=i ){
            for( k=j; k<=N+2; ++k ){
                t=d[i][k];
                d[i][k]=d[p][k];
                d[p][k]=t;    
            }    
        } 
        for( k=i+1; k<=N; ++k ){
            if( d[k][j] ){ 
                for( t=j; t<=N+2; ++t ){
                    d[k][t]^=d[i][t]; 
                }    
            } 
        } 
    
    }  
    int flag1=0, flag2=0, ans, ans1, ans2; 
    for( int p=i; p<=N; ++p ){
        if( d[p][N+1]&& (!d[p][N+2]) )flag1=1;
        if( (!d[p][N+2])&& d[p][N+2] )flag2=1;
        if( d[p][N+1]&& d[p][N+2] )    {
            puts( "Impossible" );
            return;    
        } 
    } 
    if( ! flag1 ){
        ans1=slove( N+1 ); 
    } 
    if( !flag2 ){
        ans2=slove( N+2 );    
    } 
    if( !flag1 && ! flag2 ){
        ans=min( ans1, ans2 );
        printf( "%d\n", ans );
        return;    
    } 
    if( (!flag1) && flag2 ){
        printf( "%d\n", ans1 );
        return;    
    }
    if( flag1 && (!flag2) ){
        printf( "%d\n", ans2 );
        return;    
    } 
}
int main( )
{
    char s[4][4];
    while( scanf( "%s", s[0])!= EOF ){
        for(int i=1; i<4; ++ i )
            scanf( "%s", s[i] );
        memset( d, 0, sizeof d );
        for( int i=1; i<=16; ++ i )
            d[i][i]=1;
        for( int i=0; i<4;++ i ){
            for(int j=0; j<4; ++ j){
                if( s[i][j]=='b')d[i*4+j+1][N+1]=1;
                if( s[i][j]=='w')d[i*4+j+1][N+2]=1;
                if( j>0 )d[4*i+j][4*i+1+j]=1; 
                if( i>0 )d[4*(i-1)+j+1][4*i+1+j]=1;
                if( j<3 )d[4*i+j+2][4*i+1+j]=1;
                if( i<3 )d[4*(i+1)+j+1][4*i+1+j]=1; 
            }
        }
        
        Gauss( );
    }
    return 0;
}