OFFSET
1,2
COMMENTS
Let a(1) has r digits x_1, ..., x_r. The count of the digits 1, ..., 9 in a(1) is y_1, ..., y_9; y_i >= 0. For the digits 1, ..., 9 there exist constants c_1, ..., c_9; c_j irrational. Let k = (c_1*y_1 + ... + c_9*y_9). Then the number of digits of a(n) is approximately equal to k * (1.290554)^n and the sum of the digits of a(n) is approximately equal to 3.719315 * k * (1.290554)^n. This sequence has a(1) = 1, thus c_1*y_1 = 0.670173 * 1 and the number of digits of a(n) is approximately equal to 0.670173 * (1.290554)^n. The sum of digits of a(n) is approximately equal to 2.492585 * (1.290554)^n. This is based on empirical observations.
EXAMPLE
a(1) = 1;
a(2) = 3*1 = 3;
a(3) = 3*3 = 9;
a(4) = 3*9 = 27;
a(5) = 2/2 concatenated with 3*7 = 121;
a(6) = 3*1 concatenated with 2/2 concatenated with 3*1 = 313;
and so on.
PROG
(PARI) replace_digits(n) = my(d=digits(n), e=[]); for(k=1, #d, if(d[k]%2==0, e=concat(e, d[k]/2), my(ee=digits(3*d[k])); for(r=1, #ee, e=concat(e, ee[r])))); subst(Pol(e), x, 10)
terms(n) = my(x=1, i=0); while(i < n, print1(x, ", "); i++; x=replace_digits(x))
/* Print initial 17 terms as follows: */
terms(17) \\ Felix Fröhlich, Mar 05 2019
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Ctibor O. Zizka, Mar 05 2019
EXTENSIONS
More terms from Felix Fröhlich, Mar 05 2019
STATUS
approved