BZOJ5366 [Lydsy1805月赛]代码派对

最新推荐文章于 2021-10-28 00:47:23 发布

原创最新推荐文章于 2021-10-28 00:47:23 发布 · 492 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#BZOJ #差分

蒟蒻zxl的Blog专栏同时被 3 个专栏收录

460 篇文章

订阅专栏

BZOJ

281 篇文章

订阅专栏

其他---技巧题/数学题/STL

34 篇文章

订阅专栏

该博客详细解析了BZOJ5366题目的解题思路，主要涉及差分技术。博主通过讨论如何唯一表示一个方案，引入了s1到s4四个矩阵的概念，并分析了每个坐标的矩形贡献。最后，博主提到使用瞎差分方法计算这四个矩阵，给出了扫面线算法的解决方案。

差分乱搞

题目传送门

%%%法老学长秒了这道题，我只会 $1000^2Tn^3$ 暴力。

考虑一个方案怎么被唯一表示。

三个矩形的交一定也是一个矩形。那么我们不妨设这个矩形的左下角为 $(x,y)$ ，然后记这样四个东西：

$s1[i][j]$ 表示包含但左下角不为 $(x,y)$ 的矩形个数。
$s2[i][j]$ 表示左下角为 $(x,y)$ 的矩形个数。
$s3[i][j]$ 表示 $(x,y)$ 在下边界但不是左下角的矩形个数。
$s4[i][j]$ 表示 $(x,y)$ 在左边界但不是左下角的矩形个数。

分类讨论每个坐标的贡献：

有左下角为 $(x,y)$ 的矩形： $C_{s2[x][y]}^3+C_{s2[x][y]}^2\times s1[x][y]+C_{s1[x][y]}^2\times s2[i][j]$
没有左下角为 $(x,y)$ 的矩形，此时 $(x,y)$ 一定在一个矩形的左边界和另一个矩形的下边界上： $C_{s3[i][j]}^2\times s4[i][j]+C_{s4[i][j]}^2\times s3[i][j]+s3[i][j]\times s4[i][j]\times s1[i][j]$

扫一遍过去就好了。

至于这四个东西怎么算，瞎差分就好啦。

代码：

#include<vector>
#include<cctype>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1005
#define M 100005
#define F inline
using namespace std;
struct node{ int p,x; };
int t,n,s1[N][N],s2[N][N],s3[N][N],s4[N][N];
int xa[M],xb[M],ya[M],yb[M],cf[N];
vector <node> a[N];
F char readc(){
    static char buf[100000],*l=buf,*r=buf;
    if (l==r) r=(l=buf)+fread(buf,1,100000,stdin);
    return l==r?EOF:*l++;
}
F int _read(){
    int x=0; char ch=readc();
    while (!isdigit(ch)) ch=readc();
    while (isdigit(ch)) x=(x<<3)+(x<<1)+(ch^48),ch=readc();
    return x;
}
#define C(x,y) (1ll*x*(x-1)/2*y+1ll*y*(y-1)/2*x)
#define P(x,y,w) a[x].push_back((node){y,w})
int main(){
    for (t=_read();t;t--){
        n=_read();
        for (int i=1;i<=1000;i++){
            for (int j=1;j<=1000;j++)
                s1[i][j]=s2[i][j]=s3[i][j]=s4[i][j]=0;
            a[i].clear(),cf[i]=0;
        }
        for (int i=1;i<=n;i++){
            xa[i]=_read(),ya[i]=_read(),xb[i]=_read(),yb[i]=_read();
            if (xa[i]>xb[i]) swap(xa[i],xb[i]);
            if (ya[i]>yb[i]) swap(ya[i],yb[i]);
            P(xa[i],ya[i],1),P(xa[i],yb[i]+1,-1);
            P(xb[i]+1,ya[i],-1),P(xb[i]+1,yb[i]+1,1);
        }
        for (int i=1;i<=1000;i++){
            for (int j=0,bj=a[i].size();j<bj;j++)
                cf[a[i][j].p]+=a[i][j].x;
            for (int j=1;j<=1000;j++)
                s1[i][j]=s1[i][j-1]+cf[j];
        }
        for (int i=1;i<=n;i++)
            s2[xb[i]][ya[i]]++,s1[xb[i]][ya[i]]--;
        for (int i=1;i<=n;i++){
            if (ya[i]<yb[i])
                s3[xb[i]][ya[i]+1]++,s3[xb[i]][yb[i]+1]--;
            if (xa[i]<xb[i])
                s4[xa[i]][ya[i]]++,s4[xb[i]][ya[i]]--;
        }
        for (int i=1;i<=1000;i++)
            for (int j=1;j<=1000;j++)
                s3[i][j]+=s3[i][j-1],s4[i][j]+=s4[i-1][j];
        long long ans=0;
        for (int i=1;i<=1000;i++)
            for (int j=1;j<=1000;j++){
                ans+=1ll*s2[i][j]*(s2[i][j]-1)*(s2[i][j]-2)/6;
                ans+=C(s2[i][j],s1[i][j])+C(s3[i][j],s4[i][j]);
                ans+=1ll*s3[i][j]*s4[i][j]*(s1[i][j]-s3[i][j]-s4[i][j]);
            }
        printf("%lld\n",ans);
    }
    return 0;
}