%I#23 2014年12月23日00:36:28
%S 1,1,2,2,3,3,5,4,6,5,7,6,11,7,10,8,13,9,17,10,14,11,19,12,15,13,18,14,
%电话23,15,29,16,22,17,21,18,31,19,26,20,37,21,41,22,30,23,43,24,35,25,34,
%U 26,47,27,33,28,38,29,53,30,59,31,42,32,39,33,61,34,46,35,67,36,71,37,50,38,55,39,73,40,54,41,79,42
%N将最小素数的一个实例向较小素数移动一步:a(1)=1,a(2n)=N,对于奇数>1:a(N)=(N/素数(s))*素数(s-1),其中s=A055396(N),最小素数除以N的指数。
%C在A056239(n)次迭代后,从任意n迭代为a(n)、a(a(n。
%C偶数平分得到自然数A000027,从第三项开始的奇数平分是A129128:2,3,5,6,7,11,10。。。
%H Antti Karttunen,n的表,n=1..10000的a(n)</a>
%F a(1)=1;对于n>1:a(n)=A008578(A055396(n))*A032742(n)。[与A064989的类似公式进行比较。]
%F其他身份。对于所有n>=1:
%F a(2n)=n。
%F如果n是奇数,则A001222(a(n))=A001225(n)。
%F如果n是偶数,则A001222(a(n))=A001225(n)-1。
%t a252461[n_Integer]:=块[{a008578,a032742,a055396,a},
%t a008578[x_]:=如果[x==1,1,素数[x-1]];
%t a032742[x_]:=如果[x==1,1,除数[x][[2]]];
%t a055396[x_]:=PrimePi[FactorInteger[x][[1]][[1];
%ta[1]=1;
%ta[x_]:=a008578[a055396[x]]*a032742[x];
%t数组[a,n]];a252461[84](*迈克尔·德弗里格,2014年12月21日*)
%o(方案)(定义(A252461 n)(如果(=1 n)n(*(A008578(A055396 n))(A032742 n)))
%Y变体:A252462、A252463。
%Y等分:A000027,A129128。
%Y参考A001222、A008578、A032742、A055396、A056239、A064989。
%K nonn公司
%氧1,3
%2014年12月20日,安蒂·卡图宁
|