OFFSET
1,2
COMMENTS
Of the first 30 terms, all except a(1) and a(24) are primes.
LINKS
EXAMPLE
a(6) = 19 because 19 = 1 + (2 + 2 * 2 * 2 * 2 ), and 19 cannot be built with five 1 and 2's.
a(10) = 107 because 107 = 1 + (2 + 2 * 2 * (2 + 2 * 2 * (2 + 2 * 2 ))), and 107 cannot be built with nine 1 and 2's.
MAPLE
N:= 10^6: # for terms <= N
M[1]:= {1, 2}: T[1]:= M[1]: A:= 1:
for n from 2 do
M[n]:= `union`(seq({seq(seq(x+y, x = select(`<=`, M[i], N-y)), y=M[n-i])}, i=1..n/2),
seq({seq(seq(x*y, x = select(`<=`, M[i], N/y)), y=M[n-i])}, i=1..n/2)) minus T[n-1];
T[n]:= T[n-1] union M[n];
if M[n] = {} then break fi;
A:= A, min(M[n]);
od:
A; # Robert Israel, Jun 09 2025
MATHEMATICA
num=1500;
b=Array[99999&, num]; a={};
b[[1]]=b[[2]]=1;
r=1;
Do[Do[s=b[[k]]+b[[n/k]]; If[s<b[[n]], b[[n]]=s], {k, Select[Divisors[n], 1<#^2<=n&]}];
Do[s=b[[k]]+b[[n-k]]; If[s<b[[n]], b[[n]]=s], {k, n/2}]; If[b[[n]]==r, Print[{r, n}]; AppendTo[a, n]; r++], {n, num}]; a
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Zhining Yang, Jun 02 2025
EXTENSIONS
a(28)-a(32) from Hongyang Cao, Jun 10 2025
STATUS
approved
