結果
問題 | No.1596 Distance Sum in 2D Plane |
ユーザー |
![]() |
提出日時 | 2023-04-18 22:55:13 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 188 ms / 2,000 ms |
コード長 | 1,525 bytes |
コンパイル時間 | 207 ms |
コンパイル使用メモリ | 82,688 KB |
実行使用メモリ | 101,504 KB |
最終ジャッジ日時 | 2024-10-13 22:48:55 |
合計ジャッジ時間 | 3,928 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
from sys import stdininput=lambda :stdin.readline()[:-1]mod=10**9+7table_size=2*10**5fac=[1]*(table_size+1)finv=[1]*(table_size+1)for i in range(2,table_size+1):fac[i]=fac[i-1]*i%modfinv[table_size]=pow(fac[table_size],mod-2,mod)for i in range(table_size-1,-1,-1):finv[i]=finv[i+1]*(i+1)%moddef binom(n,k):if n<0 or k<0:return 0if k>n:return 0global table_size,fac,finvif n>table_size:fac+=[0]*(n-table_size)finv+=[0]*(n-table_size)for i in range(table_size+1,n+1):fac[i]=fac[i-1]*i%modfinv[n]=inv(fac[n])for i in range(n-1,table_size,-1):finv[i]=finv[i+1]*(i+1)%modtable_size=nreturn (fac[n]*finv[k]%mod)*finv[n-k]%moddef fpow(x,k):res=1while k:if k&1:res=res*x%modx=x*x%modk>>=1return resdef inv(a):if a<table_size:return fac[a-1]*finv[a]%modreturn fpow(a,mod-2)def modint_to_frac(a):a%=modif a==0:return '0/1'for X in range(1,10000):for Y in range(1,X+1):Z=Y*inv(X)%modif a==Z:return str(Y)+'/'+str(X)if a==mod-Z:return '-'+str(Y)+'/'+str(X)Z=X*inv(Y)%modif a==Z:return str(X)+'/'+str(Y)if a==mod-Z:return '-'+str(X)+'/'+str(Y)n,m=map(int,input().split())ans=0for i in range(m):t,x,y=map(int,input().split())if t==1:ans-=binom(x+y,x)*binom(2*n-(x+y+1),n-(x+1))else:ans-=binom(x+y,x)*binom(2*n-(x+y+1),n-(y+1))ans%=modans+=binom(2*n,n)*2*nprint(ans%mod)