这是一道NOIp07年的原题,题目本身并不难。
题目看上去很熟悉,第一次看完题目后往贪心的方面去想的,设计了两种贪心策略:1、每次从两端选取最小的数字;2、从后向前倒推,使最后一次取到的数字最大。两种贪心是不同的,而且都是错的,竞赛原题中给的样例数据都过不去。但是如果不会DP的话,感觉上第二种贪心更易导致最优解(只是感觉)。
动态规划:分析出实际上行与行之间是互不影响的,就是说对每行的决策可以单独进行。给定一个区间[i,j],每次可能的决策是选择a[i]或者a[j],这样一个区间被分成了更小的一部分,考虑动态规划。容易写出状态转移方程:d[i][j]=max(d[i+1][j]+w[i],d[i][j-1]+w[j]),其中w[i]可以很容易推出来,不再给出。
此题需要注意的地方:1、需要预处理2^n,因为要被多次用到;2、要用高精度,高精度一般用得比较少,一开始我是把高精度结果按返回值处理,出了错;后来改成传记录结果的地址,在高精度函数中直接修改数值,第8、9两个点超时;分析了原因之后,发现是高精度时传递加数和乘数耗了太多时,就把加数、乘数、结果全部地址传递,vijos上全部数据9ms。导致我对于一直没有太多注意的高精度感慨万千。
以下是我的代码:
#include<stdio.h>
#define size 81
#define max(a,b) (a>b?a:b)
typedef struct
{
int len;
long s[30];
}high;
int n,m,a[size][size]={0};
high d[120][120],Two_[size];
high zero={1,{0}},one={1,{1}},two={1,{2}};
void add(high *c,high *a,high *b)
{// 高精度加上高精度
int l,i;
l=(a->len>b->len?a->len:b->len);
for(i=0;i<=l;i++)
c->s[i]=0;
for(i=0;i<l;i++)
{
c->s[i]+=a->s[i]+b->s[i];
if(c->s[i]>=100000)
{
c->s[i+1]++;
c->s[i]%=100000;
}
}
c->len=l;
if(c->s[c->len]!=0) c->len++;
}
void mul(high *c,high *a,long b)
{// 高精度乘上单精度
int i;
for(i=0;i<=a->len;i++)
c->s[i]=0;
for(i=0;i<a->len;i++)
{
c->s[i]+=a->s[i]*b;
if(c->s[i]>=100000)
{
c->s[i+1]+=c->s[i]/100000;
c->s[i]%=100000;
}
}
c->len=a->len;
if(c->s[c->len]!=0) c->len++;
}
int high_max(high *a,high *b)
{// 比较单精度数
long i;
if(a->len>b->len) return 1;
if(a->len<b->len) return 0;
for(i=a->len-1;i>=0;i--)
if(a->s[i]>b->s[i])
return 1;
else if(a->s[i]<b->s[i])
return 0;
return 1;
}
void print(high *a)
{// 输出高精度
int i;
printf("%ld",a->s[a->len-1]);
for(i=a->len-2;i>=0;i--)
{
if(a->s[i]<10000) printf("0");
if(a->s[i]<1000) printf("0");
if(a->s[i]<100) printf("0");
if(a->s[i]<10) printf("0");
printf("%ld",a->s[i]);
}
printf("\n");
}
void init()
{// 打开文件 读入 预处理2^n 初始化表
int i,j;
scanf("%ld%ld",&n,&m);
for(i=1;i<=n;i++)
for(j=1;j<=m;j++)
scanf("%ld",&a[i][j]);
Two_[0]=one;
Two_[1]=two;
for(i=2;i<=m;i++)
mul(&Two_[i],&Two_[i-1],2);
for(i=0;i<=m;i++)
for(j=0;j<=m;j++)
d[i][j]=zero;
}
high dp(int nn)
{// 对第 nn 行动态规划
int i,j;
high t1,t2,t3,t4;
for(j=0;j<=m-1;j++)
for(i=1;i<=m-j;i++)
{
mul(&t1,&Two_[m-j],a[nn][i]);
add(&t2,&d[i+1][i+j],&t1);
mul(&t3,&Two_[m-j],a[nn][i+j]);
add(&t4,&d[i][i+j-1],&t3);
if(high_max(&t2,&t4))
d[i][i+j]=t2;
else d[i][i+j]=t4;
}
return d[1][m];
}
void end()
{
int i;
high t1,t2,ans=zero;
for(i=1;i<=n;i++)
{
t1=ans;
t2=dp(i);
add(&ans,&t1,&t2);
}
print(&ans);
}
int main()
{
init();
end();
return 0;
}
posted on 2010-01-06 19:43
lee1r 阅读(774)
评论(0) 编辑 收藏 引用 所属分类:
题目分类:动态规划