hdu4427 Math Magic

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int MOD=1000000007;
int dp[2][1010][1010];
int num[1000],LCM[1010][1010];
int gcd(int a,int b){
    if(b==0) return a;
    return gcd(b,a%b);
}
int lcm(int a,int b){
    return a/gcd(a,b)*b;
}
int main(){
    int n,m,k;
    int i,j,t,p;
    for(i=1;i<=1000;i++)
        for(j=1;j<=1000;j++)
            LCM[i][j]=lcm(i,j);
    while(scanf("%d%d%d",&n,&m,&k)!=EOF){
        int cnt=0;
        for(i=1;i<=m;i++)
            if(m%i==0) num[cnt++]=i;
        int now=0;
        for(i=0;i<=n;i++)
            for(j=0;j<cnt;j++)
                dp[now][i][num[j]]=0;
        dp[now][0][1]=1;
        for(t=1;t<=k;t++){
            now^=1;
            for(i=0;i<=n;i++)
                for(j=0;j<cnt;j++)
                    dp[now][i][num[j]]=0;
            for(i=t-1;i<=n;i++)
                for(j=0;j<cnt;j++){
                    if(dp[now^1][i][num[j]]==0) continue;
                    for(p=0;p<cnt;p++){
                        int x=i+num[p];
                        int y=LCM[num[j]][num[p]];
                        if(x>n||m%y!=0) continue;
                        dp[now][x][y]+=dp[now^1][i][num[j]];
                        dp[now][x][y]%=MOD;
                    }
                }
        }
        printf("%d\n",dp[now][n][m]);
    }
    return 0;
}