链接: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;
}