Putting $x_k=a_k^{\frac{1}{n}}$, one can reformulate the inequality to
be shown as
$$
P_n : (1+x_1^n)(1+x_2^n)\ldots (1+x_n^n) \geq (1+x_1x_2\ldots x_n)^n
\ (\text{for any nonnegative} \ x_1,x_2,\ldots ,x_n)
$$
Let us show $P_n$ by induction on $n$. We have that $P_1$ is an equality,
$P_2$ is easily equivalent to $(x_1-x_2)^2 \geq 0$. Suppose now that $n\geq 3$ and that $P_{n-1}$ holds. Let $x_1,x_2, \ldots x_n$ be nonnegative
numbers. We want to show that
$$
(1+x_1^n)(1+x_2^n)\ldots (1+x_n^n) \geq (1+x_1x_2\ldots x_n)^n \tag{1}
$$
If we put $b=(x_2x_3\ldots x_n)^{\frac{n}{n-1}}$, then applying $P_{n-1}$ to the
$n-1$ numbers $x_2,x_3, \ldots x_n$, we see that it suffices to show (1) when
all the $x_i$ are equal to $b$. In other words, we may assume without loss
of generality that $x_2=x_3=\ldots =x_n$, in which case (1) becomes
$$
(1+x_1^n)(1+x_2^n)^{n-1} \geq (1+x_1x_2^{n-1})^n \tag{2}
$$
If we put $y=x_2^n$ and $z=\frac{x_1}{x_2}$, (2) can be rewritten as
$$
(1+yz^n)(1+y)^{n-1} \geq (1+yz)^n \tag{3}
$$
So if we put $\delta=(1+yz^n)(1+y)^{n-1} - (1+yz)^n$, we must show that
$\delta$ is nonnegative. But expanding the Newton binomials leads to
$$
\delta=\sum_{k=1}^{n-1}n\binom{n-1}{k-1}y^k \bigg(\frac{1-z^k}{k}-
\frac{1-z^n}{n}\bigg) \tag{4}
$$
So all we need to show now is that the sequence $(\frac{1-z^n}{n})$ is decreasing in
$n \geq 1$, for any $z>0$. The algebraic equality
$$
\frac{1-z^{n}}{n}-\frac{1-z^{n+1}}{n+1}=
\frac{(z-1)^2}{n(n+1)}\bigg(
\sum_{k=0}^{n-1} (k+1)z^k
\bigg) \tag{5}
$$
finishes the proof.