OFFSET
1,2
COMMENTS
All terms must be odd. - Harvey P. Dale, Jul 29 2015
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..3294
EXAMPLE
137 belongs since 51375 is a multiple 137 (137*375).
MAPLE
a:=proc(n) local nn: nn:=convert(n, base, 10): if type((5+10*n+5*10^(nops(nn)+1))/n, integer)=true then n else fi end: seq(a(n), n=1..10000); # Emeric Deutsch, Feb 28 2006
MATHEMATICA
f[k_, d_] := Flatten@Table[Select[Divisors[k*(10^(i + 1) + 1)], IntegerLength[ # ] == i &], {i, d}]; f[5, 10] (* Ray Chandler, May 11 2007 *)
PROG
(Python)
from sympy import isprime
from itertools import count, islice
def agen(): # generator of terms
yield from [1, 5]
for k in count(2):
t = 5*(10**(k+1) + 1)
yield from (t//i for i in range(500, 50, -1) if t%i == 0)
print(list(islice(agen(), 33))) # Michael S. Branicky, Mar 26 2023
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Giovanni Resta, Feb 15 2006
EXTENSIONS
a(31) and beyond from Michael S. Branicky, Mar 26 2023
STATUS
approved