交代置換

Question

この投稿は私の前の投稿の拡張です：小さい要素と大きい要素の順列

私は交代置換に取り組んでいます。これが私の論理です。

function Factorial(n) { 
    
    var res=1; 
          
    for (var i = 2; i <= n; i++) 
        res = res * i; 
    return res; 
} 
    
    
let n = 4;
let A = [];
let C = [];
let a = Factorial(n);
for(let i=0; i<=n;i++) {
    A[i] = 0;
}
A[1] = 1;
for(let k=0; k<n; k++) {
    let b = Factorial(k)*Factorial(n-k);
        
    A[k] = a/b * A[k]*A[n-k]/2;
}
console.log(A);

prints [0, 0, 0, 0]

以前の投稿のように、入力n = 4に対してA [n + 1] = 5を期待しています。しかし、私はすべてゼロになっています。この問題を修正する方法。

Answer 1

参照しているアンドレの式に合計（∑）を実装していません。

また、浮動小数点の精度の制限にすぐにぶつかるので、階乗の計算は避けます。代わりに、パスカルの三角形を使用できます。

これは、Aシーケンスの無限ジェネレーターの実装です。デモでは、A [0]からA [20]までの値が出力されます。

function * andre() { // generator
    let pascal = [1];
    let a = [1, 1];
    yield a[0];
    
    let n = 1;
    while (true) {
        yield a[n];
        // Build the next row in Pascal's Triangle
        pascal[n] = 1;
        for (let i = n - 1; i > 0; i--) {
            pascal[i] += pascal[i-1];
        }
        // Apply André's formula
        let sum = 0;
        for (let k = 0; k <= n; k++) {
            sum += pascal[k] * a[k] * a[n-k]
        }
        a[++n] = sum / 2;
    }
}

// demo: display A[0] up to A[20]
let i = 0;
for (let a of andre()) {
    console.log(`A[${i++}] = ${a}`);
    if (i > 20) break;
}

このシーケンスは急速に増加するため、すぐに高い精度が必要になります。あなたが必要な場合はA [i]のをより高い値のために、私は、JavaScriptを使用してのBigIntデータ型を。