結果
問題 | No.1035 Color Box |
ユーザー |
![]() |
提出日時 | 2020-12-14 05:06:04 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 97 ms / 2,000 ms |
コード長 | 1,147 bytes |
コンパイル時間 | 137 ms |
コンパイル使用メモリ | 81,920 KB |
実行使用メモリ | 76,676 KB |
最終ジャッジ日時 | 2024-09-20 00:25:43 |
合計ジャッジ時間 | 3,272 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 36 |
ソースコード
import syssys.setrecursionlimit(10**7)def I(): return int(sys.stdin.readline().rstrip())def MI(): return map(int,sys.stdin.readline().rstrip().split())def LI(): return list(map(int,sys.stdin.readline().rstrip().split()))def LI2(): return list(map(int,sys.stdin.readline().rstrip()))def S(): return sys.stdin.readline().rstrip()def LS(): return list(sys.stdin.readline().rstrip().split())def LS2(): return list(sys.stdin.readline().rstrip())def factorial(N,p): # 0!~N!(mod p) のリストres = [0]*(N+1)res[0] = 1for i in range(1,N+1):res[i] = (res[i-1]*i) % preturn resdef factorial_inverse(N,p): # 0!~N!(mod p) の逆元のリストfac = factorial(N,p)res = [0]*(N+1)a = pow(fac[-1],p-2,p)for i in range(N,-1,-1):res[i] = aa *= ia %= preturn fac,resN,M = MI()mod = 10**9+7fac,fac_inv = factorial_inverse(M,mod)def nCr(n,r):return (fac[n]*fac_inv[r]*fac_inv[n-r]) % modans = 0for i in range(M+1):a = nCr(M,i)*pow(M-i,N,mod)a %= modif i % 2 == 0:ans += aelse:ans -= aans %= modprint(ans)