2014年第五届蓝桥杯C/C++本科B组省赛 地宫取宝

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

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

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

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

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

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

【数据格式】

    输入一行3个整数，用空格分开：n m k (1<=n,m<=50, 1<=k<=12)

    接下来有 n 行数据，每行有 m 个整数 Ci (0<=Ci<=12)代表这个格子上的宝物的价值

    要求输出一个整数，表示正好取k个宝贝的行动方案数。该数字可能很大，输出它对 1000000007 取模的结果。

例如，输入：
2 2 2
1 2
2 1
程序应该输出：
2

再例如，输入：
2 3 2
1 2 3
2 1 5
程序应该输出：
14

资源约定：
峰值内存消耗 < 256M
CPU消耗  < 1000ms

请严格按要求输出，不要画蛇添足地打印类似：“请您输入...” 的多余内容。

所有代码放在同一个源文件中，调试通过后，拷贝提交该源码。

注意: main函数需要返回0
注意: 只使用ANSI C/ANSI C++ 标准，不要调用依赖于编译环境或操作系统的特殊函数。
注意: 所有依赖的函数必须明确地在源文件中 #include <xxx>， 不能通过工程设置而省略常用头文件。

提交时，注意选择所期望的编译器类型。
 

思路：利用四维dp记录每种情况的种类数，第一维记录位置的横坐标，第二维记录位置的纵坐标，第三维记录取了几件宝物，第四维记录取出的宝物的最大价值。之后找到状态转移方程。因为宝物的价值可以为0，所以先将每件宝物价值加1.

状态转移方程：

没有取当前位置的宝物时 ：dp[i][j][p][q]=(dp[i-1][j][p][q]+dp[i][j-1][p][q])%mod;

取当前位置的宝物时：dp[i][j][p][s[i][j]]=(dp[i][j][p][s[i][j]]+dp[i-1][j][p-1][q]+dp[i][j-1][p-1][q])%mod;   （q<s[i][j]）

AC代码：

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<map>
#include<queue>
#include<stack>
#include<string>
#include<iostream>
#include<algorithm>
using namespace std;
#define mod 1000000007
int s[60][60];
long long dp[60][60][15][15];
int main()
{
    int n,m,k;
    scanf("%d%d%d",&n,&m,&k);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
            {
                scanf("%d",&s[i][j]);
                s[i][j]+=1;
            }
    }
    memset(dp,0,sizeof(dp));
    dp[1][1][1][s[1][1]]=1;
    dp[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 p=0;p<=k;p++)
                for(int q=0;q<=13;q++)
                    dp[i][j][p][q]=(dp[i-1][j][p][q]+dp[i][j-1][p][q])%mod;//不取宝物的状态
            for(int p=1;p<=k;p++)
            {
                for(int q=0;q<s[i][j];q++)
                {
                    dp[i][j][p][s[i][j]]=(dp[i][j][p][s[i][j]]+dp[i-1][j][p-1][q]+dp[i][j-1][p-1][q])%mod;//取宝物的状态
                }
            }
        }
    }
    long long sum=0;
    for(int i=0;i<=13;i++)
    {
        sum=(sum+dp[n][m][k][i])%mod;//求出口有k件物品的种类
    }
    printf("%lld\n",sum);
}