login
A387519
Indices of record high points in A386487.
4
1, 3, 4, 6, 7, 12, 19, 32, 58, 77, 145, 208, 237, 321, 369, 480, 838, 1290, 1533, 3031, 4003, 4965, 6900, 7653, 13365, 18980, 23022, 28225, 36034, 54189, 74152, 82124, 129792, 274543, 419570, 730075, 1679322, 1844372, 2469066, 3367329, 5991137, 7002436, 7690250, 18043250, 29261251
OFFSET
1,2
LINKS
MATHEMATICA
r = -1; s = Import["https://oeis.org/A386482/b386482.txt", "Data"][[All, -1]]; Reap[Do[If[# > r, r = #; Sow[n]] &[s[[n]] - n], {n, Length[s]}] ][[-1, 1]] (* Michael De Vlieger, Sep 01 2025 *)
PROG
(Python)
from math import gcd
from itertools import count, islice
def A387519_gen(): # generator of terms
m, mset, c, a = 2, {2}, 3, 0
yield 1
for j in count(3):
for i in range(m-1, c-1, -1):
if i not in mset and gcd(i, m)>1:
if i-j>a:
yield j
a = i-j
break
else:
for i in count(max(c, m+1)):
if i not in mset and gcd(i, m)>1:
if i-j>a:
yield j
a = i-j
break
m = i
mset.add(i)
while c in mset:
mset.remove(c)
c += 1
A387519_list = list(islice(A387519_gen(), 20)) # Chai Wah Wu, Sep 10 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Sep 01 2025
EXTENSIONS
More terms from Michael De Vlieger, Sep 01 2025
a(40)-a(45) from Chai Wah Wu, Sep 11 2025
STATUS
approved