#include<stdio.h>
#include<algorithm>
#include<iostream>
#include<string.h>
using namespace std;
int a[230][230],x[230],free_x[230],free_num,r,c;
int Gauss()
{
int k,col,max_r,i,j;
for(free_num=k=col=0;k<r&&col<c;col++) //枚举行,列从0开始
{
max_r=k;
for(i=k+1;i<r;i++)
if(abs(a[i][col])>abs(a[max_r][col]) )max_r=i; //找到这一列里面绝对值最大的
if(!a[max_r][col])
{
free_x[free_num++]=col;
continue; //最大的都是0,那么这一列不用考虑了
}
if(max_r!=k)
for(i=col;i<=c;i++) swap(a[max_r][i],a[k][i]); //否则交换
for(i=k+1;i<r;i++)
if(a[i][col]!=0)
for(j=col;j<=c;j++) a[i][j]^=a[k][j];
k++;
}
for(i=k;i<r;i++)
if(a[i][col]) return -1;
if(k<c) return c-k;
for(i=c-1;i>=0;i--)
{
x[i]=a[i][c];
for(j=i+1;j<c;j++)
x[i]^=(a[i][j]&&x[j]);
}
return 0;
}
int main()
{
int i,j,T,n,v,ans;
char s[16];
scanf("%d",&T);
while(T--)
{
memset(a,0,sizeof(a));
memset(x,0,sizeof(x));
scanf("%d",&n);
for(i=0;i<n;i++)
{
scanf("%s",s);
for(j=0;j<n;j++)
if(s[j]=='y') a[i*n+j][n*n]=0;
else a[i*n+j][n*n]=1;
}
c=r=n*n;
for(i=0;i<n;i++)
for(j=0;j<n;j++)
{
v=i*n+j;
a[v][v]=1;
if(i>0) a[v-n][v]=1;
if(i<n-1) a[v+n][v]=1;
if(j>0) a[v-1][v]=1;
if(j<n-1) a[v+1][v]=1;
}
v=Gauss();
if(v<0)
{
printf("inf\n");
continue;
}
if(v==0)
{
ans=0;
for(i=0;i<n*n;i++)
if(x[i]) ans++;
printf("%d\n",ans);
}
else
{
ans=0x7fffffff;
int tot=1<<v;
for(i=0;i<tot;i++)
{
int cnt=0;
for(j=0;j<v;j++)
if(i&(1<<j))
{
cnt++;
x[free_x[j]]=1;
}
else x[free_x[j]]=0;
for(j=c-v-1;j>=0;j--)
{
int p;
for(p=j;p<c;p++)
if(a[j][p]) break;
x[p]=a[j][c];
for(int l=p+1;l<c;l++)
if(a[j][l]) x[p]^=x[l];
if(x[p]) cnt++;
}
ans=min(ans,cnt);
}
printf("%d\n",ans);
}
}
}
高斯消元枚举自由元——poj1222
最新推荐文章于 2022-11-25 00:02:07 发布