A285318 a(n) = A048675(A285317(n)).

18, 36, 19, 8194, 37, 2147483648, 18014398509481984, 8195, 2056, 18014398509481985, 67108868, 16400, 2057, 158456325028528675187087900672, 2097184, 67108869, 5986310706507378352962293074805895248510699696029696, 11972621413014756705924586149611790497021399392059392, 16401, 295147905179352825858, 158456325028528675187087900673, 34359738376

Offset

1,1

Links

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

Formula

a(n) = A048675(A285317(n)).

Prog

(PARI)
A019565(n) = {my(j, v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ This function from M. F. Hasler
isA285317(n) = (issquarefree(n) & (A019565(n) < n));
A048675(n) = my(f = factor(n)); sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2; \\ Michel Marcus, Oct 10 2016
n=0; k=1; while(k <= 130, n=n+1; if(isA285317(n), write("b285318.txt", k, " ", A048675(n)); k=k+1));
(Scheme) (define (A285318 n) (A048675 (A285317 n)))

Crossrefs

Cf. A048675, A285317, A285319, A285332.

Sequence in context: A215137 A160844 A256878 * A040306 A346137 A303282

Adjacent sequences: A285315 A285316 A285317 * A285319 A285320 A285321

Keyword

nonn

Author

Antti Karttunen, Apr 18 2017

Status

approved