bzoj3437

斜率优化

先把原方程拆开(拆到没有括号)

然后运用结合律

用前缀和维护

数据范围只有一百万!

#include<cstdio>
#include<cstdlib>
#include<cstring>
long long a[1000005];
long long b[1000005];
long long f[1000005];
int       q[1000005];
long long s[1000005];
long long sc[1000005];
double Y(int j)
{
 return (double)(f[j]+sc[j]);   
}
double X(int j)
{
 return (double)s[j];   
}
double slop(int j1,int j2)
{
 return (Y(j2)-Y(j1))/(X(j2)-X(j1));    
}
int main()
{
 int n;
 scanf("%d",&n);
 for (int i=1;i<=n;i++)
 scanf("%lld",&a[i]);
 for (int i=1;i<=n;i++)
 {
  scanf("%lld",&b[i]);
  s[i]=s[i-1]+b[i];
  sc[i]=sc[i-1]+b[i]*i; 
 }
 q[1]=0;
 int l=1,r=1;
 for (int i=1;i<=n;i++)
 {
  while(l<r&&slop(q[l],q[l+1])<=i)l++;
  int j=q[l];
  f[i]=f[j]+i*(s[i-1]-s[j])-sc[i-1]+sc[j]+a[i];
  while(l<r&&slop(q[r],q[r-1])>slop(q[r],i))r--;
  q[++r]=i; 
 }
 printf("%lld",f[n]);
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值