%I #13 May 10 2024 08:50:00
%S 1,2,1,3,2,1,5,2,3,1,7,2,1,11,2,3,5,4,3,7,2,13,1,17,2,15,1,19,2,1,23,
%T 2,9,1,5,7,6,1,29,2,5,11,3,13,4,11,7,1,31,2,5,13,6,1,37,2,7,3,17,4,15,
%U 1,41,2,1,43,2,33,1,47,2,35,3,19,4,3,23,4,5
%N a(n) = GCD(A280864(n), A280864(n+1)).
%H Rémy Sigrist, <a href="/A370628/a370628.gp.txt">PARI program</a>
%F GCD(a(n), a(n+1)) = 1.
%F A007947(a(n)) = A280738(n).
%F a(n) = 1 iff n belongs to A280774.
%e a(89) = GCD(A280864(89), A280864(90)) = GCD(90, 135) = 45.
%t terms = 100;
%t rad[n_] := Times @@ FactorInteger[n][[All, 1]];
%t A280864 = Reap[present = 0; p = 1; pp = 1;
%t Do[forbidden = GCD[p, pp]; mandatory = p/forbidden;
%t a = mandatory;
%t While[BitGet[present, a] > 0 || GCD[forbidden, a] > 1,
%t a += mandatory];
%t Sow[a];
%t present += 2^a; pp = p; p = rad[a],
%t {terms}]][[2, 1]];
%t Clear[a];
%t a[n_] := GCD[A280864[[n]], A280864[[n + 1]]];
%t Table[a[n], {n, 1, terms - 1}] (* _Jean-François Alcover_, May 10 2024, adapted from _Rémy Sigrist_'s PARI program *)
%o (PARI) \\ See Links section.
%Y Cf. A007947, A280738, A280774, A280864.
%K nonn
%O 1,2
%A _Rémy Sigrist_, May 01 2024