POJ 3185 The Water Bowls 高斯消元+枚举变元

转载请注明出处，谢谢 http://blog.csdn.net/ACM_cxlove?viewmode=contents by---cxlove

一如既往，继续高斯，这题和之前的略有不同，需要考虑多组解的情况。

当出现自由变元时，枚举所有解，取最优解。开始尝试用位运算，后来果断搜索吧。

不过这题貌似有问题，详见 http://blog.csdn.net/acm_cxlove/article/details/7424358

/*
ID:cxlove
PROB:poj 3185 
DATA:2012.4.1
HINT:高斯消元+枚举变元
*/
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int a[25][25],ans[25],var,var0,cnt;
int x[25];
void debug(){
	for(int i=0;i<20;i++){
		for(int j=0;j<20;j++)
			printf("%d ",a[i][j]);
		printf("%d\n",a[i][20]);
	}
}
int ta[25][25];
void dfs(int v){
	if(v==20){
		int temp=0;
		for(int i=0;i<20;i++)
			x[i]=ans[i];
		for(int i=0;i<20;i++)
			for(int j=0;j<=20;j++)
				ta[i][j]=a[i][j];
		for(int i=var0-1;i>=0;i--){
			for(int j=i+1;j<20;j++)
				ta[i][20]^=(x[j]&&ta[i][j]);
			x[i]=ta[i][20];
		}
		for(int i=0;i<20;i++)
			if(x[i])
				temp++;
		cnt=min(cnt,temp);
		return ;
	}
	ans[v]=0;
	dfs(v+1);
	ans[v]=1;
	dfs(v+1);
}
void gauss(){
	int i,j;
	for(i=0,j=0;i<20&&j<20;j++){
		int k=i;
		for(;k<20;k++)
			if(a[k][j])
				break;
		if(a[k][j]){
			for(int r=0;r<=20;r++)
				swap(a[i][r],a[k][r]);
			for(int k=0;k<20;k++)
				if(k!=i&&a[k][j]){
					for(int r=0;r<=20;r++)
						a[k][r]^=a[i][r];
				}
				i++;
		}
	}
	if(i==20){
		cnt=0;
		for(int r=0;r<20;r++)
			if(a[r][20])
				cnt++;
		printf("%d\n",cnt);
		return ;
	}
	var0=i;
	cnt=100;
	dfs(var0);
	printf("%d\n",cnt);
}
int main(){
	memset(a,0,sizeof(a));
	for(int i=0;i<20;i++)
		scanf("%d",&a[i][20]);
	for(int i=0;i<20;i++){
		a[i][i]=1;
		a[max(0,i-1)][i]=1;
		a[min(19,i+1)][i]=1;
	}
	gauss();
	return 0;
}