• zoj3229 有源汇上下界最大流


    题意:有一个人每天给妹子拍照,每个妹子有最少拍照数,每天有最大拍照数,每天只能给某些特定的妹子拍照,求最大拍照数

    题解:很容易看出来的有源汇上下界最大流,对于有源汇 的上下界最大流,我们按照无源汇的操作连边,然后从汇点连一条边到源点,容量无穷大,下界为0,从超级源ss到超级汇tt跑一遍最大流,求得的流量就是可行流 的流量,此时求出来的只是可行流并不是最大流,我们还需要从原来的s到t跑一边最大流,最后的最大流就是可行流 的流量+新增广的s-t流量

    ps:这题样例一直过不了,但是我发现很多博客代码的样例输出结果和我的一样= =。还有这题已经挂了,判不出来,一直是Judge Internal Error,所以不知道代码对不对,先放着。

    #include<bits/stdc++.h>
    #include<ext/rope>
    #define fi first
    #define se second
    #define mp make_pair
    #define pb push_back
    #define pii pair<int,int>
    #define C 0.5772156649
    #define pi acos(-1.0)
    #define ll long long
    #define mod 1000000007
    #define ls l,m,rt<<1
    #define rs m+1,r,rt<<1|1
    
    using namespace std;
    using namespace __gnu_cxx;
    
    const double g=10.0,eps=1e-7;
    const int N=2000+10,maxn=400000+10,inf=0x3f3f3f;
    
    struct edge{
        int from,to,Next,c,low;
    }e[maxn<<1];
    int dis[N];
    int cnt,head[N];
    int in[N],out[N];
    int ansid[maxn];
    void add(int u,int v,int c,int low)
    {
        out[u]+=low;
        in[v]+=low;
        e[cnt].from=u;
        e[cnt].to=v;
        e[cnt].low=low;
        e[cnt].c=c;
        e[cnt].Next=head[u];
        head[u]=cnt++;
        e[cnt].from=v;
        e[cnt].to=u;
        e[cnt].low=low;
        e[cnt].c=0;
        e[cnt].Next=head[v];
        head[v]=cnt++;
    }
    void init()
    {
        cnt=0;
        memset(head,-1,sizeof head);
        memset(in,0,sizeof in);
        memset(out,0,sizeof out);
    }
    bool bfs(int s,int t)
    {
        memset(dis,-1,sizeof dis);
        dis[s]=0;
        queue<int>q;
        q.push(s);
        while(!q.empty())
        {
            int x=q.front();
            q.pop();
            if(x==t)return 1;
            for(int i=head[x];~i;i=e[i].Next)
            {
                int te=e[i].to;
                if(dis[te]==-1&&e[i].c>0)
                {
                    dis[te]=dis[x]+1;
                    q.push(te);
                }
            }
        }
        return 0;
    }
    int dfs(int x,int mx,int t)
    {
        if(x==t)return mx;
        int flow=0;
        for(int i=head[x];~i;i=e[i].Next)
        {
            int te=e[i].to,f;
            if(e[i].c>0&&dis[te]==dis[x]+1&&(f=dfs(te,min(mx-flow,e[i].c),t)))
            {
                e[i].c-=f;
                e[i^1].c+=f;
                flow+=f;
            }
        }
        if(!flow)dis[x]=-2;
        return flow;
    }
    int maxflow(int s,int t)
    {
        int ans=0,f;
        while(bfs(s,t))
        {
            while((f=dfs(s,inf,t)))ans+=f;
        }
        return ans;
    }
    int main()
    {
        //ios::sync_with_stdio(false);
        //cin.tie(0);
        int n,m;
        while(cin>>n>>m)
        {
            init();
            int s=n+m+1,t=n+m+2;//原图
            int ss=n+m+3,tt=n+m+4;//超级源超级汇
            for(int i=1;i<=m;i++)
            {
                int a;
                cin>>a;
                add(n+i,t,inf-a,a);
            }
            int res=0;
            for(int i=1;i<=n;i++)
            {
                int k,d;
                cin>>k>>d;
                add(s,i,d,0);
                while(k--)
                {
                    int a,b,c;
                    cin>>a>>b>>c;
                    ansid[++res]=cnt;
                    add(i,n+a+1,c-b,b);
                }
            }
            add(t,s,inf,0);
            int sum=0;
            for(int i=1;i<=n+m+2;i++)
            {
                if(in[i]>out[i])sum+=in[i]-out[i],add(ss,i,in[i]-out[i],0);
                else add(i,tt,out[i]-in[i],0);
            }
            if(sum!=maxflow(ss,tt))cout<<-1<<endl;
            else
            {
                cout<<maxflow(s,t)<<endl;
                for(int i=1;i<=res;i++)
                    cout<<e[ansid[i]^1].c+e[ansid[i]].low<<endl;
            }
            cout<<endl;
        }
        return 0;
    }
    /*******************
    
    ********************/
    View Code

     update:Tle,cin->scanf,ac

    #include<bits/stdc++.h>
    #include<ext/rope>
    #define fi first
    #define se second
    #define mp make_pair
    #define pb push_back
    #define pii pair<int,int>
    #define C 0.5772156649
    #define pi acos(-1.0)
    #define ll long long
    #define mod 1000000007
    #define ls l,m,rt<<1
    #define rs m+1,r,rt<<1|1
    
    using namespace std;
    using namespace __gnu_cxx;
    
    const double g=10.0,eps=1e-7;
    const int N=2000+10,maxn=400000+10,inf=0x3f3f3f;
    
    struct edge{
        int from,to,Next,c,low;
    }e[maxn<<1];
    int dis[N];
    int cnt,head[N];
    int in[N],out[N];
    int ansid[maxn];
    void add(int u,int v,int c,int low)
    {
        out[u]+=low;
        in[v]+=low;
        e[cnt].from=u;
        e[cnt].to=v;
        e[cnt].low=low;
        e[cnt].c=c;
        e[cnt].Next=head[u];
        head[u]=cnt++;
        e[cnt].from=v;
        e[cnt].to=u;
        e[cnt].low=low;
        e[cnt].c=0;
        e[cnt].Next=head[v];
        head[v]=cnt++;
    }
    void init()
    {
        cnt=0;
        memset(head,-1,sizeof head);
        memset(in,0,sizeof in);
        memset(out,0,sizeof out);
    }
    bool bfs(int s,int t)
    {
        memset(dis,-1,sizeof dis);
        dis[s]=0;
        queue<int>q;
        q.push(s);
        while(!q.empty())
        {
            int x=q.front();
            q.pop();
            if(x==t)return 1;
            for(int i=head[x];~i;i=e[i].Next)
            {
                int te=e[i].to;
                if(dis[te]==-1&&e[i].c>0)
                {
                    dis[te]=dis[x]+1;
                    q.push(te);
                }
            }
        }
        return 0;
    }
    int dfs(int x,int mx,int t)
    {
        if(x==t)return mx;
        int flow=0;
        for(int i=head[x];~i;i=e[i].Next)
        {
            int te=e[i].to,f;
            if(e[i].c>0&&dis[te]==dis[x]+1&&(f=dfs(te,min(mx-flow,e[i].c),t)))
            {
                e[i].c-=f;
                e[i^1].c+=f;
                flow+=f;
            }
        }
        if(!flow)dis[x]=-2;
        return flow;
    }
    int maxflow(int s,int t)
    {
        int ans=0,f;
        while(bfs(s,t))
        {
            while((f=dfs(s,inf,t)))ans+=f;
        }
        return ans;
    }
    int main()
    {
        int n,m;
        while(~scanf("%d%d",&n,&m))
        {
            init();
            int s=n+m+1,t=n+m+2;//原图
            int ss=n+m+3,tt=n+m+4;//超级源超级汇
            for(int i=1;i<=m;i++)
            {
                int a;scanf("%d",&a);
                add(n+i,t,inf-a,a);
            }
            int res=0;
            for(int i=1;i<=n;i++)
            {
                int k,d;scanf("%d%d",&k,&d);
                add(s,i,d,0);
                while(k--)
                {
                    int a,b,c;scanf("%d%d%d",&a,&b,&c);
                    ansid[++res]=cnt;
                    add(i,n+a+1,c-b,b);
                }
            }
            add(t,s,inf,0);
            int sum=0;
            for(int i=1;i<=n+m+2;i++)
            {
                if(in[i]>out[i])sum+=in[i]-out[i],add(ss,i,in[i]-out[i],0);
                else add(i,tt,out[i]-in[i],0);
            }
            if(sum!=maxflow(ss,tt))puts("-1");
            else
            {
                printf("%d
    ",maxflow(s,t));
                for(int i=1;i<=res;i++)
                    printf("%d
    ",e[ansid[i]^1].c+e[ansid[i]].low);
            }
            puts("");
        }
        return 0;
    }
    /*******************
    
    ********************/
    View Code
  • 相关阅读:
    CF 552(div 3) E Two Teams 线段树,模拟链表
    单词接龙
    书的复制 动规,贪心
    C++文本处理&造轮子
    【NOIP2015】【Luogu2661】信息传递(有向图最小环)
    平时上机练习的注意点(NOIP2018)
    【SHOI2009】【BZOJ2028】会场预约(线段树染色)
    【AHOI2009】【BZOJ1798】Seq 维护序列seq(线段树模板,易错提醒)
    2018 “百度之星”程序设计大赛
    2018 “百度之星”程序设计大赛
  • 原文地址:https://www.cnblogs.com/acjiumeng/p/7805893.html
Copyright © 2020-2023  润新知