题目背景

小明过生日的时候，爸爸送给他一副乌龟棋当作礼物。

题目描述

乌龟棋的棋盘是一行N个格子，每个格子上一个分数（非负整数）。棋盘第1格是唯一的起点，第N格是终点，游戏要求玩家控制一个乌龟棋子从起点出发走到终点。

乌龟棋中M张爬行卡片，分成4种不同的类型（M张卡片中不一定包含所有4种类型的卡片，见样例），每种类型的卡片上分别标有1、2、3、4四个数字之一，表示使用这种卡片后，乌龟棋子将向前爬行相应的格子数。游戏中，玩家每次需要从所有的爬行卡片中选择一张之前没有使用过的爬行卡片，控制乌龟棋子前进相应的格子数，每张卡片只能使用一次。

游戏中，乌龟棋子自动获得起点格子的分数，并且在后续的爬行中每到达一个格子，就得到该格子相应的分数。玩家最终游戏得分就是乌龟棋子从起点到终点过程中到过的所有格子的分数总和。

很明显，用不同的爬行卡片使用顺序会使得最终游戏的得分不同，小明想要找到一种卡片使用顺序使得最终游戏得分最多。

现在，告诉你棋盘上每个格子的分数和所有的爬行卡片，你能告诉小明，他最多能得到多少分吗？

输入输出格式

输入格式：

输入文件的每行中两个数之间用一个空格隔开。

第1行2个正整数N和M，分别表示棋盘格子数和爬行卡片数。

第2行N个非负整数，a1a2……aN，其中ai表示棋盘第i个格子上的分数。

第3行M个整数，b1b2……bM，表示M张爬行卡片上的数字。

输入数据保证到达终点时刚好用光M张爬行卡片。

输出格式：

输出只有1行，1个整数，表示小明最多能得到的分数。

输入输出样例

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

73

说明

每个测试点1s

小明使用爬行卡片顺序为1，1，3，1，2，得到的分数为6+10+14+8+18+17=73。注意，由于起点是1，所以自动获得第1格的分数6。

对于30%的数据有1≤N≤30，1≤M≤12。

对于50%的数据有1≤N≤120，1≤M≤50，且4种爬行卡片，每种卡片的张数不会超过20。

对于100%的数据有1≤N≤350，1≤M≤120，且4种爬行卡片，每种卡片的张数不会超过40；0≤ai≤100，1≤i≤N；1≤bi≤4，1≤i≤M。

        不难发现由于题目要求，只要1,2,3,4卡片的使用数目被规定，那么当前格子也是确定的，那么我们不妨用dp[i][j][k][l]表示1用i张，2用j张，3用k张，4用l张时的最优解。那么dp[i][j][k][l] = max(dp[i-1][j][k][l], dp[i][j-1][k][l], dp[i][j][k-1][l], dp[i][j][k][l-1]) + num[1 + i + j * 2 + k * 3 + l * 4]/*注意是从1开始1为起点*/

        PS:终于在DP中感到了register常数优化的威力，真的服气，优化 30%的时间

#include<cstdio>
using namespace std;
int dp[41][41][41][41];
int num[351], cnt[5];
inline int max(int a, int b) {return a > b ? a : b;}
inline int read() {
    int ret = 0;
    char c = getchar();
    while(c > '9' || c < '0') c = getchar();
    while(c <= '9' && c >= '0') ret = ret * 10 + c - '0', c = getchar();
    return ret;
}
int main() {
    int n, m, a, pos;
    n = read(), m = read();
    for(register int i = 1; i <= n; ++i)
        num[i] = read();
    for(register int i = 1; i <= m; ++i) {
        a = read();
        ++cnt[a];
    }
    for(register int i = 0; i <= cnt[1]; ++i) 
        for(register int j = 0; j <= cnt[2]; ++j)
            for(register int k = 0; k <= cnt[3]; ++k)  
                for(register int l = 0; l <= cnt[4]; ++l) {
                    pos = 1 + i + j * 2 + k * 3 + l * 4;
                    if(pos > n) continue;
                    if(i > 0) dp[i][j][k][l] = max(dp[i][j][k][l], dp[i-1][j][k][l]);
                    if(j > 0) dp[i][j][k][l] = max(dp[i][j][k][l], dp[i][j-1][k][l]);
                    if(k > 0) dp[i][j][k][l] = max(dp[i][j][k][l], dp[i][j][k-1][l]);
                    if(l > 0) dp[i][j][k][l] = max(dp[i][j][k][l], dp[i][j][k][l-1]);
                    dp[i][j][k][l] += num[pos];
                }
    printf("%d", dp[cnt[1]][cnt[2]][cnt[3]][cnt[4]]);
    return 0;
}