cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 10 results.

A007504 Sum of the first n primes.

Original entry on oeis.org

0, 2, 5, 10, 17, 28, 41, 58, 77, 100, 129, 160, 197, 238, 281, 328, 381, 440, 501, 568, 639, 712, 791, 874, 963, 1060, 1161, 1264, 1371, 1480, 1593, 1720, 1851, 1988, 2127, 2276, 2427, 2584, 2747, 2914, 3087, 3266, 3447, 3638, 3831, 4028, 4227, 4438, 4661, 4888
Offset: 0

Views

Author

Keywords

Comments

It appears that a(n)^2 - a(n-1)^2 = A034960(n). - Gary Detlefs, Dec 20 2011
This is true. Proof: By definition we have A034960(n) = Sum_{k = (a(n-1)+1)..a(n)} (2*k-1). Since Sum_{k = 1..n} (2*k-1) = n^2, it follows A034960(n) = a(n)^2 - a(n-1)^2, for n > 1. - Hieronymus Fischer, Sep 27 2012 [formulas above adjusted to changed offset of A034960 - Hieronymus Fischer, Oct 14 2012]
Row sums of the triangle in A037126. - Reinhard Zumkeller, Oct 01 2012
Ramanujan noticed the apparent identity between the prime parts partition numbers A000607 and the expansion of Sum_{k >= 0} x^a(k)/((1-x)...(1-x^k)), cf. A046676. See A192541 for the difference between the two. - M. F. Hasler, Mar 05 2014
For n > 0: row 1 in A254858. - Reinhard Zumkeller, Feb 08 2015
a(n) is the smallest number that can be partitioned into n distinct primes. - Alonso del Arte, May 30 2017
For a(n) < m < a(n+1), n > 0, at least one m is a perfect square.
Proof: For n = 1, 2, ..., 6, the proposition is clear. For n > 6, a(n) < ((prime(n) - 1)/2)^2, set (k - 1)^2 <= a(n) < k^2 < ((prime(n) + 1)/2)^2, then k^2 < (k - 1)^2 + prime(n) <= a(n) + prime(n) = a(n+1), so m = k^2 is this perfect square. - Jinyuan Wang, Oct 04 2018
For n >= 5 we have a(n) < ((prime(n)+1)/2)^2. This can be shown by noting that ((prime(n)+1)/2)^2 - ((prime(n-1)+1)/2)^2 - prime(n) = (prime(n)+prime(n-1))*(prime(n)-prime(n-1)-2)/4 >= 0. - Jianing Song, Nov 13 2022
Washington gives an oscillation formula for |a(n) - pi(n^2)|, see links. - Charles R Greathouse IV, Dec 07 2022

References

  • E. Bach and J. Shallit, §2.7 in Algorithmic Number Theory, Vol. 1: Efficient Algorithms, MIT Press, Cambridge, MA, 1996.
  • H. L. Nelson, "Prime Sums", J. Rec. Math., 14 (1981), 205-206.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

See A122989 for the value of Sum_{n >= 1} 1/a(n).

Programs

  • GAP
    P:=Filtered([1..250],IsPrime);;
    a:=Concatenation([0],List([1..Length(P)],i->Sum([1..i],k->P[k]))); # Muniru A Asiru, Oct 07 2018
    
  • Haskell
    a007504 n = a007504_list !! n
    a007504_list = scanl (+) 0 a000040_list
    -- Reinhard Zumkeller, Oct 01 2014, Oct 03 2011
    
  • Magma
    [0] cat [&+[ NthPrime(k): k in [1..n]]: n in [1..50]]; // Bruno Berselli, Apr 11 2011 (adapted by Vincenzo Librandi, Nov 27 2015 after Hasler's change on Mar 05 2014)
    
  • Maple
    s1:=[2]; for n from 2 to 1000 do s1:=[op(s1),s1[n-1]+ithprime(n)]; od: s1;
    A007504 := proc(n)
        add(ithprime(i), i=1..n) ;
    end proc: # R. J. Mathar, Sep 20 2015
  • Mathematica
    Accumulate[Prime[Range[100]]] (* Zak Seidov, Apr 10 2011 *)
    primeRunSum = 0; Table[primeRunSum = primeRunSum + Prime[k], {k, 100}] (* Zak Seidov, Apr 16 2011 *)
  • PARI
    A007504(n) = sum(k=1,n,prime(k)) \\ Michael B. Porter, Feb 26 2010
    
  • PARI
    a(n) = vecsum(primes(n)); \\ Michel Marcus, Feb 06 2021
    
  • Python
    from itertools import accumulate, count, islice
    from sympy import prime
    def A007504_gen(): return accumulate(prime(n) if n > 0 else 0 for n in count(0))
    A007504_list = list(islice(A007504_gen(),20)) # Chai Wah Wu, Feb 23 2022

