墓地雕塑

[题目描述]

在一个周长为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;
	
}