动态规划_DAG模型

DAG(有向无环图)上的动态规划是学习动态规划的基础。

• 最长路及其字典序
• 固定终点的最长路和最短路

【嵌套矩形问题】

有n个矩形，每个矩形可以用a,b来描述，表示长和宽。矩形X(a,b)可以嵌套在矩形Y(c,d)中当且仅当a < c,b < d或者b < c,a < d（相当于旋转X90度）。例（1,5）可以嵌套在（6,2）内，但不能嵌套在（3,4）中。你的任务是选出尽可能多的矩形排成一行，使得除最后一个外，每一个矩形都可以嵌套在下一个矩形内。

得动态转移方程：d(i)=max{d(j)+1 | (i,j)∈E}
其中E为边集。求解：

//现假设邻接矩阵已求出,存放在矩阵G中。int dp(int i) {    int& ans = d[i];//引用,方便对于d[i][j][k][l][m][n]的书写    if(ans > 0) return ans;//若已赋值,无需再次计算    ans = 1;    for(int i=0; j<=n; j++)        if(G[i][j])/*以构建DAG*/ ans = max(ans, dp(j)+1);    return ans;//打印输出void print_ans(int i) {    printf("%d", i);    for(int j=1; j<=n; j++)        if(G[i][j] && d[i]==d[j]+1) {            print_ans(j);            break;        }}

---

【硬币问题(完全背包问题)】

n种硬币，面值分别为V1,V2,…Vn,每种都有无限多。给定非负整数S，可以选用多少个硬币，使得面值之和恰好为S？输出硬币数目的最小值和最大值。

//最长路的代码 | 【初始化】mamset(d, -1, sizeof(d)); | (有瑕疵)int dp(int S) {    int& ans = d[S];    if(ans>=0) return ans;    ans = 0;    for(int i=0; i<=n; i++)        if(S>=V[i]) ans = max(ans, dp(S-V[i])+1);    return ans;}//【BUG】S不一定可到达0

//【修正】int dp(int S) {    int& ans = d[S];    if(ans!=-1) return ans;//已赋值    ans = -(1<<30);    for(int i=0; i<=n; i++)        if(S>=V[i]) ans = max(ans, dp(S-V[i])+1);    return ans;}//【TIPS】用特殊值(如-1)表示"未算过",则必须将其和其他特殊值(如无解)//区分开。求最大值时最好将初值设为"无穷小"。

//【可读性优化】用vis数组记录访问状态。用空间代价增强代码可读性，减少代码出错可能性。int dp(int S) {    if(vis[S]) return d[S];    vis[S] = 1;    int &ans = d[S];    ans = -(1<<30);    for(int i=0; i<=n; i++)        if(S>=V[i]) ans = max(ans, dp(S-V[i])+1);    return ans;}

---

【递推】注意计算顺序和边界处理

minv[0] = maxn[0] = 0;for(int i=1; i<=S; i++) {    minv[i] = INF;    maxv[i] = -INF;}for(int i=1; i<=S; i++)    for(int j=1; j<=n; j++)        if(i>=V[j]) {            minv[i] = min(minv[i], minv[i-V[j]]+1);            maxv[i] = max(maxv[i], maxv[i-V[j]]+1);        }printf("%d %d\n", minv[S], maxv[S]);

输出字典序最小的方案(状态的可逆)

//【递归打印】void print_ans(int* d, int S) {    for(int i=1; i<=n; i++)        if(S>=V[i] && d[S]==d[S-V[i]]+1) {            printf("%d ", i);            print_ans(d,  S-V[i]);            break;        }   }

//【递推打印】void print_ans(int* d, int S) {    while(S) {        printf("%d ", d[S]);        S -= V[d[S]];    }}for(int i=1; i<=S; i++)    for(int j=1; i<=n; j++)        if(i>=V[j]) {            if(min[i]>min[i-V[j]]+1) {                min[i]=min[i-V[j]+1;//更新min,寻找最小的min                min_coin[i]=j;//选择第j种硬币            }            if(max[i]<max[i-V[j]]+1) {                max[i]=max[i-V[j]+1;                max_coin[i]=j;            }        }print_ans(min_coin, S);print_ans(max_coin, S);

---

【对于某些不可逆的问题】

• “填表法”。对于每个状态i，找到f(i)依赖的所有状态。
• “刷表法”。对于每个状态i，更新f(i)所影响到的状态。