[网络流24题]最小路径覆盖问题(二分图匹配)

题目

传送门

题解

原理:最小路径覆盖数=点数-匹配数;
很好理解:两个点可以匹配,相当于是放到并查集里面,并查集的个数就是路径的条数;
但是我对建图的理解一直就是错误的;
我原来跑了二分图染色分成了两个集合,在里面跑二分图匹配;错误之处在于我错误的没有拆点;
于是我又认为每个点只能走一次,所以就每个点向所拆的点之间连了容量为1的边;错误之处在于电和点之间无法确立联系,没有调整的作用
正确的建图方法是:每个点都拆点 i 和i+n 但是不相连;S连向i 容量1; i+n连向T 容量1;如果两个点相连x连向 y+n (体会一下这样建图的巧妙);可以说两个点具有隐形的联系;
以 1–>2–>3为例
[网络流24题]最小路径覆盖问题(二分图匹配)_第1张图片
但是建图有难倒了我;nxt[]和pre[]分别记录点的后继和前驱;递归输出;

代码

#include
#include
#include
#include
#include
using namespace std;
const int maxn=1000001;
const int inf=1e9;
queue <int> q;
int n,m,k,maxflow,a[maxn],S,T;
struct Edge{
    int next,to,dis;
}edge[maxn<<1];
int num_edge=-1,head[maxn],cur[maxn],deep[maxn];
int nxt[maxn],pre[maxn];//分别表示点i的前一个点和后一个点是什么 

void add_edge(int from,int to,int dis)
{
    edge[++num_edge].next=head[from];
    edge[num_edge].dis=dis;
    edge[num_edge].to=to;
    head[from]=num_edge;
}
void add(int x,int y,int z) {add_edge(x,y,z); add_edge(y,x,0);}

bool bfs(int s,int t)
{
    memset(deep,0x7f,sizeof(deep));
    for (int i=0; i<=t; i++) cur[i]=head[i];
    while (!q.empty()) q.pop();
    q.push(s); deep[s]=0;
    while (!q.empty())
    {
        int now=q.front(); q.pop();
        for (int i=head[now]; i!=-1; i=edge[i].next)
        {
            int to=edge[i].to;
            if (deep[to]>inf && edge[i].dis)
            {
                deep[to]=deep[now]+1;
                q.push(to);
            }
        }
    }
    return deep[t]int dfs(int now,int t,int limit)
{
    if (now==t || !limit) return limit;
    int flow=0,f;
    for (int i=cur[now]; i!=-1; i=edge[i].next)
    {
        int to=edge[i].to; cur[now]=i;
        if (deep[to]==deep[now]+1 && (f=dfs(to,t,min(limit,edge[i].dis))))
        {
            flow+=f;
            limit-=f;
            edge[i].dis-=f;
            edge[i^1].dis+=f;
            if (!limit && now!=S && to!=T) nxt[now]=to-n; pre[to-n]=now;
            if (!limit) break;
        }
    }
    return flow;
}

void Dinic(int s,int t)
{
    while (bfs(s,t))
        maxflow+=dfs(s,t,inf);
}
bool vis[maxn];
void shuchu(int x)
{
    if (!x) return;//到达了源点 
    if (pre[x]!=x) shuchu(pre[x]);
    vis[x]=1;
    printf("%d ",x);
}

int main()
{
    memset(head,-1,sizeof(head));
    scanf("%d%d",&n,&m);
    S=0; T=2*n+1;
    for (int i=1; i<=n; i++) nxt[i]=pre[i]=i;
    for (int i=1; i<=n; i++) {add(S,i,1); add(i+n,T,1);}
    for (int i=1; i<=m; i++)
    {
        int x,y; scanf("%d%d",&x,&y);
        add(x,y+n,1);
    }
//  for (int i=0; i<=num_edge; i++) printf("%d: %d %d %d\n",i,edge[i^1].to,edge[i].to,edge[i].dis);
    Dinic(S,T);
//  for (int i=1; i<=n; i++) printf("%d %d\n",nxt[i],pre[i]);
    for (int i=n; i>=1; i--)
        if (nxt[i]==i && !vis[i]) {shuchu(i); printf("\n");}
    printf("%d",n-maxflow);
    return 0;
}

总结

输出方法

你可能感兴趣的:(网络流)