login
Primes p = prime(n) such that p + sum-of-digits(p) +- 1 = prime(n+1).
1

%I #6 Nov 21 2013 12:48:16

%S 2,3,23,61,131,151,331,601,661,1013,1033,1103,1123,1231,1237,1259,

%T 1307,1321,1811,2131,2621,2861,3301,4021,4159,4373,4463,4733,4759,

%U 4801,5059,5101,6151,6229,6397,6737,7079,7369,7433,8191,9109,10181,10691,11119

%N Primes p = prime(n) such that p + sum-of-digits(p) +- 1 = prime(n+1).

%H Harvey P. Dale, <a href="/A090180/b090180.txt">Table of n, a(n) for n = 1..1000</a>

%e a(2)=23. S(d)=5. 23+5=28. 28+1=29, the next prime in sequence. a(3)=61. S(d)=7. 61+7=68. 68-1=67, the next prime in sequence.

%t ppQ[n_]:=Module[{sidn=Total[IntegerDigits[n[[1]]]]},n[[1]]+sidn+1 == n[[2]]||n[[1]]+sidn-1==n[[2]]]; Transpose[Select[Partition[Prime[ Range[1500]],2,1],ppQ]][[1]] (* _Harvey P. Dale_, Mar 19 2012 *)

%Y Cf. A089323.

%K easy,nonn,base

%O 1,1

%A _Enoch Haga_, Jan 19 2004