OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
MAPLE
filter:= proc(n) local t; uses NumberTheory;
t:= tau(n);
ormap(d -> tau(n/d+d) = t, Divisors(n))
end proc:
select(filter, [$1..300]); # Robert Israel, Jan 08 2026
MATHEMATICA
Select[Range[230], Sum[Boole[DivisorSigma[0, #]==DivisorSigma[0, #/d+d]], {d, Divisors[#]}]>0 &] (* Stefano Spezia, Jan 08 2026 *)
PROG
(Magma) [k: k in [1..230] | not #[d: d in Divisors(k) | #Divisors(k) eq #Divisors((k div d) + d)] eq 0];
(PARI) isok(k) = #select(x->(x==numdiv(k)), apply(x->numdiv(k/x+x), divisors(k))); \\ Michel Marcus, Jan 08 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Jan 07 2026
STATUS
approved
