OFFSET
1,1
LINKS
N. J. A. Sloane, Table of n, a(n) for n = 1..641 (Based on John P. Linderman's computation of a million terms of A084937)
PROG
(AWK, using a-file for A084937)
awk ' BEGIN {s = 0}
NR%3 == 0 { s = $2 }
NR%3 == 1 { t = $2
if (s > t)
print( (NR-1)/3 )
} ' a084937.txt | awk '{ print NR, $1}' >b249689.txt
(Python)
from __future__ import division
from fractions import gcd
A249689_list, l1, l2, s, b = [], 2, 1, 3, set()
for n in range(3, 10**4):
....i = s
....while True:
........if not i in b and gcd(i, l1) == 1 and gcd(i, l2) == 1:
............l2, l1 = l1, i
............b.add(i)
............if l2 > l1 and n % 3 == 1:
................A249689_list.append((n-1)//3)
............while s in b:
................b.remove(s)
................s += 1
............break
........i += 1 # Chai Wah Wu, Dec 12 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Nov 12 2014
STATUS
approved