Silverman

silv.ml

let silv n =
	let u = make_vect (2*n+1) 0 in
	let k = ref 2 in let j = ref 3 in let p = ref 3 in
		u.(1)<- 1; u.(2)<-2; u.(3)<-2;
		while u.(n)==0 do
			for i=1 to !k do u.(!p+i) <- !j; done;
			p:=!p+(!k); j:=!j+1; k:=u.(!j);
		done;
	u;;