bzoj4171 or 省队集训day3 chess: Rhl的游戏

【题目描述】

RHL最近迷上一个小游戏：Flip it。游戏的规则很简单，在一个N*M的格子上，有一些格子是黑色，有一些是白色。每选择一个格子按一次，格子以及周围边相邻的格子都会翻转颜色（边相邻指至少与该格子有一条公共边的格子），黑变白，白变黑。

RHL希望把所有格子都变成白色的。不幸的是，有一些格子坏掉了，无法被按下。这时，它可以完成游戏吗？

【输入格式】

第一行一个整数T，表示T组数据。

每组数据开始于三个整数n,m,k，分别表示格子的高度和宽度、坏掉格子的个数。接下来的n行，每行一个长度m的字符串，表示格子状态为’B’或‘W’。最后k行，每行两个整数Xi,Yi(1≤Xi≤n,1≤Yi≤m)，表示坏掉的格子。

【输出格式】

对于每组数据，先输出一行Case #i: (1≤i≤T)

如果可以成功，输出YES，否则输出NO。

【样例输入】

2

3 3 0

WBW

BBB

WBW

3 3 2

WBW

BBB

WBW

2 2

3 2

【样例输出】

Case #1:

YES

Case #2:

NO

【数据范围】

30%,n,m,k<=10

100%,n,m,k<=256,T<=10
http://www.cnblogs.com/chenyushuo/p/4685182.html
和这个类似的设个xor方程组，对于不能按的方块，直接将它定为0即可

code：

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#include<algorithm>
#define maxn 260
using namespace std;
char ch,s[maxn];
int T,n,m,N,M,k,x,y;
unsigned int c[maxn][maxn][maxn>>5],a[maxn<<1][maxn>>5];
bool col[maxn][maxn],ok,d[maxn][maxn],b[maxn<<1];
inline void read(int &x){
for (ok=0,ch=getchar();!isdigit(ch);ch=getchar()) if (ch=='-') ok=1;
for (x=0;isdigit(ch);x=x*10+ch-'0',ch=getchar());
if (ok) x=-x;
}
bool gauss(){
int i,j,k,p,q;
for (i=0,k=1;i<N;i++){
for (p=(1<<(i&31)),j=k;j<=M&&!(a[j][i>>5]&p);j++);
if (j<=M){
for (q=(i>>5);q<=((N-1)>>5);q++) swap(a[k][q],a[j][q]);
swap(b[k],b[j]);
for (j=j+1;j<=M;j++)
if (a[j][i>>5]&p){
for (q=(i>>5);q<=((N-1)>>5);q++) a[j][q]^=a[k][q];
b[j]^=b[k];
}
k++;
}
}
for (;k<=M;k++) if (b[k]) return false;
return true;
}
int main(){
read(T);
for (int t=1;t<=T;t++){
read(m),read(n),read(k),N=n,M=n;
for (int i=1;i<=m;i++){
scanf("%s",s+1);
for (int j=1;j<=n;j++) col[i][j]=(s[j]=='B');
}
for (int i=1;i<=n;i++) c[1][i][(i-1)>>5]=(1<<((i-1)&31));
for (int i=2;i<=m;i++)
for (int j=1;j<=n;j++){
for (int k=0;k<=((n-1)>>5);k++)
c[i][j][k]=c[i-1][j-1][k]^c[i-1][j][k]^c[i-1][j+1][k]^c[i-2][j][k];
d[i][j]=d[i-1][j-1]^d[i-1][j]^d[i-1][j+1]^d[i-2][j]^col[i-1][j];
}
for (int i=1;i<=n;i++){
for (int j=0;j<=((n-1)>>5);j++)
a[i][j]=c[m][i][j]^c[m][i-1][j]^c[m][i+1][j]^c[m-1][i][j];
b[i]=col[m][i]^d[m][i-1]^d[m][i]^d[m][i+1]^d[m-1][i];
}
while (k--){
read(x),read(y),++M;
for (int i=0;i<=((n-1)>>5);i++) a[M][i]=c[x][y][i];
b[M]=d[x][y];
}
printf("Case #%d:\n",t);
if (gauss()) puts("YES");
else puts("NO");
}
return 0;
}