結果

問題 No.1596 Distance Sum in 2D Plane
ユーザー moharan627
提出日時 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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

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