OFFSET
1,6
COMMENTS
Also, for n>1, a(n) equals the index of the class of n relatively to the algorithm described in A275246 (i.e., if a(n)=k, then n is of class P_k).
From Rémy Sigrist, Dec 21 2016: (Start)
- b_p is injective: b_p(n)=b_p(m) implies p*n=p*m or gcd(p*n,p*m)=1; as p>1, gcd(p*n,p*m)>1, so p*n=p*m and n=m.
- b_p is surjective: by contradiction: let k be the least number such that b_p(n) never equals k; we have a set of k terms (i_1,...,i_k) such that b_p(i_j) = j-1 for any j between 1 and k; let l be the least value such that p^l > max({1, i_1,...,i_k}). Then, by definition of a, a(p^l)=k, and b_p(p^(l-1))=k, which is a contraction.
(End)
LINKS
Rémy Sigrist, Table of n, a(n) for n = 1..10000
FORMULA
a(2*n) = n-1 for any n>0.
a(n)=0 iff n belongs to A008578.
a(n)=1 iff n belongs to A001248.
a(n)=2 iff n belongs to A089581.
a(n)=3 iff n belongs to A275246.
a(n)=4 iff n belongs to A275248.
a(n)=5 iff n belongs to A275249.
a(n)=6 iff n belongs to A275251.
a(n)=7 iff n belongs to A275252.
a(n)=8 iff n belongs to A275253.
PROG
(PARI) g = vector(76, i, 1); for (n=1, #g, a = 0; while (gcd(g[a+1], n)>1, a++); g[a+1] *= n; print1 (a ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Rémy Sigrist, Dec 06 2016
STATUS
approved