dp[0][0] = 1
for i in [1, R*S]
 for j in [1, R]
  if i ≦ j*S
   dp[i][j] = dp[i-1][j-1] + dp[i-1][j]
return dp[R*S][R] * factorial(R)