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

来源：互联网发布：马布里最强数据编辑：程序博客网时间：2024/05/29 18:12

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

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

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

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

/*ID:cxlovePROB:poj 3185 DATA:2012.4.1HINT:高斯消元+枚举变元*/#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;}