OFFSET
1,2
COMMENTS
Previous name was: Least integers so ascending descending base exponent transforms all semiprime.
Semiprime analogy to A113320. The sequence is probably infinite, but it is hard to characterize the asymptotic cost of adding an n-th term. The ascending descending base exponent transform of semiprimes is A113173.
The sequence is infinite because a(n) is the minimum k such that a(1)^k + k^a(1) + Sum_{i=2..n-1} a(i)^a(n-i+1) is semiprime, and since a(1)=1 this is equal to 1+k+T where T does not depend on k, thus k is the smallest positive value that makes 1+k+T semiprime, which exists because semiprimes are infinite. - Giovanni Resta, Jan 03 2020
FORMULA
a(1) = 1. For n>1, a(n) = min {k>0: a(1)^k + k^a(1) + Sum_{i=2..n-1} a(i)^a(n-i+1) is in A001358}.
EXAMPLE
a(1) = 1 by definition.
a(2) = 3 because 3 is the min x such that 1^x + x^1 is semiprime, i.e., 1^3 + 3^1 = 4 = 2*2.
a(3) = 5 because 1^5 + 3^3 + 5^1 = 33 = 3 * 11 is semiprime.
a(4) = 2 because 1^2 + 3^5 + 5^3 + 2^1 = 371 = 7 * 53.
a(5) = 4 because 1^4 + 3^2 + 5^5 + 2^3 + 4^1 = 3147 = 3 * 1049.
a(6) = 2 because 1^2 + 3^4 + 5^2 + 2^5 + 4^3 + 2^1 = 205 = 5 * 41.
a(7) = 2 because 1^2 + 3^2 + 5^4 + 2^2 + 4^5 + 2^3 + 2^1 = 1673 = 7 * 239.
a(8) = 4 because 1^4 + 3^2 + 5^2 + 2^4 + 4^2 + 2^5 + 2^3 + 4^1 = 111 = 3 * 37.
MATHEMATICA
semipQ[n_] := PrimeOmega[n] == 2; inve[w_] := Total[w^Reverse[w]]; a[1] = 1; a[n_] := a[n] = Block[{k = 0}, While[! semipQ[ inve@ Append[ Array[a, n - 1], ++k]]]; k]; Array[a, 81] (* Giovanni Resta, Jun 13 2016 *)
PROG
(PARI) lista(n)={my(a=vector(n)); a[1]=1; print1(1, ", "); for(n=2, #a, my(t=sum(i=2, n-1, a[i]^a[n-i+1])); my(k=1); while(2!=bigomega(t+1+k), k++); a[n]=k; print1(k, ", "))} \\ Andrew Howroyd, Jan 03 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Jonathan Vos Post, Jan 08 2006
EXTENSIONS
Corrected and extended by Giovanni Resta, Jun 13 2016
STATUS
approved