combine(x,y) =					\
	my (f, pk);					\
	if (y>1,					\
		f = factor(y);			\
		for (i=1, #f[,1],		\
			pk = f[i,1]^f[i,2];	\
			if (x % pk==0,		\
				x = x/pk,		\
				x = x*pk		\
			)					\
		)						\
	);							\
	return (x)

a = 1
for (n=1, 3365, a=combine(a,n); print(n " " a))

quit