OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
sPracQ[n_] := sPracQ[n] = Module[{sd = Sort[Map[DivisorSigma[1, #] - # &, Divisors[n]]], sm, ac}, sm = Plus @@ sd; ac = Accumulate[sd]; AllTrue[Range[Length[sd] - 1], sd[[# + 1]] <= 1 + ac[[#]] &]];
Select[Partition[Range[4300], 3, 1], AllTrue[#, sPracQ] &][[;; , 1]]
PROG
(PARI) issprac(k) = {my(sd = vecsort(apply(x -> sigma(x) - x, divisors(k))), sm = vecsum(sd), s = 0); for(i = 1, #sd-1, s += sd[i]; if(sd[i+1] > 1 + s, return(0))); 1; }
isok(k) = issprac(k) && issprac(k+1) && issprac(k+2);
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Feb 16 2026
STATUS
approved
