W = 63

H = vector(2*W+1)
V = vector(2*W+1)

{
	z=(W+1)*(1+I);
	n=0;
	dir=1;
	for (w=1, oo,
		for (t=1, 2,
			for (i=1, w,
				n++;

				v = valuation(1+bitor(V[real(z)], H[imag(z)]), 2);
				V[real(z)] += 2^v;
				H[imag(z)] += 2^v;

				print (n-1 " " v);

				if (n==(2*W+1)^2,
					break (3);
				);

				z+=dir;
			);
			dir*=I;
		);
	);
}

quit