HDOJ 4345 Permutation（分组背包 + 数论）

题意：

给你一个数N（1<=N<=1000）,求这么N有多少个不同的旋转长度。旋转长度是指，一个数最少经过多少步可以回到原来的数。

例如 N = 6 时，假如123、45、6分别为三个旋转周期，则旋转顺序为：123456，312546, 231456, 123546, 312456, 231546，123456，旋转长度为6。

当然你也可以12，34，56旋转周期，旋转长度为2。

思路：

即将一个数分成若干份，这若干份不同的最小公倍数的总和，5+4 即5+2^2，最小公倍数为20是一种选择方案

第一组为：2, 2^2, 2^3, ......

第二组为：3, 3^2, 3^3, ......

以此类推，每一组最多只能选择一种解决方案

dp[n] 为填充容量不大于n时，解决发难的总和，若最终选择结果不等于n，则相当于填充了1

其实还可以理解成第0组：1, 11, 111...... 

对于初始化就有，用第0组物品填充，于是dp[0] ~dp[n] = 1

 

#include <cstdio>

#include <cstdlib>

#include <cstring>

#include <cmath>

#include <vector>

using namespace std;



const int MAXN = 1001;

__int64 dp[1005];

vector<int> v[200];



bool isprime(int n)

{

    for (int j = 2; j <= sqrt(n*1.0); ++j)

        if (n % j == 0)

            return false;

    return true;

}



int main()

{

    int c = 0;

    for (int i = 2; i < MAXN; ++i)

    {

        if (isprime(i)) 

        {

            v[c].clear();

            for (int j = i; j < MAXN; j *= i)

               v[c].push_back(j); 

            ++c;

        }

    }

    int n;

    while (scanf("%d", &n) != EOF)

    {

        for (int i = 0; i <= n; ++i)

            dp[i] = 1;

        for (int i = 0; i < c; ++i)

            for (int j = n; j >= 1; --j)

                for (int k = 0; k < v[i].size() && v[i][k] <= j; ++k)

                    dp[j] += dp[j-v[i][k]];

        printf("%I64d\n", dp[n]);

    }

    return 0;

}