OFFSET
1,1
COMMENTS
As pointed out by Chris Thompson on the SeqFan List, the terms of the sequence are of the form q+1, where p/q are increasingly better approximations of log_10(Pi), and apparently (for all known terms*) convergents of the continued fraction. At least this easily yields upper bounds for the terms. [*not true for n=13, 16, 27, 39, ...] - M. F. Hasler, Aug 31 2017
LINKS
Joseph Myers, Table of n, a(n) for n = 1..1000
Chris Thompson, in reply to F. Fröhlich, Re: Just a quick (but hard?) funny sequence idea, SeqFan list, Aug. 2017
EXAMPLE
For n = 3: The decimal expansion of Pi starts 3.1415926... and the decimal expansion of Pi^878 starts 31447652... Since 878 is the smallest k > 1 such that the initial three digits in the decimal expansions of Pi and Pi^k are equal, a(3) = 878.
PROG
(PARI) {c=contfrac(LP=log(Pi)/log(10)); n=0; for(i=1, #c=contfracpnqn(c, #c)[2, ], while(floor(10^(frac((c[i]+1)*LP)+n))==Pi\.1^n, n++; print1(c[i]+1, ", ")))} \\ Yields a non-optimal upper bound for n=13, 16, 27, ... \\ M. F. Hasler, Aug 31 2017
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Felix Fröhlich, Aug 30 2017
EXTENSIONS
a(3)-a(6) provided by Harvey P. Dale via SeqFan on Aug 26 2017
a(7)-a(8) provided by Daniel Suteu via SeqFan on Aug 26 2017
a(9)-a(10) from Hans Havermann, Aug 31 2017
a(11) onwards from Joseph Myers, Sep 01 2017
STATUS
approved