【洛谷 P4360】 [CEOI2004]锯木厂选址（斜率优化）

最新推荐文章于 2021-04-29 15:13:33 发布

weixin_33690963

最新推荐文章于 2021-04-29 15:13:33 发布

阅读量130

点赞数

原文链接：http://www.cnblogs.com/Qihoo360/p/10327950.html

版权

题目链接
一开始我的\(dp\)方程列错了，其实也不能说列错了，毕竟我交上去还是把暴力的分都拿到了，只是和题解的不一样，然后搞半天没搞出来去看题解，又看不懂，对不上，原来状态设置不一样~~自闭了~~。
\(f[i]=all-sum[j]*dis[j]-(sum[i]-sum[j])*dis[i]\)
\(f[i]=all-sum[j]*dis[j]-sum[i]*dis[i]+sum[j]*dis[i]\)
\(sum[j]*dis[j]=dis[i]*sum[j]-sum[i]*dis[i]+all-f[i]\)

#include <cstdio>
const int MAXN = 20010;
#define ll long long
inline ll min(const ll a, const ll b){
    return a < b ? a : b;
}
int n; 
ll all, ans = 1e17, sum[MAXN], dis[MAXN];
int w[MAXN], d[MAXN];
int q[MAXN], head, tail;
inline double k(int i, int j){
    return ((double)sum[i] * dis[i] - sum[j] * dis[j]) / (sum[i] - sum[j]);
}
int main(){
    scanf("%d", &n);
    for(int i = 1; i <= n; ++i)
       scanf("%d%d", &w[i], &d[i]);
    for(int i = n; i; --i)
       all += w[i] * (dis[i] = dis[i + 1] + d[i]);
    for(int i = 1; i <= n; ++i)
       sum[i] = sum[i - 1] + w[i];
    for(int i = 1; i <= n; ++i){
       while(head < tail && k(q[head], q[head + 1]) > dis[i]) ++head;
       int j = q[head];
       ans = min(ans, all - sum[j] * dis[j] - (sum[i] - sum[j]) * dis[i]);
       while(head < tail && k(q[tail - 1], q[tail]) <= k(q[tail], i)) --tail;
       q[++tail] = i;
    }
    printf("%lld\n", ans);
    return 0;
}

转载于:https://www.cnblogs.com/Qihoo360/p/10327950.html