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

本文链接：https://blog.csdn.net/elijahqi/article/details/79124197

(http://www.elijahqi.win/2017/12/15/loj-6002-%E3%80%8C%E7%BD%91%E7%BB%9C%E6%B5%81-24-%E9%A2%98%E3%80%8D%E6%9C%80%E5%B0%8F%E8%B7%AF%E5%BE%84%E8%A6%86%E7%9B%96/%20%E2%80%8E)
题目描述

«问题描述：

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

由@zhouyonglong提供SPJ

首先建图然后从源点向每个节点的入点都连权值为1的边从每个点的出点都向汇点连权值为1的边然后如果这两点之间原来有边那么从起点的出点向终点的入点连1的边跑dinic完成最大匹配那么最后答案就是节点数-最大匹配数

我们考虑假如这两个点可以匹配那我是不是可以用一条路径给覆盖初始我每个节点都需要路径来覆盖现在有一个匹配我就可以减少1 于是答案就是这样

关于输出路径来说我应该寻找每个点然后看他反向边的残余流量如果是0 那么那个点就在我的路径上可以画图感受一下我比较傻给不出特别严谨的证明qwq

#include<queue>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define inf 0x3f3f3f3f
#define N 400
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0;char ch=gc();
    while(ch<'0'||ch>'9') ch=gc();
    while(ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=gc();}
    return x;
}
struct node{
    int y,next,z;
}data[15000];
int num=1,h[N],T,level[N],n,m;
bool visit[N];
inline void insert1(int x,int y,int z){
    data[++num].y=y;data[num].next=h[x];h[x]=num;data[num].z=z;
    data[++num].y=x;data[num].next=h[y];h[y]=num;data[num].z=0;
}
inline bool bfs(){
    queue<int> q;q.push(0);memset(level,0,sizeof(level));level[0]=1;
    while(!q.empty()){
        int x=q.front();q.pop();
        for (int i=h[x];i;i=data[i].next){
            int y=data[i].y,z=data[i].z;
            if (level[y]||!z) continue;level[y]=level[x]+1;q.push(y);if(y==T) return 1;
        }
    }return 0;
}
inline int dfs(int x,int s){
    if (x==T) return s;int ss=s;
    for (int i=h[x];i;i=data[i].next){
        int y=data[i].y,z=data[i].z;
        if (level[y]==level[x]+1&&z){
            int xx=dfs(y,min(z,s));if (!xx) level[y]=0;
            s-=xx;data[i].z-=xx;data[i^1].z+=xx;if (!s) return ss;
        }
    }return ss-s;
}
int main(){
    freopen("loj.in","r",stdin);
    freopen("loj.out","w",stdout);
    n=read();m=read();T=n*2+1;
    for (int i=1;i<=n;++i) insert1(0,i,1),insert1(i+n,T,1);
    for (int i=1;i<=m;++i){int x=read(),y=read();insert1(x,y+n,1);}int ans=0;
    while(bfs()) ans+=dfs(0,inf);
    for (int i=1;i<=n;++i){
        if (visit[i]) continue;int now=i;bool flag=0;
        do{
            printf("%d ",now);visit[now]=1;flag=0;
            for (int i=h[now];i;i=data[i].next)
                if (!data[i].z&&data[i].y>n&&data[i].y<=2*n){now=data[i].y-n;flag=1;break;}
        }while(flag);puts("");
    }
    printf("%d",n-ans);
    return 0;
}