polygon

描述：

在一个周长为10000的圆上等距分布着n个点，即这n个点是一个正n边形的顶点。现在要另加m个点到圆上，新加的m个点可以任意选择位置（可以与原有的点重合）。然后将这n+m个点中的一些点延圆周移动，最终使n+m个点均匀分布，即在一个正n+m边形的顶点上。输出最小总移动距离。

输入：

输入两个整数 n, m。 (2≤n≤1000, 1≤m≤1000).

输出：

输出最小总移动距离，保留4位小数。

输入样例：

sample input #1 2 1 sample input #2 2 3 sample input #3 3 1 sample input #4 10 10

输出样例：

sample output #1 1666.6667 sample output #2 1000.0 sample output #3 1666.6667 sample output #4 0.0 图对应前3个样例

#include <iostream>#include <cmath>#include <iomanip>using namespace std;int n, m;double Min(double a, double b){return (a < b ? a : b);}double FindSimilarOrEqual(double mod[], double ori[]){int begin = 0, end = n + m - 1, mid;double sum = 0.0;for (int i = 1; i < n; ++i){while (begin <= end){mid = (begin + end) / 2;if (mod[mid] == ori[i])break;else if (mod[mid] > ori[i])end = mid - 1;elsebegin = mid + 1;}if (begin > end)sum += Min(fabs(mod[begin] - ori[i]), fabs(mod[end] - ori[i]));begin = 0;end = n + m - 1;}return sum;}int main(void){int i;double ori_part, mod_part;cin >> n >> m;double *ori = new double[n];double *mod = new double[n + m];ori_part = 10000 / (double)n;mod_part = 10000 / (double)(n + m);for (i = 0; i < n; ++i)ori[i] = i * ori_part;for (i = 0; i < n + m; ++i)mod[i] = i * mod_part;cout << setiosflags(ios::fixed) << setprecision(4);cout << FindSimilarOrEqual(mod, ori) << endl;delete []ori;delete []mod;return 0;}