%I #11 Jan 25 2024 07:56:44
%S 1,2,4,12,3,9,45,5,25,50,8,56,7,49,147,21,42,6,30,10,70,14,154,11,121,
%T 242,16,208,13,169,507,27,297,33,66,18,234,24,408,17,289,578,32,608,
%U 19,361,1083,57,114,36,828,23,529,1058,46,230,20,220,22,286,26,130
%N a(n) = n for n <= 2; for n > 2, rad(a(m)) != rad(a(n)), m < n, implies a(n+1) is the least novel k such that rad(k) | a(n), rad(a(m)) = rad(a(n)), m < n, implies a(n+1) = q*a(n), where q is the least novel prime such that gcd(q, a(n)) = 1 and rad(a(m)) != rad(q*a(n)), m < n.
%C Let rad(n) = A007947(n).
%C In other words, starting with 1 and 2, thereafter, if a(n) is such that rad(a(m)) != rad(a(n)), m < n, then a(n+1) is the smallest k != a(m) such that rad(k) | a(n).
%C If there is a term a(m) such that rad(a(m)) = rad(a(n)), then a(n+1) is the smallest k != a(m) obtainable from a(n) by multiplying a(n) by a nondivisor prime q, and the result must be such that there is no prior term a(m) with rad(a(m)) = rad(a(n)).
%C Primes appear in order as a consequence of first condition.
%C If a(n) = p, then a(n+1) = p^2 for prime p.
%C The second condition implies a(m) != q*a(n), m < n.
%C Conjectures:
%C 1. Sequence is a permutation of natural numbers.
%C 2. Primorials outside of {1, 2, 30} enter late. a(18) = 6, a(19) = 30, a(20476) = 210, a(20477) = 2310. Thereafter many long-missing k such that 30 | k and 210 | k enter the sequence.
%H Michael De Vlieger, <a href="/A365000/b365000.txt">Table of n, a(n) for n = 1..10000</a>
%H Michael De Vlieger, <a href="/A365000/a365000.png">Log log scatterplot of a(n)</a>, n = 1..100000, with a color code showing primes in red, other prime powers in gold, squarefree composites in green, and numbers neither squarefree nor prime powers in blue, accentuating terms that are powerful but not prime powers appear in magenta.
%H Michael De Vlieger, <a href="/A365000/a365000_1.png">Log log scatterplot of a(n)</a>, n = 1..1024, with the same color function as above, except accentuating powerful numbers that are not prime powers in light blue.
%e a(3) = 4 since rad(a(2)) = rad(2) = 2, and rad(1) = 1, and the smallest missing number k such that rad(k) | 2 is 4.
%e a(4) = 12 since rad(a(3)) = rad(4) = 2, rad(2) = 2, and thus since q = 3 is the smallest prime that does not divide a(3) = 4, a(4) = q * a(3) = 3 * 4 = 12.
%e a(5) = 3 since rad(a(4)) = rad(12) = 6, and no a(m), m < 4 has rad(a(m)) = 6. Thus we select the smallest number missing number k such that rad(k) | 12, so k = 3.
%e a(6) = 9 since a(5) is squarefree and thus rad(a(5)) = 3; no a(m), m < 5, has rad(a(m)) = 3. We select the smallest missing number k such that rad(k) | 3, and that is 9.
%e a(7) = 45 since a number a(m), m < 6, with rad(a(m)) = 3 has appeared. Therefore we find the smallest prime q that does not divide 9, but also is such that no a(m), m < n, rad(a(m)) | q*a(6) appears in the sequence. Since numbers a(m) with rad(a(m)) = 6 already appear, we are forced to q = 5, hence a(7) = 5 * a(6) = 5 * 9 = 45, etc.
%t nn = 120; c[_] := False; r[_] := False;
%t rad[x_] := rad[x] = Times @@ FactorInteger[x][[All, 1]];
%t Array[Set[{a[#], c[#]}, {#, True}] &, 2];
%t j = 2; r[1] = True; u = 3;
%t Do[
%t If[r[#],
%t k = 2;
%t While[Or[! CoprimeQ[#, k], r[rad[k j]]], k = NextPrime[k]];
%t k *= j,
%t If[PrimeQ[#],
%t k = Floor@ Log[#, u]; While[c[#^k], k++]; k = #^k,
%t k = u; While[Or[c[k], ! Divisible[#, rad[k]]], k++]]] &[
%t rad[j]];
%t Set[{a[n], c[k], r[rad[j]], j}, {k, True, True, k}];
%t If[k == u, While[c[u], u++]], {n, 3, nn}];
%t Array[a, nn]
%Y Cf. A006881, A007947.
%K nonn
%O 1,2
%A _Michael De Vlieger_ and _David James Sycamore_, Sep 12 2023