結果

問題 No.109 N! mod M
コンテスト
ユーザー cologne
提出日時 2022-02-02 17:12:30
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 174 ms / 5,000 ms
コード長 884 bytes
コンパイル時間 203 ms
コンパイル使用メモリ 82,560 KB
実行使用メモリ 75,264 KB
最終ジャッジ日時 2024-06-11 09:29:25
合計ジャッジ時間 1,334 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 9
権限があれば一括ダウンロードができます

ソースコード

diff # 

import itertools
import sys


def input(): return sys.stdin.readline().rstrip()


def isPrime(N):
    for i in itertools.count(2):
        if i*i > N:
            return True
        if N % i == 0:
            return False


def fact(N, M):
    # M - N <= LIM
    LIM = 10**5

    if N >= M:
        # trivial case
        return 0
    if M <= 2*LIM:
        # small case
        a = 1 % M
        for i in range(1, N+1):
            a = a * i % M
        return a
    elif isPrime(M):
        # prime case, (M-1)! == -1 (mod M)
        a = M-1
        for i in range(N+1, M):
            a = a * i % M
        return pow(a, -1, M)
    else:
        # M = pq, q <= M/p <= M/2 <= M/2 + (M/2 - LIM) <= N
        return 0


def main():
    T = int(input())
    for i in range(T):
        N, M = map(int, input().split())
        print(fact(N, M))


if __name__ == '__main__':
    main()
0