Galaxy

Problem's Link:   http://acm.hdu.edu.cn/showproblem.php?pid=5073


Mean:

在一条数轴上,有n颗卫星,现在你可以改变k颗卫星的位置,使得剩下的n-k颗卫星到某个点(不固定)的距离的平方和最小。

抽象成数学语言后等价于:数轴上有n个点,现在去掉k个点,使得剩下的n-k个点的方差最小,求方差*n的值。

analyse:

一道让人很容易想偏的数学题。首先说一下我的思路:

1)我们最终的目的是让这n-k个点尽量的集中,所以去掉的这k个点必须是位于两边的点(想不通的请自行补脑);

2)剩下的事情就是枚举两边的数量了,但是在枚举这一步,怎样才不超时呢?咳咳,这题的关键来了。

我们可以先来推一下公式:

设Fn为这n个数的方差,d为这n个数的平均数,那么:

Fn=[(x1-d)^2+(x2-d)^2+......(xn-d)^2]/n;

 =[sum(xi^2)+n*d*d-2*d*sum(xi)]/n;

根据这个公式来枚举前后个数就简单多了,详见代码。

Time complexity: O(n)

Source code: 

//  Memory   Time
// 1347K 0MS
// by : Snarl_jsb
// 2014-11-16-22.59
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<vector>
#include<queue>
#include<stack>
#include<map>
#include<string>
#include<climits>
#include<cmath>
#define LL long long
using namespace std;
#define N 50000+10
double a[N],sum1[N],sum2[N];
int main()
{
int t;
cin>>t;
while(t--)
{
int n,k;
cin>>n>>k;
memset(sum1,0,sizeof sum1);
memset(sum2,0,sizeof sum2);
double tmp1,tmp2;
tmp1=tmp2=0.0;
for(int i=1;i<=n;++i)
scanf("%lf",&a[i]);
sort(a+1,a+1+n);
for(int i=1;i<=n;++i)
{
sum1[i]=sum1[i-1]+a[i];
sum2[i]=sum2[i-1]+a[i]*a[i];
}
if(n==k)
{
puts("0.0000000000000");
continue;
}
int m=n-k; /**< 需要选的人数 */
double d;
double tmp;
int sta,en;
double res=1000000000000000000.0;
// cout<<setprecision(10)<<res<<endl;
for(int i=1;i+m-1<=n;++i)
{
sta=i;
en=i+m-1;
d=(sum1[en]-sum1[sta-1])/m;
tmp=(sum2[en]-sum2[sta-1])+m*d*d-2*d*(sum1[en]-sum1[sta-1]);
if(tmp<res)
{
res=tmp;
}
}
printf("%.9lf\n",res);
}
return 0;
}

  

最新文章

  1. OpenCv编程
  2. [问题2014S11] 复旦高等代数II(13级)每周一题(第十一教学周)
  3. Linux_Centos中搭建nexus私服
  4. Python-面向对象编程
  5. AsyncTask的使用方法和理解
  6. Amoeba:开源的分布式数据库Porxy解决方案
  7. DevExpress 中 WaitForm 使用
  8. 浅析jQuery框架与构造对象
  9. 写一个自己定义进度颜色和圆形转动的ProgressBar(具体介绍)
  10. Uber司机手机终端问答篇
  11. CentOS6.5解压缩文件.tar.gz .war .zip
  12. 【转】C\C++代码优化的27个建议
  13. 浏览器缓存相关HTTP头部字段
  14. HDOJ 6508 Problem I. Spell Boost (01背包/DP)
  15. javascript基础修炼(2)——What&#39;s this(上)
  16. Asynchronous Programming
  17. 第一章 HTML+CSS(上)
  18. ansible系列7-mysql_user模块
  19. Bootstrap风格button
  20. 对Django框架架构和Request/Response处理流程的分析(转)

热门文章

  1. Uploadify v3.2.1 上传图片并预览
  2. js获取gridview模板列中textbox行列的值
  3. 在ps中画两个同心圆并且把两个同心圆进行任意角度切割
  4. Android开发(二十八)——基础功能函数
  5. Linux之重定向
  6. 解决“com.android.dex.DexIndexOverflowException: method ID not in [0, 0xffff]: 65536”问题(l转)
  7. JVM 参数翻译汉化解释
  8. ci配置smarty手记
  9. Android Weak Handler:可以避免内存泄漏的Handler库
  10. 8个经典炫酷的HTML5 Canvas动画欣赏