zoj 3662 dp (递推k个数 组成的和为n,最小公倍数为m的所有可能)

#include<cstdio>
#include<cstring>
#define MOD 1e9+7
int dp[2][1010][1002];
int lo[1100][1100],L[110];
int gcd(int a,int b)
{
	if(b==0)
	return a;
	else
	return gcd(b,a%b);
}
int lcm(int a,int b)
{
	return a/gcd(a,b)*b;
} 
int main()
{
	int n,m,k,s;
	for(int i=1;i<=1000;i++)
	for(int j=1;j<=1000;j++)
	lo[i][j]=lcm(i,j);
	while(~scanf("%d%d%d",&n,&m,&k))
	{
		int tot=0;
		for(int i=1;i<=m;i++)
		{
			if(m%i==0)
			L[tot++]=i;
		}
		memset(dp[0],0,sizeof(dp[0]));
		dp[0][0][1]=1;
		s=0;
		for(int i=0;i<k;i++)
		{
			s^=1;
			memset(dp[s],0,sizeof(dp[s]));
			for(int j=i;j<=n;j++)
			{
				for(int p=0;p<tot;p++)
				{
					if(dp[s^1][j][L[p]])
					{
						for(int q=0;q<tot;q++)
						{
							int t=lo[L[p]][L[q]];
							if(t==0||m%t!=0||L[q]+j>n)
							continue;
							dp[s][j+L[q]][t]+=dp[s^1][j][L[p]];
							if(dp[s][j+L[q]][t]>=MOD)
							dp[s][j+L[q]][t]-=MOD;
						}
					}
				}
			}
		}
		printf("%d\n",dp[s][n][m]);
	}
}