牛客寒假算法基础集训营1-E.小a的轰炸游戏（差分思想解题）

链接：https://ac.nowcoder.com/acm/contest/317/E
来源：牛客网
 

题目描述

小a正在玩一款即时战略游戏，现在他要用航空母舰对敌方阵地进行轰炸
地方阵地可以看做是n×mn×m的矩形
航空母舰总共会派出qq架飞机。
飞机有两种，第一种飞机会轰炸以(xi,yi)(xi,yi)为中心，对角线长为lili的正菱形(也就是两条对角线分别于xx轴 yy轴平行的正方形)，而第二种飞机只会轰炸正菱形的上半部分(包括第xixi行)
(具体看样例解释)
现在小a想知道所有格子被轰炸次数的异或和
注意：不保证被轰炸的格子一定在矩形范围内，若越界请忽略

输入描述:

第一行三个整数n,m,qn,m,q，分别表示矩阵的长/宽/询问次数
接下来qq行，每行四个整数opt,x,y,lopt,x,y,l，表示飞机类型，轰炸的坐标，以及对角线长度
保证ll为奇数！

输出描述:

一个整数，表示所有格子被轰炸次数的异或和

示例1

输入

复制

4 5 4
1 2 2 1
1 3 3 5
1 3 2 3
2 2 4 3

输出

复制

2

说明

每次的操作矩阵即操作后的矩阵的值如下

0 0 0 0 0

0 1 0 0 0
0 0 0 0 0
0 0 0 0 0

0 0 1 0 0
0 2 1 1 0
1 1 1 1 1
0 1 1 1 0

0 0 1 0 0
0 3 1 1 0
2 2 2 1 1
0 2 1 1 0

0 0 1 1 0
0 3 2 2 1
2 2 2 1 1
0 2 1 1 0
最后把所有元素异或后为2

备注:

1⩽n,m⩽10001⩽n,m⩽1000

1⩽q⩽5∗1051⩽q⩽5∗105

保证opt=1/2,1⩽x,y,l⩽max(N,M)opt=1/2,1⩽x,y,l⩽max(N,M)

读入文件过大，请使用较快的读入方式

 

思路：这个题目的话，基本上就是一个比较巧妙的差分方法，因为他是菱形的，形状虽然不像是正方形长方形那样我们运算的时候规则，但是也算是有规律可循的，

首先说一下差分吧，其实就是前缀和的思想，举个例子：

    像这样，我们通过修改a数组中的数字，就可以用sum数组来表示其中从L到R加入某个值的效果。这基本上就是差分了，

.

不过这个题目是二维的差分，其实都差不多，

再来一个图片

这里的a,b,c,d数组就是用来维护这个形状的图形的

详情看一下代码吧，这个图我觉得已经够清楚的了

代码：

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
using namespace std;

const int maxn=3010,L=1010;
int a[maxn][maxn],b[maxn][maxn];
int c[maxn][maxn],d[maxn][maxn];
int n,m,q;
void uptri(int x,int y,int l)
{
    a[x-l/2][y]++,b[x-l/2][y+1]--;
    a[x+1][y-l/2-1]--,b[x+1][y+l/2+2]++;
}
void downtri(int x,int y,int l)
{
    c[x+1][y-l/2+1]++,d[x+1][y+l/2]--;
    c[x+l/2+1][y+1]--,d[x+l/2+1][y]++;
}
int main()
{
    scanf("%d%d%d",&n,&m,&q);
    for(int i=1; i<=q; i++)
    {
        int opt,x,y,l;
        scanf("%d%d%d%d",&opt,&x,&y,&l);
        x+=L,y+=L;
        uptri(x,y,l);
        if(opt==1)
            downtri(x,y,l);
    }
    int ans=0;
    int flag=0;
    for(int i=1; i<=2*L; i++)
    {
        int t=0;
        for(int j=1; j<=2*L+m; j++)
        {
            t+=a[i][j]+b[i][j]+c[i][j]+d[i][j];
            if(i>=L+1&&i<=L+n&&j>=L+1&&j<=L+m)
            {
                ans^=t;
                //cout<<t<<' ';
            }
            a[i+1][j-1]+=a[i][j];
            b[i+1][j+1]+=b[i][j];
            c[i+1][j+1]+=c[i][j];
            d[i+1][j-1]+=d[i][j];
        }
        //if(i>=L+1&&i<=L+n)cout<<endl;
    }
    cout<<ans<<endl;
}