思路
set+并查集+玄学乱搞……
一道对STL的灵活运用题……
这个我也是看了Claris大神的博客才写出来的(%%%)
好了废话不多说,我们进入正题。
这道题其实就是求一个点为右上角能围住多少点(废话)
按照官方题解的思路,首先要用一个玄学扫描线维护。
既然要扫描线,那么肯定首先排序啊。
按y轴从大到小排序,扫描从上往下。
好了现在我们抛开被围的点,只注意栅栏的顶点。
首先,一个栅栏的顶点向左和下发出射线。
向下发出的射线在x轴上就是它本身的位置,只要在set中保留这个点就相当于保存了这条射线。
那么我们现在只考虑向左发出的射线。
若有下面两个点:
这里如果A加入的时间比D早,那么最终形成的就是这样:
如果A比D晚,则是这样:
可以看出,在上面的第二种情况下,A点要在D加入时从set中删去。(见代码中的“注释1”)
好了,这一种情况我们就解决了。
但是样例中3号点和4号点的情况呢?
可以看出,样例中是3号点包含4号点。
统计时,3号点的答案要加上4号点的答案。
但是,如果是1和3的情况,虽然3也包含1,但是并不需要把3号点加上1号点的答案。(“注释5”)
那么我们可以建立一个并查集。
在每次加入一个点的时候,记录一个在它右侧且比它的时间大的点;(见代码中的“注释2”)
最后倒序枚举每一个点,用并查集维护,这个点的答案就是这个点所属并查集的数量之和。(“注释3”)
然后将这个点与插入时在他右侧的点,加入一个并查集中。
最后一个问题:上面斜体字有啥用?
“倒序枚举”,是由于在上面统计时我们已经把所有的点加入到整个图中。
而现在需要分开统计,倒序删去每一条边,就相当于一个拓扑排序。(“注释4”)
代码
//STL有啥不好啦,代码量短,如果手写splay还要各种操作,STL多简洁明了啊
#include <cstdio>
#include <set>
#include <utility>
#include <algorithm>
const int maxn=300000;
struct point
{
int x,y,t;
};
typedef std::set<std::pair<int,int> > set_pair;
typedef std::set<std::pair<int,int> >::iterator set_pair_iterator;
set_pair s;
int cnt[maxn+10],n,m,par[maxn+10],ans[maxn+10];
point p[(maxn<<1)+10];
struct union_find_set
{
int fa[maxn+10];
int find(int x)
{
return fa[x]?fa[x]=find(fa[x]):x;
}
inline int merge(int x,int y)
{
x=find(x);
y=find(y);
if(x!=y)
{
fa[x]=y;
cnt[y]+=cnt[x];
}
return 0;
}
};
union_find_set f;
bool cmp(const point &a,const point &b)
{
if(a.y==b.y)
{
return a.x<b.x;
}
return a.y>b.y;
}
int main()
{
scanf("%d",&n);
for(register int i=1; i<=n; ++i)
{
scanf("%d%d",&p[i].x,&p[i].y);
p[i].x=(p[i].x<<1)-1;
p[i].y=(p[i].y<<1)-1;
p[i].t=0;
}
scanf("%d",&m);
for(register int i=1; i<=m; ++i)
{
scanf("%d%d",&p[n+i].x,&p[n+i].y);
p[n+i].x<<=1;
p[n+i].y<<=1;
p[n+i].t=i;
}
std::sort(p+1,p+n+m+1,cmp);
for(register int i=1; i<=n+m; ++i)
{
std::pair<int,int> q=std::make_pair(p[i].x,p[i].t);
if(p[i].t)
{
s.insert(q);
set_pair_iterator now=s.find(q);
if((++now)!=s.end())
{
par[p[i].t]=now->second;//注释2:记录在它右侧的点
}
while(1)//注释1:这个while循环就是在进行这个步骤
{
now=s.find(q);
if(now==s.begin())
{
break;
}
if((--now)->second<p[i].t)
{
break;
}
s.erase(now);
}
}
else
{
set_pair_iterator it=s.lower_bound(q);
if(it!=s.end())
{
++cnt[it->second];
}
}
}
for(register int i=m; i; --i)//注释3:这个for循环就是在进行这样的步骤
//注释4:注意倒序枚举
{
ans[i]=cnt[f.find(i)];
if(par[i])
{
f.merge(i,par[i]);
//注释5:这里并没有特殊处理出上面的情况
//因为如果par的时间大于i的时间,那么加入同一并查集了并不影响结果
}
}
for(register int i=1; i<=m; ++i)
{
printf("%d\n",ans[i]);
}
return 0;
}