OFFSET
1,1
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..10000
FORMULA
MAPLE
b:= proc(n) option remember; local k;
for k from 1+`if`(n=1, 3, b(n-1))
while isprime(k) do od; k
end:
a:= n-> mul((i[1]+1)^i[2], i=ifactors(b(n))[2]):
seq(a(n), n=1..100); # Alois P. Heinz, Oct 23 2014
MATHEMATICA
f[n_] := Module[{p, e},
Product[{p, e} = pe; (p+1)^e, {pe, FactorInteger[n]}]];
f /@ Select[Range[1000], CompositeQ] (* Jean-François Alcover, Jul 11 2021 *)
PROG
(PARI) lista(nn) = {forcomposite(c=2, nn, f = factor(c); print1(prod(i=1, #f~, (f[i, 1]+1)^f[i, 2]), ", "); ); } \\ Michel Marcus, Jun 26 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Ben Paul Thurston, May 06 2007
EXTENSIONS
More terms from Michel Marcus, Jun 26 2013
STATUS
approved