OFFSET
1,1
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..130
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));
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Apr 18 2017
STATUS
approved