HDU 4345 数学+dp

这道题的背景是一个置换，1个置换的循环是它子置换的最小公倍数，所以这道题可以转换成1个纯数论问题，将n划分为若干数，这些数组成的最小公倍数有多少种？

DP弱啊。。。竟然不会，5555

后来看了几篇题解（包括zz1215的）。大概我有个地方没有想明白刚开始

dp[i][x]表示当前和为i,且包含最大质因子不超过素数表第t个素数的情况

状态转移时dp[i][x] = dp[i-x][t-1]+dp[i-x^2][t-1]+....dp[i-x^k][t-1]

用到唯一分解定理，以质因子去构造状态，这个都想到了，但是状态转移的时候，只考虑质因子的幂，不考虑两个质因子的不同幂所组成的和，比如6（2+4），12=（3+9），这种情况，为什吗？

是因为当你有i^p时，再多加一个质因子i^p1(p1<p)，结果不变，可以用i^p1个1来代替。

那么如何处理1的问题，算最后结果的时候把dp[1][m]+...+dp[n][m]，算和就是了。因为dp的时候，没有考虑1，不用担心重复问题。

这道题很推荐，很需要反复思考总结，还用到了置换的背景，赞一个~ 附代码：

#include <iostream>
#include <cstdio>

#define N 1005
#define ss(a) scanf("%d",&a)
#define cl(a) memset(a,0,sizeof(a))     
#define ll __int64

using namespace std;

ll dp[N][N],sum[N];
int pri[N],b[N];

void prime(int &m)
{
    int i,j;
    for (i=2;i<N;i++)
        if (!pri[i])
        {
            b[++m]=i;
            for (j=i*2;j<N;j+=i) pri[j]=1; 
        }
}       

int main()
{
    int n,i,j,m=0;
    prime(m);
    dp[0][0]=1;
    for (j=1;j<=m;j++)
    {
        dp[0][j]=dp[0][j-1];
        for (i=2;i<N;i++)
        {
            int x=b[j];
            dp[i][j]+=dp[i][j-1];
            while (i>=x)
            {
                dp[i][j]+=dp[i-x][j-1];
                x*=b[j];
            } 
        }      
    }
    sum[1]=1;
    for (i=2;i<N;i++) sum[i]=sum[i-1]+dp[i][m];
    while (ss(n)!=EOF)
    {
        printf("%I64d\n",sum[n]);
    }
    return 0;
}