zoukankan      html  css  js  c++  java
  • [BZOJ2424/Luogu2517][HAOI2010]订货

    题目链接:

    BZOJ2424

    Luogu2517

    费用流大水题。

    源点与每月连边,容量(+infty)(可以生产无数件商品),费用(d_i)(商品单位成本)。

    每月与下一月连边,容量(S)(仓库容量),费用(m)(单位物品存储费用)。

    每月与汇点连边,容量(U_i)(需求量),费用(0)

    #include <queue>
    #include <cstdio>
    #include <cstring>
    
    inline int Min(int a,int b){return a<b?a:b;}
    
    int n,m,s,St,Ed,MaxFlow,MinCost;
    int Head[60],Next[405],To[405],Val[405],Cos[405],En=1;
    int Dis[60],Pre[60],Ref[60],Inq[60];
    
    inline void Add(int x,int y,int z,int w)
    {
    	Next[++En]=Head[x],To[Head[x]=En]=y,Val[En]=z,Cos[En]=+w;
    	Next[++En]=Head[y],To[Head[y]=En]=x,Val[En]=0,Cos[En]=-w;
    }
    
    bool SPFA()
    {
    	std::queue<int> q;
    	memset(Dis,0x3f,sizeof Dis);
    	q.push(St),Dis[St]=0,Ref[St]=1<<30;
    	for(int x,y;!q.empty();q.pop(),Inq[x]=0)
    		for(int i=Head[x=q.front()];i;i=Next[i])
    			if(Val[i]&&Dis[y=To[i]]>Dis[x]+Cos[i])
    			{
    				Dis[y]=Dis[x]+Cos[Pre[y]=i];
    				Ref[y]=Min(Ref[x],Val[i]);
    				if(!Inq[y])Inq[y]=1,q.push(y);
    			}
    	if(Dis[Ed]==0x3f3f3f3f)return false;
    	MaxFlow+=Ref[Ed],MinCost+=Ref[Ed]*Dis[Ed];
    	for(int x=Ed,i;x!=St;x=To[i^1])
    		Val[i=Pre[x]]-=Ref[Ed],Val[i^1]+=Ref[Ed];
    	return true;
    }
    
    int main()
    {
    	scanf("%d%d%d",&n,&m,&s),St=n+1,Ed=St+1;
    	for(int i=1,u;i<=n;++i)
    	{
    		scanf("%d",&u);
    		Add(i,Ed,u,0);
    		if(i!=n)Add(i,i+1,s,m);
    	}
    	for(int i=1,d;i<=n;++i)
    	{
    		scanf("%d",&d);
    		Add(St,i,0x3f3f3f3f,d);
    	}
    	while(SPFA());
    	printf("%d
    ",MinCost);
    	return 0;
    }
    
  • 相关阅读:
    文件处理
    三元运算
    python 第二章 二进制运算、字符编码、数据类型
    python 第一章 基础语法
    刷题-力扣-611. 有效三角形的个数
    刷题-力扣-581. 最短无序连续子数组
    刷题-力扣-LCP 07. 传递信息
    刷题-力扣-1833. 雪糕的最大数量
    刷题-力扣-168. Excel表列名称
    刷题-力扣-337. 打家劫舍 III
  • 原文地址:https://www.cnblogs.com/LanrTabe/p/10210571.html
Copyright © 2011-2022 走看看