第九届蓝桥杯省赛C++A组 付账问题（贪心）

标题：付账问题

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

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

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

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

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

第一行包含两个整数 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

请严格按要求输出，不要画蛇添足地打印类似：“请您输入...” 的多余内容。

注意：
main函数需要返回0;
只使用ANSI C/ANSI C++ 标准;
不要调用依赖于编译环境或操作系统的特殊函数。
所有依赖的函数必须明确地在源文件中 #include <xxx>
不能通过工程设置而省略常用头文件。

提交程序时，注意选择所期望的语言类型和编译器类型。
 

贪心

贪心算法的定义：
贪心算法是指在对问题求解时，总是做出在当前看来是最好的选择。也就是说，不从整体最优上加以考虑，只做出在某种意义上的局部最优解。贪心算法不是对所有问题都能得到整体最优解，关键是贪心策略的选择，选择的贪心策略必须具备无后效性，即某个状态以前的过程不会影响以后的状态，只与当前状态有关。
解题的一般步骤是：
1.建立数学模型来描述问题；
2.把求解的问题分成若干个子问题；
3.对每一子问题求解，得到子问题的局部最优解；
4.把子问题的局部最优解合成原来问题的一个解。
如果大家比较了解动态规划，就会发现它们之间的相似之处。最优解问题大部分都可以拆分成一个个的子问题，把解空间的遍历视作对子问题树的遍历，则以某种形式对树整个的遍历一遍就可以求出最优解，大部分情况下这是不可行的。贪心算法和动态规划本质上是对子问题树的一种修剪，两种算法要求问题都具有的一个性质就是子问题最优性(组成最优解的每一个子问题的解，对于这个子问题本身肯定也是最优的)。动态规划方法代表了这一类问题的一般解法，我们自底向上构造子问题的解，对每一个子树的根，求出下面每一个叶子的值，并且以其中的最优值作为自身的值，其它的值舍弃。而贪心算法是动态规划方法的一个特例，可以证明每一个子树的根的值不取决于下面叶子的值，而只取决于当前问题的状况。换句话说，不需要知道一个节点所有子树的情况，就可以求出这个节点的值。由于贪心算法的这个特性，它对解空间树的遍历不需要自底向上，而只需要自根开始，选择最优的路，一直走到底就可以了
 

思路：

这是典型的贪心问题。要求求最小标准差，就是要让每个人付出的钱尽可能接近平均值。那我们可以先求出初始的平均值S/n，如果有人付不起这个平均值，就把他所有的钱给出，用总钱数减去这个值，剩下代付的钱由剩下的人平均分，更新平均值，重复上面的操作，直到所有钱都付完为止。这里可以先将所有人的钱从小到大排序，这样只要有人付得起当前的平均值，后面的人也肯定能付的齐。就是每个人该付的钱与初始平均值的差的平方，是每个人的方差，最后累加起来在平均，就是标准差。

#include<iostream>
#include<algorithm>
#include<vector>

using namespace std;
int n, S;
double a[500005];
int main()
{
	cin >> n >> S;
	for (int i = 0; i < n; i++)
		cin >> a[i];
	sort(a,a+n);//由小到大排序
	double avg =(double) S * 1.0 / n;//计算平均值
	double newavg = avg;//新的平均值
	double ans = 0;//方差
	for (int i = 0; i < n; i++)
	{
		if (a[i] < newavg)//如果小于平均值，全部给出
		{
			ans += (double)pow(a[i] - avg, 2);//pow（）是求幂的函数，就是把a[i]的当前方差加入到ans中
			//pow也可以换成(a[i]-avg)*(a[i]-avg);
			S -= a[i];//在总钱数中减去a[i]的所有钱
			newavg = S * 1.0 / (n - i - 1);//更新平均值，剩下的钱由剩下的人平摊

		}
		else//因为经过排序，一个可以支付，剩下的都可以支付了
		{
			ans += (newavg-avg)*(newavg- avg);
			S -= newavg;

		}
	}
	double newans;//标准差
	newans = sqrt(ans / n);
	cout << newans << endl;
	system("pause");
	return 0;
}