网络流24题3最小路径覆盖问题（洛谷 P2764 最小路径覆盖问题）

题目描述

«问题描述：

给定有向图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);
}