|
|
A279119
|
|
Lexicographically earliest sequence such that, for any distinct i and j, a(i)=a(j) implies gcd(i, j)=1.
|
|
3
|
|
|
0, 0, 0, 1, 0, 2, 0, 3, 1, 4, 0, 5, 0, 6, 3, 7, 0, 8, 0, 9, 4, 10, 0, 11, 1, 12, 6, 13, 0, 14, 0, 15, 7, 16, 2, 17, 0, 18, 9, 19, 0, 20, 0, 21, 10, 22, 0, 23, 1, 24, 12, 25, 0, 26, 5, 27, 13, 28, 0, 29, 0, 30, 15, 31, 6, 32, 0, 33, 16, 34, 0, 35, 0, 36, 18, 37
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
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).
For any prime p, the sequence b_p(n)=a(p*n) is a bijection from A000027 to A001477:
- 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
|
|
|
FORMULA
|
a(2*n) = n-1 for any n>0.
|
|
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
|
|
|
STATUS
|
approved
|
|
|
|