A369394 AGM transform of the primes.

0, 1, 190, 29761, 9991618, 3349024561, 1787557230622, 1073002497284641, 913569251212186570, 1211439486817121619201, 1701996355944048723430570, 3350440495714062711027347281, 7769260076569386601943106748798, 18992268581018658446853739996365841, 54445901270324824915088660223022735282

Offset: 1

Author

Hugo Pfoertner, Jan 24 2024

Keywords

• nonn
• easy

Comments

See A368366 for the definition of the AGM transform.

Links

• Paolo Xausa, Table of n, a(n) for n = 1..195

Crossrefs

Cf. A368366.

Programs

Mathematica

A369394[n_] := With[{p = Prime[Range[n]]}, Total[p]^n - n^n*Apply[Times, p]];
Array[A369394, 15] (* Paolo Xausa, Jan 29 2024 *)

PARI

a369394(n) = {my(v=primes(n)); vecsum(v)^n - n^n*vecprod(v)};

Python

from sympy import prime, primorial
def A369394(n): return sum(prime(i) for i in range(1,n+1))**n-n**n*primorial(n) # Chai Wah Wu, Jan 25 2024