A307104 a(n) is the number which, when concatenated with A003226(n), the n-th automorphic number, gives (A003226(n))^2.

0, 0, 2, 3, 6, 57, 141, 390, 8790, 82128, 11963, 793212, 835571, 5054322, 1661682, 75880433, 45322418, 619541169, 319375992, 6745157241, 3317093849, 66891312600, 843114912509, 9837094694375, 16065496578813, 35901922360062, 67557477392256, 547721051611007

Offset

1,3

Comments

Let na and nb represent the indices of the preceding and next A003226(n)'s beginning with a 9, and where (na - nb) >= 3 (note that the first such 'zone' begins with an exception for which the index A003226(na) = 1). Then for na < n < nb and such that n == (na + 1) mod 2, it appears that A003226(n) - a(n) = A003226(n+1) - a(n+1) = k.

In such cases, it also appears that a(n)*a(n+1) = k^2 - k.

Links

Table of n, a(n) for n=1..28.

Formula

a(n) = A003226(n)*(A003226(n) - 1) / 10^A055642(A003226(n)).

Example

For n=4, A003226(4)=6, (A003226(4))^2=36. So a(4)=3.
For n=13, A003226(13)=2890625, (A003226(13))^2=8355712890625. So a(13)=835571.

Prog

(PARI) auto(n) = {n<3 & return(n-1); my(i=10, j=10, b=5, c=6, a=b); for( k=4, n, while(b<=a, b=b^2%i*=10); while(c<=a, c=(2-c)*c%j*=10); a=min(b, c)); a; } \\ A003226
a(n) = {my(m = auto(n), dm = digits(m), dm2 = digits(m^2)); fromdigits(vector(#dm2 - #dm, k, dm2[k])); } \\ Michel Marcus, May 18 2019

Crossrefs

Cf. A003226, A055642.

Sequence in context: A277374 A062309 A018390 * A180615 A018409 A226933

Adjacent sequences: A307101 A307102 A307103 * A307105 A307106 A307107

Keyword

nonn,base

Author

Christopher Hohl, Mar 24 2019

Status

approved