最小费用流

原创 2012年03月25日 22:20:23

#include <cstdio>
#include <cstring>
#define maxn 61
#define maxv (maxn*maxn*2 + 1)
#define maxe (maxv*5)
#define oo 2147483647
#define min(a,b) ((a>(b)?(b):(a)))
#define maxq maxe

using namespace std;

struct edge_type
{
    int x,y,flow,cost,next,op;
};

int n,m;
int a[maxn + 1][maxn + 1];

int total,totalv,source,sink;
int first[maxv + 1];
edge_type g[maxe + 1];

int dist[maxv + 1],fa[maxv + 1];
bool mark[maxv + 1];
int q[maxq + 1];

bool spfa()
{
    memset(mark,0,sizeof(mark));
    memset(fa,0,sizeof(fa)); //
    for (int i=0;i<=totalv;i++)
        dist[i] = -1;
    dist[source] = 0;
    int head,tail;
    head = 1; tail = 2;
    q[1] = source;
    while (head!=tail){
        int b = q[head];
        mark[b] = 0;
        for (int temp=first[b];temp;temp=g[temp].next){
            if (!g[temp].flow) continue;
            int y = g[temp].y;
            int tt;
            if ((tt = dist[b] + g[temp].cost)>dist[y]){
                dist[y] = tt;
                fa[y] = temp;
                if (!mark[y]){
                    mark[y] = 1;
                    q[tail ++] = y;
                    if (tail>maxq) tail = 1;
                }
            }
        }
        head ++;
        if (head>maxq) head = 1;
    }
    return dist[sink]>=0;
}

int mincost_maxflow()
{
    int total_flow,total_cost,cur_flow,cur_cost;
    total_flow = total_cost = 0;
    while (spfa()){/*
        for (int i=1;i<=totalv;i++){
            printf("%d ",dist[i]);
        }
        printf("\n");*/
        cur_flow = oo; cur_cost = 0;
        for (int temp=fa[sink];temp;temp=fa[g[temp].x]){
            cur_flow = min(cur_flow,g[temp].flow);
            cur_cost += g[temp].cost;
        }
        for (int temp=fa[sink];temp;temp=fa[g[temp].x]){
            g[temp].flow -= cur_flow;
            g[g[temp].op].flow += cur_flow;
        }
//        printf("%d\n",cur_flow);
        total_flow += cur_flow;
        total_cost += cur_cost*cur_flow;
    }
//    printf("%d\n",total_flow);
    return total_cost;
}

void add(int x,int y,int flow,int cost)
{
    g[++total].x = x; g[total].y = y;
    g[total].flow = flow; g[total].cost = cost;
    g[total].op = total + 1;
    g[total].next = first[x];
    first[x] = total;

    g[++total].x = y; g[total].y = x;
    g[total].flow = 0; g[total].cost = -cost;
    g[total].op = total - 1;
    g[total].next = first[y];
    first[y] = total;
}

int main()
{
    while (scanf("%d%d",&n,&m)==2){
        for (int i=1;i<=n;i++){
            for (int j=1;j<=n;j++)
                scanf("%d",&a[i][j]);
        }
        total = 0;
        memset(first,0,sizeof(first));
        source = n*n*2 + 1;
        totalv = source;
        sink = n*n*2;
        add(source,1,m,0);
        for (int i=1;i<=n;i++){
            for (int j=1;j<=n;j++){
                add((i - 1)*n + j,(i - 1)*n + j + n*n,1,a[i][j]);
                add((i - 1)*n + j,(i - 1)*n + j + n*n,oo,0);
                if (i<n)
                    add((i - 1)*n + j + n*n,i*n + j,oo,0);
                if (j<n)
                    add((i - 1)*n + j + n*n,(i - 1)*n + j + 1,oo,0);
            }
        }

        int ans = mincost_maxflow();
        printf("%d\n",ans);
    }
    return 0;
}



