OFFSET
1,1
COMMENTS
This sequence contains all positive terms of A029951. So the sequence is infinite. - Altug Alkan, Oct 07 2016
LINKS
Paolo P. Lava, Table of n, a(n) for n = 1..1000
EXAMPLE
R(132) = 231, (132 * 231) / (132 + 231) = 30492 / 363 = 84.
MAPLE
R:=proc(w) local x, y, z; x:=w; y:=0; for z from 1 to ilog10(x)+1 do y:=10*y+(x mod 10); x:=trunc(x/10); od; y; end: P:= proc(q) local n; for n from 1 to q do if type(n*R(n)/(n+R(n)), integer) then print(n); fi; od; end: P(10^4);
MATHEMATICA
Select[Range@ 770, Function[r, Mod[# r, # + r] == 0]@ FromDigits@ Reverse@ IntegerDigits@ # &] (* Michael De Vlieger, Oct 14 2016 *)
PROG
(PARI) R(n) = eval(concat(Vecrev(Str(n))));
isok(n) = n*R(n) % (n+R(n)) == 0; \\ Michel Marcus, Oct 15 2016
CROSSREFS
KEYWORD
nonn,easy,base
AUTHOR
Paolo P. Lava, Oct 07 2016
STATUS
approved