0

1500行目の合計を見つけることについてパスカルの三角形について別の質問をしました。人々がとても早く答えてくれてとてもうれしいですが、残念ながら後で気づきました。1500行目に個々の番号が必要です。

ここで、パスカルの三角形の任意の数を計算する簡単な方法を見つけましたが、コードで数式を使用しようとすると、プログラムが起動時にクラッシュします。

#include"stdio.h"
int factorial(int);

int main()
{
    int i=0;
    for(i=0; i<1501; i++)
    {
       printf("%d \n" , (factorial(1500)/factorial(1500-i))/factorial(i) );
    }
}
int factorial(int x)
{
    if(x<2)
    return 1;
    else
    {
        return x*factorial(x-1);
    }
}
4

2 に答える 2

3

結果にはbigintegerライブラリが必要です。

binom(1500,750) = 722462892448930217028116073228485295944376343040546523665632913653613596406381727723198055046187955623124069093412562720869577867557610868874271130828359513282184417776042792372322074253393127328396528996120053749558122610911178218582669317535346728464707661495135518682519172221470420360910320792434869988224466647627642393919250205687942318888922893189087379790541907686956429837978631252775258630376332505697937034877619012586751274240109457111424

通常のdouble範囲を超えています。

しかし、そのようなタイプが利用可能であると、関係

binom(n,k+1) = binom(n,k)*(n-k)/(k+1)

シンプルで比較的効率的な実装を可能にします。

bigint value = 1;
int numerator = 1500 + 1, denominator = 0;
for(; denominator <= 1500; --numerator, ++denominator, value = value*numerator/denominator)
{
    output(value);
}

outputライブラリが提供する出力メソッドはどこにありますか。

value = value*numerator/denominatorライブラリがbigintegersとnormalを乗算/除算するためのオーバーロードを提供しない場合は、パーツを適合させる必要がありますint。

于 2012-09-29T15:05:03.770 に答える
3

代わりに階乗に線形アルゴリズムを使用して、再帰によるスタック オーバーフローを回避します。

int factorial(int x) {
    if (x < 2) return 1;

    int result = x;

    while (--x) {
        result *= x;
    }

    return result;
}
于 2012-09-29T13:57:39.067 に答える