#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]);
}
}
zoj 3662 dp (递推k个数 组成的和为n,最小公倍数为m的所有可能)
最新推荐文章于 2020-03-28 11:12:40 发布