【题解】「NOIP2010」乌龟棋（DP）

题面

【题目描述】
小明过生日的时候，爸爸送给他一副乌龟棋当作礼物。
乌龟棋的棋盘是一行$N$个格子，每个格子上一个分数（非负整数）。棋盘第$1$格是唯一 的起点，第$N$格是终点，游戏要求玩家控制一个乌龟棋子从起点出发走到终点。
乌龟棋中M张爬行卡片，分成$4$种不同的类型（M张卡片中不一定包含所有4种类型 的卡片，见样例），每种类型的卡片上分别标有$1、2、3、4$四个数字之一，表示使用这种卡 片后，乌龟棋子将向前爬行相应的格子数。游戏中，玩家每次需要从所有的爬行卡片中选择 一张之前没有使用过的爬行卡片，控制乌龟棋子前进相应的格子数，每张卡片只能使用一次。 游戏中，乌龟棋子自动获得起点格子的分数，并且在后续的爬行中每到达一个格子，就得到 该格子相应的分数。玩家最终游戏得分就是乌龟棋子从起点到终点过程中到过的所有格子的 分数总和。 很明显，用不同的爬行卡片使用顺序会使得最终游戏的得分不同，小明想要找到一种卡 片使用顺序使得最终游戏得分最多。 现在，告诉你棋盘上每个格子的分数和所有的爬行卡片，你能告诉小明，他最多能得到 多少分吗？
【输入】
输入的每行中两个数之间用一个空格隔开。
第$1$行$2$个正整数$N$和$M$，分别表示棋盘格子数和爬行卡片数。
第$2$行$N$个非负整数，$a_1,a_2\dots \dots a_N$，其中$a_i$表示棋盘第$i$个格子上的分数。
第$3$行$M$个整数，$b_1，b_2，\dots \dots ，b_M$，表示M张爬行卡片上的数字。 输入数据保证到达终点时刚好用光$M$张爬行卡片
【输出】
【样例输入】

9 5
6 10 14 2 8 8 18 5 17
1 3 1 2 1

【样例输出】

73

【数据范围】
说明/提示
小明使用爬行卡片顺序为$1,1,3,1,2$,得到的分数为$6+10+14+8+18+17=73$。注意，由于起点是$1$，所以自动获得第$1$格的分数$6$。
对于$30\%$的数据有$1\le N\le 30,1\le M\le 12$。
对于$50\%$的数据有$1\le N\le 120,1\le M\le 50$，且$4$种爬行卡片，每种卡片的张数不会超过$20$。
对于$100\%$的数据有$1\le N\le 350,1\le M\le 120$，且$4$种爬行卡片，每种卡片的张数不会超过$40$；$0\le a_i\le 100,1\le i\le N,1\le b_i\le 4,1\le i\le M$。

算法分析

四维DP。
为了表示乌龟其中走到当前的位置，设当前使用了$i$张卡片$1$，j张卡片$2$，$k$张卡片$3$，$t$张卡片$4$。当前到达的格子数就是$s=1+i+2\ast j+3\ast k+4\ast t$，得到第$s$格子的分数。

设状态：
$f[i][j][k][t]$——表示使用了$i$张卡片$1$，$j$张卡片$2$，$k$张卡片$3$，$t$张卡片$4$得到的最高分数。
枚举使用的不同类型卡片的数量，找出最优解。

状态转移方程：
考虑到达当前位置时，使用的最后一张卡片类型：
最后使用的是卡片$1：f[i][j][k][t]=max(f[i][j][k][t],f[i-1][j][k][t]+score[s]);$
最后使用的是卡片$2：f[i][j][k][t]=max(f[i][j][k][t],f[i][j-1][k][t]+score[s]);$
最后使用的是卡片$3：f[i][j][k][t]=max(f[i][j][k][t],f[i][j][k-1][t]+score[s]);$
最后使用的是卡片$4：f[i][j][k][t]=max(f[i][j][k][t],f[i][j][k][t-1]+score[s]);$
最后的答案就是$f[a][b][c][d]，a,b,c,d$分别表示卡片$1，2，3，4$的总数。
时间复杂度：$O(40^4)$

参考程序

#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
int score[360],f[42][42][42][42];
int main(){
    int n,m,x;
    int a=0,b=0,c=0,d=0;
    cin>>n>>m;
    for(int i=1;i<=n;i++)
        cin>>score[i];
    for(int i=1;i<=m;i++)
    {
        cin>>x;
        if(x==1) a++;
        if(x==2) b++;
        if(x==3) c++;
        if(x==4) d++;
    }
    for(int i=0;i<=a;i++)
        for(int j=0;j<=b;j++)
            for(int k=0;k<=c;k++)
                for(int t=0;t<=d;t++)
                {
                    int s=i+j*2+k*3+t*4+1;
                    if(i) f[i][j][k][t]=max(f[i][j][k][t],f[i-1][j][k][t]+score[s]);
                    if(j) f[i][j][k][t]=max(f[i][j][k][t],f[i][j-1][k][t]+score[s]);
                    if(k) f[i][j][k][t]=max(f[i][j][k][t],f[i][j][k-1][t]+score[s]);
                    if(t) f[i][j][k][t]=max(f[i][j][k][t],f[i][j][k][t-1]+score[s]);
                }
    cout<<f[a][b][c][d]+score[1];
    return 0;
}