login
Oblong numbers that are the sum of 2 successive primes.
1

%I #27 Apr 03 2023 10:36:13

%S 12,30,42,90,210,240,462,600,702,930,1482,1560,1722,2352,2862,2970,

%T 6162,6480,6642,7656,8010,8556,10920,13572,13806,14280,14762,15006,

%U 15750,16002,21462,22350,22650,23562,24492,25122,27060,27390,29070,29412,34410,34782

%N Oblong numbers that are the sum of 2 successive primes.

%C Is this sequence infinite?

%C Includes all n*(n+1) for which n*(n+1)/2 - 2 and n*(n+1)/2 + 2 are prime. The generalized Bunyakovsky conjecture implies there are infinitely many of these. - _Robert Israel_, Feb 11 2018

%H Chai Wah Wu, <a href="/A298077/b298077.txt">Table of n, a(n) for n = 1..10000</a>

%H G. L. Honaker, Jr. and Chris Caldwell, <a href="https://t5k.org/curios/page.php?curio_id=31795">Prime Curios!</a>

%e a(4)=90 because 90 is oblong (i.e., 9*10) and the sum of 2 successive primes (i.e., 43+47).

%p filter:= proc(n) not isprime(n/2) and prevprime(n/2)+nextprime(n/2) = n end proc:

%p select(filter, [seq(n*(n+1), n=2..200)]); # _Robert Israel_, Feb 11 2018

%t Select[Total /@ Partition[Prime@ Range[2^11], 2, 1], IntegerQ@ Sqrt[4 # + 1] &] (* _Michael De Vlieger_, Jan 11 2018 *)

%o (PARI) isok(n) = my(p = 2); forprime(q=3, n, if (p+q==n, return (1)); p = q);

%o lista(nn) = {for (n=1, nn, m = n*(n+1); if (isok(m), print1(m, ", ")););} \\ _Michel Marcus_, Jan 13 2018

%o (Python)

%o from __future__ import division

%o from sympy import prevprime,nextprime,isprime

%o A298077_list = [n*(n+1) for n in range(3,10**4) if prevprime(n*(n+1)//2) + nextprime(n*(n+1)//2) == n*(n+1)] # _Chai Wah Wu_, Feb 11 2018

%Y Intersection of A001043 and A002378.

%Y Cf. A111163, A258044.

%K nonn

%O 1,1

%A _G. L. Honaker, Jr._, Jan 11 2018