結果

問題 No.1127 変形パスカルの三角形
ユーザー yuppe19 😺yuppe19 😺
提出日時 2020-07-26 12:06:26
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
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
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 30 ms
10,752 KB
testcase_01 AC 478 ms
37,248 KB
testcase_02 AC 327 ms
29,312 KB
testcase_03 AC 249 ms
24,448 KB
testcase_04 AC 82 ms
13,824 KB
testcase_05 AC 229 ms
22,912 KB
testcase_06 AC 454 ms
37,376 KB
testcase_07 AC 153 ms
18,432 KB
testcase_08 AC 151 ms
17,920 KB
testcase_09 AC 408 ms
32,512 KB
testcase_10 AC 395 ms
32,896 KB
testcase_11 AC 331 ms
29,312 KB
testcase_12 AC 308 ms
27,904 KB
testcase_13 AC 286 ms
26,752 KB
testcase_14 AC 237 ms
23,808 KB
testcase_15 AC 145 ms
17,536 KB
testcase_16 AC 260 ms
25,344 KB
testcase_17 AC 148 ms
17,792 KB
testcase_18 AC 275 ms
26,240 KB
testcase_19 AC 333 ms
29,568 KB
testcase_20 AC 345 ms
30,080 KB
testcase_21 AC 352 ms
30,592 KB
testcase_22 AC 146 ms
17,664 KB
testcase_23 AC 393 ms
32,896 KB
testcase_24 AC 259 ms
24,448 KB
testcase_25 AC 262 ms
25,344 KB
testcase_26 AC 298 ms
26,624 KB
testcase_27 AC 332 ms
29,440 KB
testcase_28 AC 326 ms
29,056 KB
testcase_29 AC 137 ms
17,152 KB
testcase_30 AC 100 ms
15,104 KB
testcase_31 AC 341 ms
30,592 KB
権限があれば一括ダウンロードができます

ソースコード

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