第九届-A10-付账问题

标题：付账问题

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

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

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

标准差的介绍：标准差是多个数与它们平均数差值的平方平均数，一般用于刻画这些数之间的“偏差有多大”。形式化地说，设第 i 个人付的钱为 bi 元，那么标准差为 : 

【输入格式】
从标准输入读入数据。

第一行包含两个整数 n、S；
第二行包含 n 个非负整数 a1, ..., an。

【输出格式】
输出到标准输出。

输出最小的标准差，四舍五入保留 4 位小数。
保证正确答案在加上或减去 10^ 9 后不会导致四舍五入的结果发生变化。

【样例1输入】
5 2333
666 666 666 666 666

【样例输出】
0.0000

【样例解释】
每个人都出 2333/5 元，标准差为 0。

再比如：
【样例输入】
10 30
2 1 4 7 4 8 3 6 4 7

【样例输出】
0.7928

【数据说明】
对于 10% 的数据，所有 ai 相等；
对于 30% 的数据，所有非 0 的 ai 相等；
对于 60% 的数据，n≤ 1000；
对于 80% 的数据，n ≤ 10^5；
对于所有数据，n ≤ 5 × 10^5, 0 ≤ ai ≤ 10^9。

资源约定：
峰值内存消耗（含虚拟机） < 256M
CPU消耗  < 1000ms

问题分析：贪心算法。
如果每一个人的钱都大于平均值，则标准差为0.0000。如果有人的钱低于平均值，则少的部分由剩下的人来平均分摊。
先求一个不变的平均费用avg = S/n,如果有人付不起avg，他们要把他们的钱全付完，不够的需要钱多于avg的人付，但是这些付不起avg的人少付的钱会抬高>avg的人的付费平均值，或许又会有人付不起这个新的avg，就是nowavg，所以循环这个过程，直到没有人付不起被抬高的avg。这样的话那些都能付得起新的avg的有钱人付的钱均相等，都是新的avg。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
using namespace std;
const int MAXN = 600000;

int a[MAXN];
	
int main(){
	int n,S;
	scanf("%d %d",&n,&S);
	for(int i = 0;i<n;i++)
		scanf("%d",&a[i]);
	double avg = (double)S/n;  //用于保存每个人应付款的平均值 
	double LS = (double)S; //保存当前剩余的人应付款的总金额 
	double nowavg = avg;   //保存当前剩余的人应付款的平均值 
	double r = 0.0;              //保存(bi - avg)^2 
	
	sort(a,a+n);//将ai升序排序，钱少的人在前面，钱多的人在后面
	
	for(int i = 0;i<n;i++){
		//如果第i个人付不起当前平均的钱，则让他把自己的钱全掏出来，不够的部分让剩下的人替他分担 
		if(a[i] < nowavg){
			r += (a[i]-avg) * (a[i]-avg);
			LS -= a[i];  //现在第i个人已经给了钱
			nowavg = LS/(n-i-1); 
		}else{
			//如果第i个人付得起当前平均的钱nowavg,那就让他付nowavg这么多 
			r += (avg-nowavg) * (avg-nowavg);
			LS -= nowavg;
		}
	} 
	
	printf("%.4lf\n",sqrt(r/n)); 
	
	return 0;	
}