OFFSET
1,1
COMMENTS
a(n) is an upper bound for A130181(n) and all the more so for A126783(n); apparently even A130181(n) < a(n)/4.
All terms are divisible by 81; the quotients a(n)/81 are in A130085.
For some n (18, 34, 35, 38, 42, 58, 59, ...) the line y = x and the graph of the staircase function y = 81*(number of digits of x^n)*(number of digits of x^(n+1)) intersect twice; this possibility has to be taken into account by the program.
LINKS
Klaus Brockhaus, Table of n, a(n) for n=1..100
EXAMPLE
Let D(n,k) = 81*(number of digits of k^n)*(number of digits of k^(n+1)).
D(2,k) > k for k = 1..4641, D(2,k) = 7776 for k = 4642..9999, D(2,k) < k for k >= 10000, hence a(2) = 7776.
D(18,k) > k for k = 1..885866, D(18,k) = 997272 for k = 885867..999999, D(18,k) = 1015335 for k = 1000000..1128837, D(18,k) < k for k >= 1128838, hence a(18) = 1015335.
PROG
(PARI) {for(n=1, 28, s=30*n; k=s; while(k<81*length(Str(k^n))*length(Str(k^(n+1))), k+=s); r=0; g=0; k-=s; b=1; while(b, p=81*length(Str(k^n))*length(Str(k^(n+1))); if(r<p, r=p; h=r; if(k>r, b=0, g=h)); k++); print1(g, ", "))}
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Klaus Brockhaus, May 20 2007
STATUS
approved