login
Primes that are the sum of the numerator and denominator of a harmonic number.
2

%I #9 May 30 2019 11:51:42

%S 2,5,17,37,197,503,9649,9901,111431,1532093,1556117,71354639,

%T 392849685203,67564865715709,67929521152309,2551756396706233,

%U 345246028268898566519,73744425850091533377950933665103,4418089000986947559043789457061782529431

%N Primes that are the sum of the numerator and denominator of a harmonic number.

%C A001008(n) + A002805(n) = A064168(n) is the sum of numerator and denominator in n-th harmonic number H(n) = 1 + 1/2 + 1/3 +...+ 1/n. A064168(n) is prime for n = {1,2,3,4,5,7,9,10,11,14,15,20,27,35,36,39,48..} = A123538.

%H Robert Israel, <a href="/A118727/b118727.txt">Table of n, a(n) for n = 1..46</a>

%e H(n) = [ 1/1, 3/2, 11/6, 25/12, 137/60, 49/20, 363/140, 761/280, 7129/2520,... ].

%e a(1) = 1 + 1 = 2,

%e a(2) = 3 + 2 = 5,

%e a(3) = 11 + 6 = 17,

%e a(4) = 25 + 12 = 37,

%e a(5) = 137 + 60 = 197,

%e a(6) = 363 + 140 = 503,

%e a(7) = 7129 + 2520 = 9649...

%p select(isprime,map(numer+denom, ListTools:-PartialSums([(1/n) $ n=1..10^3]))); # _Robert Israel_, May 30 2019

%t Select[Numerator[#]+Denominator[#]&/@HarmonicNumber[Range[100]],PrimeQ] (* _Harvey P. Dale_, Jan 29 2014 *)

%Y Cf. A064168, A001008, A002805, A123538.

%K nonn

%O 1,1

%A _Alexander Adamchuk_, May 21 2006