Edit: As Nate pointed out, the hint was wrong. I added the correct answer (from Fulton & Harris) below.
The hint I get: We have a decomposition of regular representation $R$=$U \oplus V$ where $U$ is the trivial representation. Since we have the expansion
$$ \bigwedge^k(U \oplus V) \cong \bigoplus_{a+b=k}\bigwedge^a U \otimes \bigwedge^b V $$
It suffices to compute $(\chi_{\wedge^k R},\chi_{\wedge^k R})$, and it should be $2$. Then we can deduce that $\bigwedge^k V$ is irreducible whenever $0 \leq k \leq n-1$.
I managed to compute $(\chi_{\wedge^2 R},\chi_{\wedge^2 R})$ when $n=3$. And it is indeed $2$. But I don't know how to step on. Is there any way to expand $(\chi_{\wedge^k R},\chi_{\wedge^k R})$ with respect to $U$ and $V$ (some massive computation involved I guess?). I know how to compute $\chi_{\wedge^2 R}(g)$ but how to compute $\chi_{\wedge^k R}(g)$? My guess is there is some much simpler approach for regular representation.
Any hint or solution will be sincerely appreciated!
In case people are confused with terminologies:
The inner product $(χ_{∧^kR},χ_{∧^kR})$ can be interpreted as the expectation of $(\text{Tr} \phi(g))^2$, where $\phi$ is the representation of $S_n$ on the exterior product space ${∧^k\mathbb{C}^n}$, and $g$ is a random element of $S_n$ chosen uniformly.
If we define random variables $r_a$ ($a$ runs over the $k$-element subsets of $[1, 2, ... ,n]$) as
$$r_a=\left\{ \begin{array} & 1 & \text{if the action of $g$ maps the set $a$ to $a$ and preserves the orientation} \\-1 & \text{if the action of $g$ maps the set $a$ to $a$ and reverses the orientation} \\ 0 & \text{otherwise}\end{array} \right .$$
Then $\text{Tr} \phi (g)$ is just the sum of all the $r_a$s, by the definition of exterior product space.
To compute the expectation of $(\text{Tr} \phi(g))^2$, we expand it as $$\mathbb E(\text{Tr} \phi(g) )^2= \sum_{a\in {n\choose k}} \mathbb Er_a^2+\sum_{a, b\in {n\choose k}\\{|a\cap b|=k-1} } \mathbb Er_ar_b+\sum_{a, b\in {n\choose k}\\{|a\cap b| \leq k-2} } \mathbb Er_ar_b$$
where $n \choose k $ denotes all the $k$-element subsets of $[1,2,...,n]$, and the last two sums are over ordered pairs.
It is straightforward that $\mathbb Er_a^2=\frac{k!(n-k)!}{n!}$, so the first sum evaluates to $1$.
To study the second sum, we consider the action of $g$ on $a \cap b$. For the summand to be nonzero, $g$ must map this set to itself and map the elements $a \text{\\} (a \cap b) $ and $b \text{\\} (a \cap b) $ to themselves respectively.
This happens with probability $\frac{(n-k-1)!(k-1)!}{n!}$; and whenever this happens, the value of $r_ar_b$ is $1$, whether the orientation of $a \cap b$ is reversed by $g$ or not. By counting the number of such ordered pairs $(a,b)$, i.e. $\frac{n!}{(n-k-1)!(k-1)!}$, we find out that the second term also evaluates to $1$.
Every summand in the third term equals $0$. To see this, we may assume that $b_1$ and $b_2$ are two elements in $B\text{\\}A$, and let $\tau$ be the element in $S_n$ that exchanges $b_1$ and $b_2$ and fix everything else.
The expectation $\mathbb E r_ar_b$ can be written as $\frac{1}{n!}\sum _{g\in S_n} r_a(g)r_b(g)$, just to emphasize that $r_a$ and $r_b$ depends on $g$.
But $\sum _{g\in S_n} r_a(g)r_b(g)=\sum _{g\in S_n} r_a(\tau g)r_b(\tau g)=-\sum _{g\in S_n} r_a(g)r_b(g)$, as $\tau$ reverses the orientation on $b$ but has no effect on $a$ if either of the summands are nonzero.
So the third sum evaluates to $0$ and we are done.