P10966 K-Anonymous Sequence

题目链接

就相当于是把一个序列分成若干组,每组都有至少 k1 个数字,花费就是这组的数字和 sum,再减去最小值 min 乘以这个组的 cnt,也就是 sum(min×cnt)

贪心的考虑,假设说有三个数字 ai,ak,aj,(i<k<j),那么从 aj 转移到 ai 的代价一定比从 akai 的代价更大一些。

所以我们分成的组一定是连续的一段,那么就可以设计转移方程了。

fi 表示考虑前 i 个数的花费。

fi=min{fj+[sumisumj]aj+1(ij)}

其中 ijk,此时算法复杂度为 O(n2),无法通过,考虑优化。

首先上式中的 sum(i) 和是一个定值,提取并把里头的式子分解出来。

fi=sumi+minfjsumj(ij)×aj+1

考虑两个点 j,k(j<k)

k 优于 j 时,有:

fjsumj(ij)×aj+1fksumk(ik)×ak+1

化简:

fjfk+sumksumj+j×aj+1k×ak+1×aj+1ak+1i

这里要变号,因为 aj+1ak+1 小于等于0

这种斜率的形式可以用单调队列来维护。

X(j,k)=fjfk+sumksumj+j×aj+1k×ak+1×aj+1ak+1

队首维护应该很明白了,如果 X(j,k)i,就删除 j,当然由于可能除数为 0,所以我们用乘法来维护,所以又得变一次号,不要忘了。

队尾维护:求得是最小值,所以维护一个下凸壳。

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int N=5e5+5;
int T,n,k,a[N],sum[N],f[N],q[N];
int qy(int x){
    return f[x]-sum[x]+1ll*x*a[x+1];
}
int qx(int x){
    return a[x+1];
}
void solve(){
    cin>>n>>k;
    for(int i=1;i<=n;i++){
        cin>>a[i];
        sum[i]=sum[i-1]+a[i];
    }
    int l=1,r=1;
    q[1]=0;
    for(int i=1;i<=n;i++){
        while(l<r&&qy(q[l])-qy(q[l+1])>=i*(qx(q[l])-qx(q[l+1])))++l;
//      f[i]=qy(q[l])+sum[i]-i*qx(q[l]);
        f[i]=f[q[l]]+sum[i]-sum[q[l]]-(i-q[l])*a[q[l]+1];
        int z=i-k+1;
        if(z>=k){
            while(l<r&& (qy(q[r-1]) - qy(q[r])) * (qx(q[r])-qx(z))>= (qy(q[r])-qy(z)) * (qx(q[r-1])-qx(q[r])) )--r;
            q[++r]=z;
        }
    }
    cout<<f[n]<<"\n";
}
signed main(){
    cin>>T;
    while(T--)solve();
    return 0;
}
posted @   All_Unluck_Beginning  阅读(5)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· TypeScript + Deepseek 打造卜卦网站:技术与玄学的结合
· 阿里巴巴 QwQ-32B真的超越了 DeepSeek R-1吗?
· 如何调用 DeepSeek 的自然语言处理 API 接口并集成到在线客服系统
· 【译】Visual Studio 中新的强大生产力特性
· 2025年我用 Compose 写了一个 Todo App
点击右上角即可分享
微信分享提示