cp's OEIS Frontend

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.

A048521 Primes expressible as the sum of an integer plus its digit sum.

This page as a plain text file.
%I A048521 #19 Oct 15 2019 11:56:50
%S A048521 2,11,13,17,19,23,29,37,41,43,47,59,61,67,71,73,79,83,89,101,103,107,
%T A048521 109,113,127,131,137,139,149,151,157,163,167,173,179,181,191,193,197,
%U A048521 199,223,227,229,239,241,251,257,263,269,271,281,283,293,307,311,313
%N A048521 Primes expressible as the sum of an integer plus its digit sum.
%H A048521 Reinhard Zumkeller, <a href="/A048521/b048521.txt">Table of n, a(n) for n = 1..10000</a>
%H A048521 <a href="/index/Coi#Colombian">Index entries for Colombian or self numbers and related sequences</a>
%F A048521 A107740(A049084(a(n))) > 0.
%e A048521 a(24) = prime 113 which is 106 + (1+0+6) (or 97 + (9+7)).
%t A048521 t={};Do[p=Prime[n];c=0;i=1;While[i<p&&c!=1,If[i+Total[IntegerDigits[i]] ==p,c=1;AppendTo[t,p]];i++],{n,65}];t (* _Jayanta Basu_, May 03 2013 *)
%t A048521 Union[Select[Table[n+Total[IntegerDigits[n]],{n,400}],PrimeQ]] (* _Harvey P. Dale_, Jul 14 2014 *)
%o A048521 (Haskell)
%o A048521 a048521 n = a048521_list !! (n-1)
%o A048521 a048521_list = map a000040 $ filter ((> 0) . a107740) [1..]
%o A048521 -- _Reinhard Zumkeller_, Sep 27 2014
%Y A048521 Cf. A007953, A047791, A048520.
%Y A048521 Cf. A006378, A062028, A107741.
%Y A048521 Cf. A000040, A107740, A049084.
%K A048521 nonn,base
%O A048521 1,1
%A A048521 _Patrick De Geest_, May 15 1999
%E A048521 Formula and also offset corrected by _Reinhard Zumkeller_, Sep 27 2014