A083814 Duplicate of A073805.
1, 2, 4, 6, 10, 16, 18, 20, 24, 25, 26, 28, 36, 42, 50, 52, 58, 61, 62, 63, 66, 68, 70, 80, 81
Offset: 1
This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
a(1) = 2 = 1 * 1 + 1. a(3) = 11 = 10 * 1 + 1. a(13) = 977 = 16 * 61 + 1.
N:= 6: # for terms <= 10^N where N is even S:= {}: for x from 1 to 10^(N/2)-1 do if x mod 10 = 0 then next fi; r:= rev(x); if r < x then next fi; v:= x*r; for i from 0 do w:= 10^i*v+1; if w > 10^N then break fi; if isprime(w) then S:= S union {w} fi; od od: sort(convert(S,list));
1 is a term because 10^1+1 = 11 and its digit reversal is 11, which is prime. 2 is a term because 10^1+1 + 10^2+2 = 113 and its digit reversal is 311, a prime. 3 is not a term because 10^1+1 + 10^2+2 + 10^3+3 = 1116 and R(1116) = 6111, not prime.
for(n=1,400,my(s=fromdigits(Vecrev(digits(sum(k=1,n,10^k+k)))));if(ispseudoprime(s),print1(n", ")));
Comments