问题 B: 任(duty)
时间限制: 2 Sec 内存限制: 512 MB
题目描述
liu_runda退役之后就失去梦想开始咸鱼生活了…
Bilibili夏日画板活动中,所有人都可以在一块画板上进行像素画创作.UOJ群有一群无聊的人决定在画板上创作一个50*50的UOJ的LOGO.如下图.
这块画板实际上是很大的矩形网格.一个网格是一像素.
一个人每三分钟才能画一个像素.所以liu_runda的咸鱼生活非常无聊.
郭神表示他实在是看不下去liu_rudna这只颓狗了,于是随手出了一道神题,liu_runda不会做,于是给出到联考里了.
在画板上有一片黑白相间的矩形区域满足这样的性质:如果认为相同颜色的方块可以在上下左右四个方向连通,那么任意两个黑色方块要么不连通,要么连通但之间只有一条简单路径(不重复经过同一个格子的路径).
这个矩形区域有N行M列,从上到下依次为第1,2,3…N-1,N行,从左到右依次为第1,2,3…M-1,M列.
每次郭神会询问这片矩形区域内的一个子矩形.在只考虑这个子矩形内的像素时(即从子矩形内部不能和子矩形之外的像素相连通),问这个子矩形内的黑色方块组成了多少连通块.
如果不能完成这个任务,liu_runda就会被郭神批判一番…
【输入格式】
第一行三个整数N,M,Q,表示矩形区域有N行M列,有Q组询问.
接下来N行,每行一个长为M的01字符串.0表示白色,1表示黑色.第i行第j个字符表示第i行j列的颜色,
接下来Q行,每行4个整数x1,y1,x2,y2,(x1<=x2,y1<=y2)表示选出的矩形区域的两个对角.即选出一个左上角为第x1行第y1列,右下角为第x2行第y2列,包含x2-x1+1行,y2-y1+1列的区域.
【输出格式】
Q行,第i行一个整数ans表示第i组询问的答案.
【样例输入1】
3 4 4
1101
0110
1101
1 1 3 4
1 1 3 1
2 2 3 4
1 2 2 4
【样例输出1】
3
2
2
2
【样例输入2】
5 5 6
11010
01110
10101
11101
01010
1 1 5 5
1 2 4 5
2 3 3 4
3 3 3 3
3 1 3 5
1 1 3 4
【样例输出2】
3
2
1
1
3
2
【数据范围】
对于第1,2个测试点,Q=1
对于第3,4个测试点,N=1
对于第5,6,7个测试点,N=2
对于第8个测试点,N,M<=1000
对于第9个测试点,N,M<=1500
对于全部测试点,1<=N,M<=2000,1<=Q<=200000,1<=x1<=x2<=N,1<=y1<=y2<=M,保证任意两个黑色像素之间最多只有一条简单路径.
其实题面里还有张很无聊的图，懒得粘了。。
因为保证联通块里不存在环，只能由一条边互相连通，而同一联通块满足点数-边数==1；所以求出区间里所有点和所有边，相减即为答案。
因此我们可以维护前缀和。但有些细节。
举个样例（为了方便说我把0，1改成了点的标号）
1 2 3 4
5 6 7 8
9 0 10 11
比如我们要求7，8，10，11的，那么要把3-7，4-8，6-7，0-10的边都剪掉，因此会多剪掉7-10和7-8，所以得额外维护单列的值，把它加回去，具体就不细说了。
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#define N 2005
#define inf 1000000000
using namespace std;
int read()
{
int sum=0,f=1;char x=getchar();
while(x<'0'||x>'9'){if(x=='-')f=-1;x=getchar();}
while(x>='0'&&x<='9'){sum=(sum<<1)+(sum<<3)+x-'0';x=getchar();}
return sum*f;
}
int n,m,q;
int a[2005][2005],dian[2005][2005],bian[2005][2005],h[2005][2005],z[2005][2005];
int l1,l2,r1,r2;
void init()
{
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
{
dian[i][j]=dian[i][j-1]+a[i][j];
bian[i][j]=bian[i][j-1];
if(a[i][j]==1&&a[i][j-1]==1)bian[i][j]++;
h[i][j]=bian[i][j];
}
for(int j=1;j<=m;j++)
for(int i=1;i<=n;i++)
{
dian[i][j]+=dian[i-1][j];
z[i][j]=z[i-1][j];
if(a[i][j]==1&&a[i-1][j]==1)z[i][j]++;
}
for(int i=1;i<=n;i++)
{
int s=0;
for(int j=1;j<=m;j++)
{
if(a[i][j]==1&&a[i-1][j]==1)s++;
bian[i][j]+=bian[i-1][j]+s;
}
}
/* for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
cout<<z[i][j]<<" ";
cout<<endl;
}*/
}
int main()
{
n=read();m=read();q=read();
char s[2005];
for(int i=1;i<=n;i++)
{
scanf("%s",s+1);
for(int j=1;j<=m;j++)
a[i][j]=s[j]-'0';
}
init();int ans2,ans1;
while(q--)
{
l1=read();r1=read();l2=read();r2=read();
ans1=dian[l2][r2]-dian[l2][r1-1]-dian[l1-1][r2]+dian[l1-1][r1-1];
ans2=bian[l2][r2]-bian[l2][r1]-bian[l1][r2]+bian[l1][r1];
ans2+=z[l2][r1]-z[l1][r1]+h[l1][r2]-h[l1][r1];
//cout<<ans1<<" "<<ans2<<endl;
printf("%d
",ans1-ans2);
}
}