网络流24题之运输问题

其实是水题只不过在noip前一个月时我问整个机房也没人能回答上来。

现在可以随随便便切过去了也是很爽的。

最小费用最大流，按题目要求建图即可。

#include<bits/stdc++.h>
using namespace std;
const int N=250005,inf=1e9;
int d[N],head[N],cnt=-1,n,m,w[105][105],cost,f[N],a[N],b[N],s,t;
bool v[N];
struct node{
    int w,to,nex,c,f;
}e[1000005];
void add(int x,int y,int w,int c)
{
    e[++cnt].to=y;e[cnt].nex=head[x];head[x]=cnt;e[cnt].f=x;e[cnt].w=w;e[cnt].c=c;
    e[++cnt].to=x;e[cnt].nex=head[y];head[y]=cnt;e[cnt].f=y;e[cnt].w=0;e[cnt].c=-c;
}
queue<int>q;
bool spfa()
{
    memset(d,0x3f,sizeof(d));
    memset(v,0,sizeof(v));
    memset(f,-1,sizeof(f));
    v[s]=1;d[s]=0;q.push(s);
    while(!q.empty())
    {
        int x=q.front();q.pop();v[x]=0;
        for(int i=head[x];i!=-1;i=e[i].nex)
        {
            if(!e[i].w||d[e[i].to]<=d[x]+e[i].c)continue;
            int y=e[i].to;
            d[y]=d[x]+e[i].c;f[y]=i;
            if(!v[y])
            {
                q.push(y);v[y]=1;
            }
        }
    }
    if(d[t]>inf)return 0;
    int flow=inf;
    for(int i=f[t];i!=-1;i=f[e[i].f])
    flow=min(flow,e[i].w);
    for(int i=f[t];i!=-1;i=f[e[i].f])
    e[i].w-=flow,e[i^1].w+=flow,cost+=flow*e[i].c;
    return 1;
}
int main()
{
    scanf("%d%d",&m,&n);
    s=0;t=n+m+1;memset(head,-1,sizeof(head));
    for(int i=1;i<=m;++i)scanf("%d",&a[i]);
    for(int i=1;i<=n;++i)scanf("%d",&b[i]);
    for(int i=1;i<=m;++i)
    {
        for(int j=1;j<=n;++j)
        {
            scanf("%d",&w[i][j]);
            add(i,j+m,inf,w[i][j]);
        }
    }
    for(int i=1;i<=m;++i)add(s,i,a[i],0);
    for(int i=1;i<=n;++i)add(i+m,t,b[i],0);
    while(spfa());
    printf("%d\n",cost);
    memset(head,-1,sizeof(head));
    cnt=-1;cost=0;
    for(int i=1;i<=m;++i)
    {
        for(int j=1;j<=n;++j)
        {
            add(i,j+m,inf,-w[i][j]);
        }
    }
    for(int i=1;i<=m;++i)add(s,i,a[i],0);
    for(int i=1;i<=n;++i)add(i+m,t,b[i],0);
    while(spfa());
    printf("%d\n",-cost);
    return 0;
}