Problem Statement

Nuske has a grid with N rows and M columns of squares. The rows are numbered 1 through N from top to bottom, and the columns are numbered 1 through M from left to right. Each square in the grid is painted in either blue or white. If Si,j is 1, the square at the i-th row and j-th column is blue; if Si,j is 0, the square is white. For every pair of two blue square a and b, there is at most one path that starts from a, repeatedly proceeds to an adjacent (side by side) blue square and finally reachesb, without traversing the same square more than once.

Phantom Thnook, Nuske‘s eternal rival, gives Q queries to Nuske. The i-th query consists of four integers xi,1, yi,1, xi,2 and yi,2 and asks him the following: when the rectangular region of the grid bounded by (and including) the xi,1-th row, xi,2-th row, yi,1-th column and yi,2-th column is cut out, how many connected components consisting of blue squares there are in the region?

Process all the queries.

Constraints

• 1≤N,M≤2000
• 1≤Q≤200000
• Si,j is either 0 or 1.
• Si,j satisfies the condition explained in the statement.
• 1≤xi,1≤xi,2≤N(1≤i≤Q)
• 1≤yi,1≤yi,2≤M(1≤i≤Q)

Sample Input 1

3 4 41101011011011 1 3 41 1 3 12 2 3 41 2 2 4

Sample Output 1

3222

In the first query, the whole grid is specified. There are three components consisting of blue squares, and thus 3 should be printed.

In the second query, the region within the red frame is specified. There are two components consisting of blue squares, and thus 2 should be printed. Note that squares that belong to the same component in the original grid may belong to different components.

Sample Input 2

5 5 611010011101010111101010101 1 5 51 2 4 52 3 3 43 3 3 33 1 3 51 1 3 4

Sample Output 2

321132

 题意: n*m (n,m<=2000) 的地图，q (q<=200000)次询问，地图每个格子非蓝即白，上下相连，或者左右相连算作连续，还有，每两个蓝色格子，如果连通，保证只有一条确定路径相连，对于每次询问，求范围矩形蓝色连通块个数。

//题解：这题，题意都读了很久。。。常规求连通的算法会超时，因为Q很大。对于每两个连通蓝色格子，只有一条路径连通很关键，每个连通块就变成了一颗树，，所以，即求森林个数，树的话有个性质，点数减边数等于1，所以，即求范围内点数减边数，查询时间复杂度为 O(1)

思路明白，代码比较好写

 1 #include <iostream> 2 #include <stdio.h> 3 #include <string.h> 4 using namespace std; 5 #define LL long long 6 #define INF 0x3f3f3f3f 7 #define MX 2005 8  9 int n,m,q;10 int ans;11 char mp[MX][MX];12 int d[MX][MX];  //点13 int x[MX][MX];  //线14 int up[MX][MX]; //上相连15 int lef[MX][MX];//左相连16 17 void slv()18 {19     memset(x,0,sizeof(x));20     memset(d,0,sizeof(d));21     memset(up,0,sizeof(up));22     memset(lef,0,sizeof(lef));23     for (int i=1;i<=n;i++)24     {25         for (int j=1;j<=m;j++)26         {27             d[i][j]=d[i-1][j]+d[i][j-1]-d[i-1][j-1];28             x[i][j]=x[i-1][j]+x[i][j-1]-x[i-1][j-1];29             up[i][j]=up[i][j-1];30             lef[i][j]=lef[i-1][j];31             if (mp[i][j]==‘1‘)32             {33                 d[i][j]++;34                 if (i>1&&mp[i-1][j]==‘1‘) x[i][j]++;35                 if (j>1&&mp[i][j-1]==‘1‘) x[i][j]++;36                 if (i>1&&mp[i-1][j]==‘1‘) up[i][j]++;37                 if (j>1&&mp[i][j-1]==‘1‘) lef[i][j]++;38             }39         }40     }41 }42 43 int main()44 {45     scanf("%d%d%d",&n,&m,&q);46 47     for (int i=1;i<=n;i++)48         scanf("%s",mp[i]+1);49 50     slv();51 52     while (q--)53     {54         int x1,y1,x2,y2;55         scanf("%d%d%d%d",&x1,&y1,&x2,&y2);56         if (x1>x2) swap(x1,x2);57         if (y1>y2) swap(y1,y2);58         int ans = d[x2][y2]-d[x2][y1-1]-d[x1-1][y2]+d[x1-1][y1-1];59 60         ans -= x[x2][y2]-x[x2][y1-1]-x[x1-1][y2]+x[x1-1][y1-1];61         ans += up[x1][y2]-up[x1][y1-1];62         ans += lef[x2][y1]-lef[x1-1][y1];63 64         cout<<ans<<endl;65     }66     return 0;67 }

View Code