login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

A342667
Largest prime in the numerator of ratio A341528(n)/A341529(n) (when presented in its lowest terms).
4
1, 2, 3, 13, 5, 2, 7, 2, 31, 5, 11, 13, 13, 7, 3, 11, 17, 31, 19, 13, 7, 11, 23, 2, 19, 13, 13, 13, 29, 2, 31, 13, 11, 17, 5, 31, 37, 19, 13, 5, 41, 7, 43, 11, 31, 23, 47, 11, 7, 19, 17, 13, 53, 13, 11, 7, 19, 29, 59, 13, 61, 31, 31, 1093, 13, 11, 67, 17, 23, 5, 71, 31, 73, 37, 19, 19, 7, 13, 79, 11, 71, 41, 83, 13
OFFSET
1,2
COMMENTS
Equally, largest prime in the denominator of ratio A342661(A003961(n)) / A342662(A003961(n)).
FORMULA
a(n) = A006530(A341526(n)).
a(n) = A006530(A342664(A003961(n))).
PROG
(PARI)
A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
A006530(n) = if(n>1, vecmax(factor(n)[, 1]), 1);
A341526(n) = { my(s=A003961(n)); numerator((sigma(s)*n)/(sigma(n)*s)); };
KEYWORD
nonn
AUTHOR
Antti Karttunen, Mar 23 2021
STATUS
approved