factor(Prod((n*n+1),n=1..17))