hdu4427Math Magic

4427

dp[i][j][k] i为K位的最小公倍数 j为k位的和 k以滚动数组的形式

这题最棒的是 有一个强有力的剪枝 组成公倍数m的肯定都是M的质因子 这样1000里面最多就30多个 复杂度可过了

#include <iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<stdlib.h>
#include<queue>
using namespace std;
#define mod 1000000007
int dp[1010][1010][2];
int q[110][1010],f[1010],lc[1010][1010];
int p[1010];
int gcd(int a,int b)
{
    return b==0?a:gcd(b,a%b);
}
int main()
{
    int n,m,k,i,j,g;
    for(i = 1 ; i <= 1000 ; i++)
        for(j = 1; j <= 1000 ; j++)
        lc[i][j] = i*j/gcd(i,j);
    while(scanf("%d%d%d",&n,&m,&k)!=EOF)
    {
        int o;
        int to = 0;
        for(i =1; i <= m ; i++)
        if(m%i==0)
        {
            to++;
            p[to] = i;
        }
        memset(dp,0,sizeof(dp));
        for(i = 1 ; i <= min(n,m) ; i++)
        {
            dp[i][i][1] = 1;
        }
        for(i = 2; i <= k; i++)
        {
            for(o = 1; o <= to ; o++)
                for(g = 1; g <= n ; g++)
                dp[p[o]][g][i%2] = 0;
                for(j = 1 ; j <= to ; j++)
                {
                    for(g = i-1; g <= n ; g++)
                    {
                        if(dp[p[j]][g][(i-1)%2]==0)
                        continue;
                        for(o = 1 ; o <= to ; o++)
                        {
                            int x = lc[p[j]][p[o]];
                            if(g+p[o]>n)
                            break;
                            if(x>m)
                            continue;
                            dp[x][g+p[o]][i%2] = (dp[x][g+p[o]][i%2]+dp[p[j]][g][(i-1)%2])%mod;
                        }
                    }
                }
        }
        printf("%d
",dp[m][n][k%2]);
    }
    return 0;
}