Formula

a(n) ~ n^2 * log(n) / 2. - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 24 2001 (see Bach & Shallit 1996)
a(n) = A014284(n+1) - 1. - Jaroslav Krizek, Aug 19 2009
a(n+1) - a(n) = A000040(n+1). - Jaroslav Krizek, Aug 19 2009
a(A051838(n)) = A002110(A051838(n)) / A116536(n). - Reinhard Zumkeller, Oct 03 2011
a(n) = min(A068873(n), A073619(n)) for n > 1. - Jonathan Sondow, Jul 10 2012
a(n) = A033286(n) - A152535(n). - Omar E. Pol, Aug 09 2012
For n >= 3, a(n) >= (n-1)^2 * (log(n-1) - 1/2)/2 and a(n) <= n*(n+1)*(log(n) + log(log(n))+ 1)/2. Thus a(n) = n^2 * log(n) / 2 + O(n^2*log(log(n))). It is more precise than in Fares's comment. - Vladimir Shevelev, Aug 01 2013
a(n) = (n^2/2)*(log n + log log n - 3/2 + (log log n - 3)/log n + (2 (log log n)^2 - 14 log log n + 27)/(4 log^2 n) + O((log log n/log n)^3)) [Sinha]. - Charles R Greathouse IV, Jun 11 2015
G.f: (x*b(x))/(1-x), where b(x) is the g.f. of A000040. - Mario C. Enriquez, Dec 10 2016
a(n) = A008472(A002110(n)), for n > 0. - Michel Marcus, Jul 16 2020

Extensions

More terms from Stefan Steinerberger, Apr 11 2006
a(0) = 0 prepended by M. F. Hasler, Mar 05 2014

A002522 a(n) = n^2 + 1.

Original entry on oeis.org

1, 2, 5, 10, 17, 26, 37, 50, 65, 82, 101, 122, 145, 170, 197, 226, 257, 290, 325, 362, 401, 442, 485, 530, 577, 626, 677, 730, 785, 842, 901, 962, 1025, 1090, 1157, 1226, 1297, 1370, 1445, 1522, 1601, 1682, 1765, 1850, 1937, 2026, 2117, 2210, 2305, 2402, 2501
Offset: 0

Views

Author

Keywords

Comments

