login
A337339
Denominator of (1+sigma(s)) / ((s+1)/2), where s is the square of n prime-shifted once (s = A003961(n)^2 = A003961(n^2)).
10
1, 5, 13, 41, 25, 113, 61, 365, 313, 221, 85, 1013, 145, 109, 613, 3281, 181, 2813, 265, 1985, 1513, 761, 421, 9113, 1201, 1301, 7813, 377, 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, 8821, 6613, 865
OFFSET
1,2
COMMENTS
All terms are members of A007310, because all terms of A337336 and A337337 are.
No 1's after the initial one at a(1) => No quasiperfect numbers. See comments in A336700 & A337342.
If any quasiperfect numbers qp exist, they must occur also in A325311.
Question: Is there any reliable lower bound for this sequence? See A337340, A337341.
Duplicate values are rare, but at least two cases exist: a(21) = a(74) = 1513 and a(253) = a(554) = 71065. - Antti Karttunen, Jan 03 2024
FORMULA
a(n) = A337336(n) / A337337(n) = A048673(n^2) / gcd(A048673(n^2), A336844(n^2)).
a(n) = A337336(n) / gcd(A337336(n), 1+A003973(n^2))).
PROG
(PARI)
A003961(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
A337339(n) = { my(s=(A003961(n)^2), u=(s+1)/2); (u/gcd(1+sigma(s), u)); };
\\ Or alternatively as:
A337339(n) = { my(s=A003961(n^2)); denominator((1+sigma(s))/((s+1)/2)); };
KEYWORD
nonn,frac
AUTHOR
Antti Karttunen, Aug 24 2020
STATUS
approved