BZOJ1182 Croatian2009 PLAHTE

提示:
1. 可以尝试算算每个矩形覆盖面积的二次差分

详细题解在代码后:

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <string>
#include <vector>
#include <deque>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <algorithm>

using namespace std;
const int maxt = 1e6+1e2;
typedef long long ll;

int re() {
    int n = 0, ch = getchar(); bool flag = false;
    while(!isdigit(ch)) flag |= ch == '-', ch = getchar();
    while(isdigit(ch)) n = n * 10 + ch - '0', ch = getchar();
    return flag ? -n : n;
}

int x1 , y1 , x2 , y2;

ll square(int t)
{
    ll x = min(t , x2) - max(-t , x1)+1;
    ll y = min(t , y2) - max(-t , y1)+1;
    return (x<=0 || y<=0) ? 0 : x*y;
}

int a[maxt] , a2[maxt] , a4[maxt] , b[maxt];

void modify(int t1 , int t2)
{
    ll p0 = square(t1-1);
    ll p1 = square(t1  );
    ll p2 = square(t1+1);

    a[t1] += p1-p0;
    a[t2] -= p1-p0;

    if(t2>t1+1)
    {
        ll v = p2-2*p1+p0;
        a[t2] -= v*(t2-t1);
        if(v == 2) a2[t1]++ , --a2[t2];
        if(v == 4) a4[t1]++ , --a4[t2];
    }
}


int main(int argc, char *argv[]) {

    int n;
    int q[20];
    cin>>n;

    for(int i=1;i<=n;i++)
    {
        x1 = re(); y1 = re(); x2 = re(); y2 = re();

        for(int d=-1 , cnt=0;d<=1;d++)
        {
            q[cnt++] = abs(x1+d);
            q[cnt++] = abs(x2+d);
            q[cnt++] = abs(y1+d);
            q[cnt++] = abs(y2+d);
        }

        sort(q, q+12);
        for(int i=1;i<12;i++) modify(q[i-1], q[i]);
    }

    cin>>n;
    int Mx = 0;
    for(int i=1;i<=n;i++) 
    {
        x1 = re();
        b[x1] = 1;
        Mx = max(Mx , x1);
    }

    ll res = 0 , add=0 , k2=0 , k4=0;
    for(int i=0;i<=Mx;i++)
    {
        res+= add+=a[i];
        add+= 2*(k2+=a2[i])+4*(k4+=a4[i]);
        if(b[i]) printf("%lld\n" , res);
    }


    return 0;
}

首先每个矩形的关键时刻最多12个 , 我并不确定哪些是真正的变化时刻 , 但不要紧 , 都记录一下不会错。

于是我把矩形分成若干时刻 , 在这些时刻之间矩形的面积增长是一定的。 于是查分记录一下就好。

你可能感兴趣的:(差分,Croatian)