An n X n nonnegative matrix A is primitive (see A070322) iff every element of A^k is > 0 for some power k. If A is primitive then the power which should have all positive entries is <= n^2 - 2n + 2 (Wielandt).
a(n) = Phi_4(n), where Phi_k is the k-th cyclotomic polynomial.
As the positive solution to x=2n+1/x is x=n+sqrt(a(n)), the continued fraction expansion of sqrt(a(n)) is {n; 2n, 2n, 2n, 2n, ...}. - Benoit Cloitre, Dec 07 2001
a(n) is one less than the arithmetic mean of its neighbors: a(n) = (a(n-1) + a(n+1))/2 - 1. E.g., 2 = (1+5)/2 - 1, 5 = (2+10)/2 - 1. - Amarnath Murthy, Jul 29 2003
Equivalently, the continued fraction expansion of sqrt(a(n)) is (n;2n,2n,2n,...). - Franz Vrabec, Jan 23 2006
Number of {12,1*2*,21}-avoiding signed permutations in the hyperoctahedral group.
The number of squares of side 1 which can be drawn without lifting the pencil, starting at one corner of an n X n grid and never visiting an edge twice is n^2-2n+2. - Sébastien Dumortier, Jun 16 2005
Also, numbers m such that m^3 - m^2 is a square, (n*(1 + n^2))^2. - Zak Seidov
1 + 2/2 + 2/5 + 2/10 + ... = Pi*coth Pi [Jolley], see A113319. - Gary W. Adamson, Dec 21 2006
For n >= 1, a(n-1) is the minimal number of choices from an n-set such that at least one particular element has been chosen at least n times or each of the n elements has been chosen at least once. Some games define "matches" this way; e.g., in the classic Parker Brothers, now Hasbro, board game Risk, a(2)=5 is the number of cards of three available types (suits) required to guarantee at least one match of three different types or of three of the same type (ignoring any jokers or wildcards). - Rick L. Shepherd, Nov 18 2007
Positive X values of solutions to the equation X^3 + (X - 1)^2 + X - 2 = Y^2. To prove that X = n^2 + 1: Y^2 = X^3 + (X - 1)^2 + X - 2 = X^3 + X^2 - X - 1 = (X - 1)(X^2 + 2X + 1) = (X - 1)*(X + 1)^2 it means: (X - 1) must be a perfect square, so X = n^2 + 1 and Y = n(n^2 + 2). - Mohamed Bouhamida, Nov 29 2007
{a(k): 0 <= k < 4} = divisors of 10. - Reinhard Zumkeller, Jun 17 2009
Appears in A054413 and A086902 in relation to sequences related to the numerators and denominators of continued fractions convergents to sqrt((2*n)^2/4 + 1), n=1, 2, 3, ... . - Johannes W. Meijer, Jun 12 2010
For n > 0, continued fraction [n,n] = n/a(n); e.g., [5,5] = 5/26. - Gary W. Adamson, Jul 15 2010
The only real solution of the form f(x) = A*x^p with negative p which satisfies f^(m)(x) = f^[-1](x), x >= 0, m >= 1, with f^(m) the m-th derivative and f^[-1] the compositional inverse of f, is obtained for m=2*n, p=p(n)= -(sqrt(a(n))-n) and A=A(n)=(fallfac(p(n),2*n))^(-p(n)/(p(n)+1)), with fallfac(x,k):=Product_{j=0..k-1} (x-j) (falling factorials). See the T. Koshy reference, pp. 263-4 (there are also two solutions for positive p, see the corresponding comment in A087475). - Wolfdieter Lang, Oct 21 2010
n + sqrt(a(n)) = [2*n;2*n,2*n,...] with the regular continued fraction with period 1. This is the even case. For the general case see A087475 with the Schroeder reference and comments. For the odd case see A078370.
a(n-1) counts configurations of non-attacking bishops on a 2 X n strip [Chaiken et al., Ann. Combin. 14 (2010) 419]. - R. J. Mathar, Jun 16 2011
Also numbers k such that 4*k-4 is a square. Hence this sequence is the union of A053755 and A069894. - Arkadiusz Wesolowski, Aug 02 2011
a(n) is also the Moore lower bound on the order, A191595(n), of an (n,5)-cage. - Jason Kimberley, Oct 17 2011
Left edge of the triangle in A195437: a(n+1) = A195437(n,0). - Reinhard Zumkeller, Nov 23 2011
If h (5,17,37,65,101,...) is prime is relatively prime to 6, then h^2-1 is divisible by 24. - Vincenzo Librandi, Apr 14 2014
The identity (4*n^2+2)^2 - (n^2+1)*(4*n)^2 = 4 can be written as A005899(n)^2 - a(n)*A008586(n)^2 = 4. - Vincenzo Librandi, Jun 15 2014
a(n) is also the number of permutations simultaneously avoiding 213 and 321 in the classical sense which can be realized as labels on an increasing strict binary tree with 2n-1 nodes. See A245904 for more information on increasing strict binary trees. - Manda Riehl, Aug 07 2014
a(n-1) is the maximum number of stages in the Gale-Shapley algorithm for finding a stable matching between two sets of n elements given an ordering of preferences for each element (see Gura et al.). - Melvin Peralta, Feb 07 2016
Because of Fermat's little theorem, a(n) is never divisible by 3. - Altug Alkan, Apr 08 2016
For n > 0, if a(n) points are placed inside an n X n square, it will always be the case that at least two of the points will be a distance of sqrt(2) units apart or less. - Melvin Peralta, Jan 21 2017
Also the limit as q->1^- of the unimodal polynomial (1-q^(n*k+1))/(1-q) after making the simplification k=n. The unimodal polynomial is from O'Hara's proof of unimodality of q-binomials after making the restriction to partitions of size <= 1. See G_1(n,k) from arXiv:1711.11252. As the size restriction s increases, G_s->G_infinity=G: the q-binomials. Then substituting k=n and q=1 yields the central binomial coefficients: A000984. - Bryan T. Ek, Apr 11 2018
a(n) is the smallest number congruent to both 1 (mod n) and 2 (mod n+1). - David James Sycamore, Apr 04 2019
a(n) is the number of permutations of 1,2,...,n+1 with exactly one reduced decomposition. - Richard Stanley, Dec 22 2022
From Klaus Purath, Apr 03 2025: (Start)
The odd prime factors of these terms are always of the form 4*k + 1.
All a(n) = D satisfy the Pell equation (k*x)^2 - D*y^2 = -1. The values for k and the solutions x, y can be calculated using the following algorithm: k = n, x(0) = 1, x(1) = 4*D - 1, y(0) = 1, y(1) = 4*D - 3. The two recurrences are of the form (4*D - 2, -1). The solutions x, y of the Pell equations for n = {1 ... 14} are in OEIS.
It follows from the above that this sequence is a subsequence of A031396. (End)

