OFFSET
0,1
COMMENTS
For n>=1, a(n) = 2^n*A014210(n). - Michel Marcus, Mar 07 2013
EXAMPLE
945=3^3*5*7 least primitive abundant number with no factor 2, a(0)=945
6=2*3 perfect number, a(1)=6
20=2^2*5 primitive abundant number < 28=2^2*7 perfect number, a(2)=20
88=2^3*11
272=2^4*17
1184=2^5*37
4288=2^6*67
16768=2^7*131
65792=2^8*257
266752=2^9*521
MATHEMATICA
{945}~Join~Array[2^#*Prime[1 + PrimePi[2^#]] &, 22] (* Michael De Vlieger, Oct 29 2023 *)
PROG
(PARI) isprab(v) = {my(sig = sigma(v)); if (sig < 2*v, return (0)); if (sig == 2*v, return (1)); fordiv (v, d, if ((d != v) && (sigma(d)>=2*d), return (0)); ); return (1); }
a(n) = {my(p = 2^n, k = 3); while (1, v = p * k; if (isprab(v), return (v)); k += 2; ); } \\ Michel Marcus, Mar 07 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Pierre CAMI, Jan 06 2008
EXTENSIONS
More terms from Michel Marcus, Mar 07 2013
STATUS
approved