poj1681 Painter's Problem<高斯消元>

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

#include <stdio.h> 
#include <string.h>
#include <algorithm>
#include <cmath>
using namespace std;
int d[230][230], N, M;
char s[16][16]; 
void solve( int n)
{
    int x[230],m[20], t, ans=1000, temp;
    t=M-n;
    for( int i=1; i<=pow( 2, t ); ++i ){
        memset( x, 0, sizeof x  );
        for(int j=0; j<t; ++j){
            m[j]=(i&(1<<j));
            if(m[j]){
                x[n+j]=1;
            }    
        }
        temp=0;
        for( int p=n+1;p>=1;--p ){
            for(int k=i+1;k<=M;++k){
                    x[p]^=( x[k]&&d[p][k] );
                
            }
            x[p]^=d[p][M+1];
            if( x[p])
                temp++;
        } 
        ans=min( ans, temp );
    }
    printf( "%d\n", ans );
}
void Gauss(  )
{
    int i=1, j, k , p, t;
    for(j=1; j<=M; ++j){
        for(p=i; p<=M; ++p){
            if(d[p][j])break;
        }
        if( p>M )continue;
        if(p!=i){
            for(k=j;k<=M+1; ++k){
                t=d[p][k],d[p][k]=d[i][k],d[i][k]=t;
            }    
        }
        for(p=i+1; p<=M; ++p){
            if( d[p][j] ){
                for( k=j; k<=M+1; ++k ){
                    d[p][k]^=d[i][k];
                }
            }
        }
        ++i;
    }
    for( p=i; p<=M; ++p ){
        if( d[p][M+1] ){
            puts("inf");
            return;
        }
    }    
    solve( i-1 );
}
int main( )
{
    int T; 
    scanf("%d", &T );
    while( T --) {
        scanf( "%d", &N );
        M=N*N;
        for( int i=0; i<N; ++ i ){
            scanf( "%s", s[i] ); 
        } 
        memset( d, 0, sizeof d );
        for( int i=1; i<=M; ++ i )
            d[i][i]=1;
        for( int i=0; i<N; ++i ){
            for( int j=0;j<N; ++j ){
                if( s[i][j]=='w' ) d[i*N+j+1][M+1]=1;
                
                if( i>0 )d[(i-1)*N+j+1][i*N+j+1]=1;
                if( j>0 )d[i*N+j][i*N+j+1]=1;
                if( i<N-1 )d[(i+1)*N+j+1][i*N+j+1]=1;
                if( j<N-1 )d[i*N+j+2][i*N+j+1]=1;        
            } 
        }
        Gauss( ); /*
        for( int i=1; i<=M; ++i ){
            for(int j=1;j<=M; ++ j ){
                printf( "%d ", d[i][j] );
            }
            puts( "" );
        }*/
    } 
    return 0;    
}