bzoj 2535 航空管制 —— 贪心+拓扑序

Posted zinn

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj 2535 航空管制 —— 贪心+拓扑序相关的知识,希望对你有一定的参考价值。

题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2535

这个题,如果正着考虑,也就是先考虑放在前面的再考虑放在后面的,决策时会有矛盾;
也就是,如果要求 pos[a] < pos[b],则先考虑放 a,因为许多点放在 a 后面,所以 a 尽量往前放可以给它们留出空位;
但又有限制最晚起飞时间,那么先考虑的 a 应该尽量放在靠近它最晚起飞时间的地方,以免后面的点无法满足起飞时间的限制;
于是尽量往前放和尽量往后放形成了矛盾,无法决策;
所以需要转化一个限制,不妨对于 pos[a] < pos[b],先决策 b;
于是两个限制合起来的要求就都变成尽量往后放了(真神奇);
所以我们就这样倒着做,一定能得到合法的序列;
然后考虑每个点的可能最前位置,其实就是做的过程中把这个点以及它限制的点都暂时去掉,让其他点放好,那么不能再放的时候就意味着只有放上这个点才能继续,这时这个点按刚才做法找到的位置就是它的最前位置(其他点已经尽量把它的可行位置挤到前面了);
注意放点位置是最晚起飞时间和限制它的点的起飞时间取 min 的!
注意答案是起飞序列!而不是每个点的起飞位置!
如果用 set 找某位置前面第一个空位会很慢,还是并查集好。
代码如下:

技术分享图片
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<set>
using namespace std;
int const xn=2005,xm=10005;
int n,m,hd[xn],ct,to[xm],nxt[xm],deg[xn],tmp[xn],ans[xn],d[xn],mn[xn];
struct N{
  int v,id;
  N(int v=0,int i=0):v(v),id(i) {}
};
queue<int>q;
set<int>st;
set<int>::iterator it;
int rd()
{
  int ret=0,f=1; char ch=getchar();
  while(ch<0||ch>9){if(ch==-)f=0; ch=getchar();}
  while(ch>=0&&ch<=9)ret=ret*10+ch-0,ch=getchar();
  return f?ret:-ret;
}
void add(int x,int y){to[++ct]=y; nxt[ct]=hd[x]; hd[x]=ct;}
void topo(int nw)
{
  for(int i=1;i<=n;i++)
    {
      mn[i]=min(n,d[i]); tmp[i]=deg[i];
      if(!tmp[i]&&i!=nw)q.push(i);
    }
  while(q.size())
    {
      int x=q.front(); q.pop();
      it=st.lower_bound(n-mn[x]+1);
      int tim=*it;
      //ans[x]=n-tim+1; //!!!
      ans[n-tim+1]=x;
      st.erase(tim);
      for(int i=hd[x],u;i;i=nxt[i])
    {
      if((u=to[i])==nw)continue;
      tmp[u]--; mn[u]=min(mn[u],mn[x]);//!!
      if(!tmp[u])q.push(u);
    }
    }
}
int solve(int x)
{
  st.clear();
  for(int i=1;i<=n;i++)st.insert(i);
  topo(x);
  it=st.lower_bound(n-d[x]+1);
  return n-*it+1;
}
int main()
{
  n=rd(); m=rd();
  for(int i=1;i<=n;i++)d[i]=rd(),st.insert(i);
  for(int i=1,x,y;i<=m;i++)x=rd(),y=rd(),add(y,x),deg[x]++;
  topo(0);
  for(int i=1;i<=n;i++)printf("%d ",ans[i]); puts("");
  for(int i=1;i<=n;i++)printf("%d ",solve(i)); puts("");
  return 0;
}
set
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
int const xn=2005,xm=10005;
int n,m,hd[xn],ct,to[xm],nxt[xm],deg[xn],tmp[xn],ans[xn],d[xn],mn[xn];
int fa[xn],q[xn];
int rd()
{
  int ret=0,f=1; char ch=getchar();
  while(ch<0||ch>9){if(ch==-)f=0; ch=getchar();}
  while(ch>=0&&ch<=9)ret=ret*10+ch-0,ch=getchar();
  return f?ret:-ret;
}
void add(int x,int y){to[++ct]=y; nxt[ct]=hd[x]; hd[x]=ct;}
int find(int x){return fa[x]==x?x:fa[x]=find(fa[x]);}
void topo(int nw)
{
  int cnt=0;
  for(int i=1;i<=n;i++)
    {
      mn[i]=min(n,d[i]); tmp[i]=deg[i]; fa[i]=i;
      if(!tmp[i]&&i!=nw)q[++cnt]=i;
    }
  while(cnt)
    {
      int x=q[cnt--]; int tim=find(mn[x]);
      ans[tim]=x; fa[tim]=tim-1;
      for(int i=hd[x],u;i;i=nxt[i])
    {
      if((u=to[i])==nw)continue;
      tmp[u]--; mn[u]=min(mn[u],mn[x]);//!!
      if(!tmp[u])q[++cnt]=u;
    }
    }
}
int solve(int x)
{
  topo(x); return find(d[x]);
}
int main()
{
  n=rd(); m=rd();
  for(int i=1;i<=n;i++)d[i]=rd();
  for(int i=1,x,y;i<=m;i++)x=rd(),y=rd(),add(y,x),deg[x]++;
  topo(0);
  for(int i=1;i<=n;i++)printf("%d ",ans[i]); puts("");
  for(int i=1;i<=n;i++)printf("%d ",solve(i)); puts("");
  return 0;
}

 












以上是关于bzoj 2535 航空管制 —— 贪心+拓扑序的主要内容,如果未能解决你的问题,请参考以下文章

BZOJ 2109 航空管制(拓扑排序+贪心)

Bzoj2535 [Noi2010]Plane 航空管制2

BZOJ 2535: [Noi2010]Plane 航空管制2

BZOJ 2535 Plane 航空管制2

BZOJ2109 NOI2010 Plane 航空管制 拓扑排序

bzoj 2109: [Noi2010]Plane 航空管制