【网络流24题】运输问题

（网络流24题大多需要spj，所以需要一个有spj的oj，本系列代码均在www.oj.swust.edu.cn测试通过）
非常裸的一道费用流，只不过要跑一遍最大流一遍最小流，所以这里就不赘述了，看代码吧。

#include<cstdio>
#include<cstdlib>
#include<ctime>
#include<cmath>
#include<cstring>
#include<string>
#include<iostream>
#include<iomanip>
#include<algorithm>
using namespace std;
#define INF 100000000
struct bian
{
    int l,r,f,v;
}a[1000000];
int fir[1000000];
int nex[1000000];
int tot=1;
void add_edge(int l,int r,int f,int v)
{
    a[++tot].l=l;
    a[tot].r=r;
    a[tot].f=f;
    a[tot].v=v;
    nex[tot]=fir[l];
    fir[l]=tot;
}
int n,m;
int v1[1001];
int v2[1001];
int val[1001][1001];
int s=0,t=9999;
void jianbian()
{
    memset(fir,0,sizeof(fir));
    memset(nex,0,sizeof(nex));
    tot=1;
    for(int i=1;i<=n;i++)
    {
        add_edge(s,i,v1[i],0);
        add_edge(i,s,0,0);
    }
    for(int i=1;i<=m;i++)
    {
        add_edge(n+i,t,v2[i],0);
        add_edge(t,n+i,0,0);
    }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            add_edge(i,n+j,INF,val[i][j]);
            add_edge(n+j,i,0,-val[i][j]);
        }
}
int dis[10000];
bool pd[10000];
int fro[10000];
bool spfamin()
{
    static int dui[1000000];
    memset(dis,0x1f,sizeof(dis));
    int top=1,my_final=2;
    dui[top]=s;
    dis[s]=0;
    pd[s]=true;
    while(top<my_final)
    {
        int u=dui[top++];
        for(int o=fir[u];o;o=nex[o])
        {
            if(dis[a[o].r]>dis[u]+a[o].v && a[o].f)
            {
                dis[a[o].r]=dis[u]+a[o].v;
                fro[a[o].r]=o;
                if(!pd[a[o].r]) pd[a[o].r]=true,dui[my_final++]=a[o].r;
            }
        }
        pd[u]=false;
    }
    if(dis[t]==0x1f1f1f1f) return false;
    return true;
}
bool spfamax()
{
    static int dui[1000000];
    memset(dis,0x8f,sizeof(dis));
    int top=1,my_final=2;
    dui[top]=s;
    dis[s]=0;
    pd[s]=true;
    while(top<my_final)
    {
        int u=dui[top++];
        for(int o=fir[u];o;o=nex[o])
        {
            if(dis[a[o].r]<dis[u]+a[o].v && a[o].f)
            {
                dis[a[o].r]=dis[u]+a[o].v;
                fro[a[o].r]=o;
                if(!pd[a[o].r]) pd[a[o].r]=true,dui[my_final++]=a[o].r;
            }
        }
        pd[u]=false;
    }
    if(dis[t]==0x8f8f8f8f) return false;
    return true;
}
int cost;
void add_flow()
{
    int mid=t;
    int temp=2147483647;
    while(mid!=s)
    {
        temp=min(temp,a[fro[mid]].f);
        mid=a[fro[mid]^1].r;
    }
    cost+=temp*dis[t];
    mid=t;
    while(mid!=s)
    {
        a[fro[mid]].f-=temp;
        a[fro[mid]^1].f+=temp;
        mid=a[fro[mid]^1].r;
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&v1[i]);
    for(int i=1;i<=m;i++) scanf("%d",&v2[i]);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            scanf("%d",&val[i][j]);
    jianbian();
    while(spfamin()) add_flow();
    cout<<cost<<endl;
    jianbian();
    cost=0;
    while(spfamax()) add_flow();
    cout<<cost;
    return 0;
}