|
MATHEMATICA
|
(* computes a(n), 1 <= n <= 5 *)
sfRun[n_]:=Module[{i=n}, While[SquareFreeQ[i], i++]; i-n]
nsfRun[n_]:=Module[{i=n}, While[!SquareFreeQ[i], i++]; i-n]
a270344[{low_, high_}]:=Module[{i, next, r, s, list=Table[{}, 5]}, For[i=low, i<=high, i+=next, r=nsfRun[i]; If[r==0, next=sfRun[i], s=nsfRun[i+r+1]; If[r==s, If[list[[r]]=={}, list[[r]]={i, i+r, i+2r}]; next=2r+2, next=r+1]]]; list]
Map[#[[2]]&, a270344[{0, 10000000}]]]] (* data *)
|