OFFSET
1,1
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
FORMULA
Lim_{n->infinity} a(n)/a(n-1) = 1.
EXAMPLE
16 = 4*4, 36 = 6*6, 56 = 4*14, 96 = 4*24 = 6*16, 136 = 4*34, 156 = 6*26, ...
MATHEMATICA
a={}; For[n=0, n<=150, n++, For[k=0, k<=n, k++, If[Mod[10*n+6, 10*k+4]==0 && Mod[(10*n+6)/(10*k+4), 10]==4 && 10*n+6>Max[a] || Mod[10*n+6, 10*k+6]==0 && Mod[(10*n+6)/(10*k+6), 10]==6 && 10*n+6>Max[a], AppendTo[a, 10*n+6]]]]; a
tisdQ[n_]:=AnyTrue[{Mod[#, 10], Mod[n/#, 10]}&/@Divisors[n], #[[1]] == #[[2]]&]; Select[10 Range[150]+6, tisdQ] (* Harvey P. Dale, Dec 27 2021 *)
PROG
(Python)
def aupto(lim): return sorted(set(a*b for a in range(4, lim//4+1, 10) for b in range(a, lim//a+1, 10)) | set(a*b for a in range(6, lim//6+1, 10) for b in range(a, lim//a+1, 10)))
print(aupto(1317)) # Michael S. Branicky, Sep 12 2021
(PARI) isok(m) = if ((m % 10) == 6, fordiv(m, d, if ((d % 10) == (m/d % 10), return(1)))); \\ Michel Marcus, Oct 06 2021
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Stefano Spezia, Sep 12 2021
STATUS
approved