DAG最小路径覆盖

参考博客:https://www.cnblogs.com/justPassBy/p/5369930.html

DAG的最小路径覆盖

定义:在一个有向图中,找出最少的路径,使得这些路径经过了所有的点。

最小路径覆盖分为最小不相交路径覆盖最小可相交路径覆盖

最小不相交路径覆盖:每一条路径经过的顶点各不相同。如图,其最小路径覆盖数为3。即1->3>4,2,5。

最小可相交路径覆盖:每一条路径经过的顶点可以相同。如果其最小路径覆盖数为2。即1->3->4,2->3>5。

特别的,每个点自己也可以称为是路径覆盖,只不过路径的长度是0。

DAG的最小不相交路径覆盖

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

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

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

习题http://acm.hdu.edu.cn/showproblem.php?pid=1151

#include <stdio.h>
#include <string.h>

const int maxn = 150;
int match[maxn],G[maxn][maxn],vis[maxn];
int n,m;

bool find(int u)
{
    for(int i = 1; i <= n; i++){
        if(G[u][i] && !vis[i]){
            vis[i] = 1;
            if(match[i] == -1 || find(match[i])){
                match[i] = u;
                return true;
            }
        }
    }
    return false;
}
int max_match()
{
    int ans = 0;
    memset(match,-1,sizeof(match));
    for(int i = 1; i <= n; i++){
        memset(vis,0,sizeof(vis));
        if(find(i)){
            ans++;
        }
    }
    return ans;
}
int main(void)
{
    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        memset(G,0,sizeof(G));
        for(int i = 1; i <= m; i++){
            int a,b;
            scanf("%d%d",&a,&b);
            G[a][b] = 1;
        }
        printf("%d\n",n - max_match());
    }
    return 0;
}

DAG的最小可相交路径覆盖

 

算法:先用floyd求出原图的传递闭包,即如果a到b有路径,那么就加边a->b。然后就转化成了最小不相交路径覆盖问题。

证明:为了连通两个点,某条路径可能经过其它路径的中间点。比如1->3->4,2->4->5。但是如果两个点a和b是连通的,只不过中间需要经过其它的点,那么可以在这两个点之间加边,那么a就可以直达b,不必经过中点的,那么就转化成了最小不相交路径覆盖。

习题:http://poj.org/problem?id=2594

#include <stdio.h>
#include <string.h>
const int maxn = 505;
int match[maxn],G[maxn][maxn],vis[maxn];
int n,m;
bool find(int u)
{
    for(int i = 1; i <= n; i++) {
        if(G[u][i] && !vis[i]) {
            vis[i] = 1;
            if(match[i] == -1 || find(match[i])) {
                match[i] = u;
                return true;
            }
        }
    }
    return false;
}
int max_match()
{
    int ans = 0;
    memset(match,-1,sizeof(match));
    for(int i = 1; i <= n; i++) {
        memset(vis,0,sizeof(vis));
        if(find(i)) {
            ans++;
        }
    }
    return ans;
}
void floyd()
{
    for(int k = 1; k <= n; k++) {
        for(int i = 1; i <= n; i++) {
            for(int j = 1; j <= n; j++) {
                if(G[i][k] && G[k][j]) {
                    G[i][j] = 1;
                }
            }
        }
    }
}
int main(void)
{
    int T;
    while( scanf("%d%d",&n,&m) != EOF && (n + m)) {
        memset(G,0,sizeof(G));
        for(int i = 1; i <= m; i++) {
            int a,b;
            scanf("%d%d",&a,&b);
            G[a][b] = 1;
        }
        floyd();
        printf("%d\n",n - max_match());
    }
    return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值