Examples

			G.f. = 1 + 2*x + 5*x^2 + 10*x^3 + 17*x^4 + 26*x^5 + 37*x^6 + 50*x^7 + 65*x^8 + ...
		

References

  • S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (see p. 120).
  • E. Gura and M. Maschler, Insights into Game Theory: An Alternative Mathematical Experience, Cambridge, 2008; p. 26.
  • Thomas Koshy, Fibonacci and Lucas Numbers with Applications, John Wiley and Sons, New York, 2001.

Crossrefs

Left edge of A055096.
Cf. A059100, A117950, A087475, A117951, A114949, A117619 (sequences of form n^2 + K).
a(n+1) = A101220(n, n+1, 3).
Moore lower bound on the order of a (k,g) cage: A198300 (square); rows: A000027 (k=2), A027383 (k=3), A062318 (k=4), A061547 (k=5), A198306 (k=6), A198307 (k=7), A198308 (k=8), A198309 (k=9), A198310 (k=10), A094626 (k=11); columns: A020725 (g=3), A005843 (g=4), this sequence (g=5), A051890 (g=6), A188377 (g=7). - Jason Kimberley, Oct 30 2011
Cf. A002496 (primes).
Cf. A254858.
Subsequence of A031396.

Programs

Formula

O.g.f.: (1-x+2*x^2)/((1-x)^3). - Eric Werley, Jun 27 2011
Sequences of the form a(n) = n^2 + K with offset 0 have o.g.f. (K - 2*K*x + K*x^2 + x + x^2)/(1-x)^3 and recurrence a(n) = 3*a(n-1) - 3*a(n-2) + a*(n-3). - R. J. Mathar, Apr 28 2008
For n > 0: a(n-1) = A143053(A000290(n)) - 1. - Reinhard Zumkeller, Jul 20 2008
A143053(a(n)) = A000290(n+1). - Reinhard Zumkeller, Jul 20 2008
a(n)*a(n-2) = (n-1)^4 + 4. - Reinhard Zumkeller, Feb 12 2009
a(n) = A156798(n)/A087475(n). - Reinhard Zumkeller, Feb 16 2009
From Reinhard Zumkeller, Mar 08 2010: (Start)
a(n) = A170949(A002061(n+1));
A170949(a(n)) = A132411(n+1);
A170950(a(n)) = A002061(n+1). (End)
For n > 1, a(n)^2 + (a(n) + 1)^2 + ... + (a(n) + n - 2)^2 + (a(n) + n - 1 + a(n) + n)^2 = (n+1) *(6*n^4 + 18*n^3 + 26*n^2 + 19*n + 6) / 6 = (a(n) + n)^2 + ... + (a(n) + 2*n)^2. - Charlie Marion, Jan 10 2011
From Eric Werley, Jun 27 2011: (Start)
a(n) = 2*a(n-1) - a(n-2) + 2.
a(n) = a(n-1) + 2*n - 1. (End)
a(n) = (n-1)^2 + 2(n-1) + 2 = 122 read in base n-1 (for n > 3). - Jason Kimberley, Oct 20 2011
a(n)*a(n+1) = a(n*(n+1) + 1) so a(1)*a(2) = a(3). More generally, a(n)*a(n+k) = a(n*(n+k) + 1) + k^2 - 1. - Jon Perry, Aug 01 2012
a(n) = (n!)^2* [x^n] BesselI(0, 2*sqrt(x))*(1+x). - Peter Luschny, Aug 25 2012
a(n) = A070216(n,1) for n > 0. - Reinhard Zumkeller, Nov 11 2012
E.g.f.: exp(x)*(1 + x + x^2). - Geoffrey Critzer, Aug 30 2013
a(n) = A254858(n-2,3) for n > 2. - Reinhard Zumkeller, Feb 09 2015
Sum_{n>=0} (-1)^n / a(n) = (1+Pi/sinh(Pi))/2 = 0.636014527491... = A367976 . - Vaclav Kotesovec, Feb 14 2015
Sum_{n>=0} 1/a(n) = (1 + Pi*coth(Pi))/2 = 2.076674... = A113319. - Vaclav Kotesovec, Apr 10 2016
4*a(n) = A001105(n-1) + A001105(n+1). - Bruno Berselli, Jul 03 2017
From Amiram Eldar, Jan 20 2021: (Start)
Product_{n>=0} (1 + 1/a(n)) = sqrt(2)*csch(Pi)*sinh(sqrt(2)*Pi).
Product_{n>=1} (1 - 1/a(n)) = Pi*csch(Pi). (End)
Sum_{n>=0} a(n)/n! = 3*e. - Davide Rotondo, Feb 16 2025

