zoj3662(dp)

dp还是比较好想的，但是时间还是比较坑的。

要预处理还加些优化才行 。

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <algorithm>
#include <math.h>
#include <map>
#include <queue>
#include <sstream>
#include <iostream>
using namespace std;
#define INF 0x3fffffff
#define MOD 1000000007

int n,m,k;

int dp[][][];
int save[][];
int g[];
int cnt=;
int mark[];
//bool mark1[1100][1100];

int gcd(int x,int y)
{
    int tmp;
    while(x)
    {
        tmp=x;
        x=y%x;
        y=tmp;
    }
    return y;
}

int lcm(int x,int y)
{
    return x*y/gcd(x,y);
}

int main()
{
    //freopen("C:\\Users\\Administrator\\Desktop\\in.txt","r",stdin);
    //freopen("C:\\Users\\Administrator\\Desktop\\in.txt","w",stdout);
    for(int i=;i<=;i++)
        for(int j=;j<=;j++)
        {
            save[i][j]=lcm(i,j);
        }
        //提前预处理

    while(scanf("%d%d%d",&n,&m,&k)!=EOF)
    {
        //memset(mark1,0,sizeof(mark1));
        //memset(mark,0,sizeof(mark));
        memset(dp,,sizeof(dp));
        cnt=;
        for(int i=;i<=m;i++)
            if(m%i==)
            {
                mark[i]=cnt;
                g[cnt++]=i;

            }
        dp[][][]=;
        //mark1[0][0]=1;
        int tmp,tmp1,tmp2;
        for(int i=;i<k;i++) //对于k个数
        {
            for(int j=i;j<=n;j++)
            {
                //if(mark1[i][j]==0) continue;
                for(int p=;p<cnt;p++)
                {
                    tmp2=dp[i][j][p];
                    if(tmp2==) continue;
                    for(int s=;s<cnt;s++) //选择g[s]
                    {
                        if(g[s]+j > n-(k-i)+) break;
                        tmp1=g[s];
                        tmp=mark[ save[ g[p] ][ tmp1 ] ];

                        dp[i+][ j+tmp1 ][ tmp ]=(dp[i+][ j+tmp1 ][ tmp ]+tmp2)%MOD;

                        //if(dp[i+1][ j+tmp1 ][tmp]>MOD) dp[i+1][ j+tmp1 ][tmp]-=MOD;
                    }
                }
            }
        }

        printf("%d\n",dp[k][n][cnt-]);
    }
    return ;
}