E. Plan of Lectures（给定拓扑图求输出序列满足存在给定链）

#include<bits/stdc++.h>
using namespace std;
#define pb push_back
const int M=3e5+5;
int p[M],vis[M],nex[M],in[M],rt[M],rnk[M],cmp[M];
vector<int>G[M],ans;
queue<int>que;
int main(){
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&p[i]);
    for(int u,v,i=1;i<=m;i++){
        scanf("%d%d",&u,&v);
        nex[u]=v;
        vis[v]=1;///在不是链头就标记
    }
    int tot=0;
    for(int i=1;i<=n;i++){
        if(vis[i]) continue;
        cmp[i]=++tot;
        int now=0,x=nex[i];
        rnk[i]=now;
        rt[tot]=i;
        while(x!=0&&cmp[x]==0){
            rnk[x]=++now;
            cmp[x]=tot;
            x=nex[x];
        }
    }
    for(int i=1;i<=n;i++)///check circle
        if(!cmp[i]||(cmp[p[i]]==cmp[i]&&rnk[i]<rnk[p[i]])) return puts("0"),0;

    for(int i=1;i<=n;i++){///construct new graph
        if(p[i]==0) continue;
        if(cmp[i]!=cmp[p[i]]){
            G[cmp[p[i]]].pb(cmp[i]);
            in[cmp[i]]++;
        }
    }

    ///tuopu
    for(int i=1;i<=tot;i++)
        if(in[i]==0)
            que.push(i);
    while(!que.empty()){
        int u=que.front();que.pop();
        ans.pb(u);
        for(auto v:G[u])
            if(--in[v]==0) que.push(v);
    }
    if(ans.size()==tot){
        for(auto v:ans){
            int nowi=rt[v];
            while(nowi!=0){
                printf("%d ",nowi);
                nowi=nex[nowi];
            }
        }
    }
    else puts("0");
    return 0;
}