高精度-四位压缩（结构体重载运算符）+DP（洛谷P1005-方格取数）

洛谷P1005方格取数
这是一道典型的dp题
且根据题意得行与行之间没有任何的联系，只需找到每一行的最优解相加就得出最终答案。
设dp数组为f[maxn][maxn],并将每一行的数据储存在数组a[maxn]中，
分析题意dp规则为:
f[q][w]=max(f[q-1][w]+pow(2,m-(w-q+1))*a[q-1],f[q][w+1]+pow(2,m-(w-q+1))a[w+1]);
但是提交后会发现只能通过60%的数据原因就下面：

其中1<=n,m<=80，也就意味着要计算2的80次方，而c++中longlong的范围为范围：-922 3372 0368 5477 5808 ～ 922 3372 0368 5477 5807 (92210^16) ，所以就必须要调用高精算法（其实用Java中的大数类也能实现）
高精度模板之前已经发了：
高精度模板（结构体重载运算符（+ - *））
最后一步：通过pow2()函数预先算出2的幂，并储存在数组p[maxn]中，以便调用

#include <iostream>
#include <algorithm>
#include <cmath>
#include <vector>
#include <string>
#include <cstring>
#include <cstdio>
#define maxn 86
#define mod 10000
using namespace std;
int n,m,a[maxn];
struct BI{
int l[505],len;
BI(){
    memset(l,0,sizeof(l));
    len=0;
}
void print(){
    printf("%d",l[len]);
for(int q=len-1;q>0;q--){
if(l[q]==0){
printf("0000");
continue;
}
for(int w=10;w*l[q]<mod;w*=10){
printf("0");
}
printf("%d",l[q]);
}}}f[maxn][maxn],p[maxn],sum;
BI operator +(const BI &a,const BI &b){
BI c;c.len=max(a.len,b.len);int x=0;
for(int q=1;q<=c.len;q++){
c.l[q]=a.l[q]+b.l[q]+x;
x=c.l[q]/mod;
c.l[q]%=mod;
}
if(x>0){c.l[++c.len]=x;}
return c;
}
BI operator *(const BI &a,int b){
BI c;c.len=a.len;int x=0;
for(int q=1;q<=c.len;q++){
c.l[q]=a.l[q]*b+x;
x=c.l[q]/mod;
c.l[q]%=mod;
}
while(x>0){
c.l[++c.len]=x%mod;
x/=mod;
}
return c;
}
BI  max(const BI &a,const BI &b){
BI c=a;
if(a.len!=b.len)
c=a.len>b.len?a:b;
else{
for(int q=a.len;q>0;q--){
if(a.l[q]!=b.l[q]){
c=a.l[q]>b.l[q]?a:b;
break;
}}}
return c;
}
void pow2(){
p[0].l[1]=1;p[0].len=1;
for(int q=1;q<=m+2;q++){
p[q]=p[q-1]*2;
}
}
BI ys(int a[]){
BI maxx;
maxx.len=0;
for(int q=1;q<=m;q++){
for(int w=m;w>=q;w--){
f[q][w]=max(f[q-1][w]+p[m-(w-q+1)]*a[q-1],f[q][w+1]+p[m-(w-q+1)]*a[w+1]);
}
f[q][q]=f[q][q]+p[m]*a[q];
maxx=max(maxx,f[q][q]);
}
return maxx;
}
int main()
{
    scanf("%d%d",&n,&m);
    pow2();
    int a[maxn];
    for(int q=0;q<n;q++){
    memset(f,0,sizeof(f));
    memset(a,0,sizeof(a));
    for(int w=1;w<=m;w++){
    scanf("%d",&a[w]);
    }
    sum=sum+ys(a);
    }
    sum.print();
    return 0;
}