For each k=1,2,…,n, denote bk=ak2a1a2…an and consider the following polynomial
P(x)=(x−b1)(x−b2)…(x−bn)=xn+cn−1xn−1+⋯+c1x+c0.
Based on Vieta's theorem, one can see that cn−1 is an integer, and for any k integers i1,i2,…,ik∈{1,2,…,n}, bi1bi2…bik is also an integer, which implies that all coefficients ck are integers. On the other hand, P(x) is a monic polynomial which implies that all of its roots b1,b2,…,bn are also integers. Hence,
a12a1a2…an,a22a1a2…an,…,an2a1a2…an∈Z.