标签：平方 int db CF long st Access Points define

题面传送门
首先距离的平方等于\(x\)轴平方加上\(y\)轴平方。所以\(x,y\)轴独立。
然后考虑\(x\)，如果\(x\)单调递增，那么直接令所有位置都在\(x\)上即可。
但是不会这样，如果\(x\)是单调下降的，通过简单的二次函数，我们可以知道这一段应该被赋值成同一个值，为\(\frac{\sum\limits_{i=a}^{b}{x_i}}{b-a}\)，即这一段的平均数。
但是还有一个要求就是\(x\)递增，不难发现直接维护一个单调栈合并即可。
时间复杂度\(O(n)\)
code:

#include<bits/stdc++.h>
#define I inline
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define abs(x) ((x)>0?(x):-(x))
#define re register
#define RI re int
#define ll long long
#define db double
#define lb long db
#define N 100000
#define M N*N+5
#define mod 998244353
#define Mod (mod-1)
#define eps (1e-5)
#define U unsigned int
#define it iterator
#define Gc() getchar() 
#define Me(x,y) memset(x,y,sizeof(x))
#define d(x,y) (m*x+(y))
#define R(n) (rand()*rand()%(n)+1)
#define Pc(x) putchar(x)
#define LB lower_bound
#define UB upper_bound  
using namespace std;
int n,m,k,st[N+5],H;db ToT,x[N+5],y[N+5],A[N+5],Q[N+5];
I db calc(int j,int i){return (Q[i]-Q[j])/(i-j);}
I db Solve(){
	db Ans=0;RI i;H=0;for(i=1;i<=n;i++){Q[i]=Q[i-1]+A[i];Ans+=A[i]*A[i];
		while(H&&calc(st[H-1],st[H])>calc(st[H],i)) H--;st[++H]=i;
	}for(i=1;i<=H;i++) Ans-=calc(st[i-1],st[i])*calc(st[i-1],st[i])*(st[i]-st[i-1]);return Ans;
}
int main(){
	freopen("1.in","r",stdin);
	scanf("%d",&n);for(RI i=1;i<=n;i++) scanf("%lf%lf",&x[i],&y[i]);memcpy(A,x,sizeof(x));ToT+=Solve();memcpy(A,y,sizeof(y));printf("%.8lf\n",Solve()+ToT);
}

标签：平方,int,db,CF,long,st,Access,Points,define
来源： https://www.cnblogs.com/275307894a/p/15863587.html

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