Consider an \(n\times k\) matrix of i.i.d. Bernoulli random numbers with \(p=1/2\) . The dual RSK algorithm gives a bijection of this matrix to a pair of Young tableaux of conjugate shape, which is manifestation of skew Howe \(GL_{n}\times GL_{k}\) -duality. Thus the probability measure on zero–ones matrix leads to the probability measure on Young diagrams proportional to the ratio of the dimension of \(GL_{n}\times GL_{k}\) -representation and the dimension of the exterior algebra \(\bigwedge \left( \mathbb {C}^{n}\otimes \mathbb {C}^{k}\right) \) . Similarly, by applying Proctor’s algorithm based on Berele’s modification of the Schensted insertion, we get skew Howe duality for the pairs of groups \(Sp_{2n}\times Sp_{2k}\) . In the limit when \(n,k\rightarrow \infty \) \(GL\) -case is relatively easily studied by use of free-fermionic representation for the correlation kernel. But for the symplectic groups there is no convenient free-fermionic representation. We use Christoffel transformation to obtain the semiclassical orthogonal polynomials for \(Sp_{2n}\times Sp_{2k}\) from Krawtchouk polynomials that describe \(GL_{2n}\times GL_{2k}\) case. We derive an integral representation for semiclassical polynomials. The study of the asymptotic of this integral representation gives us the description of the limit shapes and fluctuations of the random Young diagrams for symplectic groups.
