[CERC2017]Buffalo Barricades

2 篇文章 0 订阅
1 篇文章 0 订阅

思路

set+并查集+玄学乱搞……
一道对STL的灵活运用题……
这个我也是看了Claris大神的博客才写出来的(%%%)
好了废话不多说,我们进入正题。
这道题其实就是求一个点为右上角能围住多少点(废话)
按照官方题解的思路,首先要用一个玄学扫描线维护。
既然要扫描线,那么肯定首先排序啊。
按y轴从大到小排序,扫描从上往下。
好了现在我们抛开被围的点,只注意栅栏的顶点。
首先,一个栅栏的顶点向左和下发出射线。
向下发出的射线在x轴上就是它本身的位置,只要在set中保留这个点就相当于保存了这条射线。
那么我们现在只考虑向左发出的射线。
若有下面两个点:
A(1,3),D(3,2)
这里如果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;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值