本质上就是一个在一行中头尾取数情况下求出分数最大值。设d【i,j】表示该行从i到j的合并分数
则d[i,j]=2*max{d[i+1][j]+map[a][i],d[i][j-1]+map[a][j]};
由于时间比较晚了,可能不够时间写高精度,先存普通版本算法
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#define MAXN 100
int map[MAXN][MAXN];
int d[MAXN][MAXN];
int main(void)
{
int n,m;
int i,j,k;
int a,len;
int x,y;
int r=0;
scanf("%d%d",&n,&m);
for(i=1;i<=n;i++)
{
for(j=1;j<=m;j++)
{
scanf("%d",&map[i][j]);
}
}
for(a=1;a<=n;a++)
{
for(i=1;i<=m;i++)
{
d[i][i]=map[a][i]*2;
}
for(i=m-1;i>=1;i--)
{
for(j=i+1;j<=m;j++)
{
x=2*(d[i+1][j]+map[a][i]);
y=2*(d[i][j-1]+map[a][j]);
d[i][j]=x>y?x:y;
}
}
r+=d[1][m];
}
printf("%d\n",r);
return 0;
}
50分代码
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#define MAXN 100
#define MAX 100
typedef struct
{
int n[MAX];
int len;
}bigint;
bigint change(char ch[])
{
int i,j,k;
bigint a;
memset(a.n,0,sizeof(a.n));
a.len=strlen(ch);
for(i=1;i<=a.len;i++)
{
a.n[i]=ch[a.len-i]-'0';
}
return a;
}
bigint addtion(bigint a,bigint b)
{
bigint c;
memset(c.n,0,sizeof(c.n));
int i,j,k;
for(i=1;i<=a.len||i<=b.len;i++)
{
c.n[i]+=(a.n[i]+b.n[i])%10;
c.n[i+1]=(a.n[i]+b.n[i])/10;
}
c.len=i-1;
if(c.n[c.len+1]!=0)
c.len++;
return c;
}
bigint plus(bigint a,int b)
{
int i,j,k;
bigint c;
memset(c.n,0,sizeof(c.n));
for(i=1;i<=a.len;i++)
{
c.n[i]+=(a.n[i]*b)%10;
if(c.n[i]>=10)
{
c.n[i]-=10;
c.n[i+1]++;
}
c.n[i+1]+=(a.n[i]*b)/10;
}
c.len=a.len;
if(c.n[c.len+1]!=0)
c.len++;
return c;
}
void print(bigint a)
{
int i;
for(i=a.len;i>=1;i--)
{
printf("%d",a.n[i]);
}
printf("\n");
}
int compare(bigint a,bigint b)
{
int i,j;
if(a.len>b.len)
return 1;
else if(a.len<b.len)
return -1;
else
{
for(i=a.len;i>=1;i--)
{
if(a.n[i]>b.n[i])
return 1;
else if(a.n[i]<b.n[i])
return -1;
}
}
return 0;
}
bigint map[MAXN][MAXN];
bigint d[MAXN][MAXN];
int main(void)
{
char ch[MAX];
int n,m;
int i,j,k;
int a,len;
bigint x,y;
bigint r;
memset(r.n,0,sizeof(r.n));
r.len=0;
scanf("%d%d",&n,&m);
for(i=1;i<=n;i++)
{
for(j=1;j<=m;j++)
{
scanf("%s",ch);
map[i][j]=change(ch);
}
}
for(a=1;a<=n;a++)
{
for(i=1;i<=m;i++)
{
d[i][i]=plus(map[a][i],2);
}
for(i=m-1;i>=1;i--)
{
for(j=i+1;j<=m;j++)
{//大整数之间的加法、大整数与integer的乘法、大整数的比较
x=addtion(d[i+1][j],map[a][i]);
x=plus(x,2);
y=addtion(d[i][j-1],map[a][j]);
y=plus(y,2);
if(compare(x,y)==1)
{
d[i][j]=x;
}
else
d[i][j]=y;
}
}
r=addtion(r,d[1][m]);
}
print(r);
return 0;
}
C++100代码
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxdigit=30+10;
struct Bign
{
int dig,num[maxdigit];
void init()
{
dig=1;
memset(num,0,sizeof(num));
}
}f[80+10][80+10],_pow[80+10],ans;
Bign Two=(Bign){1,{0,2}};
Bign operator + (Bign A,Bign B)
{
Bign C;
memset(C.num,0,sizeof(C.num));
for(int i=1;i<=(C.dig=max(A.dig,B.dig));i++)
{
C.num[i]+=A.num[i]+B.num[i];
C.num[i+1]=C.num[i]/10;
C.num[i]%=10;
}
if(C.num[C.dig+1]>0)
C.dig++;
return C;
}
Bign operator * (Bign A,Bign B)
{
Bign C;
memset(C.num,0,sizeof(C.num));
for(int i=1;i<=A.dig;i++)
for(int j=1;j<=B.dig;j++)
C.num[i+j-1]+=A.num[i]*B.num[j];
for(int i=1;i<=(C.dig=A.dig+B.dig-1);i++)
{
C.num[i+1]+=C.num[i]/10;
C.num[i]%=10;
}
if(C.num[C.dig+1]>0)
C.dig++;
return C;
}
Bign operator * (Bign A,int B)
{
Bign C;
memset(C.num,0,sizeof(C.num));
for(int i=1;i<=A.dig;i++)
C.num[i]+=A.num[i]*B;
for(C.dig=1;;C.dig++)
{
C.num[C.dig+1]+=C.num[C.dig]/10;
C.num[C.dig]%=10;
if(C.num[C.dig+1]==0)
break;
}
return C;
}
Bign _max(Bign A,Bign B)
{
if(A.dig>B.dig)
return A;
if(A.dig<B.dig)
return B;
for(int i=A.dig;i>=1;i--)
if(A.num[i]>B.num[i])
return A;
else
if(A.num[i]<B.num[i])
return B;
return A;
}
int n,m,a[80+10];
int main()
{
scanf("%d%d",&n,&m);
_pow[0]=(Bign){1,{0,1}};
for(int i=1;i<=m;i++)
_pow[i]=_pow[i-1]*Two;
for(int r=1;r<=n;r++)
{
for(int c=1;c<=m;c++)
scanf("%d",&a[c]);
for(int i=1;i<=m;i++)
for(int j=1;j<=m;j++)
f[i][j].init();
for(int i=1;i<=m;i++)
f[i][i]=_pow[m]*a[i];
for(int i=m-1;i>=1;i--)
for(int j=i+1;j<=m;j++)
f[i][j]=_max(f[i+1][j]+_pow[m-j+i]*a[i],f[i][j-1]+_pow[m-j+i]*a[j]);
ans=ans+f[1][m];
}
for(int i=ans.dig;i>=1;i--)
printf("%d",ans.num[i]);
printf("\n");
return 0;
}