首先可以发现对于每个点
v
v
最多添加一条边指向它,假设添加了边 换成
x→y,y→v
x
→
y
,
y
→
v
不会更劣。
我们贪心选点加入拓扑序。那么考虑当前局面一个入度为
0
0
的点的集合,我们会在当前可添加边数允许的情况下,尽可能选择一个编号大的点放入拓扑序,假设我们选择点,那么对于集合中
<i
<
i
<script type="math/tex" id="MathJax-Element-4236">
j
j
,要为其安排一条边指向它,这样才一定能选到。但直接安排
i
i
连向这些点不一定是最优的,我们考虑到当前局面的所有点已经被强制限制在了之后进入拓扑序,那么安排
i
i
之后的点连向这些点也是合法的,所以我们不妨暂且设是连向了这些点,后面遇到其它合法方案再覆盖掉即可。
具体地,我们设
fi
f
i
为准备安排了
fi
f
i
连向
i
i
(表示没有安排),我们维护入度为
0
0
的点的集合,每次就要进行以下操作:
1. 假设当前还剩下条边可以安排,我们找到集合中最大的
i
i
,满足所有在集合中有的
j
j
的个数不超过。
2. 连边
fi→i
f
i
→
i
(
fi=0
f
i
=
0
则不连)。
3. 把所有的
j<i
j
<
i
的
fj
f
j
改成
i
i
,并得到新的。
4. 删除点
i
i
,并把新的入度为的点加入集合。
随便拿个数据结构维护一下即可。
代码:
#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;
}