結果

問題 No.1596 Distance Sum in 2D Plane
ユーザー moharan627moharan627
提出日時 2021-08-05 20:52:50
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 477 ms / 2,000 ms
コード長 1,440 bytes
コンパイル時間 333 ms
コンパイル使用メモリ 82,224 KB
実行使用メモリ 80,540 KB
最終ジャッジ日時 2024-09-16 15:24:29
合計ジャッジ時間 7,476 ms
ジャッジサーバーID
(参考情報)
judge2 / judge4
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 48 ms
62,400 KB
testcase_01 AC 47 ms
62,064 KB
testcase_02 AC 476 ms
80,192 KB
testcase_03 AC 477 ms
80,428 KB
testcase_04 AC 456 ms
80,184 KB
testcase_05 AC 457 ms
79,796 KB
testcase_06 AC 461 ms
80,388 KB
testcase_07 AC 462 ms
79,996 KB
testcase_08 AC 463 ms
80,028 KB
testcase_09 AC 464 ms
80,056 KB
testcase_10 AC 468 ms
79,940 KB
testcase_11 AC 436 ms
80,056 KB
testcase_12 AC 438 ms
80,532 KB
testcase_13 AC 433 ms
80,540 KB
testcase_14 AC 45 ms
62,080 KB
testcase_15 AC 46 ms
62,080 KB
testcase_16 AC 46 ms
61,792 KB
testcase_17 AC 46 ms
63,288 KB
testcase_18 AC 46 ms
63,160 KB
testcase_19 AC 46 ms
61,828 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

import sys
INF = float('inf')
#10**20,2**63,float('inf')
MOD = 10**9 + 7
MOD2 = 998244353
#from collections import defaultdict
def solve():
    def II(): return int(sys.stdin.readline())
    def LI(): return list(map(int, sys.stdin.readline().split()))
    def LC(): return list(input())
    def IC(): return [int(c) for c in input()]
    def MI(): return map(int, sys.stdin.readline().split())
    N,M = MI()
    def modinv(x):
        return pow(x, MOD - 2, MOD)
    MAX = 500050
    Fact = [0]*MAX
    Fact[0] = 1
    Fact[1] = 1
    for n in range(2,MAX):
        Fact[n] = (n*(Fact[n-1]))%MOD
    #print((Fact[N+K]*modinv(Fact[N]*Fact[K])) % MOD)
    All = Fact[2*N]*modinv(Fact[N]*Fact[N]) % MOD
    All *= (2*N)
    All %= MOD
    #print(All)
    for _ in range(M):
        T,X,Y = MI()
        Bef = Fact[X+Y]*modinv(Fact[X]*Fact[Y])%MOD
        Up = N-X
        Right = N-Y
        if T==1:
            #print(Fact[(Up-1)+Right]*modinv(Fact[Up-1]*Fact[Right]) % MOD)
            Tmp = (Fact[(Up-1)+Right]*modinv(Fact[Up-1]*Fact[Right]))%MOD
            All-= Bef*Tmp
            All %= MOD
        else:
            #print(Fact[Up+ Right-1] * modinv(Fact[Up] * Fact[Right-1]) % MOD)
            Tmp = (Fact[Up+ Right-1] * modinv(Fact[Up] * Fact[Right-1]))%MOD
            All -= Bef*Tmp
            All %= MOD
    print(All)
    return
solve()
#sys.setrecursionlimit(10 ** 6)#再帰関数ではコメントにしないこと!!
0