[luogu P3648] [APIO2014]序列分割

最新推荐文章于 2023-08-19 13:58:25 发布

a8233821

最新推荐文章于 2023-08-19 13:58:25 发布

阅读量95

点赞数

原文链接：http://www.cnblogs.com/whc200305/p/7659240.html

版权

[luogu P3648] [APIO2014]序列分割

题目描述

小H最近迷上了一个分隔序列的游戏。在这个游戏里，小H需要将一个长度为n的非负整数序列分割成k+1个非空的子序列。为了得到k+1个子序列，小H需要重复k次以下的步骤：

1.小H首先选择一个长度超过1的序列（一开始小H只有一个长度为n的序列——也就是一开始得到的整个序列）；

2.选择一个位置，并通过这个位置将这个序列分割成连续的两个非空的新序列。

每次进行上述步骤之后，小H将会得到一定的分数。这个分数为两个新序列中元素和的乘积。小H希望选择一种最佳的分割方式，使得k轮之后，小H的总得分最大。

输入输出格式

输入格式：

输入第一行包含两个整数n，k（k+1≤n）。

第二行包含n个非负整数a1，a2，...，an（0≤ai≤10^4），表示一开始小H得到的序列。

输出格式：

第一行包含一个整数，为小H可以得到的最大分数。

第二行输出

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

输入输出样例

输入样例#1：

7 3
4 1 3 4 0 2 3

输出样例#1：

108
1 3 5

说明

你可以通过下面这些操作获得

初始时你有一块 $\times (1 + 3 + 4 + 0 + 2 + 3) = 524×(1+3+4+0+2+3)=52 分。$

你现在有两块 $\times (4 + 0 + 2 + 3) = 36(1+3)×(4+0+2+3)=36 分。$

你现在有三块 $\times (2 + 3) = 20(4+0)×(2+3)=20 分。$

所以，经过这些操作后你可以获得四块

限制与约定

第一个子任务共 11 分，满足 $\leq k < n \leq 101≤k<n≤10。$

第二个子任务共 11 分，满足 $\leq k < n \leq 501≤k<n≤50。$

第三个子任务共 11 分，满足 $\leq k < n \leq 2001≤k<n≤200。$

第四个子任务共 17 分，满足 $\leq n \leq 1000, 1 \leq k \leq \min{n - 1, 200}2≤n≤1000,1≤k≤minn−1,200。$

第五个子任务共 21 分，满足 $\leq n \leq 10000, 1 \leq k \leq \min{n - 1, 200}2≤n≤10000,1≤k≤minn−1,200。$

第六个子任务共 29 分，满足 $\leq n \leq 100000, 1 \leq k \leq \min{n - 1, 200}2≤n≤100000,1≤k≤minn−1,200。$

$-s[i]<((f[j]-s[j]^2)-(f[k]-s[k]^2))/(s[j]-s[k])$

$这很像一个斜率式，我们可以将-s[i]看做点(f[j]-s[j]^2,s[j])和点(f[k]-s[k]^2,s[k])之间连线的斜率，只有满足上式（大于-s[i]），j才比k更优越。$

现在我们简单说一下如何维护这个单调队列，

现在我们想要计算第

然后有了决策点，我们就可以计算出

假设我们要添加的点为

（Hint:部分摘自：http://www.cnblogs.com/akhpl/p/6715148.html，大力推荐）

 1 %:pragma GCC optimize(2)
 2 #include<bits/stdc++.h>
 3 #define sqr(x) ((x)*(x))
 4 #define LL long long
 5 using namespace std;
 6 const int N=100005;
 7 const double inf=1e18;
 8 int n,k,l[2],r[2],roa[205][N];
 9 LL s[N],f[2][N];
10 struct point {
11     LL x,y; int i;
12     point() {}
13     point(LL _x,LL _y,int _i):x(_x),y(_y),i(_i) {}
14 }st[2][N];
15 int read() {
16     int x=0,f=1; char ch=getchar();
17     while (ch<'0'||ch>'9') f=(ch=='-')?-1:1,ch=getchar();
18     while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
19     return x*f;
20 }
21 double slope(point u,point v) {
22     return u.x==v.x?(u.y<v.y?inf:-inf):1.0*(v.y-u.y)/(v.x-u.x);
23 }
24 LL get(int y,int x,int i,LL k) {
25     while (l[i]<r[i]&&slope(st[i][l[i]],st[i][l[i]+1])>1.0*k) l[i]++;
26     roa[y][x]=st[i][l[i]].i;
27     return st[i][l[i]].y-k*st[i][l[i]].x;
28 }
29 void insert(int i,point cur) {
30     while (l[i]<r[i]&&slope(st[i][r[i]-1],st[i][r[i]])<slope(st[i][r[i]-1],cur)) r[i]--;
31     st[i][++r[i]]=cur;
32 }
33 int main() {
34     n=read(),k=read(),s[0]=0;
35     for (int i=1; i<=n; i++) s[i]=s[i-1]+read();
36     for (int j=1,now,pre; j<=k; j++) {
37         now=j&1,pre=1-now,l[pre]=1,r[pre]=0,st[++r[pre]][pre]=point(0,0,0);
38         for (int i=1; i<=n; i++) {
39             f[now][i]=get(j,i,pre,-s[i]);
40             insert(pre,point(s[i],f[pre][i]-s[i]*s[i],i));
41         }
42     }
43     printf("%lld\n",f[k&1][n]);
44     for (int i=n,j=k; j>=1; j--) i=roa[j][i],printf("%d ",i);
45     return 0;
46 }

View Code

转载于:https://www.cnblogs.com/whc200305/p/7659240.html

a8233821

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
[luogu P3648] [APIO2014]序列分割

[luogu P3648] [APIO2014]序列分割题目描述小H最近迷上了一个分隔序列的游戏。在这个游戏里，小H需要将一个长度为n的非负整数序列分割成k+1个非空的子序列。为了得到k+1个子序列，小H需要重复k次以下的步骤：1.小H首先选择一个长度超过1的序列（一开始小H只有一个长度为n的序列——也就是一开始得到的整个序列）；2.选择一个位置，并通过这个位置将这个序列...
复制链接

扫一扫