Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A243345
a(1)=1; thereafter, if n is k-th squarefree number [i.e., n = A005117(k)], a(n) = 2*a(k-1); otherwise, when n is k-th nonsquarefree number [i.e., n = A013929(k)], a(n) = 2*a(k)+1.
10
1, 2, 4, 3, 8, 6, 16, 5, 9, 12, 32, 7, 10, 18, 24, 17, 64, 13, 14, 33, 20, 36, 48, 11, 19, 34, 25, 65, 128, 26, 28, 15, 66, 40, 72, 21, 96, 22, 38, 37, 68, 50, 130, 49, 35, 256, 52, 129, 27, 29, 56, 67, 30, 41, 132, 73, 80, 144, 42, 97, 192, 44, 23, 39, 76, 74, 136, 69, 100
OFFSET
1,2
COMMENTS
Any other fixed points than 1, 2, 6, 9, 135, 147, 914, ... ?
Any other points than 4, 21, 39, 839, 4893, 12884, ... where a(n) = n-1 ?
FORMULA
a(1) = 1, and for n>1, if mu(n) = 0, a(n) = 1 + 2*a(A057627(n)), otherwise a(n) = 2*a(A013928(n)), where mu is Moebius mu function (A008683).
For all n > 1, A000035(a(n)+1) = A008966(n) = A008683(n)^2, or equally, a(n) = mu(n) + 1 modulo 2.
PROG
(Scheme, with memoizing definec-macro from Antti Karttunen's IntSeq-library)
(definec (A243345 n) (cond ((<= n 1) n) ((zero? (A008966 n)) (+ 1 (* 2 (A243345 (A057627 n))))) (else (* 2 (A243345 (A013928 n))))))
KEYWORD
nonn,look
AUTHOR
Antti Karttunen, Jun 03 2014
STATUS
approved