login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A161164
Absolute value of (the n-th twin prime minus the n-th non-twin prime).
0
1, 18, 30, 36, 40, 50, 60, 54, 58, 56, 70, 68, 70, 86, 90, 66, 70, 104, 114, 96, 112, 108, 112, 98, 112, 116, 124, 134, 138, 126, 130, 128, 132, 110, 112, 108, 114, 90, 96, 92, 94, 30, 36, 36, 46, 26, 28, 22, 20, 60, 30, 52, 44, 54, 42, 54, 50, 52, 48, 178, 164, 168, 150
OFFSET
1,2
FORMULA
a(n) = |A001097(n)-A007510(n)|.
EXAMPLE
a(1)=abs(3-2)=1, a(2)=abs(5-23)=abs(-18)=18, a(3)=abs(7-37)=abs(-30)=30, etc.
PROG
(Python)
from sympy import nextprime
def aupton(terms):
n, p, q = 1, 2, 3
alst, non_twins, twins = [], [2], [3]
while True:
p, q = q, nextprime(q)
if q - p == 2:
if p != twins[-1]: twins.append(p)
twins.append(q)
else:
if p != twins[-1]: non_twins.append(p)
if len(twins) >= n and len(non_twins) >= n:
alst.append(abs(twins[n-1] - non_twins[n-1]))
if n == terms: break
n += 1
return alst
print(aupton(63)) # Michael S. Branicky, Jan 18 2021
CROSSREFS
Sequence in context: A032612 A151700 A154367 * A092357 A141115 A377185
KEYWORD
nonn
AUTHOR
EXTENSIONS
Terms corrected by D. S. McNeil, Dec 10 2009
STATUS
approved