OFFSET
1,1
COMMENTS
The values 399..994 all seem to appear.
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..10000
MATHEMATICA
Union @@ Reap[ Do[Sow[k /. List@ToRules@ Reduce[k/x == 1 + (x - 10^n)/(10^n - 1) && 10^n <= x < 10^(n + 1) && x <= k < 2 x, {k, x}, Integers]], {n, 6}]][[2, 1]] (* Giovanni Resta, Jun 30 2017 *)
PROG
(Python 3)
from math import sqrt
def is_square(n):
root = int(sqrt(n))
return root*root == n
def find_sols(length):
count = 0
k=10**length
for n in range(k, 4*k-2):
discr= (2*k-1)*(2*k-1) - 4*(k*(k-1)-(k-1)*n)
if is_square(discr):
count+=1
b=(-(2*k-1)+sqrt(discr))/2
print(n, k+b, n/(k+b))
return count
for i in range(8):
print(find_sols(i))
CROSSREFS
KEYWORD
nonn,base
AUTHOR
James Kilfiger, Jun 15 2017
EXTENSIONS
Definition corrected by Giovanni Resta, Jun 30 2017
STATUS
approved