Extensions

Partially edited by Joerg Arndt, Mar 11 2010

A007395 Constant sequence: the all 2's sequence.

Original entry on oeis.org

2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2
Offset: 1

Views

Author

Keywords

Comments

Continued fraction for 1 + sqrt(2). - Philippe Deléham, Nov 14 2006
a(n) = A213999(n,1). - Reinhard Zumkeller, Jul 03 2012
The least witness function W(k) is defined for odd composite numbers k. The sequence W(k) does not have its own entry in the OEIS because W(k) = 2 for all k with 9 <= k < 2047; then W(2047)=3. Cf. A089105. - N. J. A. Sloane, Sep 17 2014
a(n) = A254858(n-1,1). - Reinhard Zumkeller, Feb 09 2015
a(n) = number of permutations of length n+2 having exactly one ascent such that the first element the permutation is 2. - Ran Pan, Apr 20 2015
With alternating signs, this is the sequence of determinants of the 3 X 3 matrices m with m(i,j) = Fibonacci(n+i+j-2)^2. - Michel Marcus, Dec 23 2015
For p = prime(n+2), a(n) = ord_p(H_(p-1)), where ord_p denotes the p-adic valuation and H_i = 1 + 1/2 + ... + 1/i is a harmonic sum, except for n = 1944 and n = 157504, where ord_p(H_(p-1)) = 3, and any other term of A088164 that may exist (see Conrad link). The sequence a(n) = ord_p(H_(p-1)) does not have its own entry in the OEIS. - Felix Fröhlich, Mar 16 2016
This sequence is the only infinite bounded sequence of positive integers such that a(n) = (a(n-1) + a(n-2)) / gcd(a(n-1), a(n-2)) for all n >= 2. - Bernard Schott, Dec 28 2018

