ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

CF643C Levels and Regions 题解

2021-10-09 07:31:26  阅读:152  来源: 互联网

标签:frac CF643C int 题解 sum Regions read sm double


Codeforces
Luogu

Description.

给定所有 \(n\) 个数,划分成 \(K\) 个区间。
设第 \(i\) 个数所在段的左端点是 \(l_i\),则它的贡献是下式

\[\frac{\sum_{x=l_i}^it_i}{t_i} \]

问划分成 \(K\) 个区间后所有数的贡献和最小。

Solution.

一段的贡献

\[\begin{aligned} &=\sum_{a=l}^{r}\frac1{t_a}\sum_{b=l}^at_i\\ &=\sum_{a=l}^{r}\frac1{t_a}(s_a-s_{l-1})\\ &=\sum_{a=l}^{r}\frac{s_a}{t_a}-s_{l-1}\sum_{a=l}^{r}\frac{1}{t_a}\\ &=(sm_r-sm_{l-1})-s_{l-1}\times(q_r-q_{l-1})\\ \end{aligned} \]

其中 \(s_i=S(t_i),sm_i=S(\frac{s_i}{t_i}),q_i=S(\frac1{t_i})\)
然后直接斜率优化即可

Coding.

点击查看代码
//是啊,你就是那只鬼了,所以被你碰到以后,就轮到我变成鬼了{{{
#include<bits/stdc++.h>
using namespace std;typedef long long ll;
template<typename T>inline void read(T &x)
{
	x=0;char c=getchar(),f=0;
	for(;c<48||c>57;c=getchar()) if(!(c^45)) f=1;
	for(;c>=48&&c<=57;c=getchar()) x=(x<<1)+(x<<3)+(c^48);
	f?x=-x:x;
}
template<typename T,typename...L>inline void read(T &x,L&...l) {read(x),read(l...);}//}}}
const int N=200005;int n,K,t[N],qr[N],hd,tl;
double s[N],sm[N],q[N],dp[51][N];int nwk;
inline double X(int a) {return s[a];}
inline double Y(int a) {return dp[nwk-1][a]-sm[a]+s[a]*q[a];}
inline double slop(int x,int y) {return(Y(y)-Y(x))/(X(y)-X(x));}
int main()
{
	read(n,K);for(int i=1;i<=n;i++) read(t[i]);
	for(int i=1;i<=n;i++) s[i]=s[i-1]+t[i],sm[i]=sm[i-1]+s[i]/t[i],q[i]=q[i-1]+1.0/t[i];
	for(int i=1;i<=n;i++) dp[0][i]=sm[i];
	for(int k=1;k<K;k++)
	{
		qr[hd=tl=1]=0,nwk=k;for(int i=1;i<=n;i++)
		{
			while(hd<tl&&slop(qr[hd],qr[hd+1])<q[i]) hd++;
			int j=qr[hd];dp[k][i]=dp[k-1][j]-sm[j]+s[j]*q[j]-s[j]*q[i]+sm[i];
			//printf("%d -> %d\n",j,i);
			while(hd<tl&&slop(qr[tl],i)<=slop(qr[tl-1],qr[tl])) tl--;
			qr[++tl]=i;
		}
		//for(int i=1;i<=n;i++) printf("%.5lf%c",dp[k][i],i==n?'\n':' ');
	}
	return printf("%.10lf\n",dp[K-1][n]),0;
}

标签:frac,CF643C,int,题解,sum,Regions,read,sm,double
来源: https://www.cnblogs.com/pealfrog/p/15383871.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有