%I #15 Sep 17 2023 12:08:25
%S 1,1,1,2,3,1,4,1,1,5,6,2,7,4,3,7,8,1,9,5,4,10,11,1,3,8,1,12,13,5,14,
%T 12,6,9,15,2,16,16,7,9,17,4,18,19,3,20,21,7,22,3,8,10,23,1,5,14,9,24,
%U 25,5,26,27,4,28,29,10,30,31,11,10,32,1,33,21,3,34,35,8,36,15,1,37,38,12,37,38,13,39,40,5,1
%N Lexicographically earliest infinite sequence such that a(i) = a(j) => A365465(i) = A365465(j) for all i, j >= 1, where A365465(n) = A356867(n) / gcd(n, A356867(n)), and A356867 is Sycamore's Doudna variant D(3).
%C Restricted growth sequence transform of A365465.
%C Compare to the scatter plots of A365431 (analogous sequence for Doudna variant D(2)), and also of A365393 and A365718.
%H Antti Karttunen, <a href="/A365715/b365715.txt">Table of n, a(n) for n = 1..59049</a>
%o (PARI)
%o \\ Needs also program from A356867:
%o up_to = 59049; \\ = 3^10.
%o rgs_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(om,invec[i]), my(pp = mapget(om, invec[i])); outvec[i] = outvec[pp] , mapput(om,invec[i],i); outvec[i] = u; u++ )); outvec; };
%o A365465(n) = (A356867(n)/gcd(n, A356867(n)));
%o v365715 = rgs_transform(vector(up_to,n,A365465(n)));
%o A365715(n) = v365715[n];
%Y Cf. A356867, A365465.
%Y Cf. also A365431, A365393, A365718.
%K nonn,look
%O 1,4
%A _Antti Karttunen_, Sep 17 2023