题目描述 Description###
世博期间，上海的航空客运量大大超过了平时，随之而来的航空管制也频频 发生。最近，小X就因为航空管制，连续两次在机场被延误超过了两小时。对此， 小X表示很不满意。 在这次来烟台的路上，小 X不幸又一次碰上了航空管制。于是小 X开始思考 关于航空管制的问题。 假设目前被延误航班共有 n个，编号为 1至n。机场只有一条起飞跑道，所 有的航班需按某个顺序依次起飞（称这个顺序为起飞序列）。定义一个航班的起 飞序号为该航班在起飞序列中的位置，即是第几个起飞的航班。 起飞序列还存在两类限制条件：
第一类（最晚起飞时间限制）：编号为 i的航班起飞序号不得超过 (k_i) ;
第二类（相对起飞顺序限制）：存在一些相对起飞顺序限制(a, b) ，表示 航班 a的起飞时间必须早于航班 b，即航班 a的起飞序号必须小于航班 b 的起飞序号。
小X 思考的第一个问题是，若给定以上两类限制条件，是否可以计算出一个 可行的起飞序列。第二个问题则是，在考虑两类限制条件的情况下，如何求出每 个航班在所有可行的起飞序列中的最小起飞序号。
输入描述 Input Description###
第一行包含两个正整数 n和m，n表示航班数目，m表示 第二类限制条件（相对起飞顺序限制）的数目。 第二行包含 n个正整数 (k_1, k_2, „, k_n) 。 接下来 m行，每行两个正整数 a和b，表示一对相对起飞顺序限制(a, b)， 其中$1≤a,b≤n $ , 表示航班 a必须先于航班 b起飞。
输出描述 Output Description###
包含 n个整数 (t_1, t_2, „, t_n) ，其中 ti表示航班i可能的最小起飞序 号，相邻两个整数用空格分隔。
样例输入 Sample Input###
5 5
4 5 2 5 4
1 2
3 2
5 1
3 4
3 1
样例输出 Sample Output###
3 4 1 2 1
数据范围及提示 Data Size & Hint###
(n≤2,000，m≤10,000)
之前的一些废话###
为什么我还做NOI的题，我没救了
不过反正这个博客没人看，管他呢
题解###
第二种限制显然是一个拓扑关系，第一种限制告诉你最晚出发时间
正难则反
考虑反向建图，第二种限制除了边反了没有什么区别，第二种限制变成了最早出发时间。
不考虑任何点的情况下，倒着枚举时间，每次选择k最大的，度数为0的点
考虑点i的情况下，先忽视点i，进行一次拓扑排序(注意每次入队的点不光要保证入度为0，也必须保证在该点最早时间以后)
当某个时刻没有点可以选择时，这时必须选择点i，这就是点i出现的最早时间
用堆来维护即可。
复杂度：(O(n^2))
代码###
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<queue>
using namespace std;
#define mem(a,b) memset(a,b,sizeof(a))
typedef long long LL;
typedef pair<int,int> PII;
#define X first
#define Y second
inline int read()
{
int x=0,f=1;char c=getchar();
while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}
while(isdigit(c)){x=x*10+c-'0';c=getchar();}
return x*f;
}
const int maxn=2010,maxm=10010;
struct Edge
{
int u,v,next;
Edge() {}
Edge(int _1,int _2,int _3):u(_1),v(_2),next(_3) {}
}e[maxm];
int n,m,ce=-1,first[maxn],a[maxn],x,y,s[maxn],top,r[maxn],R[maxn];
priority_queue <PII,vector<PII>,greater<PII> > Q;
void addEdge(int a,int b){e[++ce]=Edge(a,b,first[a]);first[a]=ce;r[b]++;}
int calc(int cur)
{
int now=1;
mem(s,0);top=0;while(Q.size())Q.pop();
for(int i=1;i<=n;i++)R[i]=r[i];
for(int i=1;i<=n;i++)if(R[i]==0 && i!=cur)Q.push(make_pair(a[i],i));
while(Q.size() && Q.top().X<=now)
{
s[++top]=Q.top().Y;
Q.pop();
now++;
}
while(top)
{
int x=s[top--];
for(int i=first[x];i!=-1;i=e[i].next)
{
R[e[i].v]--;
if(!R[e[i].v] && e[i].v!=cur)Q.push(make_pair(a[e[i].v],e[i].v));
}
while(Q.size() && Q.top().X<=now)
{
s[++top]=Q.top().Y;
Q.pop();
now++;
}
}
return now;
}
int main()
{
mem(first,-1);
n=read();m=read();
for(int i=1;i<=n;i++)a[i]=n-read()+1;
while(m--)x=read(),y=read(),addEdge(y,x);
for(int i=1;i<=n;i++)
{
if(i>1)putchar(' ');
printf("%d",n-calc(i)+1);
}
putchar('
');
return 0;
}
总结###
正难则反