Vijos 1790 拓扑编号（拓扑排序+堆优化）

H国有n个城市，城市与城市之间有m条单向道路，满足任何城市不能通过某条路径回到自己。
现在国王想给城市重新编号，令第i个城市的新的编号为a[i]，满足所有城市的新的编号都互不相同，并且编号为[1,n]之间的整数。国王认为一个编号方案是优美的当且仅当对于任意的两个城市i,j，如果i能够到达j，那么a[i]应当小于a[j]。
优美的编号方案有很多种，国王希望使1号城市的编号尽可能小，在此前提下，使得2号城市的编号尽可能小…依此类推。
格式
输入格式

第一行读入n,m，表示n个城市，m条有向路径。
接下来读入m行，每行两个整数:x,y
表示第x个城市到第y个城市有一条有向路径。
输出格式

输出一行:n个整数
第i个整数表示第i个城市的新编号a[i]，输出应保证是一个关于1到n的排列。
样例1
样例输入1
5 4
4 1
1 3
5 3
2 5
样例输出
2 3 5 1 4
30%的测试点满足：n <= 10, m <= 10
70%的测试点满足：n <= 1000, m <= 10000
100%的测试点满足：n <= 100000, m <= 200000
输入数据可能有重边，可能不连通，但保证是有向无环图。

思路：因为题目说，对于任意x通道y，编号x要小于y，如要正向做，发现无法确定前面的点到底是第几小，所以舍弃这种方法，既然想要由大到小不行，那么就由小到大，所以反向建边，反向拓扑，用大根堆维护，输出即可。

program df;
type point=^node;
node=record
date,ends:longint;
next:point;
end;
var i,j,n,m,x,y,z,k,t,len:longint;
path:array[0..200000] of point;
a,dis,d,f:array[0..200000] of longint;
procedure com(x,y:longint);
var i:point;
begin
i:=path[x];
new(path[x]);
path[x]^.ends:=y;
path[x]^.next:=i;
end;

procedure put(x:longint);
var i,j:longint;
begin
inc(len);
a[len]:=x;
i:=len;
while (i>1) and (a[i]>a[i div 2]) do
begin
j:=a[i]; a[i]:=a[i div 2]; a[i div 2]:=j;
i:=i div 2;
end;
end;

function get:longint;
var i,j,dd:longint;
begin
get:=a[1];
a[1]:=a[len];
dec(len);
i:=1;
while i*2<=len do
begin
if (i*2+1>len) or (a[i*2]>a[i*2+1]) then j:=i*2
else j:=i*2+1;
if a[j]>a[i] then
begin
dd:=a[i]; a[i]:=a[j]; a[j]:=dd;
i:=j;
end
else break;
end;
end;

procedure dfs;
var y,u:longint;
i:point;
begin
while len>0 do
begin
u:=get;
inc(t);
dis[t]:=u; //每次取出当前最大的
i:=path[u];
while i<>nil do
begin
y:=i^.ends;
dec(d[y]);
if d[y]=0 then begin put(y); d[y]:=maxlongint div 3; end;
i:=i^.next;
end;
end;
end;

begin

readln(n,m);
for i:=1 to m do
begin
readln(x,y);
com(y,x);
inc(d[x]);
end;
len:=0; t:=0;
for i:=1 to n do
if d[i]=0 then begin put(i); d[i]:=maxlongint div 3; end;
dfs;
for i:=1 to t do
f[dis[i]]:=n-i+1;
for i:=1 to t do
write(f[i],’ ‘);
end.