rt
我的高精挂了啊,求助各位dalao(我觉得估计是哪个高精写错了)
//NOIp2007,提高T3
//我不想写高精555
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn=100;
int n,m,matrix[maxn][maxn],tmp[500];
int f[maxn][maxn][500];
int ans[500],base[500];
int t1[500],t2[500];
void output(int a[]){
for(int i=a[0];i>=1;i--){
printf("%d",a[i]);
}
if(a[0]==0){
printf("0");
}
printf("\n");
}
void mult1(int a[],int b){
for(int i=1;i<=a[0];i++){
a[i]*=b;
}
for(int i=1;i<=a[0]+30;i++){
a[i+1]+=a[i]/10;
a[i]%=10;
}
for(int i=a[0]+30;i>=1;i--){
if(a[i]!=0){
a[0]=i;
return;
}
}
}
void mult2(int a[],int b[]){
memset(tmp,0,sizeof tmp);
for(int i=1;i<=a[0];i++){
for(int j=1;j<=b[0];j++){
tmp[i+j-1]+=a[i]*b[j];
}
}
//最多a[0]+b[0]位
for(int i=1;i<=a[0]+b[0];i++){
tmp[i+1]+=tmp[i]/10;
tmp[i]%=10;
}
for(int i=a[0]+b[0];i>=1;i--){
if(tmp[i]!=0){
memset(a,0,sizeof a);
a[0]=i;
for(int j=1;j<=i;j++){
a[j]=tmp[j];
}
return;
}
}
}
void divide(int a[]){
//a/2
memset(tmp,0,sizeof tmp);
int yushu=0;
for(int i=a[0];i>=1;i--){
yushu*=10;yushu+=a[i];
tmp[i]=yushu/2;
yushu%=2;
}
//从低位往高位
for(int i=a[0];i>=1;i--){
if(tmp[i]!=0){
a[0]=i;
for(int j=1;j<=i;j++){
a[j]=tmp[j];
}
return;
}
}
}
void add(int a[],int b[]){
memset(tmp,0,sizeof tmp);
for(int i=1;i<=a[0]||i<=b[0];i++){
tmp[i]=a[i]+b[i];
}
for(int i=1;i<=max(a[0],b[0])+1;i++){
tmp[i+1]+=tmp[i]/10;
tmp[i]%=10;
}
for(int i=max(a[0],b[0])+1;i>=1;i--){
if(tmp[i]!=0){
for(int j=1;j<=i;j++){
a[j]=tmp[j];
}
a[0]=i;
return;
}
}
}
void maxer(int a[],int b[]){
if(a[0] < b[0]){
for(int j=1;j<=b[0];j++){
a[j]=b[j];
}
a[0]=b[0];
return;
}else if(a[0]==b[0]){
//逐个比较
for(int i=a[0];i>=1;i--){
if(a[i]<b[i]){
for(int j=1;j<=b[0];j++){
a[j]=b[j];
}
a[0]=b[0];
return;
}
}
}
}
int main(){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++){
scanf("%d",&matrix[i][j]);
}
}
for(int k=1;k<=n;k++){
memset(f,0,sizeof f);
memset(base,0,sizeof base);
base[0]=base[1]=1;
for(int i=1;i<=m;i++){
mult1(base,2);
}
for(int j=1;j<=m;j++){
int tmp=matrix[k][j];
while(tmp){
f[j][j][0]++;
f[j][j][f[j][j][0]]=tmp%10;
tmp/=10;
}
mult2(f[j][j],base);
}
for(int len=2;len<=m;len++){
divide(base);
for(int i=1;i<=m;i++){
int j=i+len-1;
if(j>m)break;
//f[i][j]=max(f[i][j],f[i][j-1]+matrix[k][j]*base)
memset(t1,0,sizeof t1);
memset(t2,0,sizeof t2);
for(int x=1;x<=f[i][j-1][0];x++){
t1[x]=f[i][j-1][x];
}
t1[0]=f[i][j-1][0];
for(int x=1;x<=base[0];x++){
t2[x]=base[x];
}
t2[0]=base[0];
mult1(t2,matrix[k][j]);
//t1+t2
add(t1,t2);
//max(f[i][j],t1)
maxer(f[i][j],t1);
//f[i][j]=max(f[i][j],f[i+1][j]+matrix[k][i]*base)
memset(t1,0,sizeof t1);memset(t2,0,sizeof t2);
for(int x=1;x<=f[i+1][j][0];x++){
t1[x]=f[i+1][j][x];
}
t1[0]=f[i+1][j][0];
for(int x=1;x<=base[0];x++){
t2[x]=base[x];
}
t2[0]=base[0];
mult1(t2,matrix[k][i]);
add(t1,t2);
maxer(f[i][j],t1);
}
}
add(ans,f[1][m]);
}
output(ans);
return 0;
}