多源多汇费用流——poj2516

网络流的题好难。。感觉有点遭不住了

这题用矩阵存图，然后把k个物品，每个物品都求一次费用流

/*
多源多汇的费用流
其实是k个费用流
*/
#include<bits/stdc++.h>
using namespace std;
#define inf 0x3f3f3f3f
#define maxn 205
struct Edge{int to,nxt,w,c;}e[maxn<<];
int head[maxn],tot,n,m,s,t;
void init(){}
void add(int u,int v,int w,int c){}

int sum_need[maxn],sum_offer[maxn];
int offer[maxn][maxn],need[maxn][maxn],mp[maxn][maxn],cost[maxn][maxn],K,sum,pre[maxn];

int dis[maxn],vis[maxn];
int spfa(){
    queue<int>q;
    for(int i=;i<=t;i++){vis[i]=false;dis[i]=inf;}

    vis[s]=;dis[s]=;q.push(s);
    while(!q.empty()){
        int k=q.front();q.pop();vis[k]=false;
        for(int i=;i<=t;i++)
            if(mp[k][i] && dis[i]>dis[k]+cost[k][i]){
                dis[i]=dis[k]+cost[k][i];
                pre[i]=k;
                if(!vis[i]){
                    vis[i]=;
                    q.push(i);
                }
            }
    }

    if(dis[t]!=inf)return ;
    return ;
}
int fond(){
    int Min=inf,res=;
    while(spfa()){
        for(int i=t;i!=s;i=pre[i])
            Min=min(Min,mp[pre[i]][i]);
        for(int i=t;i!=s;i=pre[i]){
            mp[pre[i]][i]-=Min;
            mp[i][pre[i]]+=Min;
            res+=cost[pre[i]][i]*Min;
        }
    }
    return res;
}

int main(){
    while(cin>>n>>m>>K,n&&m&&K){
        sum=;
        memset(sum_need,,sizeof sum_need);
        memset(sum_offer,,sizeof sum_offer);

        for(int i=;i<=n;i++)
            for(int j=;j<=K;j++){
                scanf("%d",&need[i][j]);
                sum_need[j]+=need[i][j];
            }
        for(int i=;i<=m;i++)
            for(int j=;j<=K;j++){
                scanf("%d",&offer[i][j]);
                sum_offer[j]+=offer[i][j];
            }
        int sign=;
        for(int i=;i<=K;i++)
            if(sum_offer[i]<sum_need[i]){
                sign=;
                break;
            }
        s=;
        t=n+m+;
        for(int k=;k<=K;k++){
            memset(mp,,sizeof mp);
            memset(cost,,sizeof cost);
            //第k件物品 供应商->店主的运费
            for(int i=+m;i<=n+m;i++)
                for(int j=;j<=m;j++){
                    scanf("%d",&cost[j][i]);
                    cost[i][j]-=cost[j][i];//反向边的费用
                }

            if(sign==)continue;

            for(int i=;i<=m;i++)//s->供应商
                mp[s][i]=offer[i][k];
            for(int i=;i<=m;i++)//供应商->店主
                for(int j=m+;j<=m+n;j++)
                    mp[i][j]=offer[i][k];
             for(int i=m+;i<=m+n;i++)//店主->t
                 mp[i][t]=need[i-m][k];
            sum+=fond();
        }
        if(sign==)
            printf("-1\n");
        else printf("%d\n",sum);
    }
}