计算 fact(y) mod x

最新推荐文章于 2022-12-04 07:22:30 发布

GGN_2015

最新推荐文章于 2022-12-04 07:22:30 发布

阅读量221

点赞数

分类专栏：数学文章标签： python

本文链接：https://blog.csdn.net/ggn_2015/article/details/121025085

版权

数学专栏收录该内容

32 篇文章 2 订阅

订阅专栏

题面

已知

x = 11124440021748127159092076861405454814981575144744508857178576572929321435002942998531420985771090167262256877805902135304112271641074498386662361391760451
y = 11124440021748127159092076861405454814981575144744508857178576572929321435002942998531420985771090167262256877805902135304112271641074498386662361391661439

计算

$y!\pmod x$

的值。

做法

通过米勒拉宾质数测试，不难发现 $x$ 是质数，并且 $x - y = 99012$ 很小。根据威尔逊定理可知：

米勒拉宾素数测试

$(x-1)!\mod x =-1$

变形，记 $t = 99012$ ：

$(x-1)(x-2)\cdots(x-t+1)y!\mod x = -1$

因此：

$y!\mod x = -1\cdot\left((x-1)(x-2)\cdots(x-t+1)\right)^{-1}$

通过费马小定理计算 $\left((x-1)(x-2)\cdots(x-t+1)\right)^{-1}$ 即可。

代码

x = 11124440021748127159092076861405454814981575144744508857178576572929321435002942998531420985771090167262256877805902135304112271641074498386662361391760451
y = 11124440021748127159092076861405454814981575144744508857178576572929321435002942998531420985771090167262256877805902135304112271641074498386662361391661439

left = 1 # 从 x-1 乘到 y+1
for i in range(1, t):
    left *= (x - i)
    left %= x
    
def qpow(A, B, P): # 快速幂
    ans = 1
    mul = A % P
    while (B != 0):
        if((B & 1) == 1):
            ans *= mul
            ans %= P
        B >>= 1
        mul *= mul
        
        mul %= P
    return ans

right = qpow(left, x - 2, x)
ans = (-right % x + x)%x
print(ans)

结果

10569944080090591401315432556965818857327680380269154543273468441025963038065648915158194147019839932524599260058098616377893091051396090650574162446875193