[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=,f=;char chr=getchar();
     while(!isdigit(chr)){if(chr=='-')f=-;chr=getchar();}
     while(isdigit(chr)) {ans=(ans<<)+(ans<<)+chr-;chr=getchar();}
     return ans*f;
 }const int M = ;
 int n,m,Q,s[][],v[][],sx,sy,tx,ty,ans;
 int dx[]={,,,-},
     dy[]={,,-,};
 char a[][];
 void dfs(int x,int y){
     v[x][y]=;
     for(int i=;i<;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]=='') continue;
         dfs(fx,fy);
     }
 }
 inline void BF(){
     while(Q--){
         memset(v,,sizeof(v));ans=;
         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]=='') ans++,dfs(i,j);
         printf("%d\n",ans);
     }
 }int s1[M][M],s2[M][M],s3[M][M];
 inline void Solve(){
     for(int i=;i<=n;i++)
         for(int j=;j<=m;j++)
             a[i][j]=(a[i][j]=='');
     for(int i=;i<=n;i++)
         for(int j=;j<=m;j++)
             s1[i][j]=a[i][j]+s1[i-][j]+s1[i][j-]-s1[i-][j-],
             s2[i][j]=(a[i][j]&&a[i+][j])+s2[i-][j]+s2[i][j-]-s2[i-][j-],
             s3[i][j]=(a[i][j]&&a[i][j+])+s3[i-][j]+s3[i][j-]-s3[i-][j-];
     while(Q--){
         sx=read(),sy=read(),tx=read(),ty=read();
         int V=s1[tx][ty]-s1[sx-][ty]-s1[tx][sy-]+s1[sx-][sy-];
         int E=s2[tx-][ty]-s2[tx-][sy-]-s2[sx-][ty]+s2[sx-][sy-]
              +s3[tx][ty-]-s3[sx-][ty-]-s3[tx][sy-]+s3[sx-][sy-];
         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=;i<=n;i++)for(int j=;j<=m;j++)cin>>a[i][j];
 //    BF();
     Solve();
     return ;
 }