题意:
有n个矩形,第i个矩形的左下角是(i-1,l),右上角是(i,r),现在要你画一条平行于x轴的水平线,对于穿过这条线的矩形,保留它关于这条线轴对称的部分,就像这样:
保留黑色部分。现在问你保留部分面积最大是多少。
题解:
首先我们可以知道,这条线一定是某个矩形的中线,为什么?
假设现在有一些矩形就像这样:
这条线如果往上移x,那么1号矩形可以多出2x的面积,2号矩形会少2x的面积。但是当移到与3号矩形相交的时候,
这时候还没到第一个矩形的中线部分,那么向上移x,是否1号与3号都会加2x,2号-2x,那么所得的就是2x。这种情况直到移到1号中线为止。
所以我们大胆猜测这条线一定是某个矩形的中线(并且这是对的)
那么这时候问题来了,当一条线往上走的时候,越过中线前后的情况是不一样的,怎么办?
我们可以将一个矩形分成下与上两部分,分开做。先做下面部分:
我们按照下端点排序,依次做中线
首先中线在这,我们就知道它穿过3个矩形,那么由与下端点是递增的,我们可以用一个指针维护现在到哪个位置了。
那么我们怎么算答案?
我们可以用一个sum维护当前穿过的所有矩形的下端点的值,再用一个num维护当前穿过的矩形的数量。这样我们就可以得到当前的答案了:
当前中线高度*矩形数量-下端点的和
还需要两个数组维护每条中线所拥有的数量以及下端点的和,这样我们枚举中线的时候可以很快的知道到下个中线的时候失去了多少矩形。
比如到这里的时候我们就减去上个中线位置所拥有的矩形数量以及所有下端点的和。如此可以维护我们所需要的值。
对于右半部分处理的方法类似,不过sum维护的值有所不同。
接下来的问题就是,由与我们枚举的中线在左和右会算两次,那么我们最终需要减掉一次。
最后一个问题就是矩形的高度可能是个奇数的值,会在0.5的时候划分。我不知道double会不会有精度的问题,所以选择将所有点的高度*2,最后再/2即可
#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N=3e5+5;
struct node
{
ll l,r,id;
}e[N],line[N];
bool cmp1(node a,node b)
{
return a.l<b.l;
}
bool cmp2(node a,node b)
{
return a.r>b.r;
}
ll ans[N],del[N],a[N];
ll del_num[N],del_sum[N];
int main()
{
int n;
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
scanf("%lld%lld",&line[i].l,&line[i].r);
line[i].l*=2,line[i].r*=2;
line[i].id=(line[i].l+line[i].r)/2ll;
a[i]=line[i].id;
}
sort(a+1,a+1+n);
int all=unique(a+1,a+1+n)-a-1;
for(int i=1;i<=n;i++)
{
line[i].id=lower_bound(a+1,a+1+all,line[i].id)-a;
del_num[line[i].id]++;
del_sum[line[i].id]+=line[i].l;
if((line[i].r-line[i].l)%2==0)
del[line[i].id]+=line[i].r-line[i].l;
}
for(int i=1;i<=n;i++)
e[i].l=line[i].l,e[i].r=a[line[i].id];
sort(e+1,e+1+n,cmp1);
ll sum=0,num=0,p=1;
for(int i=1;i<=all;i++)
{
while(e[p].l<=a[i]&&p<=n)
sum+=e[p].l,num++,p++;
ans[i]+=a[i]*num-sum;
num-=del_num[i],sum-=del_sum[i];
}
sum=0,num=0,p=1;
for(int i=1;i<=all;i++)
del_sum[i]=0;
for(int i=1;i<=n;i++)
{
e[i].l=a[line[i].id],e[i].r=line[i].r;
del_sum[line[i].id]+=e[i].r;
}
sort(e+1,e+1+n,cmp2);
for(int i=all;i>=1;i--)
{
while(e[p].r>=a[i]&&p<=n)
sum+=e[p].r,num++,p++;
ans[i]+=sum-a[i]*num;
num-=del_num[i],sum-=del_sum[i];
}
ll mx=0;
for(int i=1;i<=all;i++)
mx=max(mx,(ans[i]*2-del[i])/2);
printf("%lld\n",mx);
return 0;
}