最大费用最大流和最小费用最大流一样可以用SPFA + Edmonds-Karp解决。//显然SPFA可以求最长路的

貌似先用spfa找出一条增光路,再用ek增加流

有时间再研究研究

http://rekonq.appspot.com/%3Fp%3D119001


最小费用流---ZKW算法

网络流的一些基本概念 很多同学建立过网络流模型做题目, 也学过了各种算法, 但是对于基本的概念反而说不清楚. 虽然不同的模型在具体叫法上可能不相同, 但是不同叫法对应的思想是一致的. 下面的讨论力求...
  • little_nai
  • little_nai
  • 2017年04月04日 19:36
  • 3206

使用dijkstra求解最小费用最大流网络

前言在介绍如何使用dijkstra算法求解最小费用最大流问题的时候,假设看这篇博文的读者已经知道什么是最小费用最大流问题及熟悉dijkstra单源最短路径算法。在这篇博文里面,我并不会过多强调网络拓扑...
  • flydreamforever
  • flydreamforever
  • 2017年04月08日 23:03
  • 938

最小费用最大流算法及题集

貌似网上最小费用最大流的讲解的不多。
  • y990041769
  • y990041769
  • 2014年10月14日 14:18
  • 11285

最小费用流算法不完全指南-2017华为软件精英挑战赛

前言整个比赛历时70多天,是一场持久战,期间看了不少资料,实现了一些算法,现总结如下。初赛的思路是启发式+费用流,启发式实现了遗传算法、模拟退火,费用流实现了连续最短路(spfa),zkw/prima...
  • shanchuan2012
  • shanchuan2012
  • 2017年05月17日 17:32
  • 2096

bzoj 4108: [Wf2015]Catering|带上下界最小费用可行流

终于搞懂 上下界网络流了! 注意边的数量 #include #include #include #include using namespace std; const in...
  • sxb_201
  • sxb_201
  • 2016年04月25日 21:35
  • 805

最小费用流bellman-ford算法总结

最小费用流问题是寻找流为F时的最小费用,普通的最大流中边只有流量限制,而在费用流中还加上了费用,为保证过程可逆和正确性残余网络中的负权边是原边费用的相反数,求解此类问题的方法依然是贪心,我们在残余网络...
  • liuzhan214
  • liuzhan214
  • 2016年07月30日 10:59
  • 2125

网络流 ( 最小费用最大流 )——Going Home ( HDU 1533 )

题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=1533 分析: 给出一个网格图,H代表房子,m代表人,两者数量相等,人每移动一格花费1$,求所有人...
  • FeBr2
  • FeBr2
  • 2016年07月27日 21:33
  • 478

网络流问题-最小费用最大流

继续沿着2017华为软件精英挑战赛,查阅了一些【最小费用最大流】相关的资料。今年的赛题归结为:组合优化+最小费用最大流两个子问题,给定的服务器选址下求出最小费用最大流作为底层的算法支撑会很好提升最终的...
  • Tcorpion
  • Tcorpion
  • 2017年03月22日 16:02
  • 1703

网络流(最大流和最小费用流)

几个月前学过,然而一下就忘记了,于是决定系统的复习一下。 关于网络流各路神犇早已有很好的讲解,于是我就整理一下(其实我是蒟蒻,看到的果断关掉吧) http://blog.csdn.net/leol...
  • thy0311
  • thy0311
  • 2015年07月20日 10:19
  • 2430

网络流-最小费用最大流

(在最小费用最大流中,最大流量是唯一的,但是最大流不唯一,在保证最大流的条件下,加上一些参数,确定最小费用最大流的唯一性) 最小费用流其实是线性规划的一种特殊类型。所以解决最小费用流的方式其实可以为...
  • sentimental_dog
  • sentimental_dog
  • 2016年07月12日 15:40
  • 752
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:最小费用流
举报原因:
原因补充:

(最多只允许输入30个字)