BZOJ1025 SCOI2009游戏

1025: [SCOI2009]游戏

Time Limit: 1 Sec  Memory Limit:162 MB
Submit: 1594  Solved: 1003
[Submit][Status][Discuss]

Description

windy学会了一种游戏。对于1到N这N个数字，都有唯一且不同的1到N的数字与之对应。最开始windy把数字按顺序1，2，3，……，N写一排在纸上。然后再在这一排下面写上它们对应的数字。然后又在新的一排下面写上它们对应的数字。如此反复，直到序列再次变为1，2，3，……，N。 如： 1 2 3 4 5 6 对应的关系为 1->2 2->3 3->1 4->5 5->4 6->6 windy的操作如下 1 2 3 4 5 6 2 3 1 5 4 6 3 1 2 4 5 6 1 2 3 5 4 6 2 3 1 4 5 6 3 1 2 5 4 6 1 2 3 4 5 6 这时，我们就有若干排1到N的排列，上例中有7排。现在windy想知道，对于所有可能的对应关系，有多少种可能的排数。

Input

包含一个整数，N。

Output

包含一个整数，可能的排数。

Sample Input

【输入样例一】
3
【输入样例二】
10

Sample Output

【输出样例一】
3
【输出样例二】
16

HINT

【数据规模和约定】

100%的数据，满足 1 <= N <= 1000 。

先YY一下，本题题意为：将N拆分成若干个数并计算他们的LCM，问不同的LCM个数。

设N = A1 + A2 + ... + A[t] 再将每个数A分解成Pi ^ Bi的乘积，则LCM = Pi ^ max(Bi) 的乘积（妈蛋优快云有公式编辑器么）那么对于每个Pi枚举Bi即可。

令F[i][j]表示在前i个质数内选择若干个且和为j的方案数，则：F[i][j] = F[i - 1][j] + F[i - 1][j - P[i] ^ k].

代码如下：

/**************************************************************
    Problem: 1025
    User: duyixian
    Language: C++
    Result: Accepted
    Time:20 ms
    Memory:20844 kb
****************************************************************/
 
/* 
* @Author: 逸闲
* @Date:   2015-10-01 11:49:05
* @Last Modified by:   逸闲
* @Last Modified time: 2015-10-01 16:03:23
*/
 
#include "cstdio"
#include "cstdlib"
#include "iostream"
#include "algorithm"
#include "cstring"
#include "queue"
 
using namespace std;
 
#define INF 0x3F3F3F3F
#define MAX_SIZE 5000
#define Eps 
#define Mod 
 
inline int Get_Int()
{
    int Num = 0;
    char ch;
    do
        ch = getchar();
    while(ch < '0' || ch > '9');
    do
    {
        Num = Num * 10 + ch - '0';
        ch = getchar();
    }
    while(ch >= '0' && ch <= '9');
    return Num;
}
 
int N, Total;
int Flag[MAX_SIZE], P[MAX_SIZE];
 
long long Ans;
long long DP[500][MAX_SIZE];
 
 
int main()
{
    cin >> N;
    for(int i = 2; i <= N; ++i)
        if(!Flag[i])
        {
            P[++Total] = i;
            for(int j = i + i; j <= N; j += i)
                Flag[j] = true;
        }
    DP[0][0] = 1;
    for(int i = 1; i <= Total; ++i)
        for(int j = 0; j <= N; ++j)
        {
            DP[i][j] = DP[i - 1][j];
            for(int k = P[i]; k <= j; k *= P[i])
                DP[i][j] += DP[i - 1][j - k];
        }
    for(int i = 0; i <= N; ++i)
        Ans += DP[Total][i];
    cout << Ans << endl;
    return 0;
}