Hdu 5073 Galaxy 精度问题

　　思路： 其实求解很简单直接说解法，移动K个后 上下的角动量最小，能肯定是相连的（n-k）个，至于为什么 你自己好好想想（easy）；

　　　　　对于一些等质量的质点中心在 所在位置和除以点的个数

　　　　　average=sum[l,l+(n-k)-1]/(n-k);

　　　　　 一个点的值： (pi-average)* (pi-average)

　　　　　 也就是 pi^2+avery^2 - 2*pi*average

　　　　　 多个点相加也就是 ∑pi^2+(n-k)*sum*sum/(n-k)/(n-k)   -  2*∑pi*sum/(n-k);

　　　　　　　　　　　　=   ∑pi^2+(n-k)*sum*sum/(n-k)/(n-k)   -  2*sum*sum/(n-k);

　　　　　　　　　　　　=   ∑pi^2   -  sum*sum/(n-k);

　　　　　所以要处理一下 ”普通和“，“平方和” 就好了

　　　　　但是这个卡精度：我的做法  每次比较   （n-k）*∑pi^2   -  sum*sum

　　　　　剩下最小的 再除以（n-k） 就好了

代码.cpp

 1 #include <string>
 2 #include <cstring>
 3 #include <iostream>
 4 #include <cstdio>
 5 #include <algorithm>
 6 using namespace std;
 7 typedef long long ll;
 8 int n,k,d[50005];
 9 ll sum[50005];
10 ll getsum(int i,int j)
11 {
12     if(i)return sum[j-1]-sum[i-1];
13     else return sum[j-1];
14 }
15
16 int main()
17 {
18     int t;
19     scanf("%d",&t);
20     while(t--)
21     {
22         memset(d,0,sizeof d);
23         memset(sum,0,sizeof sum);
24         scanf("%d%d",&n,&k);
25         for(int i=0;i<n;i++)
26             scanf("%d",&d[i]);
27         sort(d,d+n);
28         for(int i=0;i<n;i++)
29             if(i) sum[i]=sum[i-1]+(ll)d[i];
30             else sum[i]=(ll)d[i];
31         double messc=0.0,minv=0.0,now=0.0;
32         for(int j=0;j<=k;j++)
33         {
34             ll nows=getsum(j,n-(k-j));
35             messc=(double)nows/(double)(n-k);
36             now=0.0;
37             for(int i=j;i<=(n-1)-(k-j);i++)
38             {
39                 double dis=(double)d[i]-messc;
40                 now+=(dis*dis);
41                 if(j && now>=minv) break;
42             }
43             if(!j)minv=now;
44             else if(now<minv) minv=now;
45         }
46         printf("%.12f\n",minv);
47     }
48     return 0;
49 }