printf <- function(...) cat(sprintf(...))
 
ithDuplicatedPermutation <- function(n, m, i)
{
    c <- rep(m, n)
 
    N <- n * m
    a <- integer(N)
 
    P <- factorial(N) / factorial(m) ^ n
 
    for (j in 1:N) {
        for (k in 1:n) {
            if (!c[k]) next
            p <- P * c[k] / (N - j + 1)
            if (i <= p) {
                a[j] <- k
                c[k] <- c[k] - 1
                P <- p
                break
            }
            i <- i - p
        }
    }
 
    a
}
 
for (n in c(1, 2, 3, 123456, 234567, 369600)) {
    printf("入力: %s\n", n)
    printf("出力: [%s]\n\n", toString(ithDuplicatedPermutation(4, 3, n)))
}

cHJpbnRmIDwtIGZ1bmN0aW9uKC4uLikgY2F0KHNwcmludGYoLi4uKSkKCml0aER1cGxpY2F0ZWRQZXJtdXRhdGlvbiA8LSBmdW5jdGlvbihuLCBtLCBpKQp7CiAgICBjIDwtIHJlcChtLCBuKQoKICAgIE4gPC0gbiAqIG0KICAgIGEgPC0gaW50ZWdlcihOKQoKICAgIFAgPC0gZmFjdG9yaWFsKE4pIC8gZmFjdG9yaWFsKG0pIF4gbgoKICAgIGZvciAoaiBpbiAxOk4pIHsKICAgICAgICBmb3IgKGsgaW4gMTpuKSB7CiAgICAgICAgICAgIGlmICghY1trXSkgbmV4dAogICAgICAgICAgICBwIDwtIFAgKiBjW2tdIC8gKE4gLSBqICsgMSkKICAgICAgICAgICAgaWYgKGkgPD0gcCkgewogICAgICAgICAgICAgICAgYVtqXSA8LSBrCiAgICAgICAgICAgICAgICBjW2tdIDwtIGNba10gLSAxCiAgICAgICAgICAgICAgICBQIDwtIHAKICAgICAgICAgICAgICAgIGJyZWFrCiAgICAgICAgICAgIH0KICAgICAgICAgICAgaSA8LSBpIC0gcAogICAgICAgIH0KICAgIH0KCiAgICBhCn0KCmZvciAobiBpbiBjKDEsIDIsIDMsIDEyMzQ1NiwgMjM0NTY3LCAzNjk2MDApKSB7CiAgICBwcmludGYoIuWFpeWKmzogJXNcbiIsIG4pCiAgICBwcmludGYoIuWHuuWKmzogWyVzXVxuXG4iLCB0b1N0cmluZyhpdGhEdXBsaWNhdGVkUGVybXV0YXRpb24oNCwgMywgbikpKQp9