P1005 [NOIP2007 提高组] 矩阵取数游戏

#include<bits/stdc++.h>

#define rep(i , a , b) for(int i = a ; i <= b ; i ++)

using namespace std;

template <typename T> void Out(T x)
{
    if(x < 0) putchar('-') , x = -x;
    if(x > 9) Out(x / 10);
    putchar(x % 10 + '0');
}

const int M = 1e2 + 10;

__int128 two[M] , dp[M][M] , ans;

int a[M][M];

signed main()
{

    two[0] = 1;

    rep(i , 1 , M) two[i] = two[i - 1] * 2;

    int n , m;

    cin >> n >> m;

    for(int i = 1 ; i <= n ; i ++) for(int j = 1 ; j <= m ; j ++) cin >> a[i][j];

    for(int i = 1 ; i <= n ; i ++)
    {
        memset(dp , 0 , sizeof(dp));
        for(int len = m ; len >= 1 ; len --)
        {
            for(int l = 1 ; l + len - 1 <= m ; l ++)
            {
                int r = l + len - 1;
                __int128 c = two[m - (r - l + 1)]; 
                dp[l][r] = max(dp[l - 1][r] + c * a[i][l - 1] , dp[l][r + 1] + c * a[i][r + 1]);
            }
        }
        __int128 ma = 0;
        for(int j = 1 ; j <= m ; j ++) ma = max(ma , dp[j][j] + two[m] * a[i][j]);
        ans += ma;
    }

    Out(ans) , puts("");

    return 0;
}