結果
問題 | No.1035 Color Box |
ユーザー |
![]() |
提出日時 | 2020-04-24 22:42:56 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 428 ms / 2,000 ms |
コード長 | 1,060 bytes |
コンパイル時間 | 257 ms |
コンパイル使用メモリ | 81,964 KB |
実行使用メモリ | 247,160 KB |
最終ジャッジ日時 | 2024-10-15 03:16:04 |
合計ジャッジ時間 | 12,690 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 36 |
ソースコード
import sys mod = pow(10, 9) + 7 sys.setrecursionlimit(pow(10, 8)) def power(x, y): if y == 0: return 1 elif y == 1 : return x % mod elif y % 2 == 0 : return power(x, y//2)**2 % mod else: return power(x, (y-1)//2)**2 * x % mod def mul(a, b): return ((a % mod) * (b % mod)) % mod def div(a, b): return mul(a, power(b, mod-2)) def div2(a, b): return mul(a, modinv(b)) def modinv(a): b, u, v = mod, 1, 0 while b: t = a//b a, u = a-t*b, u-t*v a, b, u, v = b, a, v, u u %= mod return u def cmb(n, r): if ( r<0 or r>n ): return 0 r = min(r, n-r) return g1[n] * g2[r] * g2[n-r] % mod NNN = (10**6) g1 = [1, 1] g2 = [1, 1] inverse = [0, 1] for i in range( 2, NNN + 1 ): g1.append( ( g1[-1] * i ) % mod ) inverse.append( ( -inverse[mod % i] * (mod//i) ) % mod ) g2.append( (g2[-1] * inverse[-1]) % mod ) N, M = map(int, input().split()) R = 0 for i in range(1, M+1): R += (-1 if (M-i)%2 else 1) * (mul(cmb(M, i), power(i, N))) R %= mod print(R)