Problem B: qwb与矩阵

直接贴代码，NOIP真题的改遍，但一直TLE，我用的是DP，没想到还可以记忆化搜索。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<cmath>
#include<vector>
#include<map>
#include<set>
using namespace std;
const int inf=(1<<30)-1;
const int maxn=10010;
const int M=10010;
#define REP(i,n) for(int i=(0);i<(n);i++)
#define FOR(i,j,n) for(int i=(j);i<=(n);i++)
typedef long long ll;
int T;
int n,m;
int dp[22][maxn];
int a[22][maxn];
vector<int> G[M];
int main()
{
	for(int i=1;i<M;i++)
    for(int j=i;j<M;j+=i)
	G[j].push_back(i);
    scanf("%d",&T);
    while(T--)
    {
    	memset(dp,0,sizeof(dp));
        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++)
            scanf("%d",&a[i][j]);
        for(int i=0;i<22;i++) dp[i][0]=-inf;
        for(int i=0;i<10002;i++) dp[0][i]=-inf;
        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][j-1],dp[i-1][j]);
                for(int k=0;k<G[j].size();k++)
                dp[i][j]=max(dp[i][j],dp[i][G[j][k]]);
                dp[i][j]+=a[i][j];
            }
        }
        printf("%d\n",dp[n][m]);
    }
    return 0;
}