login
a(n) = prime(1+n)*A024451(n-1) - A024451(n), where A024451(n) is the numerator of Sum_{i = 1..n} 1/prime(i).
1

%I #10 Nov 19 2024 17:24:21

%S -1,0,4,94,284,9398,50692,2354158,75006408,445719928,53055577416,

%T 1055507209474,16217635288124,1663779309692618,115680647722172136,

%U 5515487010932110572,76943944440184239772,17660133817084175986164,686773689508904350332526,19312334585726976150166616,5276558856319725444255594528,245447734128317092747434820766

%N a(n) = prime(1+n)*A024451(n-1) - A024451(n), where A024451(n) is the numerator of Sum_{i = 1..n} 1/prime(i).

%C Note that for n > 1, A070826(n) < A024451(n) < A070826(1+n) and A070826(1+n) = prime(1+n)*A070826(n).

%H Antti Karttunen, <a href="/A376419/b376419.txt">Table of n, a(n) for n = 1..350</a>

%H <a href="/index/Pri#primorial_numbers">Index entries for sequences related to primorial numbers</a>

%o (PARI)

%o A024451(n) = numerator(sum(i=1, n, 1/prime(i))); \\ From A024451

%o A376419(n) = ((prime(1+n)*A024451(n-1)) - A024451(n));

%Y Cf. A000040, A024451, A070826, A077011.

%K sign

%O 1,3

%A _Antti Karttunen_, Nov 18 2024