最短路径覆盖问题

可紊 提交于 2020-03-04 07:42:32

首先诚挚感谢BYvoid大神博客里提供的好题和数据

此为线性规划与网络流24题的第三题。


«问题描述:
给定有向图G=(V,E)。设P 是G 的一个简单路(顶点不相交)的集合。如果V 中每个
顶点恰好在P 的一条路上,则称P是G 的一个路径覆盖。P 中路径可以从V 的任何一个顶
点开始,长度也是任意的,特别地,可以为0。G 的最小路径覆盖是G 的所含路径条数最少
的路径覆盖。
设计一个有效算法求一个有向无环图G 的最小路径覆盖。

«编程任务:

对于给定的给定有向无环图G,编程找出G的一个最小路径覆盖。
«数据输入:
由文件input.txt提供输入数据。文件第1 行有2个正整数n和m。n是给定有向无环图
G 的顶点数,m是G 的边数。接下来的m行,每行有2 个正整数i和j,表示一条有向边(i,j)。
«结果输出:
程序运行结束时,将最小路径覆盖输出到文件output.txt 中。从第1 行开始,每行输出
一条路径。文件的最后一行是最少路径数。

建模部分分析:

考虑到在一个满足要求的路径覆盖中:

每个顶点属于且仅属于一条路径

每个点处最多出发一条边到达另一顶点

于是将一个点i拆成两个点Xi和Yi,分别代表从点i出发和到达点i。从而建立一个二分图,若原图中存在边(i,j),则在二分图中构造边(Xi,Yi),容量设为1。

再建立源点S,汇点T,S与每个X连接一条容量为1的边,T与每个Y连接一条容量为1的边。用最大流求最大二分匹配即可。

因为每增加一个匹配数,就能多“赚”到一个点。(贪心思想)

所以二分匹配数最大时,路径数最小。

最少的路径条数即为点的总数减去最大匹配数。

#include<cstdio>
#include<cstring>
#define rep(i,a,b) for(int i=a;i<=b;++i)
using namespace std;
const int MAXN=1010;
int graph[MAXN][MAXN];
int match[MAXN];
int visit[MAXN];
int path[MAXN];
int n,m,f=0;
void init()
{
    scanf("%d%d",&n,&m);
    int x,y;
    rep(i,1,m)
    {
        scanf("%d%d",&x,&y);
        graph[x][++graph[x][0]]=y;
    }
}
bool crosspath(int k)
{
    rep(i,1,graph[k][0])
    {
        int &t=graph[k][i];                 //试试能不能让k和t配
        if(!visit[t])
        {
            visit[t]=1;                     //标记t在第i轮寻找增广路的时候已经被访问过了
            if(!match[t]||crosspath(match[t]))          //如果t单身或者t的伴可以另外找到一个伴
            {
                match[t]=k;                              //把t许配给k
                path[k]=t;
                return true;
            }
        }
    }
    return false;
}
void hungary()
{
    rep(i,1,n)
    {
        if(crosspath(i)) f++;
        memset(visit,0,sizeof(visit));
    }
}
void pout()
{
    rep(i,1,n)
    {
        int t=i;
        if(!visit[i])
        {
            while(path[t]) t=path[t];
            while(match[t])
            {
                visit[t]=1;
                printf("%d ",t);
                t=match[t];
            }
            printf("%d\n",t);
        }
    }
    printf("%d",n-f);
}
int main()
{
    freopen("path2.in","r",stdin);
    init();
    hungary();
    pout();
    return 0;
}

 

易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!