Painter's Problem POJ - 1681

题意：

涂色，每涂一个格子就会将其临近的格子都改变，黄的变成白的，白的变成黄的，

给出画板初始情况，问最少多少次操作可以使得画板全部变成黄色

思路：

消元后，自由变元的取值情况的变化会造成后边的解的变化

剩下最后n个变元，遍历每种取值情况后，回溯找到解中操作次数最少的

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
using namespace std;
const int MAXN=300;
int a[MAXN][MAXN];
char m[MAXN][MAXN];
int s[MAXN][MAXN];
int x[MAXN];
int x_free[MAXN];
int n;

void pp(int p[MAXN][MAXN],int a,int b)
{
	for(int i=0;i<a;i++)
	for(int j=0;j<b;j++)
	printf("%d%c",p[i][j]," \n"[j==b-1]);
	printf("\n\n\n");
}

void Gauss(int equ,int var) //有equ个方程，var个变元。增广矩阵行数为equ,分别为0到equ-1,列数为var+1,分别为0到var.
{
	int i,j;
    int max_r;// 当前这列绝对值最大的行.
    int col;//当前处理的列
    int k;
    int len=0,temp;
     memset(x,0,sizeof(x));
	for( col=0,k=0;k<equ&&col<var;col++,k++)  
	{
		max_r=k;
		for(int i=k;i<equ;i++)
		if(a[i][col]==1)
		{
			max_r=i;
			break; 
		} 
		
		if(k!=max_r)
		for(int i=col;i<=var;i++)
		swap(a[k][i],a[max_r][i]);
		
		if(a[k][col]==0) 
		{
			x_free[len++]=col;
			
			k--;
			continue;
		}
		
		for(int i=k+1;i<equ;i++)
		{
			if(a[i][col]!=0)
			for(int j=col;j<=var;j++)
			a[i][j]^=a[k][j];
		}
	}
//无解判断 
	for(int i=k;i<equ;i++)
		if(a[i][var]!=0 ) 
	{
		printf("inf\n");
		return ;
	}
//遍历全部解集 
	int anss,tk=n*n;
	int l=1<<(var-k),p;

	for(int j=0;j<l;j++)
	{
	p=j;
	anss=0;
	for(int i=0;i<var-k;i++)
	{
		if(p&1) 
		{
		x[x_free[i]]=1;
		anss++;
		}
	    else
	    x[x_free[i]]=0;

	}

	for (int i = k - 1; i >= 0; i--)
    {
    	temp = a[i][var];
    	for(int jj=i+1;jj<var;jj++)
    	if(a[i][jj]) temp^=x[jj];
    	x[i]=temp;
         if(temp)anss++;
    }
        
    tk=min(anss,tk);
	}
    printf("%d\n",tk);
	

}
int fx[5]={0,0,0,-1,1};
int fy[5]={0,1,-1,0,0};

void doit(int i,int j)
{
	
	int s0,dx,dy;
	int s1=i*n+j;
	a[s1][n*n]=s[i][j]^1;
	
    for(int k=0;k<5;k++)
	{
		dx=i+fx[k];
		dy=j+fy[k];
		if(dx<0||dy<0||dx>=n||dy>=n) continue;
		s0=dx*n+dy;
		a[s0][s1]=1;
	} 	
}

int main()
{
	int T,tt;
	cin>>T;
	while(T--)
	{
		memset(a,0,sizeof(a));
		scanf("%d",&n);
	
		for(int i=0;i<n;i++)
		{
			getchar();
			scanf("%s",m[i]);	
	    }
		
		for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
	    if(m[i][j]=='y') s[i][j]=1;
	    else s[i][j]=0;
		

		for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
	    doit(i,j);
	    
	    Gauss(n*n,n*n);		
		
	}
	
	return 0;
}