[题目描述]
在一个周长为10000的圆上等距分布着n个雕塑。现在又有m个雕塑加入(位置可以随便放),希望所有n+m个雕塑在圆周上均匀分布。这就要移动其中一些原有的雕塑。要求n个雕塑移动的总距离最小。
[输入格式]
包含多组数据,每组仅一行,包含两个整数,n (2<=n<=1000),m(1<=m<=1000),即原始雕塑数和新加雕塑数。
[输出格式]
输出仅一行,精确到10^-4
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;
const double NR=10000.0;
int n,m;
double ans;
double a[10001];//n个
double b[10001];//n+m个
double jdz(double n)//double的abs
{
if(n>=0) return n;
else return -n;
}
double mins(double a,double b)//double的min
{
if(a>=b) return b;
else return a;
}
int main()
{
while(scanf("%d%d",&n,&m)!=-1) // 按要求读入
{
//清空数组和变量
ans=0;
memset(a,0,sizeof(a));
memset(b,0,sizeof(b));
a[1]=0;
for(int i=2;i<=n;i++)
{
a[i]=a[i-1]+NR/n;//算出a[i]与a[1]的圆周距离
}
b[1]=0;
for(int i=2;i<=n+m;i++)
{
b[i]=b[i-1]+NR/(n+m);//算出b[i]与b[1]的圆周距离
}
for(int i=1;i<=n;i++)
{
double min_s=99999999.9;
for(int j=1;j<=n+m;j++)
{
min_s=mins(min_s,jdz(a[i]-b[j]));//贪心算法算差
}
ans+=min_s;
}
//按要求输出
int p;
if(ans==(int)ans) p=1;
else if(ans*10==(int)ans*10) p=2;
else if(ans*100==(int)ans*100) p=3;
else p=4;
switch(p)
{
case 1:printf("%.1f\n",ans); break;
case 2:printf("%.2f\n",ans); break;
case 3:printf("%.3f\n",ans); break;
case 4:printf("%.4f\n",ans); break;
}
}
return 0;
}