E - Minimum Cost（费用流）

题意：有N个供应商，M个店主，K种物品。每个供应商对每种物品的的供应量已知，每个店主对每种物品的需求量的已知，从不同的供应商运送不同的货物到不同的店主手上需要不同的花费，又已知从供应商Mj送第kind种货物的单位数量到店主Ni手上所需的单位花费。

问：供应是否满足需求？如果满足，最小运费是多少？

解决方法是用费用流。

个人觉得网络流都是比较暴力的算法，就是两个收索结合，都是比较暴力的。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#include <iostream>
#define INF 0x3f3f3f3f
using namespace std;

//建图所用
const int maxn=1e4+100;
const int maxe=1e4+1000;
const int maxv=1e4+100;
struct edge
{
    int to,w,next,c;//w是费用，c是残量网络的流量
}e[maxe<<1];
int head[maxv<<1],depth[maxv],cnt;
void init()
{
    memset(head,-1,sizeof(head));
    cnt=-1;
}
void add_edge(int u,int v,int w,int c)
{
    e[++cnt].to=v;
    e[cnt].w=w;
    e[cnt].c=c;
    e[cnt].next=head[u];
    head[u]=cnt;
}
void _add(int u,int v,int w,int c)
{
    add_edge(u,v,w,c);
    add_edge(v,u,-w,0);
}

//费用流算法
int S,T;
int vis[maxn],dis[maxn],pre_v[maxn],pre_e[maxn];
bool spfa()//类似于EK算法，通过不停的找增广路,即找一条距离最小的路径
{
    memset(vis,0,sizeof(vis));
    queue<int> que;
    memset(dis,0x3f,sizeof(dis));//初始距离为无穷大
    dis[S]=0;
    que.push(S);
    vis[S]=1;

    while(!que.empty())//bfs
    {
        int u=que.front();que.pop();
        vis[u]=0;
        for(int i=head[u];i!=-1;i=e[i].next)
        {
            int v=e[i].to;
            if(e[i].c>0 && dis[v]>dis[u]+e[i].w)//w表示费用，c表示残量网络的流量
            {
                dis[v]=dis[u]+e[i].w;//距离指标是费用
                pre_e[v]=i;//记录前驱边，即记录增广路
                pre_v[v]=u;//记录前驱点
                if(!vis[v])//bfs
                {
                    vis[v]=1;
                    que.push(v);
                }
            }
        }
    }
    if(dis[T]<INF) return true;
    return false;
}
int maxflow=0;
int augment()//调整残量网络并返回增广路的费用
{
    int flow=INF;//flow最终为当前增广路的流量
    int u=T;
    while(u!=S)//从汇点到源点逆着bfs中得到的路径走
    {
        flow=min(flow,e[pre_e[u]].c);//最大流是增广路流量的最小值
        u=pre_v[u];
    }
    u=T;
    while(u!=S)
    {
        e[pre_e[u]].c-=flow;//将增广路上的流量都减去最大流
        e[pre_e[u]^1].c+=flow;
        u=pre_v[u];
    }
    maxflow+=flow;
    return flow*dis[T];//dis[T]表示T到汇点的距离，距离指标是费用
}
int cost()
{
    int ans=0;
    while(spfa())//不停的找增广路
    {
        ans+=augment();//累加增广路的费用
    }
    return ans;
}
int N,M,K;
int a[100][100],b[100][100];
int main()
{
    while(~scanf("%d%d%d",&N,&M,&K))
    {

        if(N==0 && M==0 && K==0) break;
        S=0,T=200;
        int ans=0,sum=0,tot=0;
        for(int i=1;i<=N;i++)//第i号店进第j种物品的数量
            for(int j=1;j<=K;j++)
            {
                scanf("%d",&a[i][j]);
                sum+=a[i][j];
            }

        for(int i=1;i<=M;i++)//第i号仓库储备第j种物品的数量
            for(int j=1;j<=K;j++)
            {
                scanf("%d",&b[i][j]);
            }

        for(int k=1;k<=K;k++)//将一个非常稠密大图分成k个较稀疏的小图以减少时间
        {
            init();
            maxflow=0;
            for(int i=1;i<=N;i++) _add(S,i,0,a[i][k]);
            for(int i=1;i<=M;i++) _add(i+N,T,0,b[i][k]);
            for(int i=1;i<=N;i++)
              for(int j=1;j<=M;j++)
            {
                int x;
                scanf("%d",&x);
                _add(i,j+N,x,INF);
            }
            ans+=cost();//累加费用
            tot+=maxflow;//累加最大流
        }
        if(tot==sum)
            printf("%d\n",ans);
        else
            printf("-1\n");
    }
    return 0;
}