|
%I #28 Jun 01 2017 10:18:39
%S 55410683,102291263,141430363,226383163,280064453,457433213,531290533,
%T 542418463,555695713,582949903,629444003,664652203,665813153,
%U 777809113,802919653,852404053,887653633,894328243,898734673,979048313,993517643
%N Initial prime in set of 4 consecutive primes with common gap 48.
%H Chai Wah Wu, <a href="/A067388/b067388.txt">Table of n, a(n) for n = 1..1000</a> (terms 1..102 from Zak Seidov)
%H <a href="/index/Pri#primes_AP">Index entries for sequences related to primes in arithmetic progressions</a>
%t Transpose[Select[Partition[Prime[Range[10000000]], 4, 1], Union[Differences[#]]=={48}&]][[1]] (* _Vincenzo Librandi_, Jun 21 2015 *)
%o (Python)
%o from sympy import isprime, nextprime
%o A067388_list, p = [], 2
%o q, r, s = p+48, p+96, p+144
%o while s <= 10**10:
%o np = nextprime(p)
%o if np == q and isprime(r) and isprime(s) and nextprime(q) == r and nextprime(r) == s:
%o A067388_list.append(p)
%o p, q, r, s = np, np+48, np+96, np+144 # _Chai Wah Wu_, Jun 01 2017
%Y Analogous sequences (with differences in square brackets): A033451[6], A033447[12], A033448[18], A052242[24], A052243[30], A058252[36], A058323[42], this sequence[48].
%K nonn
%O 1,1
%A _Labos Elemer_, Jan 21 2002
%E a(7)-a(21) from _Donovan Johnson_, Sep 05 2008
|
