%I#35 2021年4月28日20:41:04
%S 1,2,2,6,2,10,2,42,6,10,2146,2,10,101806,2146,10,10,223226,
%电话:6,10,42146,2314,23263442,10,10,1042814,2,10,1023226,23146,
%U 146,10,2542731938,6146,10146,22326,1023226,10,10,2141578,2141610650056950806,10
%N a(N)=1+和{d|N,d<N}a(d)^2。
%C a(n)仅依赖于n的素数签名(参见公式)_伯纳德·肖特,2021年4月24日
%H Alois P.Heinz,n的表格,n=1..3000时的a(n)</a>
%H<a href=“/index/Pri#prime_signature”>与素数签名相关的序列的索引条目</a>
%系数:x/(1-x)+和{n>=1}a(n)^2*x^(2*n)/(1-x^n)。
%对于p素数,F a(p^k)=A007018(k)。
%F From _Bernard Schott,2021年4月24日:(开始)
%签名[1,1]的F a(A006881(n))=10。
%签名[2,1]的F a(A054753(n))=146。
%F a(A007304(n))=314用于签名[1,1,1]。
%签名[3,1]的F a(A065036(n))=23226。
%签名[2,2]的F a(A085986(n))=42814。
%签名[2,1,1]的F a(A085987(n))=141578。(结束)
%p a:=proc(n)选项记忆;
%p1+加(a(d)^2,d=numtheory[除数](n)减去{n})
%p端:
%p序列(a(n),n=1..65);#_阿洛伊斯·海因茨,2021年4月17日
%ta[n_]:=a[n]=1+和[如果[d<n,a[d]^2,0],{d,除数[n]}];表[a[n],{n,65}]
%o(Python)
%o从functools导入lru_cache
%o来自sympy导入除数
%o@lru_cache(maxsize=无)
%o def A343511(n):返回1+sum(A343511(d)**2用于除数(n)中的d,如果d<n)#_Chai Wah Wu_,2021年4月17日
%o(PARI)列表a(nn)={my(va=向量(nn
%Y参见A025487、A006881、A007018、A007304、A067824、A082588、A333120。
%K nonn公司
%O 1,2号机组
%A _Ilya Gutkovskiy,2021年4月17日
|