zoukankan      html  css  js  c++  java
  • HDU 3667

    http://acm.hdu.edu.cn/showproblem.php?pid=3667

    最小费用最大流

    本题流量和费用不是线性关系,fee=a*flow*flow,所以常规套模板spfa无法得到最小费用

    观察到每条边流量上限只有5,则可以把一条流量为f的边拆成f条流量为1的边,每条边费用是a*(2*i-1)(1<=i<=f)

    #include <iostream>
    #include <cstring>
    #include <cstdio>
    #include <queue>
    using namespace std ; 
    const int INF=0xfffffff ;
    struct node{
        int s,t,cap,cost,nxt ;
    }e[200005] ;
    int sumflow ;
    int n,m,k,cnt,head[1005],vis[1005],dis[1005],pre[1005] ;
    void add(int s,int t,int cap,int cost)
    {
        e[cnt].s=s ;e[cnt].t=t ;e[cnt].cap=cap ;e[cnt].cost=cost ;e[cnt].nxt=head[s] ;head[s]=cnt++ ;
        e[cnt].s=t ;e[cnt].t=s ;e[cnt].cap=0 ;e[cnt].cost=-cost ;e[cnt].nxt=head[t] ;head[t]=cnt++ ;
    }
    int spfa(int s,int t,int N)
    {
        for(int i=0 ;i<=N ;i++)
            dis[i]=INF ;
        dis[s]=0 ;
        memset(vis,0,sizeof(vis)) ;
        memset(pre,-1,sizeof(pre)) ;
        vis[s]=1 ;
        queue <int> q ;
        q.push(s) ;
        while(!q.empty())
        {
            int u=q.front() ;
            q.pop() ;
            vis[u]=0 ;
            for(int i=head[u] ;i!=-1 ;i=e[i].nxt)
            {
                int tt=e[i].t ;
                if(e[i].cap && dis[tt]>dis[u]+e[i].cost)
                {
                    dis[tt]=dis[u]+e[i].cost ;
                    pre[tt]=i ;
                    if(!vis[tt])
                    {
                        vis[tt]=1 ;
                        q.push(tt) ;
                    }
                }
            }
        }
        if(dis[t]==INF)return 0 ;
        return 1 ;
    }
    int MCMF(int s,int t,int N)
    {
        int flow,minflow,mincost ;
        mincost=flow=0 ;
        while(spfa(s,t,N))                                                
        {
            minflow=INF ;
            for(int i=pre[t] ;i!=-1 ;i=pre[e[i].s])
                minflow=min(minflow,e[i].cap) ;
            flow+=minflow ;
            for(int i=pre[t] ;i!=-1 ;i=pre[e[i].s])
            {
                e[i].cap-=minflow ;
                e[i^1].cap+=minflow ;
            }
            mincost+=dis[t]*minflow ;
        }
        sumflow=flow ;//最大流 
        if(sumflow<k)return -1 ;
        return mincost ;
    }
    int main()
    {
        while(~scanf("%d%d%d",&n,&m,&k))
        {
            cnt=0 ;
            memset(head,-1,sizeof(head)) ;
            while(m--)
            {
                int u,v,a,c ;
                scanf("%d%d%d%d",&u,&v,&a,&c) ;
                for(int i=1 ;i<=c ;i++)
                    add(u,v,1,a*(2*i-1)) ;//Σa*(2*i-1)=a*c*c 
            }
            add(0,1,k,0) ;
            printf("%d
    ",MCMF(0,n,n+1)) ;
        }
        return 0 ;
    }
    View Code
  • 相关阅读:
    mysql笔记
    ssh学习笔记
    oracle数据向历史表数据迁移————procedure
    关于避免模糊查询索引时效问题
    css圣杯布局和双飞翼布局篇
    什么是hasLayout
    高度自适应的水平垂直居中布局
    MAC的GIF动图录屏软件LICECAP
    sublime text 快捷键shortcuts
    为表格合并边框的样式
  • 原文地址:https://www.cnblogs.com/xiaohongmao/p/3714633.html
Copyright © 2011-2022 走看看