結果
| 問題 | No.1102 Remnants | 
| コンテスト | |
| ユーザー |  | 
| 提出日時 | 2020-07-05 10:27:06 | 
| 言語 | PyPy3 (7.3.15) | 
| 結果 | 
                                AC
                                 
                             | 
| 実行時間 | 145 ms / 2,000 ms | 
| コード長 | 1,535 bytes | 
| コンパイル時間 | 384 ms | 
| コンパイル使用メモリ | 82,572 KB | 
| 実行使用メモリ | 120,020 KB | 
| 最終ジャッジ日時 | 2024-09-22 05:19:10 | 
| 合計ジャッジ時間 | 5,026 ms | 
| ジャッジサーバーID (参考情報) | judge1 / judge3 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 3 | 
| other | AC * 25 | 
ソースコード
import sys
sys.setrecursionlimit(10 ** 7)
rl = sys.stdin.readline
class Combination:
    def __init__(self, n: int, mod: int):
        self.mod = mod
        self.fact = [0] * (n + 1)
        self.factinv = [0] * (n + 1)
        self.inv = [0] * (n + 1)
        
        self.fact[0] = self.fact[1] = 1
        self.factinv[0] = self.factinv[1] = 1
        self.inv[1] = 1
        for i in range(2, n + 1):
            self.fact[i] = (self.fact[i - 1] * i) % mod
            self.inv[i] = (-self.inv[mod % i] * (mod // i)) % mod
            self.factinv[i] = (self.factinv[i - 1] * self.inv[i]) % mod
    
    def ncr(self, n: int, r: int):
        if r < 0 or n < r:
            return 0
        r = min(r, n - r)
        return self.fact[n] * self.factinv[r] % self.mod * self.factinv[n - r] % self.mod
    
    def nhr(self, n: int, r: int):
        return self.ncr(n + r - 1, r)
    
    def npr(self, n: int, r: int):
        if r < 0 or n < r:
            return 0
        return self.fact[n] * self.factinv[n - r] % self.mod
def solve():
    N, K = map(int, rl().split())
    A = list(map(int, rl().split()))
    MOD = 10 ** 9 + 7
    
    MAX_SIZE = 2 * 10 ** 5 + 100
    com = Combination(MAX_SIZE, MOD)
    nhk = [1]
    for i in range(1, MAX_SIZE):
        nhk.append(nhk[-1] * (K + i) % MOD * com.factinv[i] % MOD * com.fact[i - 1] % MOD)
    
    ans = 0
    for i, ai in enumerate(A):
        ans += nhk[i] * nhk[N - i - 1] % MOD * ai % MOD
        ans %= MOD
    print(ans)
if __name__ == '__main__':
    solve()
            
            
            
        