結果
問題 | No.1596 Distance Sum in 2D Plane |
ユーザー | ophhdn |
提出日時 | 2021-07-09 21:43:13 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 297 ms / 2,000 ms |
コード長 | 1,556 bytes |
コンパイル時間 | 153 ms |
コンパイル使用メモリ | 82,064 KB |
実行使用メモリ | 88,396 KB |
最終ジャッジ日時 | 2024-07-01 15:41:17 |
合計ジャッジ時間 | 5,482 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
from collections import defaultdict, deque, Counter from heapq import heappush, heappop, heapify import math import bisect import random from itertools import permutations, accumulate, combinations, product import sys import string from bisect import bisect_left, bisect_right from math import factorial, ceil, floor,atan2,sin,cos from operator import mul from functools import reduce from pprint import pprint sys.setrecursionlimit(2147483647) INF = 10 ** 13 def LI(): return list(map(int, sys.stdin.readline().split())) def I(): return int(sys.stdin.readline()) def LS(): return sys.stdin.buffer.readline().rstrip().decode('utf-8').split() def S(): return sys.stdin.buffer.readline().rstrip().decode('utf-8') def IR(n): return [I() for i in range(n)] def LIR(n): return [LI() for i in range(n)] def SR(n): return [S() for i in range(n)] def LSR(n): return [LS() for i in range(n)] def SRL(n): return [list(S()) for i in range(n)] def MSRL(n): return [[int(j) for j in list(S())] for i in range(n)] mod = 10**9+7 k,m=LI() n=2*k fac = [1] * (n + 1) inv = [1] * (n + 1) for j in range(1, n + 1): fac[j] = fac[j-1] * j % mod inv[n] = pow(fac[n], mod-2, mod) for j in range(n-1, -1, -1): inv[j] = inv[j+1] * (j+1) % mod def comb(n, r): if r > n or n < 0 or r < 0: return 0 return fac[n] * inv[n - r] * inv[r] % mod ans=comb(k*2,k)*(2*k)%mod for i in range(m): t,x,y=LI() if t==1: ans-=comb(x+y,y)*comb(k-x-1+k-y,k-y) else: ans -= comb(x + y, y) * comb(k- y - 1 +k- x,k- x) print(ans%mod)