OFFSET
1,2
MATHEMATICA
Table[Length[Select[Divisors[n], DivisorSigma[0, #^(1+n)+n]==2^PrimeNu[#^(1+n)+n]&]], {n, 45}] (* James C. McMahon, Mar 05 2025 *)
PROG
(Magma) [#[d: d in Divisors(n) | #Divisors(d^(1+n)+n) eq 2^#PrimeDivisors(d^(1+n)+n)]: n in [1..50]];
(PARI) a(n) = sumdiv(n, d, my(f=factor(d^(1 + n) + n)); numdiv(f) == 2^omega(f)); \\ Michel Marcus, Feb 15 2025
CROSSREFS
KEYWORD
nonn,new
AUTHOR
Juri-Stepan Gerasimov, Feb 15 2025
EXTENSIONS
More terms from Jinyuan Wang, Mar 09 2025
STATUS
approved