結果

問題 No.1127 変形パスカルの三角形
ユーザー yuppe19 😺
提出日時 2020-07-26 12:06:26
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 478 ms / 1,500 ms
コード長 962 bytes
コンパイル時間 94 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 37,376 KB
最終ジャッジ日時 2024-06-28 11:42:05
合計ジャッジ時間 9,655 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

#!/usr/bin/env python3
# †
M = 10**9 + 7

###
if __name__ == '__main__':
    a, b = map(int, input().split())
    n, k = map(int, input().split())
    N = n+5
    fac = [None] * N
    fac[0] = 1
    for i in range(N-1):
        fac[i+1] = fac[i] * (i+1) % M
    inv = [None] * N
    inv[N-1] = pow(fac[N-1], M-2, M)
    for i in reversed(range(N-1)):
        inv[i] = inv[i+1] * (i+1) % M

    def comb(n, k):
        return fac[n] * inv[n-k] * inv[k] % M

    def solve1(n, k):
        n -= 1
        k -= 1
        x = comb(n, k)
        y = comb(n, k-1)
        res = a * x + b * y
        res %= M
        return res

    def solve2(n):
        n -= 1
        L = [comb(n, i) for i in range(n+1)]
        xs = L + [0]
        ys = [0] + L
        res = 0
        for x, y in zip(xs, ys):
            res += a*a*x*x
            res += b*b*y*y
            res += a*b*x*y*2
            res %= M
        return res

    print(solve1(n, k))
    print(solve2(n))
0