OFFSET
1,1
FORMULA
Lim_{n->infinity} a(n)/a(n-1) = 1.
EXAMPLE
9 = 3*3, 39 = 3*13, 49 = 7*7, 69 = 3*23, 99 = 3*33, 119 = 7*17, 129 = 3*43, 159 = 3*53, 169 = 13*13, 189 = 3*63 = 7*27, ...
MATHEMATICA
a={}; For[n=0, n<=100, n++, For[k=0, k<=n, k++, If[Mod[10*n+9, 10*k+3]==0 && Mod[(10*n+9)/(10*k+3), 10]==3 && 10*n+9>Max[a] || Mod[10*n+9, 10*k+7]==0 && Mod[(10*n+9)/(10*k+7), 10]==7 && 10*n+9>Max[a], AppendTo[a, 10*n+9]]]]; a
PROG
(PARI) isok(m) = ((m%10) == 9) && sumdiv(m, d, (d % 10) == (m/d % 10)); \\ Michel Marcus, Oct 22 2021
(Python)
def aupto(lim): return sorted(set(a*b for a in range(3, lim//3+1, 10) for b in range(a, lim//a+1, 10)) | set(a*b for a in range(7, lim//7+1, 10) for b in range(a, lim//a+1, 10)))
print(aupto(999)) # Michael S. Branicky, Oct 22 2021
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Stefano Spezia, Oct 22 2021
STATUS
approved