hdu 4427

来源：互联网发布：php短信接口开发实例编辑：程序博客网时间：2024/06/05 13:02

原本以为是数论知识，比赛的时候一直往那方面想，不过后来发现有些方面不能实现，结果比赛的时候就GG了。

优化的时间比没优化的时间快了快一倍，不过没优化也能卡过

/*滚动数组。。。dp[i+1][j+v][lcm[k][v]]+=dp[i][j][k];*/#include<iostream>#include<stdio.h>#include<algorithm>using namespace std;#define Maxn 1010#define Mod 1000000007int lcm[Maxn][Maxn];int dp[2][Maxn][Maxn];int gcd(int a,int b){return b==0?a:gcd(b,a%b);}int LCM(int a,int b){return (a*b)/gcd(a,b);}void Init(){int i,j;for(i=1;i<1001;i++)for(j=1;j<1001;j++)lcm[i][j]=LCM(i,j);}int main(){int n,m,i,j,k;Init();int tt[Maxn];while(~scanf("%d%d%d",&n,&m,&k)){int tot=0;for(i=1;i<=m;i++){if(m%i==0){tt[tot++]=i;}}int now=0;for(int a=0;a<=n;a++)//初始化nowfor(int b=0;b<tot;b++)dp[now][a][tt[b]]=0;dp[now][0][1]=1;for(int t=0;t<k;t++){now^=1;for(int a=0;a<=n;a++)//初始化nowfor(int b=0;b<tot;b++)dp[now][a][tt[b]]=0;for(i=;i<=n;i++)//初始和for(j=0;j<tot;j++){//初始lcmif(dp[now^1][i][tt[j]]==0) continue;//上次状态不存在和为i，lcm为tt[j]这种情况for(int z=0;z<tot;z++){int x=i+tt[z];int tmp=lcm[tt[j]][tt[z]];if(x>n) continue;dp[now][x][tmp]+=dp[now^1][i][tt[j]];dp[now][x][tmp]%=Mod;}}}printf("%d\n",dp[now][n][m]);}return 0;}