[Noi2010]Plane 航空管制

这道题做法很简单,年鉴上讲的很复杂的样子……

首先第一问,对于每个节点来说,他肯定要比自己的后继先,在此基础上越往后越好(尽量满足限制严格的其他点)

那么可以将一个节点i的k更新为min(k[i],k[j]-1)j为i的后继

注意更新的时候要按照拓扑序来,否则有可能更新不完全//他儿子还没更新完就更新他了

//我一开始没有注意拓扑序,后面懒得改了仿照Bellman-Ford在外面再套了一层控制次数,于是超慢……

最后按照这个k排序输出即可,可以证明这样一定满足条件


重点是第二问

对于每个点i,首先DFS一下,那么他肯定要放在他这棵子树的所有点的后面  //实际上不是棵树,但还是认为是一棵树吧

记当前答案为ans ,已经放了cnt个点了,那么初始cnt=ans,然后每访问一个没放到序列中的点,cnt++

然后对于其他的点,按照第一问求得的k从小到大访问

有可能存在点的k>ans,意味着这个点一定要放在ans前面,那么ans++

还有可能存在点的k<=cnt,也就是说点i无论如何一定要放在这个点后面,否则怎样都无法满足,那么ans=k+1


然后搞定……


//Lib
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<ctime>

#include<iostream>
#include<algorithm>
#include<vector>
#include<string>
#include<queue>
using namespace std;
//Macro
#define rep(i,a,b) for(int i=a,tt=b;i<=tt;++i)
#define drep(i,a,b) for(int i=a,tt=b;i>=tt;--i)
#define erep(i,e,x) for(int i=x;i;i=e[i].next)
#define irep(i,x) for(__typedef(x.begin()) i=x.begin();i!=x.end();i++)
#define read() (strtol(ipos,&ipos,10))
#define sqr(x) ((x)*(x))
#define pb push_back
#define PS system("pause");
typedef long long ll;
typedef pair<int,int> pii;
const int oo=~0U>>1;
const double inf=1e100;
const double eps=1e-6;
string name="plane",in=".in",out=".out";
//Var
struct PLANE
{
	int idx,k;
	bool operator <(const PLANE &o)const{return k<o.k;}
}p[2008];
struct E
{
	int next,node;	
}e[20008];
int n,m,tot,cnt,pos[2008],h[2008],A[2008];
bool vis[2008];
void add(int a,int b){e[++tot].next=h[a];e[tot].node=b;h[a]=tot;}
void Init()
{
	scanf("%d%d",&n,&m);int a,b;
	rep(i,1,n)scanf("%d",&p[i].k),p[i].idx=i;
	rep(i,1,m)scanf("%d%d",&a,&b),add(a,b),add(b,a);
}
bool Deal(int u)
{
	erep(i,e,h[u])
	{
		if(i&1)
		{
			int v=e[i].node;
			p[u].k=min(p[u].k,p[v].k-1);
		}
	}
}
void Solve1()
{
	rep(j,1,n)
		rep(i,1,n)
			Deal(i);
	sort(p+1,p+1+n);
	rep(i,1,n-1)printf("%d ",p[i].idx);
	printf("%d\n",p[n].idx);
}
int Count(int u)
{
	int ret=0,v;vis[u]=true;
	erep(i,e,h[u])
		if(!(i&1)&&!vis[v=e[i].node])
			ret+=Count(v);
	return ret+1;
}
void Solve2()
{
	rep(i,1,n)
	{
		memset(vis,0,sizeof vis);
		int ans=Count(i);
		cnt=ans;
		rep(j,1,n)
		{
			if(!vis[p[j].idx])
			{
				cnt++;
				if(p[j].k<=ans)ans++;
				else if(cnt>p[j].k)ans=p[j].k+1;
			}
		}
		A[i]=ans;
	}
	rep(i,1,n-1)printf("%d ",A[i]);
	printf("%d\n",A[n]);
}
void Work()
{
	Solve1();
	Solve2();
}
int main()
{
//	freopen((name+in).c_str(),"r",stdin);
//	freopen((name+out).c_str(),"w",stdout);
	Init();
	Work();
//	PS;
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值