OFFSET
1,2
COMMENTS
a(n) is the 2nd elementary symmetric function of the first n+1 primes.
Using the identity that (x_1 + x_2 + ... + x_n)^2 - (x_1^2 + x_2^2 + ... + x_n^2) is the sum of the products taken two at a time, a(n) can be expressed with the sum of the primes and the sum of the prime squared. Since they both have asymptotic formulas, this yields an asymptotic formula for this sequence. - Timothy Varghese, May 06 2014
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
a(1) = 0, a(n+1) = prime(n+1)*(sum of first n primes) + a(n), for n > 1.
a(n) ~ (3*n^4*log^2(n) - 4*n^3*log^2(n))/24. - Timothy Varghese, May 06 2014
MAPLE
Primes:= [seq](ithprime(i), i=1..100):
(map(`^`, ListTools:-PartialSums(Primes), 2) - ListTools:-PartialSums(map(`^`, Primes, 2)))/2; # Robert Israel, Sep 24 2015
MATHEMATICA
a[1] = 0; a[n_] := a[n] = a[n-1] + Prime[n]*Total[Prime[Range[n-1]]];
Array[a, 35] (* Jean-François Alcover, Feb 28 2019 *)
PROG
(PARI) /* Extra memory allocation could be required. */
Primes=List();
forprime(x=2, prime(500000), listput(Primes, x));
/* Keep previous lines global, before a(n) */
a(n)={my(p=vector(n, j, Primes[j]), s=0); forvec(y=vector(2, i, [1, #p]), s+=(p[y[1]]*p[y[2]]), 2); s} \\ R. J. Cano, Oct 11 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved