zoj3662(01背包+预处理)

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <iostream>
using namespace std;
const int MOD=1000000007;
int gcd(int a,int b){
if(b==0)
return a;
return gcd(b,a%b);
}
int cal(int a,int b){
return a/gcd(a,b)*b;
}
int num[1005],lcm[1005][1005];
int dp[2][1005][1005];
int main(){                                     //dp[i][j][k]代表第i个物品，和为j,lcm为k的种数
int i,j,k,p,N,M,K,st,cnt,tmp;               //因为不能开一个1000*1000*100的数组，因此可以将
for(i=1;i<=1000;i++)                        //i用滚动数组表示
for(j=i;j<=1000;j++)
lcm[i][j]=lcm[j][i]=cal(i,j);               //预处理出所有lcm的值
while(scanf("%d%d%d",&N,&M,&K)!=EOF){
st=cnt=0;
for(i=1;i<=M;i++)                       //以M的因子作为物品
if(M%i==0)
num[cnt++]=i;for(i=0;i<=N;i++)
for(j=0;j<=M;j++)
dp[st][i][j]=0;
dp[st][0][1]=1;
for(k=1;k<=K;k++){
st^=1;
for(i=0;i<=N;i++)
for(j=0;j<=M;j++)
dp[st][i][j]=0;
for(i=k-1;i<=N;i++){                //前k个数和至少是k-1
for(j=0;j<cnt;j++){
if(dp[st^1][i][num[j]]==0)
continue;
for(p=0;p<cnt;p++){
tmp=lcm[num[j]][num[p]];
if(num[p]+i>N||M%tmp!=0)
continue;
dp[st][i+num[p]][tmp]=(dp[st][i+num[p]][tmp]+dp[st^1][i][num[j]])%MOD;
}
}
}
}
printf("%d\n",dp[st]
[M]);
}
return 0;
}