DP：地宫取宝（摘花生+最长上升子序列变形）

👉原题链接

题意：

X 国王有一个地宫宝库，是 n×m 个格子的矩阵，每个格子放一件宝贝，每个宝贝贴着价值标签。

地宫的入口在左上角，出口在右下角。

小明被带到地宫的入口，国王要求他只能向右或向下行走。

走过某个格子时，如果那个格子中的宝贝价值比小明手中任意宝贝价值都大，小明就可以拿起它（当然，也可以不拿）。

当小明走到出口时，如果他手中的宝贝恰好是 k 件，则这些宝贝就可以送给小明。

请你帮小明算一算，在给定的局面下，他有多少种不同的行动方案能获得这 k 件宝贝。

输入格式：
第一行 3 个整数，n,m,k，含义见题目描述。

接下来 n 行，每行有 m 个整数 Ci 用来描述宝库矩阵每个格子的宝贝价值。

输出格式：
输出一个整数，表示正好取 k 个宝贝的行动方案数。

该数字可能很大，输出它对 1000000007 取模的结果。

数据范围：
1≤n,m≤50,
1≤k≤12,
0≤Ci≤12

输入样例1：
2 2 2
1 2
2 1
输出样例1：
2

输入样例2：
2 3 2
1 2 3
2 1 5
输出样例2：
14

---

这题限制条件较多：

1. 只能向下或者向右走；👉dp数组需要记录坐标位置
2. 宝物价值需要越来越大才能拿；👉需要记录最后一次拿的宝物的价值
3. 拿满k件；👉需要记录件数
   由此得出需要四位数组来转移状态（分析时间复杂度发现很契合）
   根据闫式dp分析法：
   
   注意：dp问题一般都是从最后一步结果倒推来思考的。

#include <bits/stdc++.h>
using namespace std;

const int N=55,MOD=1000000007;
int n,m,k;
int w[N][N],f[N][N][13][14];//f表示所有从起点走到(i,j)，且已经取了k件物品，且最后一件物品价值是C的合法方案集合

int main()
{
    cin>>n>>m>>k;

    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            cin>>w[i][j];
            w[i][j]++;//因为w可以取0，那么起点处不取物品的情况价值要小于0（取-1），为了方便作数组下标，我们把w都加1
        }
    //初始化起点状态
    f[1][1][1][w[1][1]]=1;//取起点处物品
    f[1][1][0][0]=1;      //不取

    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)//枚举各个坐标
        {
            if(i==1 && j==1) continue;//起点已经初始化，不用再进行计算
            for(int u=0;u<=12;u++)//当前已取件数
                for(int v=0;v<=13;v++)//当前（最后）物品价值
                {
                    int &val=f[i][j][u][v];
                    //对于从上到下、从左到右两种路径，不取(i,j)处物品的情况
                    val=(val+f[i-1][j][u][v])%MOD;//注意，因为题意要求对1e9+7取模，三倍它会爆int，所以这里最多只能同时两个数相加
                    val=(val+f[i][j-1][u][v])%MOD;

                    if(u>0 && v == w[i][j])//当取件数大于0，且当前(最后)物品价值等于C(也就是w[i][j])时
                        for(int c=0;c<v;c++)//枚举取(i,j)处物品的情况
                        {
                            val=(val+f[i-1][j][u-1][c])%MOD;
                            val=(val+f[i][j-1][u-1][c])%MOD;
                        }
                }
        }

    int res=0;
    for(int i=1;i<=13;i++)//把从起点走到(n,m)处且取了k件的情况统计一下
        res=(res+f[n][m][k][i])%MOD;

    cout<<res<<endl;

    return 0;
}