bzoj2424: [HAOI2010]订货(费用流)

Posted 2020-10-28 CHerish_OI

2424: [HAOI2010]订货

题目：传送门 

 

---

 

 

题解：

　　 做多了最小割，做一下费用流练手

　　 其实很容易就看出来是费用流啊，伏地魔肉老师用单调队列ORZ

　　 st直接连每个月，流量无限（随便买，花钱而已），费用就是给出来的

　　 然后因为可以贮存嘛，那就第i个月连第i+1个月，流量为S（仓库容量），费用为m就OK

　　 最后就是卖出去咯，第i个月连ed流量为需要的，费用为0

　　 难受，打个spfa忘记出队列...

 　　 

---

 

 

代码：

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<cstdlib>
 4 #include<cmath>
 5 #include<algorithm>
 6 #define inf 999999999
 7 using namespace std;
 8 struct node
 9 {
10     int x,y,c,d,next,other;
11 }a[1100];int len,last[110];
12 void ins(int x,int y,int c,int d)
13 {
14     int k1,k2;
15     k1=++len;
16     a[len].x=x;a[len].y=y;a[len].c=c;a[len].d=d;
17     a[len].next=last[x];last[x]=len;
18     
19     k2=++len;
20     a[len].x=y;a[len].y=x;a[len].c=0;a[len].d=-d;
21     a[len].next=last[y];last[y]=len;
22     
23     a[k1].other=k2;
24     a[k2].other=k1;
25 }
26 int list[11000],la[11000],d[11000],st,ed,n,m,head,tail,ans,S;
27 bool v[11000];
28 bool spfa()
29 {
30     for(int i=1;i<=ed;i++)d[i]=inf;
31     memset(v,false,sizeof(v));memset(la,0,sizeof(la));
32     list[1]=st;d[st]=0;v[st]=true;head=1;tail=2;
33     while(head!=tail)
34     {
35         int x=list[head];
36         for(int k=last[x];k;k=a[k].next)
37         {
38             int y=a[k].y;
39             if(a[k].c>0 && d[y]>d[x]+a[k].d)
40             {
41                 d[y]=d[x]+a[k].d;la[y]=k;
42                 if(v[y]==false)
43                 {
44                     v[y]=true;
45                     list[tail++]=y;
46                     if(tail==ed+1)tail=1;
47                 }
48             }
49         }
50         head++;if(head==ed+1)head=1;
51         v[x]=false;
52     }
53     if(d[ed]==inf)return false;
54     return true;
55 }
56 int g_f()
57 {
58     int x,k,f=inf,ans=0;
59     x=ed;while(x!=st){k=la[x];f=min(f,a[k].c);x=a[k].x;}
60     x=ed;while(x!=st){k=la[x];ans+=f*a[k].d;a[k].c-=f;a[a[k].other].c+=f;x=a[k].x;}
61     return ans;
62 }
63 int A[110],B[110];
64 int main()
65 {
66     scanf("%d%d%d",&n,&m,&S);len=0;memset(last,0,sizeof(last));st=n+1,ed=st+1;
67     for(int i=1;i<=n;i++)scanf("%d",&A[i]);
68     for(int i=1;i<=n;i++)scanf("%d",&B[i]);
69     for(int i=1;i<=n;i++)ins(st,i,inf,B[i]);
70     for(int i=1;i<=n;i++){if(i!=n)ins(i,i+1,S,m);ins(i,ed,A[i],0);}
71     int ans=0;while(spfa())ans+=g_f();
72     printf("%d\n",ans);
73     return 0;
74 }