P8664 [蓝桥杯 2018 省 A] 付账问题

[蓝桥杯 2018 省 A] 付账问题

最小标准差

题目描述

几个人一起出去吃饭是常有的事。但在结帐的时候，常常会出现一些争执。

现在有 $n$ 个人出去吃饭，他们总共消费了 $S$ 元。其中第 $i$ 个人带了 $a_i$ 元。幸运的是，所有人带的钱的总数是足够付账的，但现在问题来了：每个人分别要出多少钱呢？

为了公平起见，我们希望在总付钱量恰好为 $S$ 的前提下，最后每个人付的钱的标准差最小。这里我们约定，每个人支付的钱数可以是任意非负实数，即可以不是 $1$ 分钱的整数倍。你需要输出最小的标准差是多少。

标准差的介绍：标准差是多个数与它们平均数差值的平方平均数，一般用于刻画这些数之间的“偏差有多大”。形式化地说，设第 $i$ 个人付的钱为 $b_i$ 元，那么标准差为 $s=\sqrt{\frac{1}{n}{\sum }_{i=1}^n(b_i-\frac{1}{n}{\sum }_{i=1}^n{b}_i)}$

输入格式

第一行包含两个整数 $n$、$S$；

第二行包含 $n$ 个非负整数 $a_1,\cdots ,a_n$。

输出格式

输出到标准输出。

输出最小的标准差，四舍五入保留 $4$ 位小数。

保证正确答案在加上或减去 $10^{-9}$ 后不会导致四舍五入的结果发生变化。

样例 #1

样例输入 #1

5 2333
666 666 666 666 666

样例输出 #1

0.0000

样例 #2

样例输入 #2

10 30
2 1 4 7 4 8 3 6 4 7

样例输出 #2

0.7928

提示

【样例解释】

1. 每个人都出 2333/5 元，标准差为 0。

【数据约定】

对于 $10\%$ 的数据，所有 $a_i$ 相等；

对于 $30\%$ 的数据，所有非 $0$ 的 $a_i$ 相等；

对于 $60\%$ 的数据，$n\le 1000$；

对于 $80\%$ 的数据，$n\le 10^5$；

对于所有数据，$n\le 5\times 10^5,0\le a_i\le 10^9$。

#include<bits/stdc++.h>
using namespace std;
const int N=5e5;//开到最大 
long long arr[N];
//不足均值的上交全部,够均值的要多给,由于这两波人支付的钱不一样,为了达到每个人支付一样多的钱,需要用到标准差(方差的平方根) ,使每个人实际给的钱离自己的额度越近越好 
//求标准差(方差的平方根),即对数组每个元素减去均值的差的取平方再相加,得到的和除总数再开方
//最小标准差,降序排列数列,然后找边界 
 
int main(){
	int n;long long s;cin>>n>>s;//人数和理论总金额 
	for(int i=0;i<n;++i)cin>>arr[i];//不同人的额度 
	sort(arr,arr+n);//降序从小到大 
	
	double avg=1.0*s/n,sum=0.0;//计算理论均值,初始化标准差里面的和值(值1) 
	for(int i=0;i<n;++i){//找到一个边界点,在该点处右边的人随便怎么用钱,左边的人钱少不足均值,支付其所有的钱 
		//最坏情况下,if会到达for的最后一个人时,其中每个人都需要支付其所有的钱 
		if(arr[i]*(n-i)<s){//找边界:边界存在于执行else时的i处,即在i这个点,左边所有元素都需要上交自己的总额度
		//!!!重点:i右边那个值,因为是降序,更右边的始终比它大,在该点处都满足了一个条件,即剩余的元素个数全部按此额度上交都可以大于等于剩余的理论总金额,更何况后面的值更大(从小到大的降序) 
			sum+=pow(arr[i]-avg,2);//偏差值的平方for和 
			s-=arr[i];//理论总金额减去这笔钱(即这个人上交全部额度),else前最终得到剩余金额 
		}else{
			//随便取,故取剩余理论金额的均值(这个均值小于等于退出if时的i额度arr[i]),而不是取僧与额度的均值,因为剩余额度的均值又会导致该剩余区间前面的值不够均值(即足以支付理论均值的人,可能支付不起剩余理论均值,他们仍需要上交所有的钱,故我们这个题才需要在降序序列中找到一个点,使该点的值乘以剩下的个数能补足缺的钱,) 
			double remavg=1.0*s/(n-i);//最终剩余金额取均值:剩余均值 
			sum+=pow(remavg-avg,2)*(n-i);//此处直接取剩余的样本全乘以剩余偏差值的平方,得到最终和式值 
			//↑此处不用for,因为是右边的样本随便取,不需要forif一个一个判断 ,直接乘以剩余样本个数即可 
			cout<<fixed<<setprecision(4)<<sqrt(sum/n);//和式除n开放得标准差(最小标准差) 
			break;
		}
	}
	
	return 0;
}