結果
問題 | No.573 a^2[i] = a[i] |
ユーザー |
|
提出日時 | 2022-02-11 14:39:52 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 283 ms / 2,000 ms |
コード長 | 1,204 bytes |
コンパイル時間 | 134 ms |
コンパイル使用メモリ | 82,172 KB |
実行使用メモリ | 94,764 KB |
最終ジャッジ日時 | 2024-06-27 08:49:22 |
合計ジャッジ時間 | 3,731 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 47 |
ソースコード
class NCK:'''最大値 :Nmod'''def __init__(self,N,mod):self.mod = modself.N = Nself.factmod = [1,1]for i in range(2,self.N+1):self.factmod.append((self.factmod[-1]*i)%self.mod)self.inv = [0,1]for i in range(2, N + 1):self.inv.append((-self.inv[self.mod % i] * (self.mod // i)) % self.mod)self.invfact = [1,1]for i in range(2, N + 1):self.invfact.append((self.invfact[-1]*self.inv[i])%self.mod)def nCk(self,a,b):if a < b:return 0return ((self.factmod[a]*self.invfact[b]%self.mod)*self.invfact[a-b])%self.moddef invnCk(self,a,b):return ((self.factmod[a-b]*self.invfact[a]%self.mod)*self.factmod[b])%self.moddef modPow(a,n,mod):#繰り返し二乗法 a**n % modif n==0:return 1if n==1:return a%modif n & 1:return (a*modPow(a,n-1,mod)) % modt = modPow(a,n>>1,mod)return (t*t)%modN = int(input())# Σ(k=1,N) NCk*K**(N-K)mod = 10**9+7ans = 0nck = NCK(N+1,mod)for i in range(1,N+1):ans += nck.nCk(N,i)*modPow(i,N-i,mod)ans %= modprint(ans)