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想知道，对于所有可能的对应关系，有多少种可能的排数。
题意是让我们求回归到原序列的置换次数，也就是对于这个置换的每一个循环大小的最小公倍数
如（123）（45）（6）
对于第一个循环回归需要3x次(x代表任意自然数)，第二个循环需要2x次，第三个循环需要x次
显然对于这样的置换答案是6次
有多少种可能的排数，也就是这个最小公倍数有多少种可能..?
看起来比较奇怪的问题...
我们用一个序列来表示最小公倍数
2 3 5 7 11...997
x1 x2 x3 x4 x5...x168
最小公倍数就是2^x1*3^x2*5^x3...997^x168（说白了就是质因数分解>w<)
对于每一个循环，它的大小质因数分解之后对于每个质因子的系数都<=上面的序列
然后原问题就是满足题意的x1,x2..x168有多少种
我们来考虑怎样才会不满足题意...
首先x1,x2..x168再小也无所谓，因为就算全为0也是满足题意的，剩下的都可以用1来补全
显然对于这些x最小的n就是2^x1+3^x2...+997^x168
因为一旦2出现了x1次，说明一定有一个循环的个数中含有2^x1这个因子，我们令其中一个等于这个
为了使n最小化，其它都不含2这个因子
另外加法比乘法代价小这个显然...
也就是说对于x序列的限制也就是2^x1+3^x2...997^x168<=n
非常眼熟...非常简单的01背包...
1 /************************************************************** 2 Problem: 1025 3 User: mjy0724 4 Language: Pascal 5 Result: Accepted 6 Time:112 ms 7 Memory:8236 kb 8 ****************************************************************/ 9 10 program bzoj1025; 11 const maxn=1010; 12 var n,ans:int64; 13 i,j,k:longint; 14 p:array[-1..maxn]of int64; 15 f:array[-1..maxn,-1..maxn]of int64; 16 vis:array[-1..maxn]of boolean; 17 18 procedure Euler; 19 var i,j:longint; 20 begin 21 fillchar(vis,sizeof(vis),true); 22 p[0]:=0; 23 for i:=2 to n do 24 begin 25 if vis[i] then 26 begin 27 inc(p[0]); 28 p[p[0]]:=i; 29 end; 30 for j:=1 to p[0] do 31 begin 32 if i*p[j]>n then break; 33 vis[i*p[j]]:=false; 34 if i mod p[j]=0 then break; 35 end; 36 end; 37 end; 38 39 begin 40 readln(n); 41 Euler; 42 fillchar(f,sizeof(f),0); 43 f[0,0]:=1; 44 for i:=1 to p[0] do 45 for j:=0 to n do 46 begin 47 f[i][j]:=f[i-1][j]; 48 k:=p[i]; 49 while k<=j do 50 begin 51 inc(f[i,j],f[i-1,j-k]); 52 k:=k*p[i]; 53 end; 54 end; 55 ans:=0; 56 for i:=0 to n do inc(ans,f[p[0],i]); 57 writeln(ans); 58 end.
