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 93 results. Next

A045345 Numbers k such that k divides sum of first k primes A007504(k).

Original entry on oeis.org

1, 23, 53, 853, 11869, 117267, 339615, 3600489, 96643287, 2664167025, 43435512311, 501169672991, 745288471601, 12255356398093, 153713440932055, 6361476515268337
Offset: 1

Views

Author

Keywords

Comments

a(10) and a(11) were found by Giovanni Resta (Nov 15 2004). He states that there are no other terms for primes p < 4011201392413. See link to Prime Puzzles, Puzzle 31 below. - Alexander Adamchuk, Aug 21 2006
a(13) > pi(2*10^13). - Donovan Johnson, Aug 23 2010
a(15) > 1.42*10^13. - Giovanni Resta, Jan 07 2020
a(16) > 1.55*10^14. - Bruce Garner, Mar 06 2021
a(17) > 6.5*10^15. - Paul W. Dyson, Sep 26 2022
Numbers k such that A090396(k) = 0. - Felix Fröhlich, May 05 2021

Examples

			23 is in the sequence because the sum of the first 23 primes is 874 and that's 23 * 38.
53 is in the sequence because the sum of the first 53 primes is 5830 and that's 53 * 110.
83 is not in the sequence because the sum of the first 83 primes is 15968, which leaves a remainder of 32 when divided by 83.
The sum of the first a(14) primes is equal to a(14)*196523412770096.
		

Crossrefs

