ludic numbers

ludic()={
a=[1];
b=[2..100000];
until(#b==0,
       a=concat(a,b[1]);
       c=ceil(#b/b[1]);
       b=setminus(b,vecextract(b,vector(c,i,b[1]*(i-1)+1)))
      );
a
}