Correct calculation of random distribution in openPlatonic

[vinniefalco]

I believe this algorithm is correct, and runs in constant average time

void openPlatonic( int amount, int* blessings )
{
    std::random_device rd;
    double const weights[] = {
        13200.0, 13200.0, 13200.0, 13200.0,
        396.0, 1.0, 1.0, 1.0, 1.0 };
    double sum = 40000.0;
    for( int i = 0; i < 8; ++ i)
    {
        int n = std::binomial_distribution<int>(
            amount, weights[i] / sum)(rd);
        blessings[i] += n;
        amount -= n;
        sum -= weights[i];
    }
}

It would need to be ported from C++ to JS, including an implementation of binomial_distribution.