login
A337336
a(n) = A048673(n^2).
3
1, 5, 13, 41, 25, 113, 61, 365, 313, 221, 85, 1013, 145, 545, 613, 3281, 181, 2813, 265, 1985, 1513, 761, 421, 9113, 1201, 1301, 7813, 4901, 481, 5513, 685, 29525, 2113, 1625, 2965, 25313, 841, 2381, 3613, 17861, 925, 13613, 1105, 6845, 15313, 3785, 1405, 82013, 7321, 10805, 4513, 11705, 1741, 70313, 4141, 44105
OFFSET
1,2
COMMENTS
All terms are odd and neither there are multiples of 3, thus only terms of A007310 occur here.
FORMULA
a(n) = A048673(A000290(n)) = (1+(A003961(n)^2))/2.
For all n>= 1, A010872(a(n)) = A040001(n).
PROG
(PARI)
A048673(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); (factorback(f)+1)/2; };
A337336(n) = A048673(n^2);
KEYWORD
nonn
AUTHOR
Antti Karttunen, Aug 24 2020
STATUS
approved