login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A243151
Least number k not divisible by 10 such that the decimal expansion of k^n contains some digit exactly n times.
1
1, 11, 36, 34, 99, 258, 391, 163, 341, 951, 867, 1692, 1114, 793, 4792, 3019, 1935, 5469, 6398, 6152, 8906, 1987, 15815, 19603, 16292, 26216, 32113, 19718, 24354, 45903, 15776, 42202, 34572, 44411, 46911, 67972, 39291, 52299, 30499, 28383, 38001, 89782, 95017, 55954
OFFSET
1,2
COMMENTS
If k were divisible by 10, all of those numbers would work for any n and the sequence would be 1, 10, 10, 10, 10, 10, 10, 10, ....
Does a(n) exist for each n? - Charles R Greathouse IV, May 31 2014
FORMULA
a(n) > 10 for all n > 1. (Proof: check up to 21, then note that 9^22 < 10^21.) Charles R Greathouse IV, May 31 2014
EXAMPLE
1^2, 2^2, 3^2, 4^2, ... 9^2 all have different digits. 11^2 = 121 has two of the same digit. So a(2) = 11.
PROG
(Python)
def c(n):
..for k in range(10**5):
....if k%10 !=0:
......count = 0
......for i in range(10):
........if str(k**n).count(str(i)) == n:
..........return k
n = 1
while n < 100:
..print(c(n))
..n+=1
(PARI) digitct(n)=my(d=digits(n)); vector(10, i, sum(j=1, #d, d[j]==i-1))
a(n)=if(n==1, return(1)); my(k=9); until(k++%10 && #select(i->i==n, digitct(k^n)), ); k \\ Charles R Greathouse IV, May 31 2014
CROSSREFS
Sequence in context: A159493 A012644 A138893 * A225130 A191292 A107280
KEYWORD
nonn,base
AUTHOR
Derek Orr, May 31 2014
STATUS
approved