题目大意:
计算一个有向图的最小点路径覆盖并打印路径。
思路分析:
1:对于一个路径覆盖,有如下概念:
①:每个顶点只属于一个路径;
②:路径上除终点外,从每个顶点出发只有一条边指向路径上的另一顶点。
2:有向无环图的最小路径覆盖问题包括两种(设G是一个有向无环图,S是G的一个路径集合):
①:最小边路径覆盖:满足对于G中所有的边E,E在S中的一条路径中出现,且只在S中的一条路径中出现,求S的最小容量;
建模方法:暂略。
②:最小点路径覆盖:满足对于G中所有的点i,i在S中的一条路径中出现,且只在S中的一条路径中出现,求S的最小容量;建模方法:有向无环图的最小路径覆盖可以转化成二分图最大匹配问题,即建立一个新图,将G中的每个点i在新图中拆成两个点i'、i'',若G中存在边<i, j>则在新图中连边<i', j''>,显然新图是一个二分图,求其最大匹配,则(N-新图最大匹配的值)就是最小点路径覆盖值。
建模分析:我们可以把每个顶点理解成两个顶点,一个是出发,一个是目标,建立二分图模型。该二分图的任何一个匹配方案,都对应了一个路径覆盖方案。如果匹配数为0,那么显然路径数=顶点数。每增加一条匹配边,那么路径覆盖数就减少一个,所以路径数=顶点数 - 匹配数。要想使路径数最少,则应最大化匹配数,所以要求二分图的最大匹配。
代码实现:
#include<cstdio>
#include<cstring>
#include<queue>
#include<iostream>
using namespace std;
const int N=310,M=12610,INF=0x3f3f3f3f;
int n,m,top,s,t,head[N],dis[N],gap[N],pre[N],cur[N],path[N],cnt,vis[N];
struct Edge{
int to,next,flow;
Edge(int _to=0,int _next=0,int _flow=0):to(_to),next(_next),flow(_flow){}
}edge[M];
void Addedge(int from,int to,int flow){
edge[top]=Edge(to,head[from],flow);
head[from]=top++;
edge[top]=Edge(from,head[to],0);
head[to]=top++;
}
void Bfs(){
queue<int> q;
memset(dis,-1,sizeof(dis));
memset(gap,0,sizeof(gap));
gap[0]=1,dis[t]=0;
q.push(t);
while(!q.empty()){
int u=q.front();
q.pop();
for(int i=head[u];i+1;i=edge[i].next){
if(dis[edge[i].to]==-1){
dis[edge[i].to]=dis[u]+1;
gap[dis[edge[i].to]]++;
q.push(edge[i].to);
}
}
}
}
int Sap(){
Bfs();
memset(pre,-1,sizeof(pre));
for(int i=0;i<=t;++i) cur[i]=head[i];
int u=s,cur_flow,max_flow=0,neck,tmp;
while(dis[s]<=t){
if(u==t){
cur_flow=INF;
for(int i=s;i!=t;i=edge[cur[i]].to){
if(cur_flow>edge[cur[i]].flow){
neck=i;
cur_flow=edge[cur[i]].flow;
}
}
for(int i=s;i!=t;i=edge[cur[i]].to){
tmp=cur[i];
edge[tmp].flow-=cur_flow;
edge[tmp^1].flow+=cur_flow;
}
max_flow+=cur_flow;
u=neck;
}
int i;
for(i=cur[u];i+1;i=edge[i].next){
if(edge[i].flow&&dis[u]==dis[edge[i].to]+1) break;
}
if(i!=-1){
cur[u]=i;
pre[edge[i].to]=u;
u=edge[i].to;
}else{
if(--gap[dis[u]]==0) break;
cur[u]=head[u];
int mindis=t;
for(i=head[u];i+1;i=edge[i].next){
if(edge[i].flow&&mindis>dis[edge[i].to])
mindis=dis[edge[i].to];
}
dis[u]=mindis+1;
gap[dis[u]]++;
if(u!=s) u=pre[u];
}
}
return max_flow;
}
void Dfs(int u){
vis[u]=1;
path[cnt++]=u;
for(int i=head[u];i+1;i=edge[i].next){
if(!vis[edge[i].to]&&edge[i].to&&!edge[i].flow) Dfs(edge[i].to-n);
}
}
int main(){
freopen("path3.in","r",stdin);
freopen("path3.out","w",stdout);
memset(head,-1,sizeof(head));
memset(vis,0,sizeof(vis));
scanf("%d%d",&n,&m);
top=s=0,t=n*2+1;
int u,v;
for(int i=1;i<=n;++i) Addedge(s,i,1),Addedge(i+n,t,1);
for(int i=0;i<m;++i){
scanf("%d%d",&u,&v);
Addedge(u,v+n,1);
}
int res=Sap();
for(int i=1;i<=n;++i){
cnt=0;
if(!vis[i]){
Dfs(i);
printf("%d",path[0]);
for(int j=1;j<cnt;++j) printf(" %d",path[j]);
printf("\n");
}
}
printf("%d\n",n-res);
}