cp's OEIS Frontend

A354075 Lexicographically earliest infinite sequence of distinct positive integers such that A(a(n+1)) is prime to A(a(n)) but not to A(a(n-1)), where A is A001414.

%I A354075 #45 Mar 22 2025 19:03:10
%S A354075 2,3,4,14,15,20,16,24,18,26,33,5,7,6,10,21,12,8,94,9,124,27,25,38,30,
%T A354075 62,32,11,35,28,36,40,39,45,42,48,44,54,46,57,86,49,74,51,13,50,22,55,
%U A354075 56,60,63,64,75,65,80,66,90,70,96,68,69,92,84,105,85,112,87
%N A354075 Lexicographically earliest infinite sequence of distinct positive integers such that A(a(n+1)) is prime to A(a(n)) but not to A(a(n-1)), where A is A001414.
%C A354075 2,3,4 is the earliest string of three consecutive numbers which satisfy the definition, therefore the sequence begins a(1)=2, a(2)=3, a(3)=4.
%C A354075 Sequence is infinite since there always exists a k which has not occurred before such that A(k) is prime to A(a(n)) but not to A(a(n-1)). Since A001414 covers N/{1} a number m can be found such that A(m)=k. Thus k can be chosen for a(n+1) unless there is a smaller number with the same property.
%C A354075 Similar to the Yellowstone sequence (A098550) in terms of coprime relations.
%C A354075 The first seven primes are in natural order but then we have ...,17,23,19,31,37,43,41,47,29,...
%C A354075 Conjectured to be a permutation of N/{0,1}.
%H A354075 Jinyuan Wang, <a href="/A354075/b354075.txt">Table of n, a(n) for n = 1..10000</a>
%e A354075 a(4)=14 because A(14)=9 is prime to A(a(3))=4 but not to A(a(2))=3, and is the smallest number not already seen in the sequence which has this property.
%o A354075 (PARI) f(n) = my(f=factor(n)); f[, 1]~*f[, 2]; \\ A001414
%o A354075 lista(nn) = {my(va = vector(nn)); va[1] = 2; va[2] = 3; for (n=3, nn, my(k=1); while ((gcd(f(va[n-1]), f(k)) != 1) || (gcd(f(va[n-2]), f(k)) == 1) || #select(x->(x==k), va), k++); va[n] = k;); va;} \\ _Michel Marcus_, Jun 12 2022
%Y A354075 Cf. A001414, A056240, A098440, A277272.
%K A354075 nonn
%O A354075 1,1
%A A354075 _David James Sycamore_, Jun 11 2022
%E A354075 Corrected and extended by _Michel Marcus_, Jun 12 2022