[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[],len;
     mat(){memset(e,,sizeof(e));len=;}
 };
 int a[][]={};
 mat b2[]={};
 mat f[][][];
 mat mul(mat x,int y){
     mat z;z.len=x.len;
     for(int i=;i<=x.len;i++){
         z.e[i]+=x.e[i]*y; z.e[i+]+=z.e[i]/;
         z.e[i]%=;
         if(z.e[i+]!=&&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=;i<=z.len;i++){
         z.e[i]+=x.e[i]+y.e[i];
         z.e[i+]+=z.e[i]/;
         z.e[i]%=;
         if(z.e[i+]!=&&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>=;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==)printf("");
     for(int i=x.len;i>=;i--){
         printf("%d",x.e[i]);
     }
     printf("\n");
 }
 int main(){
     b2[].e[]=;b2[].len=;
     for(int i=;i<=;i++){
         b2[i]=mul(b2[i-],);
     }
     scanf("%d%d",&n,&m);
     for(int i=;i<=n;i++){
         for(int j=;j<=m;j++){
             scanf("%d",&a[i][j]);
         }
     }mat z,z1,z2,z3,ans;
     for(int i=;i<=n;i++){
         for(int j=;j<=m;j++){
             for(int w=;w<=m;w++){
                 if(w>j&&w<m-j+)continue;
                 if(w<=j)f[j][w][]=cham(f[j][w][],cham(f[j-][w-][],f[j-][m-j+w+][]));
                 if(w>=m-j+) f[j][w][]=cham(f[j][w][],cham(f[j-][w+][],f[j-][j--m+w][]));
                 //put(f[j][w]);
                 f[j][w][]=plu(f[j][w][],mul(b2[j],a[i][w]));
                 f[j][w][]=plu(f[j][w][],mul(b2[j],a[i][w]));
             }
         }
         z=z3;
         for(int j=;j<=m;j++){
             z=cham(z,f[m][j][]);
             z=cham(z,f[m][j][]);
         }
         for(int j=;j<=m;j++){
             for(int w=;w<=m;w++){
                 f[j][w][]=z3;
                 f[j][w][]=z3;
             }
         }
         ans=plu(ans,z);
     }put(ans);
     return ;
 }