Let $\gamma>0$. I would like a nice way to prove that $$\int_{\begin{array}{c} 0\leq s,t\leq1\\ s+t\leq1 \end{array}}t^{\gamma-1}s^{\gamma-1}(1-t-s)^{\gamma-1}dtds=\frac{\Gamma(\gamma)^3}{\Gamma(3\gamma)}.$$
Edit/Remark: zhoraster's answer along with induction allows one to prove the generalized identity when integrating over the simplex $$\int_{\begin{array}{c} 0\leq t_{i}\leq1\\ t_{1}+\cdots+t_{k}\leq1 \end{array}}t_{1}^{x_{1}-1}\cdots t_{k}^{x_{k}-1}\left(1-\sum_{i=1}^{k}t_{i}\right)^{x_{k+1}-1}dt_{1}\cdots dt_{k}=\frac{\Gamma(x_{1})\Gamma(x_{2})\cdots\Gamma(x_{k+1})}{\Gamma(x_{1}+\cdots+x_{k+1})} $$ for any $0<x_i<1$
This integral came about while I was reading through Balog and Friedlander's paper A Hybrid Theorem of Vinogradov and Piatetski-Shapiro. Let $0<\gamma<1$, and let$$P_{\gamma}=\left\{ p\ :\ p=\lfloor n^{1/\gamma}\rfloor\ \text{for some }n\right\} $$ be the set of Piatetski-Shapiro primes. While this set of primes is very sparse, it has many nice properties and it is known that a version of the Prime number theorem holds when $\gamma$ is not too small (see this paper of Heath-Brown) $$\sum_{\begin{array}{c} p\leq x\\ p\in P_{\gamma} \end{array}}1\sim\frac{x^{\gamma}}{\log x}. $$ In the paper of Balog and Friedlander, they prove that for $\gamma$ in some range close to $1$, we have $$\frac{1}{\gamma^{3}}\sum_{\begin{array}{c} p_{1}+p_{2}+p_{3}=N\\ p_{i}\in P_{\gamma} \end{array}}p_{1}^{1-\gamma}p_{2}^{1-\gamma}p_{3}^{1-\gamma}\log p_{1}\log p_{2}\log p_{3}\sim\frac{1}{2}\mathfrak{S}(N)N^{2}, $$ where $\mathfrak{S}(N)$ is a particular singular series. Note that this implies that every sufficiently large integer can be written as the sum of three Piatetski-Shapiro primes, which proves a variant of Goldbach's Ternary conjecture for all sufficient large integers. From this they deduce without proof that $$\sum_{\begin{array}{c} p_{1}+p_{2}+p_{3}=N\\ p_{i}\in P_{\gamma} \end{array}}1\sim\frac{1}{2}\frac{\mathfrak{S}(N)N^{3\gamma-1}}{\left(\log N\right)^{3}}\frac{\gamma^{3}\Gamma(\gamma)^{3}}{\Gamma(3\gamma)}.$$ When I applied partial summation, this deduction reduced to evaluating the integral written above.
Check the signs. $$ \int_{\substack{ 0\leq s,t\leq1\\ s+t\leq1}}t^{\gamma-1}s^{\gamma-1}(1-t-s)^{\gamma-1}dt\,ds= \int_0^1 t^{\gamma-1}\int_0^{1-t}s^{\gamma-1}(1-t-s)^{\gamma-1}ds\,dt \\ = \int_0^1 t^{\gamma-1}(1-t)^{2\gamma-1}\int_0^{1}u^{\gamma-1}(1-u)^{\gamma-1}du\,dt = B(\gamma,2\gamma)B(\gamma,\gamma) = \frac{\Gamma(2\gamma)\Gamma(\gamma)^3}{\Gamma(2\gamma)\Gamma(3\gamma)} = \frac{\Gamma(\gamma)^3}{\Gamma(3\gamma)}. $$