浅析一个用PD解决的邮局选址问题

最新推荐文章于 2019-04-29 20:38:44 发布

weixin_30338461

最新推荐文章于 2019-04-29 20:38:44 发布

阅读量141

点赞数

文章标签：数据结构与算法

原文链接：http://www.cnblogs.com/norx/articles/1004564.html

版权

Preface （前言）

　　动态规划（Dynamic Programming）是一种算法设计技术，它有着相当有趣的历史。作为一种使多阶段决策过程最优的通用方法，它是在20世纪50年代由一位卓越的美国数学家Richard Bellman所发明的。因此，这个技术名字中的“programming”是计划和规划的意思，不是代表计算法中的编程。它作为一种重要的工具在应用数学中的价值被大家认同以后，起码在计算机科学的圈子里，人们不仅用它来解决特定类型的最优问题，而且最终把它作为一种通用的算法设计技术来使用。在这里，我们正是从这个角度来考虑这种技术的。

　　如果问题是由交叠的子问题所构成的，我们就可以用动态规划技术来解决它。一般来说，这样的子问题出现在对给定问题求解的递推关系中，这个递推关系中包含了相同类型的更小子问题的解。动态规划法建议，与其对交叠的子问题一次又一次地求解，还不如对每个较小的子问题只求解一次并把记录记录在表中，这样就可以从表中得出原始问题的解。

　　——《算法设计与分析基础》

Description （问题描述）

　　有一条公路经过Ｖ个村庄，每一个村庄都处在整数的坐标点上（这里假设公路拉直为Ｘ轴）．规划在这条公路上建立Ｐ个邮局，当然为了方便，这些邮局应建在某Ｐ个村庄上，但是要求让不同村庄的人到邮局要走的总路程最小．

Input

先从键盘读入两个整数Ｖ和Ｐ，然后再读入Ｖ个整数，分别表示Ｖ个村庄的坐标（坐标＞＝０）

Output

输出Ｐ个以空格分隔的整数，按坐标从小到的顺序给出Ｐ个邮局的坐标．

Sample Input

10 3

1 4 7 19 70 89 105 204 18 40

Sample Output

7 89 204

Analysis（分析）

一、问题分析

在解决本问题之前，我们有必要明确以下几点：

1、对于复杂的问题，我们有必要先从考虑简单情况出发，并由简单情况找出规律，发现思路．

2、当有V个村庄点，且只有一个邮局点时，邮局点的最佳选址位置为这V个村庄点的中位点．

3、当有V个村庄点，并且不只一个邮局点时，我们发现：此问题可以分解为多个子问题进行计算，并且邮局选址的最优解包含子问题的最优解，即这个邮局选址问题具有最佳子结构性质！

4、在我们容易地列出此问题的递归表达式（即后文即将提及的表达式1.1）之后，我们发现：在递归计算的过程中，许多子问题被重复计算了多次！

× 基于第3第4点发现，我们选择动态规划的方法来求解本题．

二、算法分析

（一）、找出递归表达式

根据从特殊到一般的思路，我们先考虑当邮局数为1时的情况，再考虑邮局数大于1的情况．

1、当邮局数P=1时，易知邮局的最佳选址位置为V个村庄的中位点．

2、当邮局数P>1时，我们可以想象每个邮局对应一个辖区，那么问题就可以被转化为如何将V个村庄划分为连续的P个辖区．为此，我们做以下定义：

（1）、Center(l,r)，表示当村庄x_l到x_r为一个辖区时，邮局所在位置．

（2）、Dis(l,r)，表示当村庄x_l到x_r为一个辖区时，该辖区内村庄到邮局的最小距离．

（3）、TotolDis(t,k)，表示将村庄x_t之后的所有村庄（包括村庄x_l）划分为k个辖区时，所有辖区的Dis(l,r) 总和的最小值，题目就转化为求TotolDis(t,k)，有

下面，我们还是通过从特殊到一般的方法，考虑如何求TotolDis(t,k)．

举特例：“在V=7,P=3时，目标是求TotolDis(t,k)．”

显而易见，TotolDis(t,1)=Dis(t,6)，t=2,3,4,5,6；也就是说TotolDis(1,2)=min{F,G,H,I,J}，是可以在原有子问题上计算出来的．
　　对于情况B、C、D、E，同理．
根据上述特殊，推广并且改造表达式（1.1），可得：

显然，这是一个典型的递归过程．

（二）、构造动态规划的程序结构

