OFFSET
1,2
COMMENTS
If k = Product_{i=1..m} p_i^e_i is in the sequence, then p_i = 1 + Sum_{j<=i} e_j. Thus the primes in the factorization are obtained as cumulative sums of the exponents plus 1.
The only prime term is 2, the only squarefree terms are 1, 2 and 6.
LINKS
Felix Huber, Table of n, a(n) for n = 1..10000
EXAMPLE
150 = 2^1*3^1*5^2 is a term since 2 - 1 = 1 (exponent of 2), 3 - 2 = 1 (exponent of 3), and 5 - 3 = 2 (exponent of 5).
21609 = 3^2*7^4 is a term since 3 - 1 = 2 (exponent of 3), and 7 - 3 = 4 (exponent of 7).
MAPLE
A392195 := proc(N)
local c, i, f, m, t;
m := evalf(ln(N));
t := table();
c := 0;
f := proc(l, v, k)
local i, p, w;
c := c + 1;
t[c] := v;
i := k;
do
p := ithprime(i);
if (p - l)*evalf(ln(p)) > m then
break
end if;
w := v*p^(p - l);
if w <= N then
f(p, w, i + 1)
end if;
i := i + 1;
end do
end proc;
f(1, 1, 1);
sort([seq(t[i], i = 1 .. c)]);
end proc:
A392195(18186318150);
CROSSREFS
KEYWORD
nonn
AUTHOR
Felix Huber, Mar 11 2026
STATUS
approved
