function [b]=flavius(s,k,n)
A=ones([n,1]);
i=k-1;
c=1;
m=0;
while norm(A,1)>s
    while (m<i | m==i)
       if A(c)==1
           m=m+1;
       end
       c=c+1;
       if c==n+1
           c=1;
       end
    end
    if c==1
        A(n)=0;
    end;
    if c>1
        A(c-1)=0;
    end
    m=0;
end
b=[];
for j=1:n
    if A(j)==1
        b=[b;j];
    end
end
endfunction
