题目描述
«问题描述:
给定有向图G=(V,E)。设P 是G 的一个简单路(顶点不相交)的集合。如果V 中每个顶点恰好在P 的一条路上,则称P是G 的一个路径覆盖。P 中路径可以从V 的任何一个顶点开始,长度也是任意的,特别地,可以为0。G 的最小路径覆盖是G 的所含路径条数最少的路径覆盖。设计一个有效算法求一个有向无环图G 的最小路径覆盖。提示:设V={1,2,…. ,n},构造网络G1=(V1,E1)如下:
每条边的容量均为1。求网络G1的( 0 x , 0 y )最大流。
«编程任务:
对于给定的给定有向无环图G,编程找出G的一个最小路径覆盖。
输入输出格式
输入格式:
件第1 行有2个正整数n和m。n是给定有向无环图G 的顶点数,m是G 的边数。接下来的m行,每行有2 个正整数i和j,表示一条有向边(i,j)。
输出格式:
从第1 行开始,每行输出一条路径。文件的最后一行是最少路径数。
输入输出样例
输入样例#1:
11 12
1 2
1 3
1 4
2 5
3 6
4 7
5 8
6 9
7 10
8 11
9 11
10 11
输出样例#1:
1 4 7 10 11
2 5 8
3 6 9
3
说明
1<=n<=150,1<=m<=6000
建模思路:
先分成两边每一边都有n个点,然后就是将有向边连起来,最后跑一遍网络流就好了。
输出路径的问题:
最大流中流量为1的边就是匹配边,先处理to[i],从vis[i]==0的点开始。
代码:
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
const int N=1e5+10;
int n,m,S,T,to[N],head[N],cur[N],dis[N];
struct node{
int v,next,cap,flow;
}e[N<<1];int tot=1;
bool vis[N];
void add(int x,int y,int z)
{
e[++tot].v=y;
e[tot].cap=z;
e[tot].next=head[x];
e[tot].flow=0;
head[x]=tot;
}
bool bfs(){
memset(dis,0,sizeof dis);
queue<int>q;
q.push(S);dis[S]=1;
while(!q.empty()){
int x=q.front();q.pop();
for(int i=head[x];i;i=e[i].next)
{
int v=e[i].v;
if(!dis[v]&&e[i].cap>e[i].flow){
dis[v]=dis[x]+1;
q.push(v);
}
}
}
return dis[T];
}
int dfs(int x,int f){
if(x==T) return f;
for(int i=head[x];i;i=e[i].next)
{
int v=e[i].v,f1;
if(dis[x]+1==dis[v]&&e[i].cap>e[i].flow)
{
if(f1=dfs(v,min(f,e[i].cap-e[i].flow)))
{
e[i].flow+=f1;e[i^1].flow-=f1;
return f1;
}
}
}
return 0;
}
int dinic(){
int ans=0;
while(bfs()) ans+=dfs(S,0x7fffffff);
return ans;
}
int main(){
scanf("%d%d",&n,&m);S=0;T=n<<1|1;
for(int i=1,x,y;i<=m;i++)
{
scanf("%d%d",&x,&y);
add(x,y+n,1);
add(y+n,x,0);
}
for(int i=1;i<=n;i++)
{
add(S,i,1);
add(i,S,0);
add(i+n,T,1);
add(T,i+n,0);
}
int ans=n-dinic();
for(int i=1;i<=n;i++)
{
for(int j=head[i];j;j=e[j].next)
{
if(e[j].flow==1)
{
to[i]=e[j].v-n;
break;
}
}
}
for(int i=1;i<=n;i++) {
if(!vis[i]){
while(i) vis[i]=1,printf("%d ",i),i=to[i];
printf("\n");
}
}
printf("%d\n",ans);
}