「动态规划专题」

区间dp

P1220 关路灯
题意：有n个点，有各自的位置且按递增顺序给出，有各自的功率。sb从第v个路灯出发，速度是1m/s。求最小电量消耗下关掉所有路灯。
见解：
从这个题来看区间dp，一般都会面临在区间的哪一头转移过来的问题。而当前的位置也有两种（左端和右端）;
把状态分为两类

1. 转移状态

2. 目前状态

   切忌一次走多步;要保证状态是连续可推的。
   解题一般步骤：

   1 . 所有当前 状态
   2 . 所有转移状态。
   3.从矩阵来看，就是从上一行推到下一行
   再看这个 题：当sb在某个点时，他面临两个选择，往左or往右。dp一般是推当前状态，所以要考虑当前状态从哪里转移过来。现在假设一个区间[i,j]的最小消耗已知，推到 下一个状态
   一.可以是sb站在i点往左走，或者网右走。
   二.可以是sb站在j点往左或者往右。
   推出的状态也可以是区间的左端点也可以是右端点。
   如果现在假设f[i][j] 是目前要推的，把他分成在i点[0],和在j点[1];就有
   f[i][j][0] = min(从i+1走到i，从j走到i，（当然还有从j-1走到j再走到i，显然劣于第二种，还有从i走到j再走回i，显然更长）)。
   f[i][j][1] = min(从j-1走到j，从i走到j，（还有两种显然傻逼的走法。。）)。

DAG上动态规划「最长，最短路」

点集配对问题

最长上升/公共子序列（LIS）（LCS）

最大连续和

货郎担 问题（TSP)

双调欧几里得旅行商问题

——货郎担问题的一种简化形式。
旅行商问题描述：平面上n个点，确定一条连接各点的最短闭合旅程。这个解的一般形式为NP的（在多项式时间内可以求出）
解法：

• （1）首先将各点按照x坐标从小到大排列，时间复杂度为O(nlgn)。
• （2）寻找子结构：定义从Pi到Pj的路径为：从Pi开始，从右到左一直到P1，然后从左到右一直到Pj。在这个路径上，会经过P1到Pmax(i,j)之间的所有点且只经过一次。在定义d(i,j)为满足这一条件的最短路径。我们只考虑i>=j的情况。同时，定义dist(i,j)为点Pi到Pj之间的直线距离。
• （3）最优解：我们需要求的是d(n,n)。 关于子问题d(i,j)的求解，分三种情况：
  A、当j < i - 1时，d(i,j) = d(i-1,j) + dist(i - 1,i)。
  由定义可知，点Pi-1一定在路径Pi-Pj上，而且又由于j < i -1,因此Pi的左边的相邻点一定是Pi-1.因此可以得出上述等式。
  B、当j = i - 1时，与Pi左相邻的那个点可能是P1到Pi-1总的任何一个。因此需要递归求出最小的那个路径：
  d(i,j) = d(i,i-1) = min{d(k,j) + dist(i,k)},其中1 <= k <= j。
  C、当j=i时，路径上最后相连的两个点可能是P1-Pi、P2-Pi…Pi-1-Pi。
  因此有：
  d(i,i) = min{d(i,1)+dist(1,i),…,d(i,i-1),dist(i-1,i)}.。
  一般解法：

double DP(tagPoint *points,int n)
{
    double b[MaxLen][MaxLen];//记录最短路径的长度
   //计算所有情况下的b[i][j]，1 <= i <= j
   //初始化
   b[1][2] = distance(points,1,2);
   for (int j = 3;j <= n;++j)
  {
     //i < j-1
     for (int i = 1;i <= j - 2;++i)
    b[i][j] = b[i][j - 1] + distance(points,j - 1,j);
    //i = j - 1,b[i][j] = min(b[k][j - 1] + distance(k,j));
    b[j - 1][j] = MaxVal;
   for (int k = 1;k <= j - 2;++k)
   {
      double temp = b[k][j - 1] + distance(points,k,j);
      b[j - 1][j] = min(temp,b[j - 1][j]);
   }
  }
   b[n][n] = b[n - 1][n] + distance(points,n - 1,n);
   return b[n][n];
}

UVA1347 Tour

大意：平面上一堆点按x递增的顺序给出，求除终起点其他点经过一次走的欧几里的距离最小。
题目价值 ：教你走最短的圈圈路径。

code

int main() {
  while(scanf("%d", &n) == 1) {
    for(int i = 1; i <= n; i++)
      scanf("%lf%lf", &x[i], &y[i]);
    for(int i = 1; i <= n; i++)
      for(int j = 1; j <= n; j++)
        dist[i][j] = sqrt((x[i]-x[j])*(x[i]-x[j]) + (y[i]-y[j])*(y[i]-y[j]));
    for(int i = n-1; i >= 2; i--)
      for(int j = 1; j < i; j++) {
        if(i == n-1) d[i][j] = dist[i][n] + dist[j][n]; 
        else d[i][j] = min(dist[i][i+1] + d[i+1][j], dist[j][i+1] + d[i+1][i]);
      }
    printf("%.2lf\n", dist[1][2] + d[2][1]);
  }
  return 0;
}

矩阵链乘 （MCM）

最优二叉排序树 （OBST）