OFFSET
1,1
COMMENTS
Includes 6*(10^k+1) for k >= 1. - Robert Israel, Jul 12 2019
LINKS
EXAMPLE
66 = k + reverse(k) for k = 15, 24, 33, 42, 51, 60; 626 = k + reverse(k) for k = 115, 214, 313, 412, 511, 610.
MAPLE
N:= 10^5:
revdigs:= proc(n) local L, i;
L:= convert(n, base, 10);
add(L[-i]*10^(i-1), i=1..nops(L))
end proc:
V:= Vector(N):
for x from 1 to N do
v:= x + revdigs(x);
if v <= N then V[v]:= V[v]+1 fi;
od:
select(t -> V[t]=6, [$1..N]); # Robert Israel, Jul 12 2019
PROG
(ARIBAS) var n, k, c, i, rev: integer; st, nst: string; end; m := 6; for n := 0 to 17500 do k := n div 2; c := 0; while k <= n and c < m + 1 do st := itoa(k); nst := ""; for i := 0 to length(st) - 1 do nst := concat(st[i], nst); end; rev := atoi(nst); if n = k + rev then inc(c); if k mod 10 <> 0 and k <> rev then inc(c); end; end; inc(k); end; if c = m then write(n, ", "); end; end;
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Klaus Brockhaus, Jun 17 2002
EXTENSIONS
Offset changed to 1 by Robert Israel, Jul 12 2019
STATUS
approved