We consider the critical p -Laplacian system { - Δ p u - λ a p | u | a - 2 u | v | b = μ 1 | u | p ∗ - 2 u + α γ p ∗ | u | α - 2 u | v | β , x ∈ Ω , - Δ p v - λ b p | u | a | v | b - 2 v = μ 2 | v | p ∗ - 2 v + β γ p ∗ | u | α | v | β - 2 v , x ∈ Ω , u , v in D 0 1 , p ( Ω ) , \left\{\begin{aligned} &\displaystyle{-}\Delta_{p}u-\frac{\lambda a}{p}\lvert u% \rvert^{a-2}u\lvert v\rvert^{b}=\mu_{1}\lvert u\rvert^{p^{\ast}-2}u+\frac{% \alpha\gamma}{p^{\ast}}\lvert u\rvert^{\alpha-2}u\lvert v\rvert^{\beta},&&% \displaystyle x\in\Omega,\\ &\displaystyle{-}\Delta_{p}v-\frac{\lambda b}{p}\lvert u\rvert^{a}\lvert v% \rvert^{b-2}v=\mu_{2}\lvert v\rvert^{p^{\ast}-2}v+\frac{\beta\gamma}{p^{\ast}}% \lvert u\rvert^{\alpha}\lvert v\rvert^{\beta-2}v,&&\displaystyle x\in\Omega,\\ &\displaystyle u,v\text{ in }D_{0}^{1,p}(\Omega),\end{aligned}\right. where Δpu:=div(|∇u|p-2∇u){\Delta_{p}u:=\operatorname{div}(\lvert\nabla u\rvert^{p-2}\nabla u)} is the p -Laplacian operator defined on D 1 , p ( ℝ N ) := { u ∈ L p ∗ ( ℝ N ) : | ∇ u | ∈ L p ( ℝ N ) } , D^{1,p}(\mathbb{R}^{N}):=\bigl{\{}u\in L^{p^{\ast}}(\mathbb{R}^{N}):\lvert% \nabla u\rvert\in L^{p}(\mathbb{R}^{N})\bigr{\}}, endowed with the norm ∥u∥D1,p:=(∫ℝN|∇u|p𝑑x)1p{{\lVert u\rVert_{D^{1,p}}:=(\int_{\mathbb{R}^{N}}\lvert\nabla u\rvert^{p}\,dx% )^{\frac{1}{p}}}}, N≥3{N\geq 3}, 1<p<N{1<p<N}, λ,μ1,μ2≥0{\lambda,\mu_{1},\mu_{2}\geq 0}, γ≠0{\gamma\neq 0}, a,b,α,β>1{a,b,\alpha,\beta>1} satisfy a+b=p{a+b=p}, α+β=p∗:=NpN-p{\alpha+\beta=p^{\ast}:=\frac{Np}{N-p}}, the critical Sobolev exponent, Ω is ℝN{\mathbb{R}^{N}} or a bounded domain in ℝN{\mathbb{R}^{N}} and D01,p(Ω){D_{0}^{1,p}(\Omega)} is the closure of C0∞(Ω){C_{0}^{\infty}(\Omega)} in D1,p(ℝN){D^{1,p}(\mathbb{R}^{N})}. Under suitable assumptions, we establish the existence and nonexistence of a positive least energy solution of this system. We also consider the existence and multiplicity of the nontrivial nonnegative solutions.