zoukankan      html  css  js  c++  java
  • luogu P3648 [APIO2014]序列分割 |单调队列优化动态规划

    题目描述

    你正在玩一个关于长度为 nnn 的非负整数序列的游戏。这个游戏中你需要把序列分成 k+1k + 1k+1 个非空的块。为了得到 k+1k + 1k+1 块,你需要重复下面的操作 kkk 次:

    选择一个有超过一个元素的块(初始时你只有一块,即整个序列)

    选择两个相邻元素把这个块从中间分开,得到两个非空的块。

    每次操作后你将获得那两个新产生的块的元素和的乘积的分数。你想要最大化最后的总得分。

    输入格式

    第一行包含两个整数 nnn 和 kkk。保证 k+1≤nk + 1 leq nk+1≤n。

    第二行包含 nnn 个非负整数 a1,a2,⋯ ,ana_1, a_2, cdots, a_na1​,a2​,⋯,an​ (0≤ai≤104)(0 leq a_i leq 10^4)(0≤ai​≤104),表示前文所述的序列。

    输出格式

    第一行输出你能获得的最大总得分。

    第二行输出 kkk 个介于 111 到 n−1n - 1n−1 之间的整数,表示为了使得总得分最大,你每次操作中分开两个块的位置。第 iii 个整数 sis_isi​ 表示第 iii 次操作将在 sis_isi​ 和 si+1s_{i + 1}si+1​ 之间把块分开。

    如果有多种方案使得总得分最大,输出任意一种方案即可。


    #include<cstdio>
    #include<cstring>
    using namespace std;
    #define ll long long
    inline int read(){
        int x=0,f=1;char ch=getchar();
        while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
        while(ch>='0'&&ch<='9'){ x=(x<<3)+(x<<1)+(ch^48); ch=getchar(); }
    	return x*f;
    }
    const int N=1e5+10;
    int a[N],q[N],pre[N][205],n,k;
    ll s[N],f[N],g[N];
    inline double shape(int i,int j){
    	if(s[i]==s[j])return -1e18;
    	return 1.0*((g[i]-s[i]*s[i])-(g[j]-s[j]*s[j]))/(s[j]-s[i]);
    }
    int main(){
    	n=read(),k=read(); for(int i=1;i<=n;i++)a[i]=read(),s[i]=s[i-1]+a[i];
    	for(int j=1;j<=k;j++){
    		int l=1,r=0; q[++r]=0;
    		for(int i=1;i<=n;i++){
    			while(l<r&&shape(q[l],q[l+1])<=s[i])l++;
    			f[i]=g[q[l]]+s[q[l]]*(s[i]-s[q[l]]);
    			pre[i][j]=q[l];
    			while(l<r&&shape(q[r-1],q[r])>=shape(q[r],i))r--;
    			q[++r]=i;
    		}
    		memcpy(g,f,sizeof(f));
    	}
    	printf("%lld
    ",f[n]);
    	for(int x=n,i=k;i>=1;i--)x=pre[x][i],printf("%d ",x);
    	printf("
    ");
    }
    
  • 相关阅读:
    [CF1106E] 边界有问题
    UOJ 67 新年的毒瘤
    BZOJ 1093 [ZJOI2007]最大半连通子图
    codeforces round#510
    codeforces round#509
    杂谈
    BZOJ 3007 [SDOI2012]拯救小云公主
    BZOJ 1799
    BZOJ 3329
    BZOJ 3209 花神的数论题
  • 原文地址:https://www.cnblogs.com/naruto-mzx/p/12146996.html
Copyright © 2011-2022 走看看