链接:http://poj.org/problem?id=1681
题意:与poj1222相似 输出变为求最小步数。
分析:由于求最小步数,当我们存在自由变元时,我们可以得到解,但是不知道需要的最小步数,所以我们枚举所有的变元。
代码:
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<string>
#include<vector>
#include<queue>
#include<cmath>
#include<stack>
#include<set>
#include<map>
#define INF 0x3f3f3f3f
#define Mn 230
#define Mm 2000005
#define mod 1000000007
#define CLR(a,b) memset((a),(b),sizeof((a)))
#define CLRS(a,b,Size) memset((a),(b),sizeof((a[0]))*(Size+1))
#define CPY(a,b) memcpy ((a), (b), sizeof((a)))
#pragma comment(linker, "/STACK:102400000,102400000")
#define ul u<<1
#define ur (u<<1)|1
using namespace std;
typedef long long ll;
int a[Mn][Mn];
int x[Mn];
int free_x[Mn];
int var,equ;
int gauss() {
int k,col,max_r,cnt=0;
for(col=0,k=0;k<equ&&col<var;k++,col++) {
max_r=k;
for(int i=k+1;i<equ;i++)
if(abs(a[i][col])>abs(a[max_r][col])) max_r=i;
if(max_r!=k)
for(int j=col;j<=var;j++) swap(a[k][j],a[max_r][j]);
if(a[k][col]==0) {
k--;
free_x[cnt++]=col;
continue;
}
for(int i=k+1;i<equ;i++) {
if(a[i][col]) {
for(int j=col;j<=var;j++) {
a[i][j]^=a[k][j];
}
}
}
}
for(int i=k;i<equ;i++) {
if(a[i][col]) return -1;
}
int res=INF;
for(int s=0;s<(1<<cnt);s++) {
int num=0;
int s0=s;
for(int i=0;i<cnt;i++) {
if(x[free_x[i]]=s0&1) num++;
s0>>=1;
}
for(int i=k-1;i>=0;i--) {
int tmp=a[i][var];
for(int j=i+1;j<var;j++) {
tmp^=(a[i][j]&&x[j]);
}
if(x[i]=tmp) num++;
}
res=min(res,num);
}
return res;
}
int n;
void init() {
CLR(a,0);
var=n*n,equ=n*n;
for(int i=0;i<n;i++) {
for(int j=0;j<n;j++) {
int k=i*n+j;
a[k][k]=1;
if(i>0) a[k-n][k]=1;
if(i<n-1) a[k+n][k]=1;
if(j>0) a[k-1][k]=1;
if(j<n-1) a[k+1][k]=1;
}
}
}
int main() {
int t;
scanf("%d",&t);
for(int cas=1;cas<=t;cas++) {
scanf("%d",&n);
init();
for(int i=0;i<n;i++) {
getchar();
for(int j=0;j<n;j++) {
char c=getchar();
if(c=='w') a[i*n+j][n*n]=1;
}
}
int ans=gauss();
if(ans==-1) printf("inf\n");
else printf("%d\n",ans);
}
return 0;
}