原题链接： 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[n][m] << endl;
    }
    return 0;
}