例题4  墓地雕塑(Graveyard, NEERC 2006, LA 3708)

在一个周长为10000的圆上等距分布着n个雕塑。现在又有m个新雕塑加入(位置可以随意放),希望所有n+m个雕塑在圆周上均匀分布。这就需要移动其中一些原有的雕塑。要求n个雕塑移动的总距离尽量小。

【输入格式】

输入包含若干组数据。每组数据仅一行,包含两个整数nm(2≤n≤1 000,1≤m
≤1 000),即原始的雕塑数量和新加的雕塑数量。输入结束标志为文件结束符(EOF)。

【输出格式】

输入仅一行,为最小总距离,精确到10-4

【样例输入】

2 1

2 3

3 1

10 10

【样例输出】

1666.6667

1000.0

1666.6667

0.0

这道题 是一边写博客 一边思考的一道题

首先 显然发现 +了雕像后 可以保证一个雕像始终不会动

这是 10个雕像  变为11个雕像的时候 显然发现因为11个雕像的时候间距小于10个雕像的时候

使得在(11)红线 前面的(10)白线 选择移动到前方更近

在(11)红线 后方的(10)白线 选择移动到后方更近

显然

在m=1的情况下  选择附近最近的搬移位置即可 且能保证不会有2个雕像搬到同一位置

考虑m=2 :

10个雕像 变为12个雕像的时候

3个红色的线就是几条类似上面11时候红色的 分界线

我们可以也发现。。  搬移到附近最近位置即可;

。。。但是只能从图片上看出

严格证明想了一段时间还是没有想出来

好吧 似乎通过严格证明再得出  搬到附近最近的位置这个结论估计是不行的。

只能通过 搬到附近最近的位置这个假设 证明出是最优解 并且证明出不会搬移到同一点(书上解法)

书上面有讲解  不过因为是直接看书得出的题解  我还是过两天自己仿书上写一遍题解吧 (今天10-10)

10-12 再次编辑:

step1:

为了使操作简便  假设长度为N+M(后面知道这一步是这道题的精髓)

step2:

原雕塑位置为pos=((N+M)/N )*i         (i=0...N-1)

step3:

根据我们上面的图的操作我们知道要找距离最近的点

(若不搬到同一处 显然是最优的,但不知道是否正确)

那么该如何找最近的点呢?

这就体现了step 1的优越性了

floor(pos+0.5)  即可

四舍五入即为一个寻找最近的整数点的方法不是吗?

新雕塑的位置全部为一个整数点;

所以是一个十分赞的方法

程序流程到此结束

证明:

为何这样是正确的?证明可行性

反证法:

假设有两个原雕塑搬到同一点 比如 2

那么 两个雕塑相距最远的可能是什么?

1.49999999999  与 2.49999999999

即使这样两者距离也不过是等于1

显然可知原雕塑距离一定要大于1的

固不可能搬到同一点(这也是当时缩放的一个好处,便于推理)

小贴士:

如果考试的时候无法证明这一点怎么办?

不如从1 for 到 10000 打表证明不存在这种情况

打表程序如下:

AC代码如下:

#include<cstdio>
#include<cmath>
using namespace std;
int main()
{
int n,m;
double pos,ans;
while(scanf("%d%d",&n,&m)!=EOF)
{
ans=0;
for(int i=0;i<n;i++)
{
pos=(double)(n+m)/n*i;
ans=ans+abs(floor(pos+0.5)-pos);
}
printf("%.4lf\n",ans*10000/(double)(n+m));
}
return 0;
}

最新文章

  1. Flume_企业中日志处理
  2. iOS-Runtime知识点整理
  3. Java—输入输出技术
  4. [译]学习IPython进行交互式计算和数据可视化(二)
  5. c++语法随身记
  6. linux epoll 开发指南-【ffrpc源码解析】
  7. OC基础--内存管理中的@property关键字以及其参数
  8. mysql数据库之基础SQL语句/语法
  9. [备忘]Visio中连接线交叉时跨线小弯的去掉方法
  10. java.lang.IllegalAccessError: tried to access method org.apache.poi.util.POILogger.log from class org.apache.poi.openxml4j.opc.ZipPackage
  11. Efficient&amp;Elegant:Java程序员入门Cpp
  12. 简单记录numpy库的某些基本功能
  13. perl 里面如何写出阅读友好的代码提示
  14. 去除inline-block出现间距的几种方法
  15. [设计模式][c++]状态切换模式
  16. 【明哥报错簿】之json转换报错---net.sf.ezmorph.bean.MorphDynaBean cannot be cast to XXXDO
  17. am335x -- led 控制
  18. 【计算几何】【bitset】Gym - 101412G - Let There Be Light
  19. 【BZOJ 2324】[ZJOI2011]营救皮卡丘 费用流
  20. 在使用Arduino中遇到的问题(无法使用中文注释、程序无法下载)

热门文章

  1. [转]Laravel 4之Eloquent ORM
  2. Coding.net代码托管空间申请与使用-安装并运行WordPress博客
  3. 关于js中select的简单操作,以及js前台计算,span简单操作
  4. Ubantu 命令
  5. Silverlight调用网站项目的Session
  6. 上拉、下拉UITableView,交互式 模态弹出(自定义弹出动画)
  7. ios 获取屏幕的属性
  8. 【HTML5】在head 设置 meta 能更方便开发
  9. android启动activity文本框不获得焦点
  10. mvn多模块开发消除重复依赖造成的打包失败