OFFSET
1,3
LINKS
R. J. Mathar, Table of n, a(n) for n = 1..171
MAPLE
isA050122 := proc(n)
local i, j, f, l ;
if n =0 then
true;
else
for i from 1 do
f := A000045(i) ;
if f > n then
return false;
else
for j from i do
l := A000032(j) ;
if f*l = n then
return true;
elif f*l > n and j> 2 then
break;
end if;
end do:
end if;
end do:
end if;
end proc:
A050122 := proc(n)
option remember;
local a;
if n = 1 then
0;
else
for a from procname(n-1)+1 do
if isA050122(a) then
return a;
end if;
end do:
end if;
end proc:
seq(A050122(n), n=1..80) ; # R. J. Mathar, Jan 04 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved