【BZOJ1096】[ZJOI2007]仓库建设【斜率优化】
2016-03-04 20:55
337 查看
斜率优化咯。
一开始不知道怎么O(1)转移,看了别人blog才发现可以这么玩。
/* Footprints In The Blood Soaked Snow */
#include <cstdio>
#include <algorithm>
using namespace std;
typedef long long LL;
const int maxn = 1000005, maxq = maxn;
int n, q[maxq];
LL dp[maxn], sum[maxn], dif[maxn], x[maxn], p[maxn], c[maxn];
inline int iread() {
int f = 1, x = 0; char ch = getchar();
for(; ch < '0' || ch > '9'; ch = getchar()) f = ch == '-' ? -1 : 1;
for(; ch >= '0' && ch <= '9'; ch = getchar()) x = x * 10 + ch - '0';
return f * x;
}
inline LL dy(int a, int b) {
return (dp[a] + dif[a]) - (dp[b] + dif[b]);
}
inline LL dx(int a, int b) {
return sum[a] - sum[b];
}
int main() {
n = iread();
for(int i = 1; i <= n; i++) x[i] = iread(), p[i] = iread(), c[i] = iread();
for(int i = 1; i <= n; i++) sum[i] = sum[i - 1] + p[i], dif[i] = dif[i - 1] + x[i] * p[i];
int h = 1, t = 0; q[++t] = 0;
for(int i = 1; i <= n; i++) {
for(; h < t && dy(q[h + 1], q[h]) < dx(q[h + 1], q[h]) * x[i]; h++);
dp[i] = dp[q[h]] + (sum[i] - sum[q[h]]) * x[i] - (dif[i] - dif[q[h]]) + c[i] ;
for(; h < t && dy(q[t], q[t - 1]) * dx(i, q[t]) >= dy(i, q[t]) * dx(q[t], q[t - 1]); t--);
q[++t] = i;
}
printf("%lld\n", dp
);
return 0;
}
一开始不知道怎么O(1)转移,看了别人blog才发现可以这么玩。
/* Footprints In The Blood Soaked Snow */
#include <cstdio>
#include <algorithm>
using namespace std;
typedef long long LL;
const int maxn = 1000005, maxq = maxn;
int n, q[maxq];
LL dp[maxn], sum[maxn], dif[maxn], x[maxn], p[maxn], c[maxn];
inline int iread() {
int f = 1, x = 0; char ch = getchar();
for(; ch < '0' || ch > '9'; ch = getchar()) f = ch == '-' ? -1 : 1;
for(; ch >= '0' && ch <= '9'; ch = getchar()) x = x * 10 + ch - '0';
return f * x;
}
inline LL dy(int a, int b) {
return (dp[a] + dif[a]) - (dp[b] + dif[b]);
}
inline LL dx(int a, int b) {
return sum[a] - sum[b];
}
int main() {
n = iread();
for(int i = 1; i <= n; i++) x[i] = iread(), p[i] = iread(), c[i] = iread();
for(int i = 1; i <= n; i++) sum[i] = sum[i - 1] + p[i], dif[i] = dif[i - 1] + x[i] * p[i];
int h = 1, t = 0; q[++t] = 0;
for(int i = 1; i <= n; i++) {
for(; h < t && dy(q[h + 1], q[h]) < dx(q[h + 1], q[h]) * x[i]; h++);
dp[i] = dp[q[h]] + (sum[i] - sum[q[h]]) * x[i] - (dif[i] - dif[q[h]]) + c[i] ;
for(; h < t && dy(q[t], q[t - 1]) * dx(i, q[t]) >= dy(i, q[t]) * dx(q[t], q[t - 1]); t--);
q[++t] = i;
}
printf("%lld\n", dp
);
return 0;
}
相关文章推荐
- [BZOJ1096][ZJOI2007][DP][斜率优化]仓库建设
- [BZOJ1492][NOI2007][斜率优化][动态凸包][DP][分治]货币兑换cash
- HDU 3045
- 【bzoj1010】【斜率优化】【HNOI2008】玩具装箱toy
- (解题报告) Uva 1616 Caravan Robbers (商队抢劫者)(上凸包+队列+结构体)
- [BZOJ1010][HNOI2008]玩具装箱toy(斜率优化)
- hdu 3507 斜率优化dp入门
- 2015湖南省队集训DAY4——hoodle
- HDU3507 [斜率优化]
- DP:斜率优化 Slope optimization
- 【BZOJ 2149】 拆迁队
- 斜率优化DP 【pascal】
- 1010: [HNOI2008]玩具装箱toy
- 【BZOJ3675】[Apio2014]序列分割【斜率优化】
- 斜率优化专题1——bzoj 1597 [Usaco2008 Mar] 土地购买 题解
- 斜率优化专题2——bzoj 1010 [HNOI2008]玩具装箱toy 题解
- 斜率优化专题3——bzoj 3156 防御准备 题解
- 斜率优化专题4——bzoj 1911: [Apio2010] 特别行动队 题解
- 斜率优化专题5——bzoj 1096 [ZJOI2007]仓库建设 题解
- BZOJ 1010: [HNOI2008]玩具装箱toy