【HAOI2010】订货

可以DP也可以是费用流，然而被我用非常简单的DP破了【开心】

原题：

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

0<=n<=50,0<=m<=10,0<=S<=10000,0<=Ui<=10000,0<=di<=100

恩紫萱的讲解说要用队列优化DP，我太弱没看懂，自己想出了一个我觉得比较妙的方法DP掉了

核心思路是在买东西的时候f[i]只与f[i-1]的最优值有关（注意是最优值）

然后要买的时候从1到S顺推，f[i]=min(f[i],f[i-1]+d[i])

恩原理我语文不好，直接丢个图吧

可以理解为f[i-1]代表了1到i-1所有的方案，但是都没有f[i]优，再往上考虑的话，由于增加的费用都是一样的，所以和f[1]到f[i-1]有关的方案不会优于f[i]相关的方案，就只考虑f[i]了

如果f[i]不是更优，f[i-1]代表了f[1]到f[i-1]的方案，只管使用即可

然后还有许多需要注意的细节（我就是因为细节问题拖了很久才A quq）

首先货物是可以进货不入库直接卖，不占用仓库空间，所以实际上仓库有M+a[i]的储存空间（每个月的需求也可能超过仓库容量），然后每个月结束的时候f中1到M不用处理，下个月直接转移即可，但是M+1到M+a[i]这一段一定要处理成正无穷，不然会影响到下个月（我就是因为这个细节卡了一下午，差点弃疗quq）

还有更多的细节问题，大家自己体验吧一。一

DP代码非常短

代码（DP）：

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #include<cmath>
 using namespace std;
 int read(){int z=,mark=;  char ch=getchar();
     while(ch<''||ch>''){if(ch=='-')mark=-;  ch=getchar();}
     while(ch>=''&&ch<=''){z=(z<<)+(z<<)+ch-'';  ch=getchar();}
     return z*mark;
 }
 const int oo=;
 int n,m,ns;  int a[],b[];
 int f[];
 int main(){//freopen("ddd.in","r",stdin);
     memset(f,,sizeof(f));
     cin>>n>>ns>>m;
     for(int i=;i<=n;i++)  a[i]=read();
     for(int i=;i<=n;i++)  b[i]=read();
     f[]=;  a[]=;
     for(int i=;i<=n;i++){
         for(int j=;j<=m;j++)  f[j]=f[j+a[i-]]+j*ns;
         for(int j=m+;j<=m+a[i-];j++)  f[j]=oo;
         for(int j=;j<=m+a[i];j++)  f[j]=min(f[j],f[j-]+b[i]);
     }
     cout<<f[a[n]]<<endl;
     return ;
 }

======================================================下面是费用流======================================================

至于费用流呐，建图也非常好想

首先利用割的思想，每个月到汇一条流量为需求，费用0的边，割掉这条边表示这个月的费用满足了

然后源到每个月一条流量正无穷，费用为这个月进货费用的边，表示物品可以随意买

每个月到下个月一条流量为仓库容量，费用为储存话费的边，表示上个月的货可以继承到下个月

至于我是怎么想出来的……直觉……我觉得还要继续发掘一下这其中潜在的思路规律

代码也很简单，标准的费用流模板：

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #include<cmath>
 using namespace std;
 int read(){int z=,mark=;  char ch=getchar();
     while(ch<''||ch>''){if(ch=='-')mark=-;  ch=getchar();}
     while(ch>=''&&ch<=''){z=(z<<)+(z<<)+ch-'';  ch=getchar();}
     return z*mark;
 }
 const int oo=;
 struct ddd{int next,y,evalue,rev,cost;}e[];  int LINK[],ltop=;
 inline void insert(int x,int y,int z,int _cost){
     e[++ltop].next=LINK[x];LINK[x]=ltop;e[ltop].y=y;e[ltop].evalue=z;e[ltop].rev=ltop+;e[ltop].cost=_cost;
     e[++ltop].next=LINK[y];LINK[y]=ltop;e[ltop].y=x;e[ltop].evalue=;e[ltop].rev=ltop-;e[ltop].cost=-_cost;
 }
 int n,m,ns;  int s,t;
 int dist[];
 int QUEUE[],head=;  bool visited[];
 int last[],last_e[];
 bool spfa(){
     memset(visited,,sizeof(visited));
     memset(dist,,sizeof(dist));
     QUEUE[head=]=s;  dist[s]=;
     for(int k=;k<=head;k++){
         for(int i=LINK[QUEUE[k]];i;i=e[i].next)
             if(e[i].evalue && dist[QUEUE[k]]+e[i].cost<dist[e[i].y]){
                 dist[e[i].y]=dist[QUEUE[k]]+e[i].cost;
                 last[e[i].y]=QUEUE[k],last_e[e[i].y]=i;
                 if(!visited[e[i].y])  QUEUE[++head]=e[i].y,visited[e[i].y]=true;
             }
         visited[QUEUE[k]]=false;
     }
     return dist[t]<oo;
 }
 int cost_flow(){
     int bowl=;
     while(spfa()){
         int min_flow=oo;
         for(int i=t;i!=s;i=last[i])  min_flow=min(min_flow,e[last_e[i]].evalue);
         for(int i=t;i!=s;i=last[i]){
             bowl+=min_flow*e[last_e[i]].cost;
             e[last_e[i]].evalue-=min_flow,e[e[last_e[i]].rev].evalue+=min_flow;
         }
     }
     return bowl;
 }
 int main(){//freopen("ddd.in","r",stdin);
     cin>>n>>ns>>m;  s=,t=n+;
     for(int i=;i<=n;i++)  insert(i,t,read(),);
     for(int i=;i<=n;i++)  insert(s,i,oo,read());
     for(int i=;i<n;i++)  insert(i,i+,m,ns);
     cout<<cost_flow()<<endl;
     return ;
 }