結果
問題 | No.1596 Distance Sum in 2D Plane |
ユーザー | shotoyoo |
提出日時 | 2021-07-09 21:37:55 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 233 ms / 2,000 ms |
コード長 | 1,255 bytes |
コンパイル時間 | 345 ms |
コンパイル使用メモリ | 86,852 KB |
実行使用メモリ | 87,452 KB |
最終ジャッジ日時 | 2023-09-14 08:03:37 |
合計ジャッジ時間 | 5,325 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge12 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 100 ms
85,084 KB |
testcase_01 | AC | 102 ms
85,156 KB |
testcase_02 | AC | 233 ms
86,748 KB |
testcase_03 | AC | 222 ms
86,752 KB |
testcase_04 | AC | 224 ms
86,760 KB |
testcase_05 | AC | 212 ms
86,524 KB |
testcase_06 | AC | 217 ms
86,648 KB |
testcase_07 | AC | 210 ms
86,744 KB |
testcase_08 | AC | 207 ms
86,472 KB |
testcase_09 | AC | 211 ms
86,396 KB |
testcase_10 | AC | 211 ms
86,492 KB |
testcase_11 | AC | 195 ms
86,840 KB |
testcase_12 | AC | 194 ms
87,452 KB |
testcase_13 | AC | 194 ms
87,112 KB |
testcase_14 | AC | 72 ms
71,088 KB |
testcase_15 | AC | 72 ms
71,088 KB |
testcase_16 | AC | 74 ms
71,096 KB |
testcase_17 | AC | 72 ms
71,300 KB |
testcase_18 | AC | 74 ms
71,084 KB |
testcase_19 | AC | 73 ms
71,212 KB |
ソースコード
import sys input = lambda : sys.stdin.readline().rstrip() sys.setrecursionlimit(2*10**5+10) write = lambda x: sys.stdout.write(x+"\n") debug = lambda x: sys.stderr.write(x+"\n") writef = lambda x: print("{:.12f}".format(x)) n,m = list(map(int, input().split())) ### 素数の逆元とCombination M = 10**9+7 # 出力の制限 N = 2*n+10 # 必要なテーブルサイズ g1 = [0] * (N+1) # 元テーブル g2 = [0] * (N+1) #逆元テーブル inverse = [0] * (N+1) #逆元テーブル計算用テーブル g1[0] = g1[1] = g2[0] = g2[1] = 1 inverse[0], inverse[1] = [0, 1] for i in range( 2, N + 1 ): g1[i] = ( g1[i-1] * i ) % M inverse[i] = ( -inverse[M % i] * (M//i) ) % M # ai+b==0 mod M <=> i==-b*a^(-1) <=> i^(-1)==-b^(-1)*aより g2[i] = (g2[i-1] * inverse[i]) % M def cmb(n, r, M=M): if ( r<0 or r>n ): return 0 r = min(r, n-r) return ((g1[n] * g2[r] % M) * g2[n-r]) % M def perm(n, r, M=M): if (r<0 or r>n): return 0 return (g1[n] * g2[n-r]) % M ans = 2*n*cmb(2*n,n)%M for i in range(m): t,x,y = map(int, input().split()) if t==1: ans -= cmb(x+y,x) * cmb(2*n - x - y - 1, n-x-1) % M else: ans -= cmb(x+y,x) * cmb(2*n - x - y - 1, n-y-1) % M ans %= M print(ans%M)