reviewed
approved
reviewed
approved
proposed
reviewed
editing
proposed
terms = 80; rad[n_] := Times @@ FactorInteger[n][[All, 1]];
A280864 = Reap[present = 0; p = 1; pp = 1; Do[forbidden = GCD[p, pp]; mandatory = p/forbidden; a = mandatory; While[BitGet[present, a] > 0 || GCD[forbidden, a] > 1, a += mandatory]; Sow[a]; present += 2^a; pp = p; p = rad[a], terms]][[2, 1]];
Clear[a];
a[1] = 1;
a[n_] := a[n] = For[b = 2a[n-1], True, b += 2, If[FreeQ[A280864[[1 ;; n-1]], b], Return[b/2]]];
Array[a, terms] (* Jean-François Alcover, Nov 23 2017, after Rémy Sigrist program for A280864 *)
approved
editing
editing
approved
editing
proposed
editing
approved