OFFSET
1,1
LINKS
David A. Corneth, Table of n, a(n) for n = 1..50
Eric Weisstein's World of Mathematics, Almost Prime.
EXAMPLE
For n = 2: 19 = 2*2 + 3*5 = 3*3 + 2*5.
MAPLE
f:= proc(n) uses priqueue; local pq, S, t, x, y, k, i, p, v, R;
initialize(pq);
insert([-2^n, 2$n], pq);
S[0]:= 1:
for i from 1 to n do S[i]:= 0 od:
do
t:= extract(pq);
x:= -t[1];
for i from n to 1 by -1 do
S[i]:= expand(S[i] + S[i-1] * y^x);
od;
if type(S[n], `+`) then
R:= select(t -> degree(t, y) < x and eval(t, y=1) = n, convert(S[n], list));
if R <> [] then return min(map(t -> degree(t, y), R)) fi;
fi;
p:= nextprime(t[-1]);
for i from n+1 to 2 by -1 while t[i] = t[-1] do
v:= x*(p/t[-1])^(n+2-i);
insert([-v, op(t[2..i-1]), p$(n+2-i)], pq)
od;
od;
end proc:
map(f, [$1..19]); # Robert Israel, Jun 10 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Sep 07 2023
EXTENSIONS
a(5)-a(19) from Robert Israel, Jun 10 2025
More terms from David A. Corneth, Jun 10 2025
STATUS
approved
