dobrodeev1978.py 文件源码

def _generate_jk(n, pm_type, j, k):
    M = fact(n) // fact(j) // fact(k) // fact(n-j-k) * 2**(j+k)
    G = fr(1, M)

    t = 1 if pm_type == 'I' else -1
    b = sqrt(fr(1, j+k) * (1 + t * (k/j * sqrt(3*(j+k)/(n+2) - 1))))
    c = sqrt(fr(1, j+k) * (1 - t * (j/k * sqrt(3*(j+k)/(n+2) - 1))))
    return G, b, c


# pylint: disable=too-many-arguments, too-many-locals