結果

問題 No.1127 変形パスカルの三角形
ユーザー yuppe19 😺yuppe19 😺
提出日時 2020-07-26 12:06:26
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 375 ms / 1,500 ms
コード長 962 bytes
コンパイル時間 84 ms
コンパイル使用メモリ 10,984 KB
実行使用メモリ 34,804 KB
最終ジャッジ日時 2023-09-10 20:35:19
合計ジャッジ時間 8,400 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 17 ms
7,916 KB
testcase_01 AC 375 ms
34,704 KB
testcase_02 AC 267 ms
26,636 KB
testcase_03 AC 199 ms
21,812 KB
testcase_04 AC 60 ms
11,504 KB
testcase_05 AC 177 ms
20,444 KB
testcase_06 AC 372 ms
34,804 KB
testcase_07 AC 119 ms
15,852 KB
testcase_08 AC 114 ms
15,612 KB
testcase_09 AC 305 ms
30,176 KB
testcase_10 AC 317 ms
30,504 KB
testcase_11 AC 264 ms
26,828 KB
testcase_12 AC 245 ms
25,472 KB
testcase_13 AC 231 ms
24,200 KB
testcase_14 AC 191 ms
21,368 KB
testcase_15 AC 105 ms
14,948 KB
testcase_16 AC 209 ms
22,596 KB
testcase_17 AC 113 ms
15,412 KB
testcase_18 AC 222 ms
23,860 KB
testcase_19 AC 268 ms
27,108 KB
testcase_20 AC 277 ms
27,760 KB
testcase_21 AC 279 ms
28,100 KB
testcase_22 AC 109 ms
15,320 KB
testcase_23 AC 312 ms
30,316 KB
testcase_24 AC 197 ms
21,956 KB
testcase_25 AC 210 ms
22,620 KB
testcase_26 AC 229 ms
24,168 KB
testcase_27 AC 263 ms
27,056 KB
testcase_28 AC 255 ms
26,544 KB
testcase_29 AC 103 ms
14,788 KB
testcase_30 AC 76 ms
12,744 KB
testcase_31 AC 284 ms
27,880 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