[联合集训6-21] 最小拓扑序 贪心

首先可以发现对于每个点$v$最多添加一条边指向它，假设添加了边$x\to v,y\to v(x>y)$ 换成$x\to y,y\to v$不会更劣。
我们贪心选点加入拓扑序。那么考虑当前局面一个入度为$0$的点的集合，我们会在当前可添加边数允许的情况下，尽可能选择一个编号大的点放入拓扑序，假设我们选择点$i$，那么对于集合中 <i < i <script type="math/tex" id="MathJax-Element-4236"> $j$，要为其安排一条边指向它，这样才一定能选到$i$。但直接安排 $i$连向这些点不一定是最优的，我们考虑到当前局面的所有点已经被强制限制在了$i$之后进入拓扑序，那么安排 $i$之后的点连向这些点也是合法的，所以我们不妨暂且设是$i$连向了这些点，后面遇到其它合法方案再覆盖掉即可。

具体地，我们设$f_i$为准备安排了$f_i$连向$i$（$f_i=0$表示没有安排），我们维护入度为$0$的点的集合，每次就要进行以下操作：
1. 假设当前还剩下$K$条边可以安排，我们找到集合中最大的$i$，满足所有在集合中有$j<i,f_j=0$的$j$的个数不超过$K$。
2. 连边$f_i\to i$（$f_i=0$则不连）。
3. 把所有的$j<i$的$f_j$改成$i$，并得到新的$K$。
4. 删除点$i$，并把新的入度为$0$的点加入集合。

随便拿个数据结构维护一下即可。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<algorithm>
#include<queue>
#define N 100010
#define fs first
#define sc second
#define mid (l+r>>1)
using namespace std;
int n,m,K,top,rd[N],tote,con[N],to[N<<1],nxt[N<<1],ord[N];
pair<int,int> ans[N];
priority_queue<int,vector<int>,greater<int> > Q;
int read()
{
    int x=0,f=1;char ch=getchar();
    for(;ch<'0'||ch>'9';ch=getchar()) if(ch=='-') f=-1;
    for(;ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0';
    return x*f;
}
void ins(int x,int y)
{
    to[++tote]=y;
    nxt[tote]=con[x];
    con[x]=tote;
}
struct tree
{
    int cvr[N<<2],num[N<<2];
    bool ex[N<<2];
    void update(int v)
    {   
        num[v]=num[v<<1]+num[v<<1|1];
        ex[v]=ex[v<<1]|ex[v<<1|1];
    }
    void pushdown(int v)
    {
        if(cvr[v])
        {
            if(ex[v<<1]) cvr[v<<1]=cvr[v],num[v<<1]=0;
            if(ex[v<<1|1]) cvr[v<<1|1]=cvr[v],num[v<<1|1]=0;
            cvr[v]=0;
        }
    }
    void build(int v,int l,int r)
    {
        if(l==r){num[v]=0;cvr[v]=0;return ;}
        build(v<<1,l,mid);
        build(v<<1|1,mid+1,r);
    }
    void pnew(int v,int l,int r,int x)
    {
        if(l==r){num[v]=1;ex[v]=1;return ;}
        pushdown(v);
        if(x<=mid) pnew(v<<1,l,mid,x);
        else pnew(v<<1|1,mid+1,r,x);
        update(v);
    }
    void mdf(int v,int l,int r,int lx,int rx,int d)
    {
        if(!ex[v]) return ;
        if(l==lx&&r==rx)
        {   
            K-=num[v];cvr[v]=d;num[v]=0;
            return ;
        }
        pushdown(v);
        if(rx<=mid) mdf(v<<1,l,mid,lx,rx,d);
        else if(lx>mid) mdf(v<<1|1,mid+1,r,lx,rx,d);
        else mdf(v<<1,l,mid,lx,mid,d),mdf(v<<1|1,mid+1,r,mid+1,rx,d);
        update(v);
    }
    int del(int v,int l,int r,int k)
    {

        if(l==r)
        {
            if(cvr[v])
            {
                ans[++top]=make_pair(cvr[v],l);
                cvr[v]=0;
            }
            ex[v]=0;num[v]=0;
            return l;
        }
        int re;
        pushdown(v);
        if(!ex[v<<1|1]||k<num[v<<1]) re=del(v<<1,l,mid,k);
        else re=del(v<<1|1,mid+1,r,k-num[v<<1]);
        update(v);
        return re;
    }
    void travel(int v,int l,int r)
    {
        if(!ex[v]) return ;
        if(l==r) 
        {
            if(cvr[v]) ans[++top]=make_pair(cvr[v],l);
        }
        pushdown(v);
        travel(v<<1,l,mid);
        travel(v<<1|1,mid+1,r);
    }
}T;
void run()
{
    memset(rd,0,sizeof(rd));
    for(int i=1;i<=tote;i++)
        rd[to[i]]++;
    for(int i=1;i<=n;i++)
        if(!rd[i]) Q.push(i);
    int tp=0;   
    while(!Q.empty())
    {
        int v=Q.top();Q.pop();
        ord[++tp]=v;
        for(int p=con[v];p;p=nxt[p])
        {
            rd[to[p]]--;
            if(!rd[to[p]]) Q.push(to[p]);   
        }
    }   
}
int main()
{
    freopen("topo.in","r",stdin);
    freopen("topo.out","w",stdout);
    n=read();m=read();K=read();
    for(int i=1;i<=m;i++)
    {
        int x=read(),y=read();
        ins(x,y);rd[y]++;
    }
    T.build(1,1,n);

    for(int i=1;i<=n;i++)
        if(!rd[i]) T.pnew(1,1,n,i);
    while(T.ex[1])
    {
        int v=T.del(1,1,n,K);
        T.mdf(1,1,n,1,v,v);
        for(int p=con[v];p;p=nxt[p])
            {
                rd[to[p]]--;
                if(!rd[to[p]]) T.pnew(1,1,n,to[p]);
            }
    }   
    for(int i=1;i<=top;i++)
        ins(ans[i].fs,ans[i].sc);
    run();
    for(int i=1;i<=n;i++)
        printf("%d ",ord[i]);
    printf("\n%d\n",top);
    for(int i=1;i<=top;i++)
        printf("%d %d\n",ans[i].fs,ans[i].sc);
    return 0;
}