\\ A067138	OR-numbral multiplication table.
U(n,k) = {
	my (v=0);
	while (n,
		if (n%2,
			v = bitor(v,k);
		);
		n \= 2;
		k *= 2;
	);
	return (v);
}

T(n,k) = {
	my (fn=factor(n), fk=factor(k),
		pin=apply(primepi, fn[,1]~),
		pik=apply(primepi, fk[,1]~),
		v=1);
	for (i=1, #fn~,
		for (j=1, #fk~,
			my (p=prime(pin[i]+pik[j]-1), x=valuation(v,p), y=U(fn[i,2], fk[j,2]));
			v *= p^(bitor(x,y)-x);
		);
	);
	return (v);
}

for (d=1, 200, for (k=1, d, print (i++ " " T(d+1-k, k))))

quit