Problem
Solution
不知道为什么网上的题解每一篇都要强调一句“此题很简单”,我觉得不是很简单啊qwqqq
按x坐标排序,然后进行分治,每次只统计跨立两个区间的矩形个数,为了避免重复统计,我们仅统计从左下角到右下角的矩形。枚举左边的点,显然我们只需要考虑y坐标比它大的点(当然,x坐标也比它大),那么我们可以这样按照y坐标从上到下进行归并排序。
由于矩形中不允许有其他的稻草人,那么考虑左边其他点对矩形的影响。假设有点p,点q,其中q在p的右下方,即满足
xp<xq,yp>yq
x
p
<
x
q
,
y
p
>
y
q
,那么我们会发现如果一个跨立的矩形如果包含p点,那么就一定会包含q点,那么加入q点后我们就可以不考虑p点的影响了,可以画个图理解一下。而且又因为我们是按照y轴从大到小加点的,所以这其实就相当于要维护一个x逐渐减小的单调栈。这个单调栈里的点就是对当前点有限制的点。
再来考虑右边的话,就会发现如果q在p的左下方,那么包含p的矩形一定会包含q,同样的我们认为加入q之后就不需要考虑p的影响了,那么右边是一个x逐渐增大的单调栈。这个单调栈里的点是可能与当前点匹配的点。
当我们去考虑当前点A能匹配的点B时,点B的y不能超过A的前一个点,而右边又是单调的,二分这个位置即可。
时间复杂度 O(nlog2n) O ( n log 2 n )
Code
#include <algorithm>
#include <cstdio>
#define rg register
using namespace std;
typedef long long ll;
const int maxn=200010;
template <typename Tp> inline void getmin(Tp &x,Tp y){if(y<x) x=y;}
template <typename Tp> inline void getmax(Tp &x,Tp y){if(y>x) x=y;}
template <typename Tp> inline void read(Tp &x)
{
x=0;int f=0;char ch=getchar();
while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
if(ch=='-') f=1,ch=getchar();
while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
if(f) x=-x;
}
struct point{
int x,y;
bool operator < (const point &t)const{return x<t.x;}
}p[maxn],tmp[maxn];
struct Stack{
int tp;point a[maxn];
void clear(){tp=0;}
point &top(){return a[tp];}
void push(point x){a[++tp]=x;}
void pop(){tp--;}
}ls,rs;
int n,tp;
ll ans;
int find(int k)
{
int res=rs.tp+1,l=1,r=rs.tp+1,m;
while(l<=r)
{
m=(l+r)>>1;
if(rs.a[m].y<k) res=m,r=m-1;
else l=m+1;
}
return res;
}
void cdq(int l,int r)
{
int temp;
if(l>=r) return ;
int m=(l+r)>>1,i=l,j=m+1;
cdq(l,m);cdq(m+1,r);
ls.clear();rs.clear();
for(;i<=m;i++)
{
while(ls.tp&&ls.top().x<p[i].x) ls.pop();
ls.push(p[i]);
while(j<=r&&p[j].y>=p[i].y)
{
while(rs.tp&&rs.top().x>p[j].x) rs.pop();
rs.push(p[j++]);
}
if(ls.tp==1) ans+=(temp=rs.tp);
else ans+=(temp=rs.tp-find(ls.a[ls.tp-1].y)+1);
}
i=l;j=m+1;tp=0;
while(i<=m&&j<=r)
{
if(p[i].y>p[j].y) tmp[++tp]=p[i++];
else tmp[++tp]=p[j++];
}
while(i<=m) tmp[++tp]=p[i++];
while(j<=r) tmp[++tp]=p[j++];
for(rg int i=1;i<=tp;i++) p[l+i-1]=tmp[i];
}
int main()
{
#ifndef ONLINE_JUDGE
freopen("in.txt","r",stdin);
#endif
read(n);
for(rg int i=1;i<=n;i++){read(p[i].x);read(p[i].y);}
sort(p+1,p+n+1);
cdq(1,n);
printf("%lld\n",ans);
return 0;
}