POJ 2516 Minimum Cost (KM最优匹配)

题意：有N家家店，每家店都对K种货物有需求；同时有M家仓库，对K钟货物有供应。对于每种货物，每个仓库送至每家店都有自己的单位费用。求满足所有店所有货物的最小费用

分析：对于每一种货物，如果总需求大于总费用，那么无解的；否则可以用KM匹配求解每一种货物的最小费用。

因为KM求的是最优匹配，所以建图的时候，建负权边即可。

将每一份需求和每一份供应都离散化，而不是把一个店的需求和每个一仓库的供应建立边。

#include<stdio.h>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
const int maxn =;
const int INF= ;

int w[maxn][maxn];
int m,n;//n左m右
int cx[maxn],cy[maxn];//顶标
bool usex[maxn],usey[maxn];//本回合使用的x，y
int link[maxn];//link[i]=x代表：在y图中的i与x相连 

int A,B;

bool dfs(int u){
    usex[u]=;
    for(int i=;i<=B;i++)
        if(!usey[i]&&cx[u]+cy[i]==w[u][i]){
            usey[i]=;
            if(link[i]==-||dfs(link[i])){
                link[i]=u;
                return ;
            }
        }
    return ;
}
int KM(){
    memset(cy,,sizeof(cy));
    memset(cx,-,sizeof(cx));
    memset(link,-,sizeof(link));
    for(int i=;i<=A;i++)
        for(int j=;j<=B;j++)
            cx[i]=max(cx[i],w[i][j]);
    for(int i=;i<=A;i++){
        while(){
            int d=INF;
            memset(usex,,sizeof(usex));
            memset(usey,,sizeof(usey));
            if(dfs(i))break;
            for(int i=;i<=A;i++)
                if(usex[i])
                    for(int j=;j<=B;j++)
                        if(!usey[j])d=min(d,cx[i]+cy[j]-w[i][j]);
            if(d==INF)return -;
            for(int i=;i<=A;i++)
                if(usex[i])cx[i]-=d;
            for(int i=;i<=B;i++)
                if(usey[i])cy[i]+=d;
        }
    }

    int ans=;
    for(int i=;i<=B;i++){
        if(~link[i]){
            ans -= w[link[i]][i];
        }
    }
    return ans;
}

int nd[maxn][maxn];
int sp[maxn][maxn];
int ct[maxn][maxn][maxn];
int cas;

int cnt1[maxn],cnt2[maxn];

int main()
{
    #ifndef ONLINE_JUDGE
        freopen("in.txt","r",stdin);
        freopen("out.txt","w",stdout);
    #endif
    int T,tmp;
    int u,v,k;
    while(scanf("%d%d%d",&n,&m,&k)==){
        if(!n && !m && !k) break;

        for(int i=;i<=n;++i){
            for(int j=;j<=k;++j){
                scanf("%d",&nd[i][j]);
            }
        }

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

        for(int i=;i<=k;++i){
            for(int j=;j<=n;++j){
                for(int t = ;t<=m;++t){
                    scanf("%d",&ct[i][j][t]);
                }
            }
        }

        bool flag = true;
        for(int cas=;cas<=k;++cas){
            int need = ,supply =;
            for(int i=;i<=n;++i)   need += nd[i][cas];
            for(int i=;i<=m;++i)   supply += sp[i][cas];
            if(need>supply){
                flag = false;
                break;
            }
        }

        if(!flag){
            puts("-1");
            continue;
        }

        int tot = ;
        for(int cas = ;cas<=k;++cas){
            A=B=;
            for(int i=;i<=n;++i){
                for(int j=;j<=nd[i][cas];++j){
                    cnt1[++A] = i;
                }
            }
            for(int i=;i<=m;++i){
                for(int j=;j<=sp[i][cas];++j){
                    cnt2[++B] = i;
                }
            }
            for(int i=;i<=A;++i){
                for(int j=;j<=B;++j){
                    w[i][j] = -ct[cas][cnt1[i]][cnt2[j]];
                }
            }
            int tmp = KM();
            tot+=tmp;
        }
        printf("%d\n",tot);
    }
    return ;
}