OFFSET
1,1
COMMENTS
MAPLE
M:=[0, 2]: isElement := proc(n) global M; local b, j, k; b:=true; for j from 1 to nops(M)-1 do for k from j+1 to nops(M) do if M[j] = n mod M[k] then b:=false; break; fi od od; return b end: for z to 1 do for n from 3 to 1000 do if isElement(n) then M:=[op(M), n] fi od od; M; map(proc(k) M[k]-M[k-1] end, [$2..nops(M)]);
MATHEMATICA
a[1] = 0; a[2] = 2; a[n_] := a[n] = Block[{k = a[n - 1] + 1, t = a[ # ] & /@ Range[n - 1]}, While[ Intersection[t, Union[ Mod[k, Rest[ t]]]] != {}, k++ ]; k]; t = Table[ a[n], {n, 71}]; Drop[t, 1] - Drop[t, -1] (* Robert G. Wilson v *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Walter Kehowski, Oct 14 2005
EXTENSIONS
More terms from Robert G. Wilson v, Oct 19 2005
STATUS
approved