题意:
把一条高速公路看成一个数轴,在此高速公路上有V个村庄。现知道每个村庄在数轴上的位置,要需要选择P个村庄在其中建立邮局,要求每个村庄到最近邮局的距离之和最小,求最小距离和。
思路:
利用动态规划算法求解。
假设数组dp[i][j]表示前j个村庄,建i个邮局的最短距离和;
数组cost[i][j]表示从村庄i到村庄j建一个邮局,所得最短距离。
当从村庄i到村庄j之间只建一个邮局时,将邮局建在(i + j)/2处所得距离和最短。
若村庄之间要建多个邮局,建设在j个村庄建i个邮局。则前j个村庄建i个邮局的最短距离和等于前k个村庄建i - 1个邮局的最短距离与第 k + 1 到第 j 个村庄建一个邮局的最短距离之和。
得转移方程为:
dp[i][j] = min(dp[i][j], dp[i - 1][k] + cost[k + 1][j]);
代码如下:
#include <iostream>
using namespace std;
#define MAX_V 310 //最大村庄数
#define MAX_P 35 //最大邮局数
#define MAXINT 1000000
void Init(int num_v,int pos[],int cost[MAX_V][MAX_V],int dp[MAX_P][MAX_V])
{
//初始化
int Ppos = 0; //邮局坐标
for(int i = 0; i < num_v; i++)
{
for(int j = i; j < num_v; j++)
{
Ppos = pos[(i + j)/2];
for(int k = i; k <= j; k++)
{
cost[i][j] = cost[i][j] + abs(pos[k] - Ppos); //i到j,建一个邮局,最短距离
}
}
}
for(int j = 0; j < num_v; j++)
{
dp[0][j] = cost[0][j];
}
}
void DP(int num_p,int num_v,int pos[],int cost[MAX_V][MAX_V],int dp[MAX_P][MAX_V])
{
//动态规划求解
for(int i = 1; i < num_p; i++)
{
for(int j = i; j < num_v; j++)
{
dp[i][j] = MAXINT;
for(int k = i; k < j; k++)
{
if(dp[i][j] > dp[i - 1][k] + cost[k + 1][j])
{
dp[i][j] = dp[i - 1][k] + cost[k + 1][j];
}
}
}
}
}
void Output(int num_p,int num_v,int dp[MAX_P][MAX_V])
{
cout<<dp[num_p - 1][num_v - 1]<<endl;
}
int main()
{
int dp[MAX_P][MAX_V] = {0}; //存储前j个村庄,建i个邮局的最短距离
int cost[MAX_V][MAX_V] = {0}; //存储存储从第i个村庄到第j个村庄建立一个邮局的最短距离
int pos[MAX_V] = {0}; //存储村庄位置
int num_v = 0; //村庄个数
int num_p = 0; //邮局个数
cin>>num_v>>num_p;
for(int i = 0; i < num_v; i++)
{
cin>>pos[i];
}
Init(num_v,pos,cost,dp);
DP(num_p,num_v,pos,cost,dp);
Output(num_p,num_v,dp);
return 0;
}