OFFSET
0,2
COMMENTS
a(n) is analogous to A033844.
FORMULA
EXAMPLE
a(7) = A005153(128) = 594.
MATHEMATICA
PracticalQ[n_] := Module[{f, p, e, prod=1, ok=True}, If[n<1||(n>1&&OddQ[n]), False, If[n==1, True, f=FactorInteger[n]; {p, e}=Transpose[f]; Do[If[p[[i]]>1+DivisorSigma[1, prod], ok=False; Break[]]; prod=prod*p[[i]]^e[[i]], {i, Length[p]}]; ok]]]; nextpractical[n1_]:=(m1=n1+1; While[!PracticalQ[m1], m1++]; m1); Table[Nest[nextpractical, 0, 2^n], {n, 0, 20}] (* using T. D. Noe's program A005153 *)
PROG
(Python)
from math import prod
from itertools import count
from sympy import factorint
def A225316(n):
if n:
c, k = 1, 1<<n
for m in count(2, 2):
f = list(factorint(m).items())
if all(f[i][0] <= 1+prod((f[j][0]**(f[j][1]+1)-1)//(f[j][0]-1) for j in range(i)) for i in range(len(f))):
c += 1
if c == k:
return m
else:
return 1 # Chai Wah Wu, Aug 04 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Frank M Jackson, May 05 2013
EXTENSIONS
a(26)-a(27) from Chai Wah Wu, Aug 07 2023
a(28)-a(29) from David A. Corneth, Aug 07 2023
More terms from David A. Corneth, Aug 14 2023
STATUS
approved