题意:
给定一个n*m的矩阵,每个矩阵元素是一个数v,现在有q次查询,每次查询分两种情况:
①格式:1 x1,y1,x2,y2,求子矩阵(x1,y1)~(x2,y2)的异或和
②格式:2 x,y,z,修改矩阵位置(x,y)中的数值为z。
范围:
0<n,m<=500,v<=10^9,q<=2*10^5。
解法:
异或和满足条件:①假设原来值为v=a^b,若加入某个数c,则v=a^b^c=v^c,②假设原来为v=a^b,若取消值b,则v=a^b^b=v^b,即对同一个数异或两次,相当于没有进行操作。
可以据此条件,可以先求出(1,1)~(x,y)的异或和,则(x1,y1)~(x2,y2)的异或和为:sum(x2,y2)^sum(x1-1,y1-1)^sum(x1-1,y2)^sum(x2,y1-1)。
为了加快运算速度,通过树状数组解决,这里用到的是二维树状数组。
这题卡常数,最开始用cout输出结果,导致了TLE,改为printf后,就过了。所以在分析题目的过程中,要认真计算时间复杂度,当时间比较紧时,用printf加快运行效率。
具体代码:
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<vector>
#include<queue>
#include<cmath>
#include<set>
using namespace std;
#define N 505
typedef long long LL;
#define aa frist
#define bb second
#define mod 1000000007
int n,m,p;
int num[N][N];
int c[N+5][N+5];
int lowbit( int x )
{
return x & (-x);
}
void modify( int x, int y, int delta )
{
int i, j;
for(i=x; i<=N; i+=lowbit(i))
{
for(j=y; j<=N; j+=lowbit(j))
{
c[i][j] ^= delta;
}
}
}
int sum( int x, int y )
{
int res = 0, i, j;
for(i=x; i>0; i-=lowbit(i))
{
for(j=y; j>0; j-=lowbit(j))
{
res ^= c[i][j];
}
}
return res;
}
int main()
{
int t;
while(scanf("%d",&t)!=EOF)
{
while(t--)
{
memset(c,0,sizeof(c));
memset(num,0,sizeof(num));
scanf("%d%d%d",&n,&m,&p);
for(int i=1; i<=n; i++)
{
for(int j=1; j<=m; j++){
scanf("%d",&num[i][j]);
modify(i,j,num[i][j]);
}
}
while(p--)
{
int f,x1,x2,y1,y2,ans=0;
scanf("%d",&f);
if(f==1)
{
scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
ans=sum(x2,y2)^sum(x1-1,y1-1)^sum(x2,y1-1)^sum(x1-1,y2);
// cout<<ans<<endl;
// if(ans) cout<<"Yes"<<endl;
// else cout<<"No"<<endl;
if(ans) printf("Yes\n");
else printf("No\n");
}
else
{
int x,y,z;
scanf("%d%d%d",&x,&y,&z);
modify(x,y,z^num[x][y]);
num[x][y]=z;
}
}
}
}
return 0;
}
/*
1
1 2 10
1 2
1 1 1 1 2
2 1 2 1
1 1 1 1 2
*/