从上到下的dp错哪了？

#include<bits/stdc++.h>
using namespace std;
#define INF 9e18
int n, k, a[110][110];
long long b[110][110], dp[2][110][5100], ans = -INF;
int main()
{
    cin >> n >> k;
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= i; j++)
        {
            cin >> a[i][j];
            b[i][j] = a[i][j] * 3;
        }
    //输入，用b数组表示3倍
    for(int i = 0; i <= k; i++)
        if(i == 0) dp[1][1][i] = a[1][1];
        else dp[1][1][i] = b[1][1];
    //初始化第一行第一列
    for(int i = 2; i <= n; i++)
        for(int j = 1; j <= i; j++)
            for(int p = k; p >= 1; p--) //p由大到小
    {
        if(j == 1)
        dp[i % 2][j][p] = max(dp[(i - 1) % 2][j][p] + a[i][j], dp[(i - 1) % 2][j][p - 1] + b[i][j]);
        //第一列时只能从上一行对应列转移
        else if(j == i)
        dp[i % 2][j][p] = max(dp[(i - 1) % 2][j - 1][p] + a[i][j], dp[(i - 1) % 2][j - 1][p - 1] + b[i][j]);
        //最后一列只能从上一行前一列转移
        else
        {
        dp[i % 2][j][p] = max(dp[(i - 1) % 2][j][p] + a[i][j], dp[i % 2][j][p]);
        dp[i % 2][j][p] = max(dp[(i - 1) % 2][j][p - 1] + b[i][j], dp[i % 2][j][p]);
        dp[i % 2][j][p] = max(dp[(i - 1) % 2][j - 1][p] + a[i][j], dp[i % 2][j][p]);
        dp[i % 2][j][p] = max(dp[(i - 1) % 2][j - 1][p - 1] + b[i][j], dp[i % 2][j][p]);
        }
        //非特殊列则有四种转移情况
    }
    for(int i = 1; i <= n; i++)
        for(int j = 0; j <= k; j++)
        ans = max(ans, dp[n % 2][i][j]);
    //ans已经初始化为负无穷，同时也枚举了k的值
    cout << ans;
    return 0;
}