結果

問題 No.1596 Distance Sum in 2D Plane
ユーザー moharan627moharan627
提出日時 2021-08-05 20:52:50
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 486 ms / 2,000 ms
コード長 1,440 bytes
コンパイル時間 294 ms
コンパイル使用メモリ 87,140 KB
実行使用メモリ 81,928 KB
最終ジャッジ日時 2023-10-14 21:41:06
合計ジャッジ時間 7,828 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 82 ms
79,512 KB
testcase_01 AC 79 ms
79,564 KB
testcase_02 AC 486 ms
81,732 KB
testcase_03 AC 474 ms
81,860 KB
testcase_04 AC 475 ms
81,304 KB
testcase_05 AC 460 ms
81,260 KB
testcase_06 AC 465 ms
81,356 KB
testcase_07 AC 462 ms
81,332 KB
testcase_08 AC 463 ms
81,200 KB
testcase_09 AC 458 ms
81,168 KB
testcase_10 AC 459 ms
81,188 KB
testcase_11 AC 445 ms
81,648 KB
testcase_12 AC 443 ms
81,928 KB
testcase_13 AC 435 ms
81,428 KB
testcase_14 AC 77 ms
79,316 KB
testcase_15 AC 78 ms
79,344 KB
testcase_16 AC 78 ms
79,480 KB
testcase_17 AC 78 ms
79,380 KB
testcase_18 AC 78 ms
79,376 KB
testcase_19 AC 77 ms
79,348 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