loj #6002. 「网络流 24 题」最小路径覆盖(二分匹配)

#6002. 「网络流 24 题」最小路径覆盖

内存限制：256 MiB 时间限制：1000 ms 标准输入输出

题目类型：传统 评测方式：Special Judge

上传者： 匿名

提交 提交记录 统计 讨论 测试数据

题目描述

给定有向图 $G=(V,E)$。设 $P$ 是 $G$ 的一个简单路（顶点不相交）的集合。如果 $V$ 中每个顶点恰好在 $P$ 的一条路上，则称 $P$ 是 $G$ 的一个路径覆盖。$P$ 中路径可以从 $V$ 的任何一个顶点开始，长度也是任意的，特别地，可以为 $0$。$G$ 的最小路径覆盖是 $G$ 的所含路径条数最少的路径覆盖。

设计一个有效算法求一个有向无环图 $G$ 的最小路径覆盖。

输入格式

第 $1$ 行有 $2$ 个正整数 $n$ 和 $m$。$n$ 是给定有向无环图 $G$ 的顶点数，$m$ 是 $G$ 的边数。
接下来的 $m$ 行，每行有 $2$ 个正整数 $u$ 和 $v$，表示一条有向边 $(i,j)$。

输出格式

从第 $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 4 7 10 11
2 5 8
3 6 9
3

数据范围与提示

$1\le n\le 200,1\le m\le 6000$

这题本来是用网络流做的，但是用二分匹配也很方便就用 匹配做了

直接拆点 连边 就是二分匹配里常见的 独立子集的模板

#include <iostream>
#include <cstdio>
#include <cstring>
#include <bits/stdc++.h>
using namespace std;
typedef long long LL ;
const int inf = 500;
const int N = 1000;
const double eps = 1e-7;
struct node
{
    int to, next;
}p[1000000];
int head[N], cnt;
int match[N], w[N], vis[N];
void init()
{
    memset(head,-1,sizeof(head));
    memset(match,-1,sizeof(match));
    memset(w,0,sizeof(w));
    cnt=0;
    return ;
}
void add(int u,int v)
{
    p[cnt].to=v,p[cnt].next=head[u];head[u]=cnt++;
    return ;
}
int dfs(int u)
{
    for(int i=head[u];i!=-1;i=p[i].next)
    {
        int v=p[i].to;
        if(vis[v]) continue;
        vis[v]=1;
        if(match[v]==-1||dfs(match[v]))
        {
            match[v]=u,w[u]=v;
            return 1;
        }
    }
    return 0;
}
int ans[N];

int main()
{
    init();
    int n, m;
    scanf("%d %d", &n, &m);
    for(int i=1;i<=m;i++)
    {
        int a, b;
        scanf("%d %d", &a, &b);
        add(n+a,b);
    }
    for(int i=n+1;i<=2*n;i++)
    {
        memset(vis,0,sizeof(vis));
        dfs(i);
    }
    int k=0;
    for(int i=1;i<=n;i++)  if(match[i]==-1) ans[++k]=i;
    for(int i=1;i<=k;i++)
    {
        printf("%d",ans[i]);
        int x=n+ans[i];
        while(w[x])
        {
            printf(" %d",w[x]);
            x=n+w[x];
        }
        puts("");
    }
    printf("%d\n",k);
    return 0;
}