#斜率优化，单调队列，动态规划#bzoj 2684 洛谷 4360 锯木场选址

最新推荐文章于 2022-02-05 15:35:36 发布

lemondinosaur

最新推荐文章于 2022-02-05 15:35:36 发布

阅读量199

点赞数

分类专栏：单调队列、斜率优化文章标签： bzoj 2684 洛谷 4360 锯木场选址

本文链接：https://blog.csdn.net/sugar_free_mint/article/details/84590042

版权

单调队列、斜率优化专栏收录该内容

30 篇文章 0 订阅

订阅专栏

题目链接

分析

设 $d p [i]$ 表示第二个锯木场修到第i个位置
那么状态转移方程是 $dp[i]=min(tot-dis[j](从j到n的距离)(后缀和)*(\sum_{k=1}^{j}w[k](前缀和))-dis[i]*(\sum_{k=j+1}^{i}w[k]))(j<i)$
然而这样肯定会超时，所以必须维护一个最小的 $j$ ，使 $d p$ 方程最优，所以说就用到了斜率优化
把 $d p$ 方程剖析，就得到
$t o t - d i s [j] * s u m [j] - d i s [i] * (s u m [i] - s u m [j])$
把已知项移出来得到
$t o t - d i s [i] * s u m [i] - d i s [i] * s u m [j] - d i s [j] * s u m [j]$
若使 $k (j < k)$ 更优秀，那么
$t o t - d i s [i] * s u m [i] - d i s [i] * s u m [j] - d i s [j] * s u m [j] > t o t - d i s [i] * s u m [i] - d i s [i] * s u m [k] - d i s [j] * s u m [k]$
最后得到 $\frac{dis[k]*sum[k]-dis[j]*sum[j]}{sum[k]-sum[j]}<dis[i]$
可以发现， $d i s$ 是单调递减的，也就是维护上凸壳，用单调队列实现

代码

#include <cstdio>
#include <cctype>
#define rr register
using namespace std;
int d[20001],w[20001],sum,n,head=1,tail=1,q[20001],ans=2147483647;
inline signed iut(){
    rr int ans=0; rr char c=getchar();
    while (!isdigit(c)) c=getchar();
    while (isdigit(c)) ans=(ans<<3)+(ans<<1)+c-48,c=getchar();
    return ans;
}
inline signed min(int a,int b){return (a<b)?a:b;}
inline double slope(int x,int y){return 1.0*(d[x]*w[x]-d[y]*w[y])/(w[x]-w[y]);}
inline signed count(int x,int y){return sum-d[y]*w[y]-d[x]*(w[x]-w[y]);}
signed main(){
    n=iut();
    for (rr int i=1;i<=n;++i) w[i]=iut(),d[i]=iut();
    for (rr int i=n-1;i;--i) d[i]+=d[i+1];
    for (rr int i=1;i<=n;++i) sum+=d[i]*w[i],w[i]+=w[i-1];
    for (rr int i=1;i<=n;++i){
        while (head<tail&&slope(q[head],q[head+1])>d[i]) ++head;//找出最优的队头
        ans=min(ans,count(i,q[head]));
        while (head<tail&&slope(q[tail-1],q[tail])<slope(q[tail],i)) --tail;//移去不优的队尾
        q[++tail]=i;
    }
    return !printf("%d",ans);
}