References

  • Titu Andreescu and Dorin Andrica, Number Theory, Birkhäuser, 2009, from 1999 Russian Mathematical Olympiad, p. 347.
  • Paulo Ribenboim, My Numbers, My Friends: Popular Lectures on Number Theory, Springer-Verlag, NY, 2000, p. 6.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

Formula

G.f.: 2/(1-x), and e.g.f.: 2*e^x. - Mohammad K. Azarian, Dec 22 2008
a(n) = A000005(A000040(n)). - Omar E. Pol, Feb 28 2018
a(n) = A002061(n) - A165900(n). - Torlach Rush, Feb 21 2019

A144396 The odd numbers greater than 1.

Original entry on oeis.org

3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 31, 33, 35, 37, 39, 41, 43, 45, 47, 49, 51, 53, 55, 57, 59, 61, 63, 65, 67, 69, 71, 73, 75, 77, 79, 81, 83, 85, 87, 89, 91, 93, 95, 97, 99, 101, 103, 105, 107, 109, 111, 113, 115, 117, 119, 121, 123, 125, 127, 129, 131, 133
Offset: 1

Views

Author

Paul Curtz, Oct 03 2008

Keywords

Comments

Last number of the n-th row of the triangle described in A142717.
If negated, these are also the values at local minima of the sequence A141620.
a(n) is the shortest leg of the n-th Pythagorean triple with consecutive longer leg and hypotenuse. The n-th such triple is given by (2n+1,2n^2+2n, 2n^2+2n+1), so that the longer legs are A046092(n) and the hypotenuses are A099776(n). - Ant King, Feb 10 2011
Numbers k such that the symmetric representation of sigma(k) has a pair of bars as its ends (cf. A237593). - Omar E. Pol, Sep 28 2018
Numbers k such that there is a prime knot with k crossings and braid index 2. (IS this true with "prime" removed?) - Charles R Greathouse IV, Feb 14 2023

Crossrefs

Complement of A004275 and of A004277.
Essentially the same as A140139, A130773, A062545, A020735, A005818.

Programs

Formula

a(n) = A005408(n+1) = A000290(n+1) - A000290(n).
G.f.: x*(3-x)/(1-x)^2. - Jaume Oliver Lafont, Aug 30 2009
a(n) = A254858(n-1,2). - Reinhard Zumkeller, Feb 09 2015

Extensions

Edited by R. J. Mathar, May 21 2009

A014148 a(n) = Sum_{m=1..n} Sum_{k=1..m} prime(k).

Original entry on oeis.org

2, 7, 17, 34, 62, 103, 161, 238, 338, 467, 627, 824, 1062, 1343, 1671, 2052, 2492, 2993, 3561, 4200, 4912, 5703, 6577, 7540, 8600, 9761, 11025, 12396, 13876, 15469, 17189, 19040, 21028, 23155, 25431, 27858, 30442, 33189, 36103, 39190, 42456, 45903
Offset: 1

Views

Author

Keywords

Comments

Previous name was: Apply partial sum operator twice to sequence of primes.
Numbers n such that a(n) is prime are listed in A122381(n) = {1, 2, 3, 6, 10, 23, 31, 46, 55, 58, 66, 70, 82, 91, 118, 131, 151, 163, 182, 187, 198, 199, ...}. Corresponding primes a(n) = a( A122381(n) ) = A122382(n) = {2, 7, 17, 103, 467, 6577, 17189, 61627, 109919, 130531, 198109, 239579, 399557, 559313, ...}. - Alexander Adamchuk, Aug 30 2006
Row 2 in A254858. - Reinhard Zumkeller, Feb 08 2015
Partial sums of A007504, n>=1. - Omar E. Pol, Nov 23 2016

Crossrefs

