杭电2571-命运 动态规划

原题链接： http://acm.hdu.edu.cn/showproblem.php?pid=2571

这题浪费了太多时间了，到现在还WA，通过了好几个样例，还是没找到错误。先记录一下学到的其他东西。

错的第一个地方，是找j的因子的时候，sqrt（j）的时候，如果j%p==0，说明找到了两个因子，但是我只写了一个。具体代码在下面。

错的第二个地方，是我对数组初始化的错误。

现在才突然想到d[10]={5}只是对第一个元素初始化为5，其他的初始化为0。所以以后还是少用这个吧。memset(d,0,sizeof(d))也行，但是注意要加头文件# include cstring。

因为memset是按字节赋值的。所以如果想初始化为最小值，可以写memset(d,128,sizeof(d))，memset(d,0x7f,sizeof(d)) . 如果初始化为最大值，可以写memset(d,127,sizeof(d))。

或者还有一种比较好的办法，memset(d,0x3f,sizeof(d))。这样是为了两个数相加时不会爆。

下面代码仍然WA：

#include <iostream>
#include <algorithm>
#include <cstring>
#include "math.h"
using namespace std;
int k[25][1005] = { 0 };
int d[25][1005];   //memset(d,128,sizeof(d))  将所有初始化为最小值
//memset(d,127,sizeof(d))  将所有初始化为最大值

int main()
{
int t;
cin >> t;
for (int i = 0; i < 1005; i++)  d[0][i] = -104;
for (int i = 0; i < 25; i++) d[i][0] = -104;
while (t--){
int n, m;
cin >> n >> m;
for (int i = 1; i <= n; i++){
for (int j = 1; j <= m; j++){
cin >> k[i][j];
}
}
d[1][1] = k[1][1];
for (int i = 1; i <= n; i++){
for (int j = 1; j <= m; j++){
int ans = d[i - 1][j];
ans = max(ans, d[i][j - 1]);
ans = max(ans, d[i][1]);
for (int p = 2; p <= sqrt(j); p++){
if (j%p == 0){
ans = max(ans, d[i][j / p]);
ans = max(ans, d[i][p]);
}
}
if (i != 1 || j != 1) //填充非d[1][1]的位置的元素
d[i][j] = ans + k[i][j];
}
}
cout << d
[m] << endl;
}
return 0;
}