在以上定义Center(l,r)、Dis(l.r)、TotalDis(t,k)的基础上，构造合理的数据结构和控制结构，自底向上的求解这个问题还是先考虑上述“V=7，P=3，村庄坐标为0～6”的特例，Center(l,r)表和Dis(l,r)表如下所示．

l"r	0	1	2	3	4	5	6	l"r	0	1	2	3	4	5	6
0	0	0	1	1	2	2	3	0	0	1	2	4	6	9	12
1		1	1	2	2	3	3	1		0	1	2	4	6	9
2			2	2	3	3	4	2			0	1	2	4	6
3				3	3	4	4	3				0	1	2	4
4					4	4	5	4					0	1	2
5						5	5	5						0	1
6							6	6							0
center[l][r] 图1-1　特例的Center(l,r)表								dis[l][r] 图1-2　特例的Dis(l,r)表

据此，我们可以很容易的编写出计算center[l][r]和dis[l][r]的代码模块．

void calculateCenter()

{
2

//--- 内存分配
4

center = new unsigned *[village]; //动态分配二维数组的第一维
6

for (unsigned i=0; i<village; i++) //动态分配二维数组的第二维
8

center[i] = new unsigned[village];
10

//--- 初始化Center(l,r)
14

for (unsigned l = 0; l < village; l++)
16

for (unsigned r = l; r < village; r++)
18

center[l][r] = xCoordinate[(r-l)/2 + l];
20

}
22

calculateCenter()函数

void calculateDis()

{
2

//--- 内存分配
4

dis = new unsigned *[village]; //动态分配二维数组的第一维
6

for (unsigned i = 0; i < village; i++) //动态分配二维数组的第二维
8

dis[i] = new unsigned[village];
10

//--- 初始化Dis(l,r)
14

for (unsigned l = 0; l < village; l++)
16

for (unsigned r = l; r < village; r++)

{
18

dis[l][r] = 0;
20

for (unsigned k = l; k <= r; k++)
22

if (center[l][r] > xCoordinate[k])
24

dis[l][r] += center[l][r] - xCoordinate[k];
26

else
28

dis[l][r] += xCoordinate[k]- center[l][r];
30

}
32

}
34

calculateDis()函数

在计算TotalDis(0,P)时，由于题目要求输出邮局的位置，所以在计算TotalDis(0,P)的过程中，只计算最小值还不够，还必须记录下取得最小值时的分辖区情况．于是，我们这么定义TotalDis二维数组．

1 struct info{
2
3 unsigned dis; //最小值
4
5 unsigned r; //标号r之前（包括r）的村庄为一个辖区
6
7};
8
9
10
11info **totalDis=NULL;
12

totalDis表的元素结构

计算过程如下所示．

void calculateTotalDis()

{
2

//--- 内存分配
4

totalDis = new info *[village]; //动态分配二维数组的第一维
6

for (unsigned i = 0; i < village; i++) //动态分配二维数组的第二维
8

totalDis[i] = new info[postoffice+1];
10

//--- 计算TotalDis(v,p+1)
14

//---- 当k=1时，根据公式(1.2)，直接计算
16

for (unsigned t = 0; t < village; t++)
18

totalDis[t][1].dis = dis[t][village-1];
20

//---- 当k=2,3,

,p时的情况
22

for (unsigned k = 2; k <= postoffice; k++)
24

for (unsigned t = 0; t < village; t++)

{
26

totalDis[t][k].dis = (unsigned)-1;
28

totalDis[t][k].r = 0;
30

for (unsigned r = t; r <= village-k; r++)

{
32

unsigned temp = dis[t][r] + totalDis[r+1][k-1].dis;
34

//---- 计算最小值
36

if (temp < totalDis[t][k].dis)

{
38

totalDis[t][k].dis = temp;
40

totalDis[t][k].r = r;
42

}
44

}
46

}
48

}
50

calculateTotalDis()函数

得到的TotalDis(t,k)表如下所示．

对于本题，还有一个地方需要注意，那就是输出．由于题目要求输出各邮局的坐标，在这里需用到递归的方法，具体输出模块如下所示．

1 void output(unsigned t, unsigned k){
2
3       if (1 == k)
4
5              cout << center[t][village-1];
6
7       else{
8
9              cout << center[t][totalDis[t][k].r] << ' ';
10
11              output(totalDis[t][k].r+1, k-1); //递归
12
13       }
14
15}
16