4.4 【付账问题】

题目描述

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

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

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

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

输入描述

输出描述

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

输入输出样例

输入：

5 2333
666 666 666 666 666

输出：

0.0000

---

最终代码

1. c/c++

#include <bits/stdc++.h>
using namespace std;
const int M=5e5;
long long a[M];
int main()
{
    int n;  long long s;   
    scanf("%d %ld",&n,&s);
    for(int i=1;i<=n;i++)  
        scanf("%ld",&a[i]);
    
    //排序，每个人带的钱从小到大
    sort(a+1,a+n+1);    
    
    
    double avg=(1.0*s)/n;  //平均值
    double sum=0.0;
    
    for(int i=1;i<=n;i++)
    {
        if(a[i]*(n+1-i)<s)
        { 
            //需要把钱全拿出的人：（1）钱不够平均数的，（2）钱够平均数，但也不是很多的
            sum+=(a[i]-avg)*(a[i]-avg);
            //更新剩余钱数
            s-=a[i];           
        }
        else //不用把钱全拿出的人：非常有钱，不管怎么平均都够
        {       
            //更新平均出钱数
            double cur_avg=1.0*s/(n+1-i);  
            
            sum += (cur_avg-avg)*(cur_avg-avg)*(n+1-i);
            break;
        }
    }
    
    printf("%.4f",sqrt(sum/n));
    return 0;
}

---

2. java

import java.io.FileNotFoundException;
import java.util.Arrays;
import java.util.Scanner;
public class Main {
    public static void main(String args[]) throws FileNotFoundException {
        int n,i;
        long S;
        double ans=0,avg;
        Scanner input=new Scanner(System.in);
        n=input.nextInt();
        S=input.nextLong();
        long a[]=new long[n];
        for(i=0;i<n;i++) {
            a[i]=input.nextLong();
        }
        Arrays.sort(a);
        avg=(double)S/n;
        for(i=0;i<n;i++) {
            if(S<=(n-i)*a[i]) {
                ans+=(n-i)*Math.pow((double)S/(n-i)-avg,2);
                break;
            }
            ans+=Math.pow(a[i]-avg,2);
            S-=a[i];
        }
        System.out.printf("%.4f\n",Math.sqrt(ans/n));
    }
}

---

3. python

from math import *
n, s = map(int,input().split())
a = list(map(int,input().split()))
a.sort()

avg = s/n
sum = 0
for i in range(n):
     if a[i]*(n-i) < s:
          sum += pow(a[i]-avg,2)
          s -= a[i]
     else:
          cur_avg = s/(n-i);      #更新平均出钱数
          sum += pow(cur_avg-avg,2)*(n-i)
          break
print('{:.4f}'.format(sqrt(sum/(n))))

---

过程理解