hdu5892Resident Evil

链接：http://acm.hdu.edu.cn/showproblem.php?pid=5892

题意：给定n和m表示有一个n*n的矩阵和m个操作，操作1：给出左上角的位置[x1,y1]和右下角的位置[x2,y2]，然后给定k表示有k对[a,b]接下来在这个给定的矩阵中每个格子中都添加b个a类物品；操作2：给定左上角位置[x1,y1]和右下角位置[x2,y2]，求所有物品在这个矩阵中的奇偶情况。

分析：首先因为物品只要50种，状态压缩和矩阵的容斥问题都很容易想到。我们把问题转化成求[1,1]~[x,y]的矩阵中的奇偶情况，但是如果自己画图看下操作1的影响就会发现，只有相隔1行或一列的那些点才会有作用。我们先简化一下这个问题，如果是一维的怎么处理呢？区间添加，区间查询。我们可以将区间添加变成单点添加，我们会发现如果在区间[x,y]上添加一个物品，那么它奇数的情况只会在x,x+2,x+4...这些间隔的位置上才有影响，看到这里我们可以想到将奇偶两种位置拆成两个数组维护就好了，而且前缀和是可以异或的。然后我们将问题回到二维上，根据[1,1]~[x,y]的x和y的奇偶我们可以将矩阵分成4份，我们只要每次维护对应矩阵就行了。如果还不懂的话就参考下代码吧。

代码：

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<bitset>
#include<math.h>
#include<vector>
#include<string>
#include<stdio.h>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=1510;
const int M=50010;
const int mod=1000000007;
const int MOD1=1000000007;
const int MOD2=1000000009;
const double EPS=0.00000001;
typedef long long ll;
const ll MOD=1000000007;
const int INF=1000000010;
const ll MAX=1ll<<55;
const double eps=1e-5;
const double inf=~0u>>1;
const double pi=acos(-1.0);
typedef double db;
typedef unsigned int uint;
typedef unsigned long long ull;
int mx;
char s[3];
ll B,t[55],f[2][2][N][N];
ll getans(int x,int y) {
    int one=x&1,two=y&1;ll ret=0;
    int i,j,X=x/2+one,Y=y/2+two;
    for (i=X;i;i-=i&-i)
        for (j=Y;j;j-=j&-j) ret^=f[one][two][i][j];
    return ret;
}
void add(int x,int y) {
    int one=x&1,two=y&1;
    int i,j,X=x/2+one,Y=y/2+two;
    for (i=X;i<=mx;i+=i&-i)
        for (j=Y;j<=mx;j+=j&-j) f[one][two][i][j]^=B;
}
int main()
{
    int i,j,k,n,m,x,y,x1,y1,x2,y2;
    for (i=0;i<=50;i++) t[i]=1ll<<i;
    while (scanf("%d%d", &n, &m)!=EOF) {
        memset(f,0,sizeof(f));mx=n/2+1;
        for (i=1;i<=m;i++) {
            scanf("%s%d%d%d%d", s, &x1, &y1, &x2, &y2);
            if (s[0]=='Q') {
                x1--;y1--;
                B=getans(x1,y1)^getans(x2,y2)^getans(x1,y2)^getans(x2,y1);
                for (j=1;j<=50;j++)
                if (B&t[j]) printf("2 ");
                else printf("1 ");
                printf("\n");
            } else {
                scanf("%d", &k);B=0;
                for (j=1;j<=k;j++) {
                    scanf("%d%d", &x, &y);
                    if (y&1) B^=t[x];
                }
                x2++;y2++;
                add(x1,y1);add(x2,y2);
                add(x1,y2);add(x2,y1);
            }
        }
    }
    return 0;
}