其他分享
首页 > 其他分享> > 【网络流24题】【洛谷P2764】最小路径覆盖问题【网络流】

【网络流24题】【洛谷P2764】最小路径覆盖问题【网络流】

作者:互联网

题目大意:

题目链接:https://www.luogu.org/problemnew/show/P2764

给定有向图G=(V,E)G=(V,E)G=(V,E)。设PPP是GGG的一个简单路(顶点不相交)的集合。如果VVV中每个定点恰好在PPP的一条路上,则称PPP是GGG的一个路径覆盖。PPP中路径可以从VVV的任何一个定点开始,长度也是任意的,特别地,可以为000。GGG的最小路径覆盖是GGG所含路径条数最少的路径覆盖.设计一个有效算法求一个GAP\text{GAP}GAP(有向无环图)GGG的最小路径覆盖。


思路:

由于题目中要求任意一个点都必须正好走一次,所以这肯定是要拆点的,把每个点拆成x1,x2x_1,x_2x1​,x2​。
对于GAPGAPGAP中每一条边(u,v)(u,v)(u,v),我们建边(u2,v)(u_2,v)(u2​,v),流量为1。这样就限制了每一条边只能走一次。
然后将源点连向x2(x[1,n])x_2(x\in [1,n])x2​(x∈[1,n]),x1(x[1,n])x_1(x\in [1,n])x1​(x∈[1,n])连向汇点,跑最大流,就得到了最多可以合并几条道路。因为流一次就相当于合并了一条道路。
那么最小路径覆盖就是 nn-n−最大流 了。
输出方案的话可以利用残余流量,递归输出方案。


代码:

#include <queue>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

const int N=100010,Inf=1e9;
int n,m,S,T,tot=1,maxflow,head[N],dep[N],cur[N];
bool vis[N];

struct edge
{
    int from,to,flow,next;
}e[N];

void add(int from,int to,int flow)
{
    e[++tot].to=to;
    e[tot].flow=flow;
    e[tot].next=head[from];
    head[from]=tot;
}

bool bfs()
{
    memcpy(cur,head,sizeof(cur));
    memset(dep,0x3f3f3f3f,sizeof(dep));
    queue<int> q;
    q.push(S);
    dep[S]=1;
    while (q.size())
    {
        int u=q.front(),v;
        q.pop();
        for (int i=head[u];~i;i=e[i].next)
        {
            v=e[i].to;
            if (dep[v]>dep[u]+1&&e[i].flow)
            {
                dep[v]=dep[u]+1;
                q.push(v);
            }
        }
    }
    return dep[T]<Inf;
}

int dfs(int x,int flow)
{
    int low=0;
    if (x==T)
    {
        maxflow+=flow;
        return flow;
    }
    int used=0;
    for (int i=cur[x];~i;i=e[i].next)
    {
        int y=e[i].to;
        cur[x]=i;
        if (dep[y]==dep[x]+1&&e[i].flow)
        {
            low=dfs(y,min(e[i].flow,flow-used));
            if (low)
            {
                used+=low;
                e[i].flow-=low;
                e[i^1].flow+=low;
                if (used==flow) break;
            }
        }
    }
    return used;
}

void dinic()
{
    while (bfs())
        dfs(S,Inf);
}

void find(int x)
{
    vis[x-n]=1;
    printf("%d ",x-n);
    for (int i=head[x];~i;i=e[i].next)
    {
        int y=e[i].to;
        if (!e[i].flow && y!=x-n && !vis[y])
        {
            find(y+n);
            return;
        }
    }
}

int main()
{
    memset(head,-1,sizeof(head));
    scanf("%d%d",&n,&m);
    for (int i=1;i<=m;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        add(x+n,y,1);
        add(y,x+n,0);
    }
    S=N-1; T=N-2;
    for (int i=1;i<=n;i++)
    {
        add(S,i+n,1);
        add(i+n,S,0);
        add(i,T,1);
        add(T,i,0);
    }
    dinic();
    vis[S]=vis[T]=1;
    for (int i=1;i<=n;i++)
        if (!vis[i])
        {
            find(i+n);
            putchar(10);
        }
    printf("%d\n",n-maxflow);
    return 0;
}

标签:24,head,洛谷,dep,路径,网络,tot,int,GGG
来源: https://blog.csdn.net/SSL_ZYC/article/details/91350229