POJ 3709 K-Anonymous Sequence

#include <cstdio>
#include <cstring>

#define maxn 5000005
#define LL long long

LL a[maxn];
LL sum[maxn];
LL dp[maxn];
int Q[maxn];

int head,tail;

LL get0(int p)
{
    return dp[p]-sum[p]+p*a[p+1];
}

LL get1(int p,int q)
{
    return get0(p)-get0(q);
}

LL get2(int p,int q)
{
    return a[p+1]-a[q+1];
}

void pop(int i)
{
    while(head+1<tail&&get1(Q[head+1],Q[head])<=i*get2(Q[head+1],Q[head]))  head++;
}

void push(int i)
{
    while(head+1<tail&&(get1(Q[tail-1],Q[tail-2])*get2(i,Q[tail-1]))>=(get1(i,Q[tail-1])*get2(Q[tail-1],Q[tail-2])))  tail--;
    Q[tail++]=i;
}

int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        int n,k;
        int i;
        scanf("%d%d",&n,&k);
        for(sum[0]=0,i=1;i<=n;i++)
        {
            scanf("%lld",&a[i]);
            sum[i]=sum[i-1]+a[i];
        }
        memset(dp,0,sizeof(dp));
        head=tail=0;
        for(int i=k;i<2*k;i++)
        {
            dp[i]=sum[i]-i*a[1];
        }
        push(k);
        for(i=2*k;i<=n;i++)
        {
            pop(i);
            dp[i]=dp[Q[head]]-sum[Q[head]]-(i-Q[head])*a[Q[head]+1]+sum[i];
            push(i-k+1);
        }
        printf("%lld\n",dp[n]);
    }
    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值