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.

A064779 Primes such that the sum of their digits and the sum of the reciprocals of their digits is also prime.

This page as a plain text file.
%I A064779 #7 Oct 03 2013 09:37:42
%S A064779 11,2441,4241,4421,12163,12613,13313,13331,16231,16363,16633,21163,
%T A064779 21613,26113,31663,32233,32323,32611,33113,33223,33311,48281,48821,
%U A064779 61231,61363,62131,62311,63211,63361,88241
%N A064779 Primes such that the sum of their digits and the sum of the reciprocals of their digits is also prime.
%C A064779 Zero, five, and seven never appear as a digit of any of the terms of this sequence. - _Harvey P. Dale_, Jul 17 2013
%H A064779 Harvey P. Dale, <a href="/A064779/b064779.txt">Table of n, a(n) for n = 1..1000</a>
%t A064779 f[ n_ ] := 1/n a[ n_ ] := Apply[ Plus, Map[ f, IntegerDigits[ n ] ] ] b[ n_ ] := Apply[ Plus, IntegerDigits[ n ] ] Select[ Range[ 100000 ], FreeQ[ IntegerDigits[ # ], 0 ] && PrimeQ[ a[ # ] ] && PrimeQ[ b[ # ] ] && PrimeQ[ # ] & ]
%t A064779 sdpQ[n_]:=Module[{idn=IntegerDigits[n]},Min[idn]>0&&And@@PrimeQ[{Total[ idn], Total[ 1/idn]}]]; Select[Prime[Range[10000]],sdpQ] (* _Harvey P. Dale_, Jul 17 2013 *)
%Y A064779 Cf. A034708.
%K A064779 base,nonn
%O A064779 1,1
%A A064779 _Santi Spadaro_, Oct 19 2001