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

把原图的每个点V拆成Vx和Vy两个点，如果有一条有向边A->B，那么就加边Ax−>By。这样就得到了一个二分图。那么最小路径覆盖=原图的结点数-新图的最大匹配数。

一开始每个点都是独立的为一条路径，总共有n条不相交路径。我们每次在二分图里找一条匹配边就相当于把两条路径合成了一条路径，也就相当于路径数减少了1。所以找到了几条匹配边，路径数就减少了多少。所以有最小路径覆盖=原图的结点数-新图的最大匹配数。

因为路径之间不能有公共点，所以加的边之间也不能有公共点，这就是匹配的定义。

转自：https://blog.csdn.net/qq_39627843/article/details/82012572

#include<bits/stdc++.h>
using namespace std;

struct node{
    int v,nx;
}edge[20000];

int tot = 1,head[500];

inline void add(int u,int v)
{
    edge[++tot].v = v;
    edge[tot].nx = head[u];
    head[u] = tot;
}

int cx[500],cy[500],vis[500],n,m,in[500];

bool dfs(int u)
{
    for (int i = head[u];i;i = edge[i].nx)
    {
        int v = edge[i].v;
        if (vis[v]) continue;
        vis[v] = 1;
        if (cy[v] == 0 ||dfs(cy[v]))
        {
            cy[v] = u;
            cx[u] = v;
            return true;
        }
    }
    return false;
}

int maxmatch()
{
    int ans = 0;
    for (int i = 1;i<=n;i++)
    {
        if (cx[i] == 0)
        {
            memset(vis,0,sizeof(vis));
            ans += dfs(i);
        }
    }
    //拓扑排序输出最小不相交路径
    memset(vis,0,sizeof(vis));
    for (int i = 1;i<=n;i++)
    {
        int tmp = i;
        if (in[i] == 0 && !vis[i])
        {
            printf("%d ",i);
            vis[i] = 1;
            while (cx[i] != 0)
            {
                for (int j = head[i];j;j = edge[j].nx)
                {
                    in[edge[j].v - 150]--;
                }
                vis[i = cx[i] - 150] = 1;
                printf("%d ",i);
            }
            printf("\n");
        }
        i = tmp;
    }
    return ans;
}

int main()
{
    scanf("%d%d",&n,&m);
    //将DAG每个点拆成两个点，0-150只有出度，150-300只有入度
    for (int i = 1,u,v;i<=m;i++)
    {
        scanf("%d%d",&u,&v);
        add(u,v+150);
        //in记录原来每个点的入度，用来后面输出最小不相交路径
        in[v]++;
    }
    cout<<n - maxmatch();
    return 0;
}