NOI: Post Office

一、题目大意：

二、解题思路

• 定义：
  • $dp[i][j]$:前$i$个村庄建$j$个邮局的最小距离.
  • $m[i][j]$:村庄$[i,j]$之间建立一个邮局的最小距离.(想象: 很明显应该建在$(i+j)/2$的邮局上)
• 目标状态: $dp[V][P]$.即前$V$个村庄建立$P$个邮局.
• 状态转移：
  • 子问题:
    • 定义：$build(s,e,j):=$在村庄$[s,e]$ 间建$j$个邮局
    • b u i l d ( 1 , i , j ) : = ⋃ k = 1 i − 1 { b u i l d ( 1 , k , j − 1 ) ∩ b u i l d ( k + 1 , i , 1 ) } build(1,i,j):=\bigcup_{k=1}^{i-1} \{build(1,k , j-1)\cap build(k+1, i, 1)\} build(1,i,j):=⋃k=1i−1​{build(1,k,j−1)∩build(k+1,i,1)}
  • 因此有: 状态转移方程: d p [ i ] [ j ] = m a x { d p [ k ] [ j − 1 ] + m [ k + 1 ] [ i ] ∣ 1 &lt; = k &lt; = i − 1 } dp[i][j] = max\{dp[k][j-1]+m[k+1][i]\quad | \quad 1&lt;=k&lt;=i-1\} dp[i][j]=max{dp[k][j−1]+m[k+1][i]∣1<=k<=i−1}（这里默认了每一个子问题下，最优情况为: 前$k$个村庄只往前$j-1$个邮局前进，如果是要往第$j$个邮局，这种情况包含在其他子问题下)因此状态转移最终结果没有问题。
  • $m[i][j]=m[i][j-1]+a[j]-a[(i+j)/2]$(想象可得出)
• 转移策略: 先求出所有的$m$, $dp[i][j]$只依赖与$dp[i^{\prime }][j-1]$，因此只需循环更新下三角矩阵。
• 初始状态:$dp[i][1]=m[1][i]，dp[i][0]=inf$

三、 代码:

#include<iostream>
using namespace std;
const int MAXV = 304;
const int MAXP = 34;
const int inf = 1 << 30;
int dp[MAXV][MAXP];
int m[MAXV][MAXV];
int a[MAXV];
int main()
{
    int V, P;
    cin >> V >> P;
    for(int i=1; i<=V; i++)
        cin >> a[i];
    for(int i=0; i<=V; i++)
        for(int j=0; j<=P; j++)
            dp[i][j] = inf;
    for(int i=1; i<=V; i++)
        m[i][i] = 0;
    for(int i=1; i<=V; i++)
        for(int j=i+1; j<=V; j++)
            m[i][j] = m[i][j-1] + a[j] - a[(i+j)/2];
    for(int i=1; i<=V; i++)
        dp[i][1] = m[1][i];

    for(int i=2; i<=V; i++)
        for(int j=1; j<=i; j++)
            for(int k=1; k<=i-1; k++)
                dp[i][j] = min(dp[i][j], dp[k][j-1]+m[k+1][i]);

    cout << dp[V][P] << endl;

    return 0;
}