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个 ， 我并不确定哪些是真正的变化时刻 ， 但不要紧 ， 都记录一下不会错。

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