Programs

  • Mathematica
    s = 0; t = {}; Do[s = s + Prime[n]; If[ Mod[s, n] == 0, AppendTo[t, n]], {n, 1000000}]; t (* Alexander Adamchuk, Aug 21 2006 *)
    nn = 4000000; With[{acpr = Accumulate[Prime[Range[nn]]]}, Select[Range[nn], Divisible[acpr[[#]], #] &]] (* Harvey P. Dale, Sep 14 2012 *)
    Select[Range[100], Mod[Sum[Prime[i], {i, #}], #] == 0 &] (* Alonso del Arte, Mar 22 2014 based on Bill McEachen's Wolfram Alpha example *)
    A007504 = Cases[Import["https://oeis.org/A007504/b007504.txt", "Table"], {, }][[All, 2]]; Select[Range[10^5], Divisible[A007504[[# + 1]], #] &] (* Robert Price, Mar 13 2020 *)
  • PARI
    s=0;n=0;forprime(p=2,1e7,s+=p;if(s%n++==0,print1(n", "))) \\ Charles R Greathouse IV, Jul 15 2011
    
  • PARI
    isok(n) = (vecsum(primes(n)) % n) == 0; \\ Michel Marcus, Nov 26 2020
    
  • Python
    from itertools import accumulate, count, islice
    from sympy import prime
    def A045345_gen(): return (i+1 for i, m in enumerate(accumulate(prime(n) for n in count(1))) if m % (i+1) == 0)
    A045345_list = list(islice(A045345_gen(),5)) # Chai Wah Wu, Feb 23 2022

Formula

Matomäki proves that a(n) >> n^(24/19). - Charles R Greathouse IV, Jun 13 2012

Extensions

More terms from Alexander Adamchuk, Aug 21 2006
a(12) from Donovan Johnson, Aug 23 2010
a(13) from Robert Price, Mar 17 2013
a(14) from Giovanni Resta, Jan 07 2020
a(15) from Bruce Garner, Mar 06 2021
a(16) from Paul W. Dyson, Sep 26 2022

A050248 Integer averages of first k primes for some k.

Original entry on oeis.org

2, 38, 110, 3066, 60020, 740282, 2340038, 29380602, 957565746, 31043311588, 569424748566, 7207204117608, 10871205353578, 196523412770096, 2665506690112870, 122498079071529726
Offset: 1

Views

Author

Keywords

Comments

a(10) and a(11) were found by Giovanni Resta (Nov 15 2004). He states that there are no other terms for primes p < 4011201392413. See link to Prime Puzzles, Puzzle 31 below.
a(13) > (sum of first pi(2*10^13) primes)/pi(2*10^13). - Donovan Johnson, Aug 23 2010
a(16) > 2688482385366706. - Bruce Garner, Mar 06 2021
a(17) > 125237452139872271. - Paul W. Dyson, Sep 26 2022

Examples

			38 is average of first k = 23 primes; 110 (k = 53); 3066 (k = 853); 60020 (k = 11869).
		

Crossrefs

Cf. A085450 (smallest m > 1 such that m divides Sum_{k=1..m} prime(k)^n).

Programs

  • Mathematica
    a=0;lst={};Do[p=Prime[n];a+=p;If[a/n==IntegerPart[a/n],AppendTo[lst,a/n]],{n,10!}];lst (* Vladimir Joseph Stephan Orlovsky, Aug 05 2009 *)
    Module[{nn=10^6,prs},prs=Prime[Range[nn]];Select[Table[Mean[Take[prs,n]],{n,nn}],IntegerQ]] (* The program generates the first 7 terms of the sequence. *) (* Harvey P. Dale, Jun 12 2024 *)
  • PARI
    s=n=0;forprime(p=2,1e9,if((s+=p)%n++==0, print1(s/n", "))) \\ Charles R Greathouse IV, Nov 07 2014

Formula

a(n) = A050247(n)/A045345(n).

Extensions

Edited by N. J. A. Sloane at the suggestion of David W. Wilson, Jun 23 2007
a(10)-a(11) from Giovanni Resta via Ray Chandler, Jul 19 2010
a(12) from Donovan Johnson, Aug 23 2010
a(13) from Robert Price, Mar 17 2013
a(14)-a(15) from Bruce Garner, Mar 06 2021
a(16) from Paul W. Dyson, Sep 26 2022

A050247 a(n) is the sum of the first A045345(n) primes.

Original entry on oeis.org

2, 874, 5830, 2615298, 712377380, 86810649294, 794712005370, 105784534314378, 92542301212047102, 82704567079549985700, 24733255676526572596026, 3612032130800990065125528, 8102184022428756416738378
Offset: 1

Views

Author

Keywords

Comments

a(10) and a(11) were found by Giovanni Resta (Nov 15 2004). He states that there are no other terms for primes p < 4011201392413. See link to Prime Puzzles, Puzzle 31 below.
a(13) > 6640510710493148698166596 (sum of first pi(2*10^13) primes). - Donovan Johnson, Aug 23 2010
a(16) > 416714769731839517991408161209 (sum of first pi(1.55*10^14) primes). - Bruce Garner, Mar 06 2021
a(17) > 814043439429001245436559390420866 (sum of first 6500000004150767 primes). - Paul W. Dyson, Sep 27 2022

Crossrefs

Cf. A085450 (smallest m > 1 such that m divides Sum_{k=1..m} prime(k)^n).

Programs

Formula

a(n) = Sum_{i=1..A045345(n)} A000040(i). - R. J. Mathar, Jan 26 2008

Extensions

a(10)-a(11) from Giovanni Resta submitted by Ray Chandler, Jul 19 2010
a(12) from Donovan Johnson, Aug 23 2010
a(13) from Robert Price, Mar 17 2013

A217599 Prime(n), where n is such that (Sum_{i=1..n} prime(i)^2) / n is an integer.

Original entry on oeis.org

2, 67, 157, 3217, 3637, 4201, 231947, 2790569, 30116309, 12021325961, 26144296151, 1380187561637, 6549419699279, 735325088697473, 1746583001138813, 68725636353488501
Offset: 1

Views

Author

Robert Price, Mar 19 2013

Keywords

Comments

a(16) > 3*10^15 if it exists. - Anders Kaseorg, Dec 02 2020
a(17) > 3.1*10^17. - Paul W. Dyson, Jan 16 2025

Examples

			a(2) = 67, because 67 is the 19th prime and the sum of the first 19 primes^2 = 24966 when divided by 19 equals 1314 which is an integer.
		

Crossrefs

Cf. A085450 = smallest m > 1 such that m divides Sum_{k=1..m} prime(k)^n.

Programs

  • Mathematica
    t = {}; sm = 0; Do[sm = sm + Prime[n]^2; If[Mod[sm, n] == 0, AppendTo[t, Prime[n]]], {n, 100000}]; t (* T. D. Noe, Mar 19 2013 *)
    k = 1; p = 2; s = 0; lst = {}; While[p < 1000000000, s = s + p^2; If[ Mod[s, k++] == 0, AppendTo[lst, p]]; p = NextPrime@ p]; lst (* Robert G. Wilson v, Mar 08 2015 *)
  • PARI
    n=s=0; forprime(p=2,1e9, if((s+=p^2)%n++==0, print1(p", "))) \\ Charles R Greathouse IV, Feb 06 2015

Extensions

a(13) from Willem Hengeveld, Nov 29 2020
a(14)-a(15) from Anders Kaseorg, Dec 02 2020
a(16) from Paul W. Dyson, Sep 03 2022

A128165 Numbers k such that k divides 1 plus the sum of the first k primes.

Original entry on oeis.org

1, 2, 6, 10, 20, 22, 28, 155, 488, 664, 992, 6162, 7840, 7975, 8793, 18961, 32422, 148220, 231625, 332198, 459121, 462932, 2115894, 8108930, 10336641, 11789731, 15500046, 23483195, 46571611, 48582404, 77033887, 105390951, 132421841, 229481560, 1224959312
Offset: 1

Views

Author

Alexander Adamchuk, Feb 22 2007

Keywords

Comments

a(44) > 4.4*10^10. - Robert Price, Dec 15 2013
a(50) > 10^14. - Bruce Garner, Jun 05 2021

Crossrefs

Cf. A085450 (smallest m > 1 such that m divides Sum_{k=1..m} prime(k)^n).

Programs

  • Mathematica
    k = 0; s = 1; p = 2; A128165 = {}; While[k < 247336000, If[Mod[s += p, ++k] == 0, AppendTo[A128165, k]; Print[{k, p}]]; p = NextPrime@ p]; A128165
  • PARI
    is(n)=sum(i=1,n,prime(i),1)%n==0 \\ Charles R Greathouse IV, Nov 07 2014
    
  • PARI
    n=0; s=1; forprime(p=2,1e9, s+=p; if(s%n++==0, print1(n", "))) \\ Charles R Greathouse IV, Nov 07 2014

Extensions

More terms from Ryan Propper, Apr 05 2007
a(34) from Robert G. Wilson v, Jan 21 2011
a(35) from Robert Price, Dec 15 2013

A217600 Sum of the squares of the first A111441(n) primes.

Original entry on oeis.org

4, 24966, 263736, 1401992410, 2040870112, 3054955450, 346739122490032, 499159078330000800, 539391065522650998496, 25318239660367402306502991202, 251882074412384639674100925616, 31734804589156174948658730855096778, 3209990334856119248883461357431048564, 3910080232300154696097509520638192488259772
Offset: 1

Views

Author

Robert Price, Mar 19 2013

Keywords

Comments

a(n) - A111441(a(n)) - 11 == 0 (mod 24) for n > 1. This is similar to the relation between A000027 and A024450. - Karl-Heinz Hofmann, Jan 11 2021

Crossrefs

Cf. A085450 (smallest m > 1 such that m divides Sum_{k=1..m} prime(k)^n).

Programs

Extensions

a(13) from Willem Hengeveld, Nov 29 2020
a(14) from Bruce Garner, Dec 16 2020
a(15) from Bruce Garner, Dec 24 2020

A128166 Numbers k such that k divides 1 + Sum_{j=1..k} prime(j)^2 = 1 + A024450(k).

Original entry on oeis.org

1, 2, 3, 4, 6, 9, 12, 13, 26, 28, 45, 66, 174, 308, 350, 366, 417, 783, 804, 3774, 5714, 7998, 17628, 17940, 63447, 67620, 83028, 137868, 216717, 297486, 425708, 659316, 674166, 883500, 1203786, 3605052, 6778607, 9516098, 19964862, 25338586, 27771732, 70980884, 91871891, 208234138, 231967260, 238066596, 289829748, 784027092, 1078515812, 33256634230
Offset: 1

Views

Author

Alexander Adamchuk, Feb 22 2007, Feb 23 2007

Keywords

Comments

a(51) > 5.3*10^10. - Robert Price, Dec 16 2013
a(67) > 7*10^13. - Bruce Garner, May 05 2021

Crossrefs

Cf. A085450 (smallest m > 1 such that m divides Sum_{k=1..m} prime(k)^n).

Programs

  • Mathematica
    s = 1; Do[s = s + Prime[n]^2; If[ Mod[s, n] == 0, Print[n]], {n, 700000}]
    (* or *)
    Select[Range[10^4], IntegerQ[(1 + Plus@@(Prime[Range[#]]^2))/#] &] (* Alonso del Arte, Jan 20 2011 *)

Extensions

More terms from Sean A. Irvine, Jan 20 2011
a(45)-a(50) from Robert Price, Dec 16 2013

A217601 Integer averages of squares of first primes.

Original entry on oeis.org

4, 1314, 7128, 3081302, 4009568, 5312966, 16834447856, 2462344442400, 289274033242208, 46671783125431818542, 221000817555367050608, 618811172463743796896678, 13954866972387224169218132, 176536110349401666017009273532, 996528450408723697487070591774
Offset: 1

Views

Author

Robert Price, Mar 19 2013

Keywords

Examples

			a(2) = 1314 is the average of squares of first 19 primes (24966/19=1314).
		

Crossrefs

Cf. A085450 (smallest m > 1 such that m divides Sum_{k=1..m} prime(k)^n).

Programs

Formula

a(n) = A217600(n)/A111441(n).

Extensions

a(13) from Karl-Heinz Hofmann, Dec 08 2020
a(14) from Karl-Heinz Hofmann, Dec 26 2020
a(15) from Karl-Heinz Hofmann, Dec 27 2020

A233523 Prime(n), where n is such that (1+sum_{i=1..n} prime(i)) / n is an integer.

Original entry on oeis.org

2, 3, 13, 29, 71, 79, 107, 907, 3491, 4967, 7853, 61223, 80051, 81547, 90901, 211811, 381629, 1990007, 3220793, 4749637, 6725027, 6784937, 34463699, 143691323, 185831033, 213609173, 285336497, 442634651, 911588849, 953122843, 1548789581, 2153787017
Offset: 1

Views

Author

Robert Price, Dec 15 2013

Keywords

Comments

a(50) > 3475385758524527. - Bruce Garner, Jun 05 2021

Examples

			a(3) = 13, because 13 is the 6th prime and the sum of the first 6 primes+1 = 42 when divided by 6 equals 7 which is an integer.
		

Crossrefs

Cf. A085450 = smallest m > 1 such that m divides Sum_{k=1..m} prime(k)^n.

Programs

  • Mathematica
    t = {}; sm = 1; Do[sm = sm + Prime[n]; If[Mod[sm, n] == 0, AppendTo[t, Prime[n]]], {n, 100000}]; t (* Derived from A217599 *)
  • PARI
    is(n)=if(!isprime(n),return(0)); my(t=primepi(n),s); forprime(p=2,n,s+=Mod(p,t)); s==0 \\ Charles R Greathouse IV, Nov 30 2013

A233862 Prime(k), where k is such that (1 + Sum_{i=1..k} prime(i)^2) / k is an integer.

Original entry on oeis.org

2, 3, 5, 7, 13, 23, 37, 41, 101, 107, 197, 317, 1033, 2029, 2357, 2473, 2879, 5987, 6173, 35437, 56369, 81769, 195691, 199457, 793187, 850027, 1062931, 1840453, 2998421, 4217771, 6200923, 9914351, 10153807, 13563889, 18878099, 60767923, 118825361, 170244929
Offset: 1

Views

Author

Robert Price, Dec 16 2013

Keywords

Comments

a(51) > 1428199016921.
a(67) > 2407033812270611. - Bruce Garner, May 05 2021

Examples

			a(5) = 13, because 13 is the 6th prime and the sum of the first 6 primes^2+1 = 378 when divided by 6 equals 63 which is an integer.
		

Crossrefs

Cf. A085450 (smallest m > 1 such that m divides Sum_{k=1..m} prime(k)^n).

Programs

  • Mathematica
    t = {}; sm = 1; Do[sm = sm + Prime[n]^2; If[Mod[sm, n] == 0, AppendTo[t, Prime[n]]], {n, 100000}]; t (* Derived from A217599 *)
    Module[{nn=9600000},Prime[#]&/@Transpose[Select[Thread[{Range[nn], 1+ Accumulate[ Prime[Range[nn]]^2]}],IntegerQ[Last[#]/First[#]]&]][[1]]] (* Harvey P. Dale, Sep 09 2014 *)
  • PARI
    is(n)=if(!isprime(n),return(0)); my(t=primepi(n),s); forprime(p=2,n,s+=Mod(p,t)^2); s==0 \\ Charles R Greathouse IV, Nov 30 2013
Showing 1-10 of 93 results. Next