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

经过昨天一战 今天开始我要努力学习dp了！

这道题的dp！是区间dp！是我从来没学过的类型 所以我那张题来感觉啥都不会哈哈

我们可以把f_i_j 看成 我们到达区间 [i,j] 时的max

那么状态转移方程就出来了 

注意这里需要用高精 看起来代码很多 其实都是模板

f_i_j 可以是左边过来 那么就要加上a[i-1]*其次幂 我们可以预处理出er【】

然后可以简单记住一个公式就是 m-j+i-1 就是当前次幂

还有就写在注释里面了

f[i][j] = Max(add(f[i - 1][j],mul(er[m - j + i - 1],a[i - 1])),add(f[i][j + 1],mul(er[m - j + i - 1],a[j + 1])));

#include 
using namespace std;
const int N = 85;
const int M = 1e4+10;
const int mod = 1e9+7;
#define int long long
#define endl '\n'
#define Endl '\n'
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return xf[N][N];
vector mul(vector A, int b) {//高精度乘法
    vector  C;
    int t = 0;
    for (int i = 0; i < A.size(); i ++) {
        t += A[i] * b;
        C.push_back(t % 10);
        t /= 10;
    }
    while (t) {
        C.push_back(t % 10);
        t /= 10;
    }
    while (C.size() > 1 && C.back() == 0) C.pop_back();
    return C;
}
vector add(vector(A), vector(B)){//高精度加法
    int t=0;
    vectorC;
    if(A.size()er[N];//预处理二次幂
void two(){
    er[0].push_back(1);
    for(int i=1;i<=80;i++)er[i]=mul(er[i-1],2);
}
vector Max(vectorA,vectorB){//取max
    if(A.size()=0;k--){
            if(A[k]==B[k])continue;
            else if(A[k]>B[k])return A;
            else if(A[k]>n>>m;
    two();
    vectorans;
    while(n--){
        memset(f,0,sizeof f);//每次都要清空啊
        for(int j=1;j<=m;j++)cin>>a[j];//输入
        for(int i=1;i<=m;i++)//状态转移
            for(int j=m;j>=i;j--)//要是>=i了就不用做了 因为再做就是重复的没有必要
                f[i][j] = Max(add(f[i - 1][j],mul(er[m - j + i - 1],a[i - 1])),add(f[i][j + 1],mul(er[m - j + i - 1],a[j + 1])));
        vectormx;//记录最大值，我们知道我们区间dp不能到长度为1的区间所以要自己一个一个遍历
        for(int i=1;i<=m;i++)mx=Max(mx,add(f[i][i],mul(er[m],a[i])));
        ans=add(ans,mx);
    }
    for(int i=ans.size()-1;i>=0;i--)cout<