login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A298313 The first of three consecutive primes the sum of which is equal to the sum of three consecutive octagonal numbers. 2

%I #14 Dec 03 2022 13:34:18

%S 12541,75521,159617,182519,271181,373091,603901,609289,851197,983819,

%T 1246757,2079997,3299081,3687421,4484737,4692497,5636171,7514477,

%U 8273437,9299831,10408577,10430921,10746557,10769281,12739037,13012487,14213621,15440531,15713959

%N The first of three consecutive primes the sum of which is equal to the sum of three consecutive octagonal numbers.

%H Chai Wah Wu, <a href="/A298313/b298313.txt">Table of n, a(n) for n = 1..10000</a> (n = 1..70 from Colin Barker)

%e 12541 is in the sequence because 12541+12547+12553 (consecutive primes) = 37641 = 12160+12545+12936 (consecutive octagonal numbers).

%t Module[{nn=5000,oct3},oct3=Total/@Partition[PolygonalNumber[8,Range[nn]],3,1];Select[ Partition[Prime[Range[PrimePi[Ceiling[Max[oct3]/3]]]],3,1],MemberQ[ oct3,Total[ #]]&]][[All,1]] (* _Harvey P. Dale_, Dec 03 2022 *)

%o (PARI) L=List(); forprime(p=2, 20000000, q=nextprime(p+1); r=nextprime(q+1); t=p+q+r; if(issquare(36*t-180, &sq) && (sq-12)%18==0, u=(sq-12)\18; listput(L, p))); Vec(L)

%o (Python)

%o from __future__ import division

%o from sympy import prevprime, nextprime

%o A298313_list, n, m = [], 1, 30

%o while len(A298313_list) < 10000:

%o k = prevprime(m//3)

%o k2 = prevprime(k)

%o k3 = nextprime(k)

%o if k2 + k + k3 == m:

%o A298313_list.append(k2)

%o elif k + k3 + nextprime(k3) == m:

%o A298313_list.append(k)

%o n += 1

%o m += 18*n + 3 # _Chai Wah Wu_, Jan 22 2018

%Y Cf. A000040, A000567, A054643, A298073, A298168, A298169, A298222, A298223, A298250, A298251, A298272, A298273, A298301, A298302, A298312.

%K nonn

%O 1,1

%A _Colin Barker_, Jan 17 2018

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 06:44 EDT 2024. Contains 371265 sequences. (Running on oeis4.)