A055265
a(n) is the smallest positive integer not already in the sequence such that a(n)+a(n-1) is prime, starting with a(1)=1.
Original entry on oeis.org
1, 2, 3, 4, 7, 6, 5, 8, 9, 10, 13, 16, 15, 14, 17, 12, 11, 18, 19, 22, 21, 20, 23, 24, 29, 30, 31, 28, 25, 34, 27, 26, 33, 38, 35, 32, 39, 40, 43, 36, 37, 42, 41, 48, 49, 52, 45, 44, 53, 50, 47, 54, 55, 46, 51, 56, 57, 70, 61, 66, 65, 62, 69, 58, 73, 64, 63, 68, 59, 72, 67, 60
Offset: 1
a(5) = 7 because 1, 2, 3 and 4 have already been used and neither 4 + 5 = 9 nor 4 + 6 = 10 are prime while 4 + 7 = 11 is prime.
Cf.
A086527 (the primes a(n)+a(n-1)).
Cf.
A070942 (n's such that a(1..n) is a permutation of (1..n)). -
Zak Seidov, Oct 19 2011
See
A282695 for deviation from identity sequence.
A073659 is a version where the partial sums must be primes.
-
import Data.List (delete)
a055265 n = a055265_list !! (n-1)
a055265_list = 1 : f 1 [2..] where
f x vs = g vs where
g (w:ws) = if a010051 (x + w) == 1
then w : f w (delete w vs) else g ws
-- Reinhard Zumkeller, Feb 14 2013
-
A055265 := proc(n)
local a,i,known ;
option remember;
if n =1 then
1;
else
for a from 1 do
known := false;
for i from 1 to n-1 do
if procname(i) = a then
known := true;
break;
end if;
end do:
if not known and isprime(procname(n-1)+a) then
return a;
end if;
end do:
end if;
end proc:
seq(A055265(n),n=1..100) ; # R. J. Mathar, Feb 25 2017
-
f[s_List] := Block[{k = 1, a = s[[ -1]]}, While[ MemberQ[s, k] || ! PrimeQ[a + k], k++ ]; Append[s, k]]; Nest[f, {1}, 71] (* Robert G. Wilson v, May 27 2009 *)
q=2000; a={1}; z=Range[2,2*q]; While[Length[z]>q-1, k=1; While[!PrimeQ[z[[k]]+Last[a]], k++]; AppendTo[a,z[[k]]]; z=Delete[z,k]]; Print[a] (*200 times faster*) (* Vladimir Joseph Stephan Orlovsky, May 03 2011 *)
-
v=[1];n=1;while(n<50,if(isprime(v[#v]+n)&&!vecsearch(vecsort(v),n), v=concat(v,n);n=0);n++);v \\ Derek Orr, Jun 01 2015
-
U=-a=1; vector(100,k, k=valuation(1+U+=1<M. F. Hasler, Feb 11 2020
A329449
For any n >= 0, exactly four sums a(n+i) + a(n+j) are prime, for 0 <= i < j <= 3: lexicographically earliest such sequence of distinct nonnegative integers.
Original entry on oeis.org
0, 1, 2, 3, 4, 9, 8, 15, 14, 5, 26, 17, 6, 11, 12, 7, 30, 29, 24, 13, 18, 19, 10, 43, 28, 31, 16, 25, 22, 21, 46, 37, 52, 27, 34, 45, 44, 39, 58, 69, 20, 51, 32, 41, 38, 35, 48, 23, 36, 53, 50, 47, 54, 59, 42, 55, 72, 65, 84, 67, 114, 79, 60, 49, 78, 71, 102, 61, 66, 91, 40, 73, 76, 33, 64, 63, 68
Offset: 0
We start with a(0) = 0, a(1) = 1, a(2) = 2, a(3) = 3, the smallest possibilities which do not lead to a contradiction. Indeed, the four sums 0 + 2, 0 + 3, 1 + 2 and 2 + 3 are prime.
Now we have 2 prime sums using {1, 2, 3}, so the next term must give two more prime when added to these. We find that a(4) = 4 is the smallest possible choice, with 1 + 4 = 5 and 3 + 4 = 7.
Then there are again 2 primes among the pairwise sums using {2, 3, 4}, so the next term must again produce two more prime sums. We find that a(5) = 9 is the smallest possibility, with 2 + 9 = 11 and 4 + 9 = 13.
a(10^4) = 9834 and all numbers up to 9834 occurred by then.
a(10^5) = 99840 and all numbers below 99777 occurred by then.
a(10^6) = 1000144 and all numbers below 999402 occurred by then.
- Eric Angelini, Prime sums from neighbouring terms, personal blog "Cinquante signes" (and post to the SeqFan list), Nov. 11, 2019.
- Eric Angelini, Prime sums from neighbouring terms [Cached copy of html file, with permission]
- Eric Angelini, Prime sums from neighbouring terms [Cached copy of pdf file, with permission]
- M. F. Hasler, Prime sums from neighboring terms, OEIS Wiki, Nov. 23, 2019.
Other sequences with N primes among pairwise sums of M consecutive terms, starting with a(o) = o, sorted by decreasing N and lowest possible M:
A329581 (N=11, M=8, o=0),
A329580 (N=10, M=8, o=0),
A329569 (N=9, M=6, o=0),
A329568 (N=9, M=6, o=1),
A329425 (N=6, M=5, o=0),
A329449 (N=4, M=4, o=0),
A329411 (N=2, M=3, o=0 or 1),
A128280 (N=1, M=2, o=0),
A055265 (N=1, M=2, o=1),
A055266 (N=0, M=2; o=1),
A253074 (N=0, M=2; o=0).
-
A329449(n, show=0, o=0, N=4, M=3, p=[], U, u=o)={for(n=o, n-1, if(show>0, print1(o", "), show<0, listput(L,o)); U+=1<<(o-u); U>>=-u+u+=valuation(U+1, 2); p=concat(if(#p>=M, p[^1], p), o); my(c=N-sum(i=2, #p, sum(j=1, i-1, isprime(p[i]+p[j])))); for(k=u, oo, bittest(U, k-u) || min(c-#[0|p<-p, isprime(p+k)], #p>=M) || [o=k, break]));show&&print([u]); o} \\ Optional args: show=1: print a(o..n-1), show=-1: append a(o..n-1) to the global list L, in both cases print [least unused number] at the end; o=1: start with a(1)=1; N, M: get N primes using M+1 consecutive terms.
A329411
Among the pairwise sums of any three consecutive terms there are exactly two prime sums: lexicographically earliest such sequence of distinct positive numbers.
Original entry on oeis.org
1, 2, 3, 4, 7, 6, 5, 8, 9, 10, 13, 16, 15, 14, 17, 12, 11, 18, 19, 22, 21, 20, 23, 24, 29, 30, 31, 28, 25, 33, 34, 26, 27, 32, 35, 36, 37, 42, 41, 38, 45, 44, 39, 40, 43, 46, 51, 50, 47, 53, 54, 48, 49, 52, 55, 57, 82, 56, 75, 62, 64, 87, 63, 76, 61, 66, 65, 71, 86, 60, 77, 67, 72, 59, 68, 69, 58, 70
Offset: 1
a(1) = 1 is the smallest possible choice; there's no restriction on the first term.
a(2) = 2 as 2 is the smallest available integer not leading to a contradiction. Note that as 1 + 2 = 3 we already have one prime sum (out of the required two) with the pair {1, 2}.
a(3) = 3 as 3 is the smallest available integer not leading to a contradiction. Since 2 + 3 = 5 we now have our two prime sums with the triplet {1, 2, 3}.
a(4) = 4 as 4 is the smallest available integer not leading to a contradiction. Since 3 + 4 = 7 we now have our two prime sums with the triplet {2, 3, 4}: they are 2 + 3 = 5 and 3 + 4 = 7.
a(5) = 7 because 5 or 6 would lead to a contradiction: indeed, both the triplets {3, 4, 5} and {3, 4, 6} will produce only one prime sum (instead of two). With a(5) = 7 we have the triplet {3, 4, 7} and the two prime sums we were looking for: 3 + 4 = 7 and 4 + 7 = 11.
And so on.
- Jean-Marc Falcoz, Table of n, a(n) for n = 1..10000
- Eric Angelini, Prime sums from neighbouring terms [Cached copy of html file, with permission]
- Eric Angelini, Prime sums from neighbouring terms [Cached copy of pdf file, with permission]
- M. F. Hasler, Prime sums from neighboring terms, OEIS wiki, Nov. 23, 2019
Cf.
A055265 (sum of two consecutive terms is always prime: differs from a(30) on).
Cf.
A329412 ..
A329416 (exactly 2 prime sums using 4, ..., 10 consecutive terms).
Cf.
A055266 (no prime sum among 2 consecutive terms),
A329405 (no prime among the pairwise sums of 3 consecutive terms).
-
a[1]=1;a[2]=2;a[n_]:=a[n]=(k=1;While[Length@Select[Plus@@@Subsets[{a[n-1],a[n-2],++k},{2}],PrimeQ]!=2||MemberQ[Array[a,n-1],k]];k);Array[a,100] (* Giorgos Kalogeropoulos, May 09 2021 *)
-
A329411(n,show=0,o=1,N=2,M=2,p=[],U,u=o)={for(n=o,n-1, show>0&& print1(o", "); show<0&& listput(L,o); U+=1<<(o-u); U>>=-u+u+=valuation(U+1,2); p=concat(if(#p>=M, p[^1], p), o); my(c=N-sum(i=2,#p, sum(j=1,i-1, isprime(p[i]+p[j])))); for(k=u,oo, bittest(U,k-u)|| min(c-#[0|p<-p, isprime(p+k)], #p>=M) ||[o=k,break]));show&&print([u]);o} \\ Optional args: show=1: print a(o..n-1), show=-1: append a(o..n-1) to the (global) list L, in both cases print [least unused number] at the end; o=0: start with a(o)=o; N, M: find N primes using M+1 consecutive terms. - M. F. Hasler, Nov 16 2019
A329416
Among the pairwise sums of any ten consecutive terms there are exactly two prime sums: lexicographically earliest such sequence of distinct positive numbers.
Original entry on oeis.org
1, 2, 3, 7, 13, 19, 23, 25, 31, 32, 17, 8, 26, 37, 43, 49, 14, 38, 55, 61, 11, 20, 35, 67, 73, 79, 57, 9, 5, 15, 21, 42, 27, 12, 33, 30, 39, 45, 47, 18, 48, 6, 51, 24, 63, 69, 72, 75, 16, 36, 54, 60, 22, 66, 10, 4, 40, 29, 28, 34, 44, 41, 46, 50, 52, 58, 64, 53, 70, 71, 59, 62, 76, 56, 82, 88, 94, 65, 100
Offset: 1
a(1) = 1 is the smallest possible choice, there's no restriction on the first term.
a(2) = 2 as 2 is the smallest available integer not leading to a contradiction. Note that as 1 + 2 = 3 we already have one prime sum (on the required two) with the 10-set {1,2,a(3),a(4),a(5),a(6),a(7),a(8),a(9),a(10)}.
a(3) = 3 as 3 is the smallest available integer not leading to a contradiction. Note that as 2 + 3 = 5 we now have the two prime sums required with the 10-set {1,2,a(3),a(4),a(5),a(6),a(7),a(8),a(9),a(10)}.
a(4) = 7 as a(4) = 4, 5 or 6 would lead to a contradiction: indeed, the 10-sets {1,2,3,4,a(5),a(6),a(7),a(8),a(9),a(10)}, {1,2,3,5,a(5),a(6),a(7),a(8),a(9),a(10)} and {1,2,3,6,a(5),a(6),a(7),a(8),a(9),a(10)} will produce more than the two required prime sums. With a(4) = 7 we have no contradiction as the 10-set {1,2,3,7,a(5),a(6),a(7),a(8),a(9),a(10)} has now two prime sums so far: 1 + 2 = 3 and 2 + 3 = 5.
a(5) = 13 as a(5) = 4, 5, 6, 8, 9, 10, 11 or 12 would again lead to a contradiction (more than 2 prime sums with the 10-set); in combination with any other term before it, a(5) = 13 will produce only composite sums.
a(6) = 19 as 19 is the smallest available integer not leading to a contradiction: indeed, the 10-set {1,2,3,7,13,19,a(7),a(8),a(9),a(10)} shows two prime sums so far: 1 + 2 = 3 and 2 + 3 = 5.
a(7) = 23 as 23 is the smallest available integer not leading to a contradiction; indeed, the 10-set {1,2,3,7,13,19,23,a(8),a(9),a(10)} shows only two prime sums so far, which are 1 + 2 = 3 and 2 + 3 = 5.
a(8) = 25 as 25 is the smallest available integer not leading to a contradiction and producing two prime sums so far with the 10-set {1,2,3,7,13,19,23,25,a(9),a(10)}; etc.
Cf.
A329333 (3 consecutive terms, exactly 1 prime sum).
Cf.
A329405 (no prime among the pairwise sums of 3 consecutive terms).
Cf.
A329406 ..
A329410 (exactly 1 prime sum using 4, ..., 10 consecutive terms).
Cf.
A329411 ..
A329415 (exactly 2 prime sums using 3, ..., 7 consecutive terms).
See also "nonnegative" variants:
A329450 (0 primes using 3 terms),
A329452 (2 primes using 4 terms),
A329453 (2 primes using 5 terms),
A329454 (3 primes using 4 terms),
A329449 (4 primes using 4 terms),
A329455 (3 primes using 5 terms),
A329456 (4 primes using 5 terms).
-
A329416(n, show=0, o=1, N=2, M=9, p=[], U, u=o)={for(n=o, n-1, show&&print1(o", "); U+=1<<(o-u); U>>=-u+u+=valuation(U+1, 2); p=concat(if(#p>=M, p[^1], p), o); my(c=N-sum(i=2, #p, sum(j=1, i-1, isprime(p[i]+p[j])))); if(#pM. F. Hasler, Nov 15 2019
A329566
For all n >= 0, exactly six sums are prime among a(n+i) + a(n+j), 0 <= i < j < 6; lexicographically earliest such sequence of distinct nonnegative numbers.
Original entry on oeis.org
0, 1, 2, 3, 4, 24, 5, 7, 6, 8, 9, 10, 11, 13, 18, 19, 16, 12, 28, 31, 17, 15, 14, 22, 26, 20, 21, 27, 23, 30, 32, 80, 41, 38, 51, 39, 62, 29, 35, 44, 34, 45, 54, 25, 49, 33, 64, 36, 37, 40, 46, 61, 47, 42, 43, 55, 66, 58, 65, 48, 72, 79, 52, 53, 59, 78, 50, 57, 60, 89, 71, 56, 68, 63, 74, 75, 76, 69, 82, 81, 67, 91, 88, 70, 100
Offset: 0
For n = 0, we consider pairwise sums of the first 6 terms a(0..5) = (0, 1, 2, 3, 4, 24): We have (a(i) + a(j), 0 <= i < j < 6) = (1; 2, 3; 3, 4, 5; 4, 5, 6, 7; 24, 25, 26, 27, 28) among which there are 6 primes, counted with repetition. This justifies taking a(0..4) = (0, ..., 4), the smallest possible choices for these first 5 terms. Since no smaller a(5) between 5 and 23 has this property, this is the start of the lexicographically earliest nonnegative sequence with this property and no duplicate terms.
Then we find that a(6) = 5 is possible, also giving 6 prime sums for n = 1, so this is the correct continuation (modulo later confirmation that the sequence can be continued without contradiction given this choice).
Next we find that a(7) = 6 is not possible, it would give only 5 prime sums using the 6 consecutive terms (2, 3, 4, 24, 5, 6). However, a(7) = 7 is a valid continuation, and so on.
Cf.
A329425 (6 primes using 5 consecutive terms).
Cf.
A329449 (4 primes using 4 consecutive terms),
A329456 (4 primes using 5 consecutive terms).
Cf.
A329454 (3 primes using 4 consecutive terms),
A329455 (3 primes using 5 consecutive terms).
Cf.
A329411 (2 primes using 3 consecutive terms),
A329452 (2 primes using 4 consecutive terms),
A329453 (2 primes using 5 consecutive terms).
-
A329566(n,show=0,o=0,N=6,M=5,p=[],U,u=o)={for(n=o,n-1, if(show>0,print1(o", "), show<0,listput(L,o)); U+=1<<(o-u); U>>=-u+u+=valuation(U+1,2); p=concat(if(#p>=M,p[^1],p),o); my(c=N-sum(i=2,#p, sum(j=1,i-1, isprime(p[i]+p[j]))));if(#p
A329563
For all n >= 1, exactly five sums are prime among a(n+i) + a(n+j), 0 <= i < j < 5; lexicographically earliest such sequence of distinct positive numbers.
Original entry on oeis.org
1, 2, 3, 4, 5, 8, 9, 14, 6, 23, 17, 7, 12, 24, 10, 13, 19, 16, 18, 25, 22, 15, 28, 21, 26, 32, 75, 20, 11, 27, 56, 30, 41, 53, 29, 38, 60, 44, 35, 113, 36, 31, 48, 61, 37, 42, 46, 33, 34, 55, 39, 40, 49, 58, 45, 43, 52, 51, 106, 57, 62, 50, 87, 47, 54, 59, 80, 66, 83, 68
Offset: 1
For n = 1, we consider pairwise sums among the first 5 terms chosen as small as possible, a(1..5) = (1, 2, 3, 4, 5). We see that we have indeed 5 primes among the sums 1+2, 1+3, 1+4, 1+5, 2+3, 2+4, 2+5, 3+4, 3+5, 4+5.
Then, to get a(6), consider first the pairwise sums among terms a(2..5), (2+3, 2+4, 2+5; 3+4, 3+5; 4+5), among which there are 3 primes, counted with multiplicity (i.e., the prime 7 is there two times). So the new term a(6) must give exactly two more prime sums with the terms a(2..5). We find that 6 or 7 would give just one more (5+6 resp. 4+7), but a(6) = 8 gives exactly two more, 3+8 and 5+8.
Cf.
A329425 (6 primes using 5 consecutive terms),
A329566 (6 primes using 6 consecutive terms).
Cf.
A329449 (4 primes using 4 consecutive terms),
A329456 (4 primes using 5 consecutive terms).
Cf.
A329454 (3 primes using 4 consecutive terms),
A329455 (3 primes using 5 consecutive terms).
Cf.
A329411 (2 primes using 3 consecutive terms),
A329452 (2 primes using 4 consecutive terms),
A329453 (2 primes using 5 consecutive terms).
-
{A329563(n,show=1,o=1,N=5,M=4,p=[],u=o,U)=for(n=o,n-1, show>0&& print1(o","); show<0&& listput(L,o); U+=1<<(o-u); U>>=-u+u+=valuation(U+1,2); p=concat(if(#p>=M,p[^1],p),o); my(c=N-sum(i=2,#p, sum(j=1,i-1, isprime(p[i]+p[j])))); if(#p
A329572
For all n >= 0, exactly 12 sums are prime among a(n+i) + a(n+j), 0 <= i < j < 7; lexicographically earliest such sequence of distinct nonnegative numbers.
Original entry on oeis.org
0, 1, 2, 5, 6, 11, 12, 17, 26, 35, 36, 47, 24, 54, 77, 7, 43, 60, 13, 30, 96, 4, 67, 97, 16, 133, 34, 3, 40, 27, 63, 100, 10, 20, 171, 9, 8, 51, 21, 22, 52, 15, 32, 38, 75, 141, 56, 41, 71, 122, 152, 45, 68, 29, 59, 14, 39, 44, 50, 23, 53, 57, 74, 107, 170, 176, 93, 134, 137, 86, 177, 65, 476, 62, 87, 92, 101
Offset: 0
-
{A329572(n,show=0,o=0,N=12,M=6,D=[3,5,4,6,5,11],p=[],u=o,U)=for(n=o+1,n, show>0&& print1(o","); show<0&& listput(L,o); U+=1<<(o-u); U>>=-u+u+=valuation(U+1,2); p=concat(if(#p>=M,p[^1],p),o); D&& D[1]==n&& [o=D[2],D=D[3..-1]]&& next; my(c=N-sum(i=2,#p, sum(j=1,i-1, isprime(p[i]+p[j])))); for(k=u,oo,bittest(U,k-u)|| min(c-#[0|p<-p,isprime(p+k)],#p>=M)|| [o=k,break]));show&&print([u]);o} \\ optional args: show=1: print a(o..n-1), show=-1: append them on global list L, in both cases print [least unused number] at the end. See the wiki page for more.
A329573
For all n >= 1, exactly 12 sums are prime among a(n+i) + a(n+j), 0 <= i < j < 7; lexicographically earliest such sequence of distinct positive numbers.
Original entry on oeis.org
1, 2, 3, 4, 9, 10, 27, 14, 20, 33, 34, 69, 39, 28, 40, 13, 19, 70, 31, 43, 180, 220, 61, 36, 66, 91, 127, 7, 12, 5, 102, 186, 11, 6, 25, 18, 55, 41, 42, 48, 65, 72, 59, 38, 125, 24, 29, 35, 54, 32, 47, 77, 164, 26, 407, 15, 116, 63, 75, 404, 416, 8, 215, 45, 56, 183, 23, 134, 206, 17, 44, 50
Offset: 1
Up to and including the 6th term, there is no constraint other than not using a term more than once, since it is impossible to have more than 12 primes as pairwise sums of 6 numbers. So one would first try to use the lexicographically smallest possible choice a(1..6) =?= (1, 2, ..., 6). But then one would have only 7 pairs (i,j) such that a(i) + a(j) is prime, 1 <= i < j <= 6. So one would need 12 - 7 = 5 more primes in {1, 2, ..., 6} + a(7), which is impossible. One can check that even a(1..5) =?= (1,...,5) does not allow one to find a(6) and a(7) in order to have 12 prime sums a(i) + a(j), 1 <= i < j <= 7. Nor is it possible to find a solution with a(5) equal to 6 or 7 or 8. One finds that a(5) = 9, and a(6) = 10, are the smallest possible choices for which a(7) can be found as to satisfy the requirement. In that case, a(7) = 27 is the smallest possible solution, which yields the 12 prime sums 1+2, 2+3, 1+4, 3+4, 2+9, 4+9, 1+10, 3+10, 9+10, 2+27, 4+27, 10+27.
Now, to satisfy the definition of the sequence for n = 2, we drop the initial 1 from the set of consecutive terms, and search for a(8) producing the same number of additional primes together with {2, 3, 4, 9, 10, 27} as did a(1) = 1, namely 3. We see that a(8) = 14 is the smallest possibility. And so on.
It seems that once a(5) and a(6) are chosen, one may always take the smallest possible choice for the next term without ever again running into difficulty. This is in strong contrast to the (exceptional) case of the variant where we require 10 prime sums among 7 consecutive terms, cf. sequence A329574.
-
{A329573(n,show=0,o=1,N=12,M=6,D=[5,9,6,10],p=[],u=o,U)=for(n=o+1,n, show>0&& print1(o","); show<0&& listput(L,o); U+=1<<(o-u); U>>=-u+u+=valuation(U+1,2); p=concat(if(#p>=M,p[^1],p),o); D&& D[1]==n&& [o=D[2],D=D[3..-1]]&& next; my(c=N-sum(i=2,#p, sum(j=1,i-1, isprime(p[i]+p[j])))); for(k=u,oo,bittest(U,k-u)|| min(c-#[0|p<-p,isprime(p+k)],#p>=M)|| [o=k,break]));show&&print([u]);o} \\ optional args: show=1: print a(o..n-1), show=-1: append them on global list L, in both cases print [least unused number] at the end. See the wiki page for more.
A329564
For all n >= 0, exactly five sums are prime among a(n+i) + a(n+j), 0 <= i < j < 5; lexicographically earliest such sequence of distinct nonnegative numbers.
Original entry on oeis.org
0, 1, 2, 3, 6, 5, 8, 11, 7, 12, 29, 18, 19, 4, 13, 9, 22, 10, 21, 14, 57, 16, 15, 17, 26, 27, 20, 23, 33, 34, 38, 45, 25, 28, 51, 46, 31, 43, 58, 30, 24, 37, 49, 35, 36, 102, 47, 42, 55, 32, 41, 48, 65, 39, 62, 44, 40, 63, 69, 50, 68, 59, 80, 71, 54, 77, 60, 53, 56, 74, 75
Offset: 0
For n = 0, we consider pairwise sums among the first 5 terms a(0..4), among which we must have 5 primes. To get a(4), consider first a(0..3) = (0, 1, 2, 3) and the pairwise sums (a(i) + a(j), 0 <= i < j <= 3) = (1; 2, 3; 3, 4, 5) among which there are 4 primes, counted with multiplicity (i.e., the prime 3 is there two times). So the additional term a(4) must give exactly one more prime sum with all of a(0..3). We find that 4 or 5 would give two more primes, but a(4) = 6 gives exactly one more, 1 + 6 = 7.
Now, for n = 1 we forget the initial 0 and consider the pairwise sums of the remaining terms {1, 2, 3, 6}. There are 3 prime sums, so the next term must give two more. The term 4 would give two more (1+4 and 3+4) primes, but thereafter we would have {2, 3, 6, 4} with only 2 prime sums and impossibility to add one term to get three more prime sums: 2+x, 6+x and 4+x can't be all prime for x > 1.
Therefore 4 isn't the next term, and we try a(5) = 5 which indeed gives the required number of primes, and also allows us to continue.
Cf.
A329425 (6 primes using 5 consecutive terms).
Cf.
A055266 &
A253074 (0 primes using 2 terms),
A329405 &
A329450 (0 primes using 3 terms),
A055265 &
A128280 (1 prime using 2 terms),
A329333,
A329406 -
A329410 (1 prime using 3, ..., 10 terms),
A329411 -
A329416 and
A329452,
A329453 (2 primes using 3, ..., 10 terms),
A329454 &
A329455 (3 primes using 4 resp. 5 terms),
A329449 &
A329456 (4 primes using 4 resp. 5 terms),
A329568 &
A329569 (9 primes using 6 terms),
A329572 &
A329573 (12 primes using 7 terms),
A329563 -
A329581: other variants.
-
{A329564(n,show=1,o=0,N=5,M=4,X=[[4,4]],p=[],u,U)=for(n=o,n-1, show>0&& print1(o","); show<0&& listput(L,o); U+=1<<(o-u); U>>=-u+u+=valuation(U+1,2); p=concat(if(#p>=M,p[^1],p),o); my(c=N-sum(i=2,#p, sum(j=1,i-1, isprime(p[i]+p[j])))); if(#p
A329576
For all n >= 1, exactly seven sums are prime among a(n+i) + a(n+j), 0 <= i < j < 6; lexicographically earliest such sequence of distinct positive numbers.
Original entry on oeis.org
1, 2, 3, 4, 5, 8, 11, 26, 15, 9, 14, 32, 17, 20, 21, 27, 10, 16, 19, 7, 12, 13, 24, 6, 23, 35, 25, 37, 18, 36, 22, 31, 61, 28, 30, 39, 40, 43, 33, 64, 38, 45, 34, 29, 63, 50, 44, 53, 42, 59, 47, 54, 48, 41, 90, 49, 55, 52, 108, 58, 46, 51, 121, 73, 78, 76, 100, 79, 81, 151, 60, 67, 112, 70, 69
Offset: 1
For n = 1, we must forbid the greedy choice for a(6) which would be 6, which leads to a dead end: there is no possibility to find a subsequent term that would give 7 prime sums together with {2, 3, 4, 5, 6}. If we take the next larger possibility, a(6) = 8, then it works for the next and all subsequent terms.
Cf.
A329425 (6 primes using 5 consecutive terms),
A329566 (6 primes using 6 consecutive terms).
Cf.
A329449 (4 primes using 4 consecutive terms),
A329456 (4 primes using 5 consecutive terms).
Cf.
A329454 (3 primes using 4 consecutive terms),
A329455 (3 primes using 5 consecutive terms).
Cf.
A329411 (2 primes using 3 consecutive terms),
A329452 (2 primes using 4 consecutive terms),
A329453 (2 primes using 5 consecutive terms).
-
{A329576(n,show=1,o=1,N=7,M=5,X=[[6,6]],p=[],u=o,U)=for(n=o+1,n, show>0&& print1(o","); show<0&& listput(L,o); U+=1<<(o-u); U>>=-u+u+=valuation(U+1,2); p=concat(if(#p>=M,p[^1],p),o); my(c=N-sum(i=2,#p, sum(j=1,i-1, isprime(p[i]+p[j])))); if(#p
Showing 1-10 of 12 results.
Comments