bzoj2424（费用流）

最新推荐文章于 2019-11-01 16:43:50 发布

zhhx2001

最新推荐文章于 2019-11-01 16:43:50 发布

阅读量355

点赞数

分类专栏：网络流

本文链接：https://blog.csdn.net/zhhx2001/article/details/52066824

版权

网络流专栏收录该内容

26 篇文章 0 订阅

订阅专栏

某公司估计市场在第i个月对某产品的需求量为Ui，已知在第i月该产品的订货单价为di，上个月月底未销完的单位产品要付存贮费用m，假定第一月月初的库存量为零，第n月月底的库存量也为零，问如何安排这n个月订购计划，才能使成本最低？每月月初订购，订购后产品立即到货，进库并供应市场，于当月被售掉则不必付存贮费。假设仓库容量为S。

还是构成一个线的建图

其实比餐巾计划简单多了。。

#include<cstdio>
#include<cstdlib>
#include<queue>
#include<algorithm>
#include<cstring>
using namespace std;
const int N=10200;
const int inf=0x3f3f3f3f;
int n,m,S;
int u[N],d[N];
struct aa
{
	int flow,cap,to,pre,w;
}edge[N*4];
int head[N],tot,dis[N],pre[N],now[N],s,t;
bool b[N];
void addedge(int x,int y,int z,int w)
{
	edge[++tot].to=y;edge[tot].w=w;edge[tot].cap=z;edge[tot].pre=head[x];head[x]=tot;
	edge[++tot].to=x;edge[tot].w=-w;edge[tot].cap=0;edge[tot].pre=head[y];head[y]=tot;
}
bool spfa()
{
	memset(dis,inf,sizeof(dis));
	memset(b,false,sizeof(b));
	memset(now,0,sizeof(now));
	memset(pre,0,sizeof(pre));
	queue<int> q;q.push(s);dis[s]=0;
	while (!q.empty())
	{
		int u=q.front();q.pop();
		for (int i=head[u];i;i=edge[i].pre)
		if (edge[i].cap>edge[i].flow&&dis[edge[i].to]>edge[i].w+dis[u])
		{
			dis[edge[i].to]=dis[u]+edge[i].w;
			now[edge[i].to]=i;pre[edge[i].to]=u;
			if (!b[edge[i].to]) {b[edge[i].to]=true;q.push(edge[i].to);}
		}
		b[u]=false;
	}
	return dis[t]!=inf;
}
int work()
{
	int ans=0;
	while(spfa())
	{
		int flow=inf;
		for (int i=t;i!=s;i=pre[i])
		flow=min(flow,edge[now[i]].cap-edge[now[i]].flow);
		ans+=flow*dis[t];
		for (int i=t;i!=s;i=pre[i])
		edge[now[i]].flow+=flow,edge[((now[i]-1)^1)+1].flow-=flow;
	}
	return ans;
}
int main()
{
	scanf("%d%d%d",&n,&m,&S);
	s=0,t=n+1;
	int x;
	for (int i=1;i<=n;i++) scanf("%d",&x),addedge(i,t,x,0);
	for (int i=1;i<=n;i++) scanf("%d",&x),addedge(s,i,inf,x);
	for (int i=1;i<n;i++) addedge(i,i+1,S,m);
	printf("%d",work()); 
	return 0;
}