login
A293219
a(n) = A000188(A292270(n)).
4
1, 1, 2, 1, 1, 5, 6, 1, 1, 9, 2, 1, 2, 11, 14, 1, 1, 1, 18, 1, 4, 5, 1, 7, 2, 1, 26, 1, 1, 29, 30, 1, 2, 33, 1, 2, 1, 1, 3, 1, 1, 41, 2, 1, 2, 1, 1, 22, 1, 1, 50, 1, 2, 53, 2, 2, 2, 7, 1, 3, 5, 4, 2, 1, 1, 65, 1, 1, 2, 69, 1, 17, 2, 1, 74, 1, 1, 7, 2, 1, 1, 81, 1, 1, 6, 7, 86, 1, 1, 89, 90, 26, 2, 1, 1, 1, 1, 1, 98, 3, 3, 6, 2, 1, 1, 105
OFFSET
0,3
LINKS
FORMULA
a(n) = A000188(A292270(n)).
PROG
(PARI)
A000188(n) = core(n, 1)[2]; \\ This function from Michel Marcus, Feb 27 2013
A000265(n) = (n >> valuation(n, 2));
A006519(n) = 2^valuation(n, 2);
A292270(n) = { my(x = n+n+1, z = ((1+x)/A006519(1+x)), m = A000265(1+x)); while(m!=1, z += ((x+m)/A006519(x+m)); m = A000265(x+m)); z; };
(Scheme) (define (A293219 n) (A000188 (A292270 n)))
CROSSREFS
Cf. A163782 (after a(1) = 1 gives the other positions where a(n) = n).
Sequence in context: A052916 A326048 A156576 * A266572 A266681 A210664
KEYWORD
nonn
AUTHOR
Antti Karttunen, Oct 07 2017
STATUS
approved