Programs

  • Haskell
    a014148 n = a014148_list !! (n-1)
    a014148_list = (iterate (scanl1 (+)) a000040_list) !! 2
    -- Reinhard Zumkeller, Feb 08 2015
  • Maple
    b:= proc(n) option remember; `if`(n<1, [0$2],
          (p-> p+[ithprime(n), p[1]])(b(n-1)))
        end:
    a:= n-> b(n+1)[2]:
    seq(a(n), n=1..42);  # Alois P. Heinz, Oct 07 2021
  • Mathematica
    Table[Sum[Sum[Prime[k],{k,1,m}],{m,1,n}],{n,1,100}] (* Alexander Adamchuk, Aug 30 2006 *)
    Accumulate[Accumulate[Prime[Range[50]]]] (* Harvey P. Dale, Dec 29 2011 *)

Formula

Convolution of the primes with the positive integers: Sum_{k=1..n} (n-k+1)*prime(k). - David Scambler, Oct 08 2006

Extensions

More terms from Alexander Adamchuk, Aug 30 2006
Name changed by Wesley Ivan Hurt, Oct 04 2021

A014150 Apply partial sum operator thrice to primes.

Original entry on oeis.org

2, 9, 26, 60, 122, 225, 386, 624, 962, 1429, 2056, 2880, 3942, 5285, 6956, 9008, 11500, 14493, 18054, 22254, 27166, 32869, 39446, 46986, 55586, 65347, 76372, 88768, 102644, 118113, 135302, 154342, 175370
Offset: 1

Views

Author

Keywords

Comments

Row 3 in A254858. - Reinhard Zumkeller, Feb 08 2015

Crossrefs

