[NOIP2019模拟赛][AT2381] Nuske vs Phantom Thnook

题目链接

评测姬好快啊（港记号?）暴力40pts变成60pts

因为题目说了保证蓝色点两两之间只有一条路径，所以肯定组成了一棵树，而对于每次询问的x1,y1,x2,y2的子矩阵中就存在着一个森林

不难知道对于一个森林，其中树（联通块）的数量为$V-E$(V为节点数，E为边数)

也就是说对于每一个询问，只要求出蓝色节点数减去边数的答案就好了

点数和边数都可以用二维前缀和求，其中边可以分横边和竖边分别记录

 

#pragma GCC optimize("Ofast")
#include<bits/stdc++.h>
using namespace std;
inline int read(){
    int ans=0,f=1;char chr=getchar();
    while(!isdigit(chr)){if(chr=='-')f=-1;chr=getchar();}
    while(isdigit(chr)) {ans=(ans<<3)+(ans<<1)+chr-48;chr=getchar();}
    return ans*f;
}const int M = 2005;
int n,m,Q,s[2005][2005],v[2005][2005],sx,sy,tx,ty,ans;
int dx[]={0,1,0,-1},
    dy[]={1,0,-1,0};
char a[2005][2005];
void dfs(int x,int y){
    v[x][y]=1;
    for(int i=0;i<4;i++){
        int fx=x+dx[i],fy=y+dy[i];
        if(v[fx][fy]||fx<sx||fx>tx||fy<sy||fy>ty||a[fx][fy]=='0') continue;
        dfs(fx,fy);
    }
}
inline void BF(){
    while(Q--){
        memset(v,0,sizeof(v));ans=0;
        sx=read(),sy=read(),tx=read(),ty=read();
        for(int i=sx;i<=tx;i++)
            for(int j=sy;j<=ty;j++)
                if(!v[i][j]&&a[i][j]=='1') ans++,dfs(i,j);
        printf("%d\n",ans);
    }
}int s1[M][M],s2[M][M],s3[M][M];
inline void Solve(){
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            a[i][j]=(a[i][j]=='1');
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            s1[i][j]=a[i][j]+s1[i-1][j]+s1[i][j-1]-s1[i-1][j-1],
            s2[i][j]=(a[i][j]&&a[i+1][j])+s2[i-1][j]+s2[i][j-1]-s2[i-1][j-1],
            s3[i][j]=(a[i][j]&&a[i][j+1])+s3[i-1][j]+s3[i][j-1]-s3[i-1][j-1];
    while(Q--){
        sx=read(),sy=read(),tx=read(),ty=read();
        int V=s1[tx][ty]-s1[sx-1][ty]-s1[tx][sy-1]+s1[sx-1][sy-1];
        int E=s2[tx-1][ty]-s2[tx-1][sy-1]-s2[sx-1][ty]+s2[sx-1][sy-1]
             +s3[tx][ty-1]-s3[sx-1][ty-1]-s3[tx][sy-1]+s3[sx-1][sy-1];
        printf("%d\n",V-E);
    }
    
}
int main(){
    freopen("wang.in","r",stdin);
    freopen("wang.out","w",stdout);
    n=read(),m=read(),Q=read();
    for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)cin>>a[i][j];
//    BF();
    Solve();
    return 0;
}