180行的高精竟然只有10pts,求助
查看原帖
180行的高精竟然只有10pts,求助
340632
Cry_For_theMoon楼主2020/8/5 13:49

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;
} 
2020/8/5 13:49
加载中...