A285326 a(0) = 0, for n > 0, a(n) = n + A006519(n).
0, 2, 4, 4, 8, 6, 8, 8, 16, 10, 12, 12, 16, 14, 16, 16, 32, 18, 20, 20, 24, 22, 24, 24, 32, 26, 28, 28, 32, 30, 32, 32, 64, 34, 36, 36, 40, 38, 40, 40, 48, 42, 44, 44, 48, 46, 48, 48, 64, 50, 52, 52, 56, 54, 56, 56, 64, 58, 60, 60, 64, 62, 64, 64, 128, 66, 68, 68, 72, 70, 72, 72, 80, 74, 76, 76, 80, 78, 80, 80, 96, 82, 84, 84, 88, 86, 88, 88, 96, 90, 92, 92
Offset: 0
Links
Crossrefs
Programs
-
Mathematica
Table[If[n>0, n + GCD[2^n, n], 0], {n, 0, 100}] (* Indranil Ghosh, Apr 20 2017 *)
-
PARI
a(n) = if(n>0, n + gcd(2^n, n), 0); \\ Indranil Ghosh, Apr 20 2017
-
PARI
A285326(n)=n+bitand(n,-n) \\ Or: {a(n)=-bitand(-n,bitneg(n))}, not faster. - M. F. Hasler, Oct 19 2019
-
Python
from sympy import gcd def a(n): return n + gcd(2**n, n) if n>0 else 0 # Indranil Ghosh, Apr 20 2017
-
Scheme
(define (A285326 n) (if (zero? n) n (+ n (A006519 n)))) (define (A285326 n) (A048675 (A065642 (A019565 n))))
Comments