BZOJ1096_仓库建设_KEY

题目传送门

一道斜率优化的题目,加深了印象。

设sum[i]=∑p[i],S[i]=∑p[i]*x[i]。

暴力方程加前缀和优化:

f[i]=min(f[j]+c[i]+(sum[i]-sum[j])*x[i]-(S[i]-S[j])};

然后变形:

f[j]+c[i]+sum[i]*x[i]-(S[i]-S[j])=x[i]*sum[j]+f[i]
y =k x +b

求最小截距,因为k单调,所以是个下凸包,单调队列维护。

//红体字在y相减时会抵消,所以无影响。

code:

复制代码
/**************************************************************
    Problem: 1096
    User: yekehe
    Language: C++
    Result: Accepted
    Time:1100 ms
    Memory:47796 kb
****************************************************************/
 
#include <cstdio>
#include <algorithm>
using namespace std;
 
typedef long long ll;
 
char tc()
{
    static char tr[100000],*A=tr,*B=tr;
    return A==B&&(B=(A=tr)+fread(tr,1,100000,stdin),A==B)?EOF:*A++;
}
 
int read()
{
    char c;while(c=tc(),c<'0'||c>'9');
    int x=c-'0';while(c=tc(),c>='0'&&c<='9')x=x*10+c-'0';
    return x;
}
 
const int MAXN=1000005;
ll N,x[MAXN],p,c[MAXN];
ll sum[MAXN],S[MAXN],f[MAXN];
ll l[MAXN],h,t;
double X(int i){return sum[i];}
double Y(int i){return f[i]+S[i];}
double get(int x,int y){return (Y(y)-Y(x))/(X(y)-X(x));}
 
int main()
{
    N=read();
    register int i,j;
        for(i=1;i<=N;i++){
            x[i]=read(),p=read(),c[i]=read();
            sum[i]=sum[i-1]+p;
            S[i]=S[i-1]+p*x[i];
        }
    h=t=0;
        for(i=1;i<=N;i++){
            while(h<t&&get(l[h],l[h+1])<x[i])h++;
            j=l[h];f[i]=f[j]+c[i]+(sum[i]-sum[j])*x[i]-(S[i]-S[j]);
            while(h<t&&get(l[t],l[t-1])>get(l[t],i))t--;
            l[++t]=i;
        }
    printf("%lld",f[N]);
    return 0;
}
复制代码

 

posted @   Cptraser  阅读(178)  评论(0编辑  收藏  举报
编辑推荐:
· DeepSeek 解答了困扰我五年的技术问题
· 为什么说在企业级应用开发中,后端往往是效率杀手?
· 用 C# 插值字符串处理器写一个 sscanf
· Java 中堆内存和栈内存上的数据分布和特点
· 开发中对象命名的一点思考
阅读排行:
· DeepSeek 解答了困扰我五年的技术问题。时代确实变了!
· PPT革命!DeepSeek+Kimi=N小时工作5分钟完成?
· What?废柴, 还在本地部署DeepSeek吗?Are you kidding?
· 赶AI大潮:在VSCode中使用DeepSeek及近百种模型的极简方法
· DeepSeek企业级部署实战指南:从服务器选型到Dify私有化落地
点击右上角即可分享
微信分享提示