%I #21 Oct 07 2023 21:36:41
%S 2,6,10,14,18,9,36,12,4,8,16,20,24,15,30,21,42,28,32,40,44,22,26,52,
%T 48,39,78,60,35,70,45,54,27,108,72,56,64,68,34,38,76,80,55,110,66,33,
%U 132,84,63,90,65,130,104,88,92,46,50,25,150,75,100,96,87,174,58,62,124,116,112,98,49,392
%N Starting with a(1) = 2, the lexicographically earliest infinite sequence of distinct positive integers such that |a(n) - a(n-1)| is a divisor of a(n)*a(n-1), where |a(n) - a(n-1)| is not a prime and greater than 1.
%C For the sequence to be infinite no term can be a prime except for a(1) = 2. One can show that if a(n) is a prime p, then the only possible value for a(n-1) is 2p or p + p^2 since, if a term is prime, the preceding term must be a multiple of that prime. However the preceding term cannot be 2p since the difference between the terms would then be prime, therefore it must be p + p^2. However the only possible value for the term after a prime p is likewise p + p^2, but that has already been used, thus allowing a term to be prime would terminate the sequence.
%H Scott R. Shannon, <a href="/A366047/b366047.txt">Table of n, a(n) for n = 1..10000</a>.
%H Michael De Vlieger, <a href="/A366047/a366047.png">Log log scatterplot of a(n)</a> n = 1..1024, showing prime powers in gold, squarefree numbers in green, and numbers neither squarefree nor prime powers in blue, highlighting powerful numbers that are not prime powers in light blue.
%e a(9) = 4 as |4 - a(8)| = |4 - 12| = 8, and 8 is a divisor of 4*12 = 48 and is not a prime. Note that |3 - 12| = 9 is a divisor of 3*12 = 36 and is not a prime, but as shown above a prime term will terminate the sequence so is not permitted.
%t nn = 120; c[_] := False; s = {2, 6};
%t f[x_] := Times @@ FactorInteger[x][[All, 1]];
%t MapIndexed[Set[{a[First[#2]], c[#1]}, {#1, True}] &, s];
%t Set[{j, u}, {s[[-1]], 4}];
%t Do[k = u;
%t While[Or[c[k], #1 < 4, PrimeQ[#1],
%t ! Divisible[j*k, #1], ! Divisible[j, #2], ! Divisible[k, #2]] & @@
%t {#, f[#]} &@ Abs[j - k], k++];
%t Set[{a[n], c[k], j}, {k, True, k}];
%t If[k == u, While[Or[c[u], PrimeQ[u]], u++]], {n, Length[s] + 1, nn}];
%t Array[a, nn] (* _Michael De Vlieger_, Sep 29 2023 *)
%Y Cf. A027750, A366111, A365984, A363576, A359799.
%K nonn
%O 1,1
%A _Scott R. Shannon_, Sep 27 2023
|