%I#11 2019年10月27日12:02:04
%S 1,1,1,2,1,3,3,6,5,8,7,14,15,21,31,39,51,69,98133177254329,
%电话:4716329021230171023703270459163848898124291725224230,
%电话:337834740566254928601301421824692562623596765052317100599979531404215
%N循环成分的数量(正整数的项链)与N相加,相邻部分(包括最后和第一部分)不可分割(无论哪种方式)。
%H Andrew Howroyd,n的表,n=1..100的a(n)</a>
%F a(n)=A328601(n)+1.-_安德鲁·霍罗伊,2019年10月27日
%e a(14)=14循环成分,相邻部分不可分割:
%e(14)
%e(3,11)(4,10)(5,9)(6,8)
%e(2,5,7)(2,7,5)(3,4,7)
%e(2,3,2,7)(2,3,1,5)(2,5,2,5)(2,5,4,3)(3,4,3,4)
%t颈部q[q_]:=数组[OrderedQ[{q,RotateRight[q,#]}]&,长度[q]-1,1,And];
%t表[Length[Select[Join@@Permutations/@Integer Partitions[n],Or[Length[#]==1,And[neckQ[#],And@@Not/@Divisible@@@Partition[#,2,1,1],And@@Not/@Divisible@@@Reverse/@Partition[#,2,1,1]]&]],{n,20}]
%o(PARI)
%o b(n,q,pred)={my(M=矩阵(n,n));对于(k=1,n,M[k,k]=pred(q,k);对于
%o序列(n)={my(v=总和(k=1,n,k*b(n,k,(i,j)->i%j<>0&j%i<>0));向量(n,n,1+sumdiv(n,d,eulerphi(d)*v[n/d])/n)}
%Y参见A000740、A008965、A059966、A167606、A285573、A303362、A304713、A316476、A318726、A318727、A318728、A318729、A328601。
%K nonn公司
%O 1,5型
%A _Gus Wiseman_,2018年9月2日
%E 2018年9月8日_Andrew Howroyd_的第a(21)条及以后条款
|