[go: up one dir, main page]

login
A285327
Row 3 of A285325: a(n) = A048675(A065642(A065642(A019565(n)))).
4
0, 3, 6, 5, 12, 7, 10, 9, 24, 11, 18, 13, 20, 15, 18, 17, 48, 19, 22, 21, 36, 23, 26, 25, 40, 27, 34, 29, 36, 31, 34, 33, 96, 35, 38, 37, 68, 39, 42, 41, 72, 43, 50, 45, 52, 47, 50, 49, 80, 51, 54, 53, 68, 55, 58, 57, 72, 59, 66, 61, 68, 63, 66, 65, 192, 67, 70, 69, 132, 71, 74, 73, 136, 75, 82, 77, 84, 79, 82, 81, 144, 83, 86, 85, 100, 87, 90, 89
OFFSET
0,2
LINKS
FORMULA
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
A048675(n) = my(f = factor(n)); sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2; \\ Michel Marcus, Oct 10 2016
A007947(n) = factorback(factorint(n)[, 1]); \\ From Andrew Lelechenko, May 09 2014
A065642(n) = { my(r=A007947(n)); if(1==n, n, n = n+r; while(A007947(n) <> r, n = n+r); n); };
(Scheme) (define (A285327 n) (A048675 (A065642 (A065642 (A019565 n)))))
CROSSREFS
Row 3 of A285325 (after the initial zero).
Sequence in context: A007479 A076535 A349102 * A239487 A246979 A246980
KEYWORD
nonn
AUTHOR
Antti Karttunen, Apr 19 2017
STATUS
approved