[JZYZOJ 1288][洛谷 1005] NOIP2007 矩阵取数 dp 高精度

https://www.luogu.org/problem/show?pid=1005

 

dp好想，高精度练手题，有点不舒服的是前后取数位置的计算，代码量太少才会写题这么慢，noip之前虽然重点放在知识点补全上但是基础还是要打扎实。

代码

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
int n,m;
struct mat{
    int e[50],len;
    mat(){memset(e,0,sizeof(e));len=0;}
};
int a[100][100]={};
mat b2[100]={};
mat f[100][100][2];
mat mul(mat x,int y){
    mat z;z.len=x.len;
    for(int i=1;i<=x.len;i++){
        z.e[i]+=x.e[i]*y; z.e[i+1]+=z.e[i]/10;
        z.e[i]%=10;
        if(z.e[i+1]!=0&&z.len==i){
            z.len++;
        }
    }
    return z;
}
mat plu(mat x,mat y){
    mat z;z.len=max(x.len,y.len);
    for(int i=1;i<=z.len;i++){
        z.e[i]+=x.e[i]+y.e[i];
        z.e[i+1]+=z.e[i]/10;
        z.e[i]%=10;
        if(z.e[i+1]!=0&&i==z.len){
            z.len++;
        }
    }
    return z;
}
mat cham(mat x,mat y){
    if(x.len>y.len)return x;
    if(x.len<y.len)return y;
    for(int i=x.len;i>=1;i--){
        if(x.e[i]>y.e[i])return x;
        if(x.e[i]<y.e[i])return y;
    }
    return x;
}
void put(mat x){
    if(x.len==0)printf("0");
    for(int i=x.len;i>=1;i--){
        printf("%d",x.e[i]);
    }
    printf("
");
}
int main(){
    b2[0].e[1]=1;b2[0].len=1;
    for(int i=1;i<=85;i++){
        b2[i]=mul(b2[i-1],2);
    }
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            scanf("%d",&a[i][j]);
        }
    }mat z,z1,z2,z3,ans;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            for(int w=1;w<=m;w++){
                if(w>j&&w<m-j+1)continue;
                if(w<=j)f[j][w][0]=cham(f[j][w][0],cham(f[j-1][w-1][0],f[j-1][m-j+w+1][1]));
                if(w>=m-j+1) f[j][w][1]=cham(f[j][w][1],cham(f[j-1][w+1][1],f[j-1][j-1-m+w][0]));
                //put(f[j][w]);
                f[j][w][1]=plu(f[j][w][1],mul(b2[j],a[i][w]));
                f[j][w][0]=plu(f[j][w][0],mul(b2[j],a[i][w]));
            }
        }
        z=z3;
        for(int j=1;j<=m;j++){
            z=cham(z,f[m][j][0]);
            z=cham(z,f[m][j][1]);
        }
        for(int j=1;j<=m;j++){
            for(int w=1;w<=m;w++){
                f[j][w][1]=z3;
                f[j][w][0]=z3;
            }
        }
        ans=plu(ans,z);
    }put(ans);
    return 0;
}