为计算 nCr 值而编写的程序 (C++)答案

【问题标题】：Program written to calculate nCr value (C++)为计算 nCr 值而编写的程序 (C++)
【发布时间】：2014-10-31 19:30:09
【问题描述】：

我已经编写了下面的代码来计算 n 个 testCount 值的 nCr。它似乎对我的价值观很有效，但黑客地球的测试用例失败了，问题链接http://www.hackerearth.com/problem/golf/minimal-combinatorial/description/ 谁能告诉我我的逻辑中潜在的谬误是什么

`#include <iostream>
 using namespace std;

 int main()
{
    int testCount; 
    int n , r;
    cin >> testCount;
    for (int i = 0 ; i < testCount ; i++)
    {

        long int a=1;
        long int b=1;
        long int c=1;
        cin >> n;
        cin >> r;
        for (int i = n ; i > 1 ; i--)
        {
             a = a * i;
        }
        for (int i = n-r; i >= 1 ; i--)
        {
            b=b*i;
        }
       for (int i=r;i >=1;i--)
       {
           c=c*i;
       }
       long int result=a/(b*c);
       cout << result<<"\n";
    }

    return 0;
}

` 简化分子和分母的情况

 `#include <iostream>
  using namespace std;

int main()
{
   int testCount; 
   int n , r;
   cin >> testCount;
  for (int i = 0 ; i < testCount ; i++)
  {

       long int numerator=1;
       long int denominator=1;
       cin >> n;
       cin >> r;
       for (int i = n ; i > r ; i--)
       {
          numerator = numerator * i;
       }
       for (int i = n-r; i >= 1 ; i--)
       {
           denominator=denominator*i;
       } 

       long int result=numerator/denominator;
       cout << result;
   }

    return 0;

} `

【问题讨论】：

你可能有溢出，你可以从分子/分母简化自己的公因数。
相关which-is-better-way-to-calculate-ncr
@Jarod42 我也这样做了，但它也失败了。我不应该溢出，因为他们已将值限制为 64 位 int 。简化版正在编辑中
最终结果适合 64 位整数，但 b * c 可能不适合。
您必须在第二个样本中使用r = max(r, n - r) 以避免其他可能的溢出。

标签： c++ binomial-coefficients

【解决方案1】：

您的两个代码都可能溢出（即使结果应该适合 64 位整数）

你可以试试：

std::uint64_t cnr(int n, int r)
{
    std::uint64_t ans = 1;
    r = std::min(r, n - r);

    for(int j = 1; j <= r; j++, n--) {
        // all following do: ans = (ans * n) / j;
        // but the 2 first cases do some simplifications
        // to reduce the possibility of overflow

        if (n % j == 0) {
            ans *= n / j; 
        } else if (ans % j == 0) {
            ans = (ans / j) * n;
        } else {
            ans = (ans * n) / j;
        }
    }
    return ans;
}

【讨论】：