OFFSET
1,2
COMMENTS
Unlike A356850 all the terms are concentrated along three straight lines. In the first 100000 terms there are ten fixed points, 1, 2, 3, ..., 27, 57, and it is likely no more exist. The sequence is conjectured to be a permutation of the positive integers.
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..16384
Michael De Vlieger, Log log scatterplot of a(n) n = 1..2^10, with records in red, local minima in blue, highlighting primes in green and other prime powers in gold, showing numbers that are neither squarefree nor prime powers in magenta.
Michael De Vlieger, Scalar scatterplot of a(n), n = 1..2^12, showing primes in red, composite prime powers in gold, composite squarefree numbers in green, products of composite prime powers in magenta, and other numbers in blue.
Michael De Vlieger, Log-log scatterplot of a(n), n = 1..2^10 using a color function showing Omega(a(n-1)).
EXAMPLE
MATHEMATICA
nn = 90; c[_] = False; Array[Set[{a[#], c[#]}, {#, True}] &, 2]; u = 3; j = 2; Do[k = u; m = Times @@ Array[a[n - #] &, PrimeOmega[j]]; While[Nand[! c[k], CoprimeQ[k, m]], k++]; Set[{a[n], c[k], j}, {k, True, k}]; If[k == u, While[c[u], u++]], {n, 3, nn}]; Array[a, nn] (* Michael De Vlieger, Nov 20 2022 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Scott R. Shannon, Nov 08 2022
STATUS
approved