HDU 2571 命运（DP）

简单dp，状态转移方程题上已给出。。。
但奇怪的一点是在找j的因子时，循环上限到了m?按理说j 应该就是前面的y*k，在y >=1时 k 应该<=j，而k >1，可能是个bug?
如果要保证j/k在迷宫内，应该有j/k>=1，j >=k

#include <iostream>
#include <stack> 
#include <cstring>
#define INF -9999999
using namespace std;

int main()
{
	int C;
	cin >> C;
	while(C--)
	{
		int n,m;
		stack<int> f;
		cin >> n >> m;
		int dp[30][2000],map[30][2000];
		for(int i = 0; i <= n; i++) dp[i][0] = INF;
		for(int j = 0; j <= m; j++) dp[0][j] = INF;
		for(int i = 1; i <= n; i++)
			for(int j = 1; j <= m; j++)
				cin >> map[i][j];
		dp[1][0]=dp[0][1] = 0;
		for(int i = 1; i <= n;	i++)
			for(int j = 1; j <= m; j++)
				{
					dp[i][j] = max(dp[i-1][j],dp[i][j-1]);
					for(int k = 2; k * k <= m;k++)//寻找因子。 
						if(j % k == 0)
						{
							f.push(k);
							f.push(j/k); 
						} 
					
					while(!f.empty())//开始把状态转移方程搞到这里面了，这样j为质数的时候就不会有更新。 
					{
						int q = f.top();f.pop();
						dp[i][j] = max(dp[i][q],dp[i][j]);
					}
					/*for(int k = 2;k <= m; k++)//这里不知道为什么写成小于j就不行 可能是个bug? 
					{
					 	if(j % k == 0) dp[i][j] = max(dp[i][j],dp[i][j/k]);
					 }*/
					dp[i][j] += map[i][j];
				}
		cout << dp[n][m] << endl;
	}
	return 0;
}