1025: [SCOI2009]游戏（完全背包求种数，易错）

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 。

题解：

问题相当于一个各项之和为N的正整数序列，求所有这样的序列中的各数的最小公倍数有多少种。

设最小公倍数为m，则m=p1^a1*p2^a2...。
而对于一个m，存在一个序列的最小公倍数为m的充要条件是p1^a1+p2^a2+....<=n。
预处理出n以内的质数后，用类似背包Dp的递推统计即可。

原（wrong）:

for(int i=1;i<=tot;i++){
int k=pri[i];
for(;k<=n;k*=pri[i])
for(int j=n;j>=k;j--){
f[j]+=f[j-k];
}
}

错因在重复，如2^2+2^1=6，本来6是不可能出现的。

(right):

for (int i=0;i<tot;i++)
    {
        memcpy(tmp,f,sizeof(f));
        for (int j=ss[i];j<=n;j*=ss[i])
        {
            for (int k=j;k<=n && k-j>=0;k++)
            {
                f[k]+=tmp[k-j];
            }
        }
    }

code：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
using namespace std;
const int MaxN=1051;


int n,tot=0;
unsigned long long f[MaxN],tmp[MaxN],ans=0;
int ss[MaxN];
bool vis[MaxN];


void treat()
{
    for (int i=2;i<=n;i++)
    {
        if (vis[i]==0)
        {
            ss[tot++]=i;
            vis[i]=1;
        }
        for (int j=0;j<tot && ss[j]*i<=n;j++)
        {
            vis[ss[j]*i]=1;
            if (i%ss[j]==0)
            {
                break;
            }
        }
    }
}


int main()
{
    cin>>n;
    treat();
    f[0]=1;
    for (int i=0;i<tot;i++)
    {
        memcpy(tmp,f,sizeof(f));
        for (int j=ss[i];j<=n;j*=ss[i])
        {
            for (int k=j;k<=n && k-j>=0;k++)
            {
                f[k]+=tmp[k-j];
            }
        }
    }
    for (int i=1;i<=n;i++)
    {
        ans+=f[i];
    }
    cout<<ans+1<<endl;
    return 0;
}