OFFSET
1,1
COMMENTS
Because 2 must be a prime factor and, among any run of consecutive integers, the number of composites cannot be greater than the number of primes, the largest allowable prime factor in any term is 13, and the largest allowable exponent is 12. Consequently, the admissible prime sets p_i are {2, 3, 5, 7, 11}, {2, 3, 5, 7, 11, 13}, and {2} together with any selection from {3, 5, 7}, excluding {2, 7}. This yields a finite sequence of 920 terms in total. The largest term is 2^4*3^6*5^8*7^9*11^10*13^12 < 10^42.
LINKS
Felix Huber, Table of n, a(n) for n = 1..920
EXAMPLE
162 = 2^1*3^4 is a term since 162 is even and 1, 2, 3, 4 are consecutive integers.
2000 = 2^4*5^3 is a term since 2000 is even and 2, 3, 4, 5 are consecutive integers.
MAPLE
with(combinat):
A392046List := proc(N)
local a, c, d, e, i, j, k, l, p, s, t;
a := [];
p := [op(select(x -> x <> [2, 7], map(x -> [2, op(x)], choose([3, 5, 7])))), [2, 3, 5, 7, 11], [2, 3, 5, 7, 11, 13]];
for i in p do
l := nops(i);
c := select(x -> not member(x, i), [seq(j, j = 2 .. i[l])]);
d := l - nops(c);
s := [[], [1], [1, i[l] + 1]][d + 1];
t := [[], [i[l] + 1], [i[l] + 1, i[l] + 2]][d + 1];
e := permute([op(s), op(c)]);
if nops(t) = 0 then
for j in e do
a := [op(a), mul(i[k]^j[k], k = 1 .. l)]
od
else
for j to nops(e) do
a:=[op(a), mul(i[k]^e[j, k], k = 1 .. l), mul(i[k]^permute([op(c), op(t)])[j, k], k = 1 .. l)]
od
fi
od;
select(x -> x <= N, sort(a))
end proc:
A392046List(10^42);
CROSSREFS
KEYWORD
nonn,easy,fini,full
AUTHOR
Felix Huber, Jan 05 2026
STATUS
approved
