p = Prime[Range[PrimePi[10^3]]]; SequenceCases[p, ({a_, b_, c_}) /; (a + c < 2 b) :> b, 36, Overlaps -> True] SequenceCases[p, ({a_, b_, c_}) /; (a + c > 2 b) :> b, 37, Overlaps -> True] p = Prime[Range[PrimePi[10^6] + 1]]; Length[Select[Partition[p, 3, 1], #[[3]] + #[[1]] < 2 #[[2]] &]] Length[Select[Partition[p, 3, 1], #[[3]] + #[[1]] > 2 #[[2]] &]] p = Prime[Range[PrimePi[10^7] + 1]]; Length[Select[Partition[p, 3, 1], #[[3]] + #[[1]] < 2 #[[2]] &]] Length[Select[Partition[p, 3, 1], #[[3]] + #[[1]] > 2 #[[2]] &]]