A376419 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, 0, 4, 94, 284, 9398, 50692, 2354158, 75006408, 445719928, 53055577416, 1055507209474, 16217635288124, 1663779309692618, 115680647722172136, 5515487010932110572, 76943944440184239772, 17660133817084175986164, 686773689508904350332526, 19312334585726976150166616, 5276558856319725444255594528, 245447734128317092747434820766

Offset: 1

Author

Antti Karttunen, Nov 18 2024

Keywords

• sign

Comments

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

Links

• Antti Karttunen, Table of n, a(n) for n = 1..350
• Index entries for sequences related to primorial numbers

Crossrefs

Cf. A000040, A024451, A070826, A077011.

Programs

PARI

A024451(n) = numerator(sum(i=1, n, 1/prime(i))); \\ From A024451
A376419(n) = ((prime(1+n)*A024451(n-1)) - A024451(n));