zoukankan      html  css  js  c++  java
  • luogu P2764 最小路径覆盖问题 网络流24

    #include <iostream>
    #include <cstdio>
    #include <cstring>
    #define N 500000
    using namespace std;
    int n,m,S,T,tmp1,tmp2,tot;
    int idx,head[N],cur[N],q[N],ans1[N],ans2[N];
    int e[N],ne[N],w[N];
    bool b[N];
    int r[N],c[N];
    bool vis[N];
    inline void add(int u,int v,int f)
    {
        e[idx]=v;
        w[idx]=f;
        ne[idx]=head[u];
        head[u]=idx++;
    }
    inline bool bfs()
    {
        int f=0,t=0;
        memset(cur,-1,sizeof(cur));
        q[t++]=S;
        cur[S]=0;
        while(f<t)
        {
            int now=q[f++];
            for(int i=head[now]; ~i; i=ne[i])
            {
                int v=e[i];
                if(cur[v]==-1&&w[i])
                {
                    cur[v]=cur[now]+1;
                    q[t++]=v;
                }
            }
        }
        if(cur[T]!=-1)
            return 1;
        return 0;
    }
    inline int dfs(int x,int f)
    {
        if(x==T)
            return f;
        int w1,used=0;
        for(int i=head[x]; ~i; i=ne[i])
        {
            int v=e[i];
            if(cur[v]==cur[x]+1&&w[i])
            {
                w1=dfs(v,min(f-used,w[i]));
                w[i]-=w1;
                w[i^1]+=w1;
                used+=w1;
                if(used==f)
                    return f;
            }
        }
        if(!used)
            cur[x]=-1;
        return used;
    }
    void dinic()
    {
        while(bfs())
            tot+=dfs(S,0x3f3f3f3f);
    }
    int main()
    {
        memset(head,-1,sizeof head);
        cin>>n>>m;
        S=0,T=n*2+1;
        for(int i=1; i<=n; i++)
        {
            add(S,i,1);
            add(i,S,0);
        }
        for(int i=1;i<=n;i++)
        {
            add(i+n,T,1);
            add(T,i+n,0);
        }
        for(int i=1; i<=m; i++)
        {
            int p,q;
            cin>>p>>q;
            add(p,q+n,1);
            add(q+n,p,0);
        }
        dinic();
        for(int i=1;i<=n;i++)
            if(!vis[i])
            {
                int now=i;
                do
                {
                    vis[now]=1;
                    cout<<now<<" ";
                    for(int j=head[now];j!=-1;j=ne[j])
                    {
                        int v=e[j];
                        if(v<=n)
                            continue;
                        if(!w[j])
                        {
                            now=v-n;
                            break;
                        }
                    }
                }while(!vis[now]);
                cout<<endl;
            }
        cout<<n-tot<<endl;
        return 0;
    }
  • 相关阅读:
    改变JupyterLab的启动路径
    副业创收的误区,为啥你总掉进坑里[转]
    消费主义陷阱[转]
    排列组合生成算法CombinationAll
    nlp跳坑基础
    通俗易懂告诉你CPU/GPU/TPU/NPU...XPU都是些什么鬼?
    Shotcuts in linux terminal for typing commands 各类终端的快捷键
    Python Why?
    2020年15.5以后的QT入坑指南
    PreparedStatement的用法
  • 原文地址:https://www.cnblogs.com/QingyuYYYYY/p/13169174.html
Copyright © 2011-2022 走看看