OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= proc(n) local W1, n1, W2; uses numtheory;
if n::odd then nops(factorset(n)) = 1 and nops(factorset((n+1)/2)) = 1
else nops(factorset(n/2)) = 1 and nops(factorset(n+1)) = 1
fi
end proc:
select(filter, [$1..1000]); # Robert Israel, Mar 12 2025
MATHEMATICA
Select[Range[600], PrimeNu[#*(#+1)/2] == 2 &] (* Amiram Eldar, Jan 12 2025 *)
PROG
(Magma) [k: k in [1..400] | #PrimeDivisors(k*(k+1) div 2) eq 2];
(PARI) isok(k) = omega(k*(k+1)/2) == 2; \\ Michel Marcus, Jan 14 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Jan 12 2025
STATUS
approved
