duty （图论题.二维前缀和）

10.7

图论题.二维前缀和的应用.
标算为:对于不包含环的图,连通块数目=点数-边数,所以利用二维前缀和进行预处理,O(1)求出矩形区域内的边数和点数.
很好写的70分算法:对每组询问都暴力求连通块数目的复杂度为O(NMQ),可以通过前7个测试点.出题人认为,在考场上,一个水平中等的选手最佳的策略是采用这个70分算法以留出时间思考第三题.
第3,4,5,6个测试点的其他做法:利用和标算同样的思路,但不需要二维前缀和,求解时较简单.
第8,9个测试点的目的是照顾常数过大的选手和数组开小的选手.
本题实际上是AtCoderGrandContest015的C题,可能会有一些选手因为做过原题而在本题快速得到高分,但估计不会很多.

#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 2010
using namespace std;

int n, m, Q, Esum[N][N], Psum[N][N];
int ln[N][N], rw[N][N];
char ss[N][N];

int main(){
    freopen( "duty.in", "r", stdin);
    freopen( "duty.out", "w", stdout);
    scanf( "%d%d%d" , &n , &m , &Q );
    for(int i=1; i<=n; i++){
        scanf("%s", ss[i] + 1);
        for(int j=1; j<=m; j++){
            ss[i][j] -= '0';
            Esum[i][j] = Esum[i-1][j] + Esum[i][j-1] - Esum[i-1][j-1];//边数 
            Psum[i][j] = Psum[i-1][j] + Psum[i][j-1] - Psum[i-1][j-1] + ss[i][j];//点数 
            ln[i][j] = ln[i][j-1];//一列穿过了多少边 
            rw[i][j] = rw[i-1][j];//一行穿过了多少边  
            if( ss[i][j] ){
                if( ss[i-1][j] ) Esum[i][j]++, ln[i][j]++;
                if( ss[i][j-1] ) Esum[i][j]++, rw[i][j]++;
            }
        }
    }
    for(int i=1, x1, y1, x2, y2; i<=Q; i++){
        scanf( "%d%d%d%d" , &x1 , &y1 , &x2 , &y2 );
        int Ptot = Psum[x2][y2] - Psum[x2][y1-1] - Psum[x1-1][y2] + Psum[x1-1][y1-1];
        int Etot = Esum[x2][y2] - Esum[x2][y1-1] - Esum[x1-1][y2] + Esum[x1-1][y1-1];
        int cc = ln[x1][y2] - ln[x1][y1-1] + rw[x2][y1] - rw[x1-1][y1];
        printf( "%d\n" , Ptot - Etot + cc );//cc是减重的边 
    }
}