Time Limit: 1 Sec Memory Limit: 162 MB
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
3
【输入样例二】
10
Sample Output
【输出样例一】
3
【输出样例二】
16
3
【输出样例二】
16
HINT
【数据规模和约定】
100%的数据,满足 1 <= N <= 1000 。
SOLUTION
样例中每一列都是由循环节组成的,比如第一列1-2-3,第二列2-3-1,第三列3-1-2,第四列4-5……那么我们可以知道,如果原数列当中的一个数为i,i所在列的循环节长度为L,那么在经过k*L次变换之后,原来i所在的那一列的数字又将变成i。若要使数列所有位都变回原装态,就要使排列的行数Lines是所有循环节长度的整数倍。
如果设循环节长度分别为L1,L2,L2,......,Ln,那么Lines=LCM(L1,L2,L3,......,Ln)。
至此,问题被转化成了:给你一个N,问你任意一坨循环节长度的LCM是N,有多少种情况。
预处理把整数转化为素数幂次的积,然后记忆化搜索。
var prime:array[1..1000]of int64; f:array[1..1000,0..1000]of int64; n,sum:int64; function p(x:int64):boolean; var i:longint; begin for i:=1 to sum do if x mod prime[i]=0 then exit(false); exit(true); end; procedure main; var i:longint; begin for i:=2 to n do if p(i) then begin inc(sum); prime[sum]:=i; end; end; function solve(step,n:int64):int64; var pow:int64; begin if step>sum then exit(1); if f[step,n]>=0 then exit(f[step,n]); solve:=0; pow:=prime[step]; while pow<=n do begin inc(solve,solve(step+1,n-pow)); pow:=pow*prime[step]; end; inc(solve,solve(step+1,n)); f[step,n]:=solve; end; procedure intt; begin assign(input,'game.in'); assign(output,'game.out'); reset(input); rewrite(output); end; procedure outt; begin close(input); close(output); end; begin intt; sum:=0; readln(n); fillchar(f,sizeof(f),char(-1)); main; writeln(solve(1,n)); outt; end.