Programs

  • Haskell
    a014150 n = a014150_list !! (n-1)
    a014150_list = (iterate (scanl1 (+)) a000040_list) !! 3
    -- Reinhard Zumkeller, Feb 08 2015
  • Mathematica
    lst={};s1=0;s2=0;s3=0;Do[s1=s1+Prime[n];s2=s2+s1;s3=s3+s2;AppendTo[lst, s3], {n, 5!}];lst (* Vladimir Joseph Stephan Orlovsky, Sep 15 2008 *)
    With[{nn=3},Nest[Accumulate[#]&,Prime[Range[50]],nn]] (* Harvey P. Dale, Feb 07 2015 *)

Extensions

Offset fixed by Reinhard Zumkeller, Feb 08 2015

A178138 Apply partial sum operator 4 times to primes.

Original entry on oeis.org

2, 11, 37, 97, 219, 444, 830, 1454, 2416, 3845, 5901, 8781, 12723, 18008, 24964, 33972, 45472, 59965, 78019, 100273, 127439, 160308, 199754, 246740, 302326, 367673, 444045, 532813, 635457, 753570, 888872, 1043214, 1218584
Offset: 1

Views

Author

Jonathan Vos Post, May 20 2010

Keywords

Comments

Unlike the results of applying the partial sum operator once (A007504), twice (A014148), or thrice (A014150) to primes, this sequence begins with 4 primes. The next prime in the sequence is a(26) = 367673.
Row 4 in A254858. - Reinhard Zumkeller, Feb 08 2015

Examples

			a(15) = 2 + 9 + 26 + 60 + 122 + 225 + 386 + 624 + 962 + 1429 + 2056 + 2880 + 3942 + 5285 + 6956 = 24964 = 2^2 x 79^2.
		

Crossrefs

Programs

  • Haskell
    a178138 n = a178138_list !! (n-1)
    a178138_list = (iterate (scanl1 (+)) a000040_list) !! 4
    -- Reinhard Zumkeller, Feb 08 2015
  • Maple
    Contribution from R. J. Mathar, Oct 19 2010: (Start)
    A007504 := proc(n) option remember; add( ithprime(i),i=1..n) ; end proc:
    A014148 := proc(n) option remember; add( A007504(i),i=1..n) ; end proc:
    A014150 := proc(n) option remember; add( A014148(i),i=1..n) ; end proc:
    A178138 := proc(n) option remember; add( A014150(i),i=1..n) ; end proc: seq(A178138(n),n=1..80) ; (End)
  • Mathematica
    Nest[Accumulate[#]&,Prime[Range[40]],4] (* Harvey P. Dale, Sep 25 2014 *)

A254784 Apply partial sum operator 5 times to primes.

Original entry on oeis.org

2, 13, 50, 147, 366, 810, 1640, 3094, 5510, 9355, 15256, 24037, 36760, 54768, 79732, 113704, 159176, 219141, 297160, 397433, 524872, 685180, 884934, 1131674, 1434000, 1801673, 2245718, 2778531, 3413988, 4167558, 5056430, 6099644, 7318228, 8735337, 10376402
Offset: 1

Views

Author

Harvey P. Dale, Feb 07 2015

Keywords

Comments

Row 5 in A254858. - Reinhard Zumkeller, Feb 08 2015

Crossrefs

Programs

  • Haskell
    a254784 n = a254784_list !! (n-1)
    a254784_list = (iterate (scanl1 (+)) a000040_list) !! 5
    -- Reinhard Zumkeller, Feb 08 2015
  • Mathematica
    With[{nn=5},Nest[Accumulate[#]&,Prime[Range[50]],nn]]

A125179 Triangle read by rows: T(n,1) = prime(n) (the n-th prime); T(n,k) = 0 for k > n; T(n,k) = T(n-1,k) + T(n-1,k-1) for 2 <= k <= n (1 <= k <= n).

Original entry on oeis.org

2, 3, 2, 5, 5, 2, 7, 10, 7, 2, 11, 17, 17, 9, 2, 13, 28, 34, 26, 11, 2, 17, 41, 62, 60, 37, 13, 2, 19, 58, 103, 122, 97, 50, 15, 2, 23, 77, 161, 225, 219, 147, 65, 17, 2, 29, 100, 238, 386, 444, 366, 212, 82, 19, 2, 31, 129, 338, 624, 830, 810, 578, 294, 101, 21, 2, 37, 160, 467
Offset: 1

Views

Author

Gary W. Adamson, Nov 22 2006

Keywords

Comments

Sum of row n = A125180(n).

Examples

			Triangle starts:
   2;
   3,  2;
   5,  5,  2;
   7, 10,  7,  2;
  11, 17, 17,  9,  2;
  13, 28, 34, 26, 11,  2;
  17, 41, 62, 60, 37, 13,  2;
		

Crossrefs

Cf. A125180 (row sums), A007442, A254858 (rows reversed).
Cf. A007504.

Programs

  • Maple
    T:=proc(n,k) if k=1 then ithprime(n) elif k>n then 0 else T(n-1,k)+T(n-1,k-1) fi end: for n from 1 to 12 do seq(T(n,k),k=1..n) od; # yields sequence in triangular form
  • Mathematica
    nmax = 11;
    row[1] = Prime[Range[nmax]];
    row[n_] := row[n] = row[n-1] // Accumulate;
    T[n_, k_] := row[n][[k]];
    Table[T[n-k+1, k], {n, 1, nmax}, {k, n, 1, -1}] // Flatten (* Jean-François Alcover, Oct 11 2021 *)

Formula

T(n,2) = A007504(n-1) (n>=2).

Extensions

Edited by N. J. A. Sloane, Dec 02 2006

A293210 a(n) = [x^n] (1/(1 - x)^n)*Sum_{k>=1} prime(k)*x^k.

Original entry on oeis.org

0, 2, 7, 26, 97, 366, 1388, 5288, 20225, 77618, 298766, 1153018, 4460072, 17287558, 67129566, 261095420, 1016994627, 3966529870, 15488964428, 60549061804, 236932924494, 927984726826, 3637661249946, 14270586372348, 56024073085546, 220089137078792, 865154426179408, 3402841810234762, 13391422390407194
Offset: 0

Views

Author

Ilya Gutkovskiy, Oct 02 2017

Keywords

Crossrefs

Programs

  • Mathematica
    Table[SeriesCoefficient[1/(1 - x)^n Sum[Prime[k] x^k, {k, 1, n}], {x, 0, n}], {n, 0, 28}]

Formula

a(n) = A254858(n,n).
Showing 1-10 of 10 results.