%I#11 2017年4月20日09:21:20
%S 0,1,2,2,3,4,3,6,4,9,6,5,8,8,4,15,8,12,5,14,10,27,8,10,6,25,12,7,16,
%电话:16,5210,16,45,10,35,16,18,5105,16,28,11462,28,81,10,21,12,20,7,
%U 154,26125,16,30,8,49,24,11,32,32,610659212420,17910,46,75,10,78,36175,20,33,20,24,63094106315,18385,32,56,17780045
%N a(N)=A048675(A285332(N))。
%C在A285332之后,该序列也可以用二叉树的形式表示:
%C 0类
%C类|
%C。。。。。。。。。。。。。。。。。。。1...................
%C 2 2个
%C 3……../\。。。。。。。。4 4......../ \........三
%C/\/\/\/\
%C/\/\/\/\
%C/\/\/\/\
%C 6 4 9 6 5 8 4
%C 15 8 12 5 14 10 27 8 10 6 25 12 7 16 5
%C等。
%H Antti Karttunen,n的表,n=0..1023的a(n)</a>
%F a(n)=A048675(A285332(n))。
%F对于所有n>=1,a(2n)=A285332(n)。
%F a(2^n)=A109162(1+n)。[树的左边。]
%F a(A000225(n))=n。[树的右边缘。]
%o(PARI)
%o A019565(n)={my(j,v);因子回复(Mat(向量(如果(n,#n=向量提取(二进制(n),“-1..1”)),j,[素数(j),n[j]])~))};\\此函数来自_M.F.Hasler_
%o A048675(n)=我的(f=系数(n));sum(k=1,#f~,f[k,2]*2^素数pi(f[k,1]))/2;\\_米歇尔·马库斯,2016年10月10日
%o A007947(n)=因子回收(因子(n)[,1]);\\发件人:Andrew Lelechenko,2014年5月9日
%o A065642(n)={my(r=A007947(n));如果(1==n,n,n=n+r;而(A007947-(n)<>r,n=n+r);n);};
%o A285332(n)={if(n≤1,n+1,if(!(n%2),A019565(A285333(n/2)),A065642(A285322(n-1)/2));};
%o A285333(n)=如果(!n,n,如果(!(n%2),A285332(n/2),A048675(A285331(n)));
%o(方案)(定义(A285333 n)(A048675(A285323 n))
%Y参见A001477、A048675、A109162、A285325、A285330、A285 332(等分)。
%K nonn,标签
%0、3
%2017年4月19日,安蒂·卡图宁
|