HDU - 2571 命运(线性dp)

最新推荐文章于 2022-11-15 11:14:21 发布

chagelo

最新推荐文章于 2022-11-15 11:14:21 发布

阅读量185

点赞数

分类专栏： ACM

本文链接：https://blog.csdn.net/UoweMee/article/details/90212523

版权

ACM 专栏收录该内容

113 篇文章 2 订阅

订阅专栏

Link

http://acm.hdu.edu.cn/showproblem.php?pid=2571

Description

穿过幽谷意味着离大魔王 $l e m o n$ 已经无限接近了！

可谁能想到， $y i f e n f e i$ 在斩杀了一些虾兵蟹将后，却再次面临命运大迷宫的考验，这是魔王 $l e m o n$ 设下的又一个机关。要知道，不论何人，若在迷宫中被困1小时以上，则必死无疑！

可怜的 $y i f e n f e i$ 为了去救 $M M$ ，义无返顾地跳进了迷宫。让我们一起帮帮执着的他吧！
yifenfei一开始在左上角，目的当然是到达右下角的大魔王所在地。迷宫的每一个格子都受到幸运女神眷恋或者痛苦魔王的诅咒，所以每个格子都对应一个值，走到那里便自动得到了对应的值。

现在规定 $y i f e n f e i$ 只能向右或者向下走，向下一次只能走一格。但是如果向右走，则每次可以走一格或者走到该行的列数是当前所在列数倍数的格子，即：如果当前格子是 $（ x, y ）$ ，下一步可以是 $（ x + 1, y ）$ ， $(x, y + 1)$ 或者 $(x, y * k)$ 其中 $k > 1$ 。

为了能够最大把握的消灭魔王 $l e m o n ， y i f e n f e i$ 希望能够在这个命运大迷宫中得到最大的幸运值。

Solution

这是一个 $d p$ 题，原型即是数字三角形，不过从三角矩阵变成了矩阵，并且点 $(i, j)$ 可以从 $(i, j - 1)$ 、 $(i - 1, j)$ 及 $(i, j * k)$ 三个点到达。

Code

#include <algorithm>
#include <iostream>
#include <cstring>
#include <string>
#include <cstdio>
#include <vector>
#include <queue>
#include <stack>
#include <cmath>
#include <set>
#include <map>

using namespace std;

#define inf -0x7f7f7f7f
#define maxn 100050
#define N 100100
#define P 2

typedef long long ll;
typedef struct {
    int u, v, next, w;
} Edge;
Edge e[10*maxn];
int cnt, head[maxn];

inline void add(int u, int v) {
    e[cnt].u = u;
    e[cnt].v = v;
//    e[cnt].w = w;
    // e[cnt].f=f;
    e[cnt].next = head[u];
    head[u] = cnt++;
    e[cnt].u = v;
    e[cnt].v = u;
//    e[cnt].w = w;
    //    e[cnt].f=-f;
    e[cnt].next = head[v];
    head[v] = cnt++;
}

inline void write(int x) {
    if (x < 0)
        putchar('-'), x = -x;
    if (x > 9)
        write(x / 10);
    putchar(x % 10 + '0');
}
inline int read() {
    int x = 0, f = 1;
    char c = getchar();
    while (c < '0' || c > '9') {
        if (c == '-')
            f = -1;
        c = getchar();
    }
    while (c >= '0' && c <= '9') {
        x = x * 10 + c - '0';
        c = getchar();
    }
    return x * f;
}
int t,n,m,dp[25][1005],a[25][1006];
int main(){
    t=read();
    while(t--){
        memset(dp,inf,sizeof(dp));
        memset(a,0,sizeof(a));
        n=read(),m=read();
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++)
                a[i][j]=read();
            dp[1][1]=a[1][1];
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++) {
                if(i!=1||j!=1)
                    dp[i][j] = max(dp[i - 1][j], dp[i][j - 1]) + a[i][j];
                for(int k=2;k<=m;k++) {
                    if(j%k==0)
                        dp[i][j] = max(dp[i][j], dp[i][j/k] + a[i][j]);
                }
            }
        cout<<dp[n][m]<<endl;
    }
}