CF724E Goods transportation 最小割 DP

照惯例CF的题不放原题链接。。。

题意:一个序列上有n个点,每个点有权值pi和si。表示这个点一开始有pi个物品,最多可以卖出si个物品,每个点都可以把物品向编号更大的点运输,但是对于i < j的任意点对(i, j)最多从i到j运c个物品。求最多能卖出多少个物品。

题解:

  如果不考虑数据范围的话,可以直接用网络流建图。s向每个点连流量为pi的边,表示一开始有pi的流量,每个点i向满足i < j的点j连流量为c的边,表示最多运送c个物品,每个点向t连流量为si的边,表示最多可以卖si个物品。

  最大流即为答案。

  但是这题数据范围过大,无法用网络流跑过,观察到这个图比较特殊,即如果我们知道了哪些点属于S集合(S在残余网络中可以到达的点),哪些点属于T集合,就可以直接算出最小割。

  即$ans = \sum_{i \in S}{s_{i}} + \sum_{i \in T}{p_{i}} + \sum_{i \in S}\sum_{j \in T}{c}$

  解释一下为什么是这样,应该是比较好理解的,从最小割这个角度来理解。

  如果一个点属于S,那么它肯定要断开与t的边,所以加上$s_{i}$,

  如果一个点属于T,那么它肯定要断开与s的边,所以加上$p_{i}$,

  如果一个点属于S,另一个点属于T,那么这两个点之间的边肯定要断开。

  因此我们设f[i][j]表示前i个点有j个属于S集的最小代价。

  那么如果这个点i我们划分到S集合里面,我们就加上$s_{i}$的代价。

  否则加上$p[i] + j \cdot c$的代价,$j \cdot c$表示这个点要与之前被划分到S集合里面的j个点都断开联系,因为边权都是c,所以代价就是$j \cdot c$.

  转移的时候注意一下j从0开始枚举。

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 #define R register int
 4 #define AC 10010
 5 #define LL long long
 6 
 7 int n;
 8 LL c, ans = 1e18;
 9 LL f[2][AC], s[AC], p[AC];
10 
11 inline int read()
12 {
13     int x = 0;char c = getchar();
14     while(c > '9' || c < '0') c = getchar();
15     while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
16     return x;
17 }
18 
19 inline void upmin(LL &a, LL b){
20     if(b < a) a = b;
21 }
22 
23 void pre()
24 {
25     n = read(), c = read();
26     for(R i = 1; i <= n; i ++) p[i] = read();
27     for(R i = 1; i <= n; i ++) s[i] = read();
28 }
29 
30 void work()
31 {
32     int now = 0;
33     for(R i = 1; i <= n; i ++)
34     {
35         now ^= 1;
36         memset(f[now], 63, sizeof(f[now]));
37         for(R j = 0; j <= i; j ++)
38             if(!j) f[now][j] = f[now ^ 1][j] + p[i];
39             else f[now][j] = min(f[now ^ 1][j] + p[i] + j * c, f[now ^ 1][j - 1] + s[i]);
40     }
41     for(R i = 0; i <= n; i ++) upmin(ans, f[now][i]);
42     printf("%lld\n", ans);
43 }
44 
45 int main()
46 {
47     freopen("in.in", "r", stdin);
48     pre();
49     work();
50     fclose(stdin);
51     return 0;
52 }
View Code

 

posted @ 2018-11-15 20:32  ww3113306  阅读(233)  评论(0编辑  收藏  举报
知识共享许可协议
本作品采用知识共享署名-非商业性使用-禁止演绎 3.0 未本地化版本许可协议进行许可。