四维dp 或者 剪枝 + dfs Codeforces Beta Round #6 (Div. 2 Only) D

http://codeforces.com/contest/6/problem/D

题目大意：有一队人，排成一列，每个人都有生命值，你每次可以攻击2~n位置的一个的人，假设每次攻击的位置为pos，那么pos位受到a点伤害，pos-1和pos+1受到b点伤害。问让所有人生命值都小于0所需要的最少操作数。

思路：最近在加强dp类型的题目，但是因为做过一点IDA*，感觉这题貌似可以用IDA*差不多的思路来做哈？（关键是我不会DP，TAT）然后我们找寻一下IDA*的可行性剪枝就好了：目前数列的maxval<=(lim - deep) * a即可，否则TLE了

这是IDA*的代码

//看看会不会爆int!数组会不会少了一维！
//取物问题一定要小心先手胜利的条件
#include <bits/stdc++.h>
using namespace std;
#define LL long long
#define ALL(a) a.begin(), a.end()
#define pb push_back
#define mk make_pair
#define fi first
#define se second
#define haha printf("haha\n")
const int maxn = ;
const int inf = 0x3f3f3f3f;
int road[maxn * maxn];
int n, a, b;
int h[maxn];

bool dfs(int pos, int deep, int lim){

    int maxcnt = -inf;
    for (int i = ; i <= n; i++) maxcnt = max(maxcnt, h[i]);
    if (deep > lim) return false;
    if (maxcnt > (lim - deep) * a) return false;//加上这个剪枝就变得超快
    if (deep == lim){
        if (maxcnt >= )return false;
        else return true;
    }
    if (pos == n) return false;

    int cnt = max(h[pos - ] / b, max(h[pos] / a, h[pos + ] / b)) + ;
    for (int cnt1 = ; cnt1 <= cnt; cnt1++){
        if (h[pos - ] >= cnt1 * b) continue;
        h[pos - ] -= b * cnt1; h[pos] -= a * cnt1; h[pos + ] -= b * cnt1;
        for (int j = deep + ; j <= deep + cnt1; j++)
            road[j] = pos;
        if (dfs(pos + , deep + cnt1, lim)) {
            for (int j = deep + ; j <= deep + cnt1; j++) road[j] = pos;
            return true;
        }
        h[pos - ] += b * cnt1; h[pos] += a * cnt1; h[pos + ] += b * cnt1;
    }
    return false;
}

int main(){
    memset(road, , sizeof(road));
    cin >> n >> a >> b;
    int beg = ;
    for (int i = ; i <= n; i++) scanf("%d", h + i);
    for (int i = ; i <=  * ; i++)
        if (dfs(, , i)){
            printf("%d\n", i);
            for (int j = ; j <= i; j++)
                printf("%d%c", road[j], j == i ? '\n' : ' ');
            break;
        }
    return ;
}

明天附上DP代码。晚安~

dp思路：

定义dp[i][j][k][z]表示目前攻击第i个人，j表示第i-1个人的生命，k表示第i个人的生命，z表示第i+1个人的生命。然后转移需要注意一下以下方面

①简单的转移：dp[i][j - b][k - a][z - b] = dp[i][j][k][z];

②如果j-b以后的生命值为0了，那么dp[i+1][k - a][z - b][h[i + 2]]=dp[i][j][k][z] + 1;

③如果j本来就为0，那么dp[i+1][k][z][h[i + 2]] = dp[i][j][k][z];

//看看会不会爆int!数组会不会少了一维！
//取物问题一定要小心先手胜利的条件
#include <bits/stdc++.h>
using namespace std;
#define LL long long
#define ALL(a) a.begin(), a.end()
#define pb push_back
#define mk make_pair
#define fi first
#define se second
#define haha printf("haha\n")
const int maxn = ;
const int inf = 0x3f3f3f3f;
struct Point{
    int i, j, k, z;
    Point(int i = , int j = , int k = , int z = ): i(i), j(j), k(k), z(z){}
}par[maxn][maxn][maxn][maxn];
int dp[maxn][maxn][maxn][maxn];
int n, a, b;
int h[maxn];

int main(){
    cin >> n >> a >> b;
    for (int i = ; i <= n; i++){
        scanf("%d", h + i); h[i]++;
    }
    memset(par, -, sizeof(par));
    memset(dp, 0x3f, sizeof(dp));
    dp[][h[]][h[]][h[]] = ;
    for (int i = ; i <= n - ; i++){
        for (int j = h[i - ]; j >= ; j--){
            for (int k = h[i]; k >= ; k--){
                for (int z = h[i + ]; z >= ; z--){
                    int nj = max(, j - b), nk = max(, k - a), nz = max(, z - b);
                    if (dp[i][nj][nk][nz] > dp[i][j][k][z] + ){
                        dp[i][nj][nk][nz] = dp[i][j][k][z] + ;
                        par[i][nj][nk][nz] = Point(i, j, k, z);
                    }
                    if (nj ==  && dp[i + ][nk][nz][h[i + ]] > dp[i][j][k][z] + ){
                        dp[i + ][nk][nz][h[i + ]] = dp[i][j][k][z] + ;
                        par[i + ][nk][nz][h[i + ]] = Point(i, j, k, z);
                    }
                    if (j ==  && dp[i + ][k][z][h[i + ]] > dp[i][j][k][z]){
                        dp[i + ][k][z][h[i + ]] = dp[i][j][k][z];
                        par[i + ][k][z][h[i + ]] = par[i][j][k][z];
                    }
                }
            }
        }
    }
    int ans = dp[n][][][];
    printf("%d\n", ans);
    Point tmp = Point(n, , , );
    while (ans--){
        tmp = par[tmp.i][tmp.j][tmp.k][tmp.z];
        printf("%d ", tmp.i);
    }
    return ;
}