A328523 a(n) = number of terms in A002182 that are of the form A002182(n)/p for some prime p|n.
0, 1, 1, 1, 2, 1, 1, 1, 1, 2, 2, 2, 2, 2, 1, 1, 2, 3, 3, 1, 1, 2, 1, 1, 1, 1, 2, 2, 2, 2, 3, 2, 2, 2, 2, 3, 2, 1, 1, 2, 3, 2, 2, 3, 2, 3, 3, 2, 2, 1, 3, 2, 1, 3, 1, 2, 2, 2, 2, 3, 2, 2, 4, 3, 1, 4, 2, 2, 2, 2, 1, 2, 2, 2, 2, 4, 3, 4, 3, 2, 3, 1, 1, 1, 1, 2, 2, 3
Offset: 1
Keywords
Programs
-
PARI
apply( {a(n)=#[0|p<-factor(n=A2182[n])[,1], setsearch(A2182,n/p)]}, [1..99]) \\ Uses global var. A2182 which must be assigned to a vector of at least n elements of that sequence. - M. F. Hasler, Jan 08 2020
Comments