

A332897


a(1) = 0, a(2) = 1, and for n > 2, a(n) = a(A332893(n)) + [n == 1 (mod 4)].


7



0, 1, 1, 1, 2, 1, 2, 1, 2, 2, 3, 1, 3, 2, 1, 1, 4, 2, 4, 2, 3, 3, 5, 1, 3, 3, 1, 2, 5, 1, 6, 1, 4, 4, 1, 2, 6, 4, 2, 2, 7, 3, 7, 3, 2, 5, 8, 1, 4, 3, 3, 3, 8, 1, 2, 2, 5, 5, 9, 1, 9, 6, 2, 1, 4, 4, 10, 4, 6, 1, 11, 2, 10, 6, 2, 4, 5, 2, 12, 2, 2, 7, 13, 3, 5, 7, 4, 3, 11, 2, 1, 5, 7, 8, 3, 1, 12, 4, 3, 3, 13, 3, 14, 3, 2
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,5


COMMENTS

Starting from x=n, iterate the map x > A332893(x) which divides even numbers by 2, and for odd n changes every 4k+1 prime in their prime factorization to 4k+3 prime and vice versa (except 3 > 2), like in A332819. a(n) counts the numbers of the form 4k+1 encountered until 1 has been reached, which is also included in the count when n > 1. This count includes also n itself when it is of the form 4k+1 (A016813) and larger than 1.


LINKS

Antti Karttunen, Table of n, a(n) for n = 1..65537


FORMULA

a(1) = 0, a(2) = 1, and for n > 2, a(n) = a(A332893(n)) + [n == 1 (mod 4)].
a(n) = A000120(A332895(n)).


PROG

(PARI) A332897(n) = if(n<=2, n1, A332897(A332893(n)) + (1==(n%4)));


CROSSREFS

Cf. A016813, A332893, A332895, A332898, A332899.
Cf. also A292375.
Sequence in context: A287820 A191373 A322873 * A349552 A026904 A057828
Adjacent sequences: A332894 A332895 A332896 * A332898 A332899 A332900


KEYWORD

nonn


AUTHOR

Antti Karttunen, Mar 04 2020


STATUS

approved



