OFFSET
1,1
COMMENTS
a(3k)= a(11k) =a(37k)=0. - Ray Chandler, Jul 30 2003
LINKS
Robert Israel, Table of n, a(n) for n = 1..3000
FORMULA
a(n) = n*A088394(n). - R. J. Mathar, Aug 06 2019
MAPLE
nextL:= proc(L)
local m, q, Lp;
for m from 1 do
if L[m] > 0 then
if m = LinearAlgebra:-Dimension(L) then return <4, 0$(m-1), 1>
else Lp:= L;
Lp[1]:= L[m]-1;
Lp[2..m]:= 0;
Lp[m+1]:= L[m+1]+1;
return Lp;
fi
fi
od;
end proc:
N:= 300:
Agenda:= remove(t -> ormap(s -> (t mod s = 0), [3, 11, 37, 101, 137, 239]), {$1..N}):
L:= <5>: x:= 5:
A:= Vector(N):
while Agenda <> {} and x < 10^40 do
x:= add(L[i]*10^(i-1), i=1..LinearAlgebra:-Dimension(L));
found, Agenda:= selectremove(t -> x mod t = 0, Agenda);
if found <> {} then
A[convert(found, list)]:= x;
fi;
L:= nextL(L);
od:
convert(A, list): # Robert Israel, Sep 05 2019
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Amarnath Murthy, Apr 01 2002
EXTENSIONS
More terms from Ray Chandler, Jul 30 2003
STATUS
approved