关于构图:
首先,这是一道最小费用最大流和最大费用最大流(就是最小费用最大流,但建边时的费用取负,就可以保证得到一个最小的费用,再取反,就得到最大费用)的题目,比较裸(类似于分配问题)。
构图:
1.首先有一个超级源点和汇点:st=0,ed=2*n+1。
2.源点向每一个仓库建一条流量为仓库的容量,费用为0的边,表示每一个仓库只能获得这么多的货物。
3.每一个零售商店向汇点建一条流量为零售商店的容量,费用为0边,表示每一个零售仓库只能存储获得这么多的货物。
4.每一个仓库向每一个零售商店建一条流量为仓库的容量,费用为x的边,表示当前这个仓库最多能向这个零售商店运输这么多货物,且费用x。
跑两边费用流(最小和最大)即可。
关于代码:
#include<cstdio>
#include<queue>
#include<cstring>
#include<queue>
#include<algorithm>
#define INF 2147483647
using namespace std;
queue<int> f;
int n,m,len,st,ed;
struct node{int x,y,c,d,next;} a[10000];
int le1[200],le2[200],b[200][200],last[10000],dis[10000],pre[10000],pos[10000],p[10000];
bool bz[10000];
void ins(int x,int y,int c,int d)
{
a[++len].x=x;a[len].y=y;a[len].c=c;a[len].d=d;a[len].next=last[x];last[x]=len;
}
void init(int x)
{
len=-1;
memset(last,-1,sizeof(last));
for(int i=1;i<=n;i++)
ins(st,i,le1[i],0),ins(i,st,0,0);
for(int i=1;i<=m;i++)
ins(i+n,ed,le2[i],0),ins(ed,i+n,0,0);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
ins(i,j+n,le1[i],b[i][j]*x),ins(j+n,i,0,-b[i][j]*x);
}
bool spfa()
{
memset(bz,true,sizeof(bz));
bz[st]=false;
memset(dis,63,sizeof(dis));
dis[st]=0;
p[st]=INF;
f.push(st);
while(!f.empty())
{
int x=f.front();
bz[x]=true;
for(int i=last[x];i>=0;i=a[i].next)
{
int y=a[i].y;
if(a[i].c>0&&dis[y]>dis[x]+a[i].d)
{
dis[y]=dis[x]+a[i].d;
pos[y]=x;
pre[y]=i;
p[y]=min(p[x],a[i].c);
if(bz[y])
{
f.push(y);
bz[y]=false;
}
}
}
f.pop();
}
return dis[ed]<1061109567;
}
int flow()
{
int ans=0;
while(spfa())
{
ans+=p[ed]*dis[ed];
for(int i=ed;i!=st;i=pos[i])
{
a[pre[i]].c-=p[ed];
a[pre[i]^1].c+=p[ed];
}
}
return ans;
}
int main()
{
scanf("%d %d",&n,&m);
st=0,ed=n+m+1;
for(int i=1;i<=n;i++)
scanf("%d",&le1[i]);
for(int i=1;i<=m;i++)
scanf("%d",&le2[i]);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
scanf("%d",&b[i][j]);
init(1);
printf("%d\n",flow());
init(-1);
printf("%d",-flow());
}