poj 2516 （最小费用最大流）

题目链接：哆啦A梦传送门

题意：题意非常难看，直接看这个题意

有n个店铺，m个供货商，k种商品。给出每个供货商的仓库里每种商品的数量、每种商品给每个店铺供货的费用，每个店铺需要的各种商品的数量，求最小费用。 

 

输入一开始是n，m，k。

然后 n行   ，每行 k列 ，第i行第j个数代表 第i个店铺 需要第j种商品多少件。

然后 m行  ，每行 k列 ， 第i行第j个数代表 第i个供货商 有第j种物品多少件。

然后是k个矩阵  ，每个n行m列，第i个矩阵的第j行第d列代表着第d个供货商给第j个店铺发第i种货物一件需要的花费。

 

最小费用最大流参考博客：

博客1

博客2

 

我们分别给k个商品建立模型求最小费用。

 

求k次最下费用最大流（在满足最短路的情况下求最大流）就行了。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>

using namespace std;

const int N=200;
const int INF=0x3f3f3f3f;

struct edge{
    int u,v,flow,cost,next;
    edge(){}
    edge(int _u,int _v,int _flow,int _cost,int _next){
        u=_u;v=_v;flow=_flow;cost=_cost;next=_next;
    }
}e[200000];

int head[N],cnt;
queue<int> que;

void init()
{
    memset(head,-1,sizeof(head));
    cnt=0;
}

void add(int x,int y,int z,int c)
{
    e[cnt]=edge(x,y,z,c,head[x]);
    head[x]=cnt++;
    e[cnt]=edge(y,x,0,-c,head[y]);
    head[y]=cnt++;
}

bool vis[N];
int dis[N],pre[N];

bool spfa(int s,int t)
{
    memset(dis,INF,sizeof(dis));
    memset(vis,0,sizeof(vis));
    memset(pre,-1,sizeof(pre));

    dis[s]=0;
    vis[s]=true;
    while(!que.empty()) que.pop();

    que.push(s);

    while(!que.empty())
    {
        int x=que.front();
        que.pop();
        vis[x]=false;

        for(int i=head[x];i!=-1;i=e[i].next){
            if(e[i].flow){
                int y=e[i].v;
                if(dis[x]+e[i].cost<dis[y]){
                    dis[y]=dis[x]+e[i].cost;
                    pre[y]=i;
                    if(!vis[y]){
                        que.push(y);
                        vis[y]=true;
                    }
                }
            }
        }

    }
     return pre[t]!=-1;

}

int costflow(int s,int t,int x)
{
    int ret=0;
    int max_flow=0;

    while(spfa(s,t))
    {
        int min_flow=INF;
        for(int i=pre[t];i!=-1;i=pre[e[i].u]){
            min_flow=min(e[i].flow,min_flow);
        }
        for(int i=pre[t];i!=-1;i=pre[e[i].u]){
            e[i].flow-=min_flow;
            e[i^1].flow+=min_flow;
        }
        ret+=dis[t]*min_flow;
        max_flow+=min_flow;
    }

    if(max_flow<x) return -1;///达不到满流
    return ret;
}

int need[100][100],supply[100][100];
int  main()
{

    int n,m,k;

    while(~scanf("%d%d%d",&n,&m,&k))
    {
        if(n==0&&m==0&&k==0) break;
        for(int i=1;i<=n;i++)
            for(int j=1;j<=k;j++)
                scanf("%d",&need[i][j]);

        for(int i=1;i<=m;i++)
            for(int j=1;j<=k;j++)
            scanf("%d",&supply[i][j]);

        int flag=1;
        int sum=0;///最小花费

        for(int x=1;x<=k;x++)///分别给每种商品计算一次最小费用最大流
        {
            init();
            int s=0,t=n+m+1;
            int ans=0;
            for(int i=1;i<=n;i++)///店铺i+m与汇点建边，权值为该店铺i+m所需此商品的需求量
            {
                add(i+m,t,need[i][x],0);
                ans+=need[i][x];
            }

            for(int i=1;i<=m;i++) ///源点向供应商i建边，权值为供应商i对此商品的存货量
                add(s,i,supply[i][x],0);

                int in;
            for(int i=1;i<=n;i++)
            {
                for(int j=1;j<=m;j++){
                    scanf("%d",&in);
                    add(j,i+m,supply[j][x],in);///供应商j与店铺i+m建边，权值为该供应商j供货给该店铺i+m的商品x所需的费用
                }
            }

            int item=costflow(s,t,ans);
//            printf("item=%d sum=%d\n",item,sum);
            if(item==-1) flag=0;
            else sum+=item;
        }
        if(!flag) puts("-1");
        else printf("%d\n",sum);
    }
    return 0;
}