|
|
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
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
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
|
|
LINKS
|
|
|
FORMULA
|
|
|
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)); };
|
|
CROSSREFS
|
Cf. A003961, A003973, A007310, A048673, A074627, A074630, A209922, A324899, A325311, A336697, A336700, A336844, A337194, A337336, A337337, A337340, A337341, A337342.
|
|
KEYWORD
|
nonn,frac
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|