結果
問題 |
No.1035 Color Box
|
ユーザー |
![]() |
提出日時 | 2020-04-25 12:18:47 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 526 ms / 2,000 ms |
コード長 | 837 bytes |
コンパイル時間 | 324 ms |
コンパイル使用メモリ | 82,432 KB |
実行使用メモリ | 246,272 KB |
最終ジャッジ日時 | 2024-11-06 22:34:10 |
合計ジャッジ時間 | 17,705 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 36 |
ソースコード
n, m = map(int, input().split()) mod = 10**9*7 def cmb1(n, r, mod): if ( r<0 or r>n ): return 0 r = min(r, n-r) return g1[n] * g2[r] * g2[n-r] % mod mod = 10**9+7 #出力の制限 N = 10**6 g1 = [1, 1] # 元テーブル g2 = [1, 1] #逆元テーブル inverse = [0, 1] #逆元テーブル計算用テーブル for i in range( 2, N + 1 ): g1.append( ( g1[-1] * i ) % mod ) inverse.append( ( -inverse[mod % i] * (mod//i) ) % mod ) g2.append( (g2[-1] * inverse[-1]) % mod ) def power(a, n, mod): bi=str(format(n,"b")) #2進数 res=1 for i in range(len(bi)): res=(res*res) %mod if bi[i]=="1": res=(res*a) %mod return res res = 0 for i in range(m+1): temp = cmb1(m, i, mod) * power(i, n, mod) if (m-i)%2 == 0: res += temp else: res -= temp print(res%mod)