A356733 Number of neighborless parts in the integer partition with Heinz number n.
0, 1, 1, 1, 1, 0, 1, 1, 1, 2, 1, 0, 1, 2, 0, 1, 1, 0, 1, 2, 2, 2, 1, 0, 1, 2, 1, 2, 1, 0, 1, 1, 2, 2, 0, 0, 1, 2, 2, 2, 1, 1, 1, 2, 0, 2, 1, 0, 1, 2, 2, 2, 1, 0, 2, 2, 2, 2, 1, 0, 1, 2, 2, 1, 2, 1, 1, 2, 2, 1, 1, 0, 1, 2, 0, 2, 0, 1, 1, 2, 1, 2, 1, 1, 2, 2, 2, 2, 1, 0, 2, 2, 2, 2, 2, 0, 1, 2, 2, 2, 1, 1, 1, 2, 0
Offset: 1
Keywords
Examples
The prime indices of 42 are {1,2,4}, of which only 4 is neighborless, so a(42) = 1.
The prime indices of 462 are {1,2,4,5}, all of which have neighbors, so a(462) = 0.
The prime indices of 1300 are {1,1,3,3,6}, with neighborless parts {1,3,6}, so a(1300) = 3.
Links
Crossrefs
Programs
-
Mathematica
primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]; Table[Length[Select[Union[primeMS[n]],!MemberQ[primeMS[n],#-1]&&!MemberQ[primeMS[n],#+1]&]],{n,100}] -
PARI
A356733(n) = if(1==n,0,my(pis=apply(primepi,factor(n)[,1])); sum(i=1, #pis, ((n%prime(pis[i]+1)) && (pis[i]==1 || (n%prime(pis[i]-1)))))); \\ Antti Karttunen, Jan 28 2025
Extensions
Data section extended to a(105) by Antti Karttunen, Jan 28 2025
Comments