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
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#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<


你可能感兴趣的:(online,judge,Hdu,数据结构,树状数组,其他,前缀和,数学,容斥)