poj2516----最小费用流（模版+构图）

【题意】

有n个商店m个供应商k种商品

给各商店的各种商品的需求量和各供应商的供应量以及供应商和商店之间的运输费用，求满足需求的最小费用

【输入】

多组数据

每组数据第一行n,m,k（均小于50）

接下来n行每行k个数，表示需求量（0到3）

接下来m行每行k个数，表示供应量（0到3）

接下来k个n*m矩阵表示商店和供应商之间的运输费用

输入由第一行三个0结束

【输出】

每组数据输出一个数，表示最小费用

//在多重循环中一定注意下标容易出错，i,j,k!!!

#include<cstdlib>
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include <algorithm>
#include<queue>
#include<set>
#define LL long long
#define inf 0x7fffffff
#define E 1e-9
#define M 100
#define N 105
using namespace std;
int nead[N][N],have[N][N];


int n,m,k;
int cost[N][N];
int flow[N][N],cap[N][N];
int d[N],p[N];
bool inq[N];
int s,t;
int num;

int getflow(int s,int t)
{
  queue<int> q;
  memset(flow,0,sizeof(flow));
  int c=0,f=0;
  for(;;)
    {
      memset(d,0x7f,sizeof(d));//
      memset(inq,0,sizeof(inq));
      d[s]=0;
      q.push(s);

      while(!q.empty())
        {
          int u=q.front();
          q.pop();
          inq[u]=0;
          for(int i=1; i<=num; i++)
            if(cap[u][i]>flow[u][i]&&d[i]>d[u]+cost[u][i])
              {
                d[i]=d[u]+cost[u][i];
                p[i]=u;
                if(!inq[i])
                  {
                    q.push(i);
                    inq[i]=1;
                  }
              }
        }


      if(d[t]==0x7f7f7f7f)0x7f7f7f7f
        break;
      int a=inf;
      for(int i=t; i!=s; i=p[i])
        {
          a=min(a,cap[p[i]][i]-flow[p[i]][i]);
        }

      for(int i=t; i!=s; i=p[i])
        {
          flow[p[i]][i]+=a;
          flow[i][p[i]]-=a;
        }
      f+=a;
      c+=d[t]*a;
    }
  return c;
}
int main()
{
#ifndef ONLINE_JUDGE
  freopen("ex.in","r",stdin);
#endif
  while(scanf("%d%d%d",&n,&m,&k))
    {
      int ans=0;
      if(!n&&!m&&!k)
        break;
      for(int i=1; i<=n; i++)
        for(int j=1; j<=k; j++)
          scanf("%d",&nead[i][j]);
      for(int i=1; i<=m; i++)
        for(int j=1; j<=k; j++)
          scanf("%d",&have[i][j]);


      int flag=0;
      for(int i=1; i<=k; i++)
        {
          int nead1=0;
          for(int j=1; j<=n; j++)
            nead1+=nead[j][i];

          int have1=0;
          for(int j=1; j<=m; j++)
            have1+=have[j][i];
          if(have1<nead1)
            {
              flag=1;
              break;
            }
        }
      for(int p=1; p<=k; p++)//在多重循环中一定注意下标容易出错，i,j,k!!!
        {
          memset(cost,0,sizeof(cost));
          memset(cap,0,sizeof(cap));
          for(int j=1; j<=n; j++)
            for(int k=1; k<=m; k++)
              {
                scanf("%d",&cost[j][k+n]);
                cost[k+n][j]=-cost[j][k+n];
                cap[j][k+n]=nead[j][p];
//                cap[k+n][j]=0;
              }
          if(flag)
            continue;

          for(int j=1; j<=n; j++)
            {
              cap[0][j]=nead[j][p];
            }
          num=n+m+1;
          for(int j=1; j<=m; j++)
            {
              cap[j+n][num]=have[j][p];
            }
          s=0;
          t=num;

          ans+=getflow(s,t);
        }
      if(flag)
        printf("-1\n");
      else
        printf("%d\n",ans);

    }
  return 0;
}