public class NegativeBinomialDistribution : DiscreteDistribution
    {
        private int trials;

        private double p;

        public double ProbabilityOfSuccess => this.p;

        public int NumberOfTrials => this.trials;

        public override double Mean => (double)this.trials * (1.0 - this.p) / this.p;

        public override double Variance => (double)this.trials * (1.0 - this.p) / (this.p * this.p);

        public override double Skewness => (2.0 - this.p) / Math.Sqrt((double)this.trials * (1.0 - this.p));

        public override double Kurtosis => (6.0 - this.p * (6.0 - this.p)) / ((double)this.trials * (1.0 - this.p));

        private void Initialize(int numberOfTrials, double probability)
        {
            if (numberOfTrials < 0)
            {
                ThrowException.ArgumentOutOfRange("numberOfTrials");
            }
            if (probability < 0.0 || probability > 1.0)
            {
                ThrowException.ArgumentOutOfRange("probability");
            }
            this.trials = numberOfTrials;
            this.p = probability;
        }

        private static int GetRandomVariateInternal(System.Random random, int n, double p)
        {
            if (p == 1.0)
            {
                return n;
            }
            if (n == 1)
            {
                return GeometricDistribution.GetRandomVariate(random, p);
            }
            double randomVariate;
            do
            {
                randomVariate = GammaDistribution.GetRandomVariate(random, n, (1.0 - p) / p);
            }
            while (randomVariate == 0.0);
            return PoissonDistribution.GetRandomVariate(random, randomVariate);
        }

        public static int GetRandomVariate(System.Random random, int numberOfTrials, double probabilityOfSuccess)
        {
            if (random == null)
            {
                ThrowException.ArgumentNull("random");
            }
            if (numberOfTrials <= 0)
            {
                ThrowException.ArgumentOutOfRange("numberOfTrials");
            }
            if (probabilityOfSuccess < 0.0 || probabilityOfSuccess > 1.0)
            {
                ThrowException.ArgumentOutOfRange("probabilityOfSuccess");
            }
            return new NegativeBinomialDistribution(numberOfTrials, probabilityOfSuccess).GetRandomVariate(random);
        }

        public NegativeBinomialDistribution(int numberOfTrials, double probabilityOfSuccess)
        {
            if (numberOfTrials <= 0)
            {
                ThrowException.ArgumentOutOfRange("numberOfTrials");
            }
            if (probabilityOfSuccess < 0.0 || probabilityOfSuccess > 1.0)
            {
                ThrowException.ArgumentOutOfRange("probabilityOfSuccess");
            }
            this.Initialize(numberOfTrials, probabilityOfSuccess);
        }

        public override double Probability(int n)
        {
            if (n < this.trials)
            {
                return 0.0;
            }
            if (this.p == 1.0)
            {
                if (n == this.trials)
                {
                    return 1.0;
                }
                return 0.0;
            }
            return Combinatorics.Combinations(n + this.trials - 1, this.trials - 1) * ElementaryFunctions.Pow(this.p, this.trials) * ElementaryFunctions.Pow(1.0 - this.p, n);
        }

        public override double DistributionFunction(int n)
        {
            return GammaFunctions.RegularizedBeta(this.trials, (double)n + 1.0, this.p);
        }

        public override int GetRandomVariate(System.Random random)
        {
            if (random == null)
            {
                ThrowException.ArgumentNull("random");
            }
            return NegativeBinomialDistribution.GetRandomVariateInternal(random, this.trials, this.p);
        }
    }

 如果对您有帮忙，非常感谢您支持一下创造者的付出！

 感谢支持技术分享，请扫码点赞支持：

技术合作交流qq:2401315930