結果
問題 | No.1596 Distance Sum in 2D Plane |
ユーザー | 沙耶花 |
提出日時 | 2021-07-09 21:35:39 |
言語 | C++17 (gcc 13.2.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 78 ms / 2,000 ms |
コード長 | 1,161 bytes |
コンパイル時間 | 3,855 ms |
コンパイル使用メモリ | 265,220 KB |
実行使用メモリ | 11,800 KB |
最終ジャッジ日時 | 2023-09-14 07:58:39 |
合計ジャッジ時間 | 6,199 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge13 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 19 ms
11,512 KB |
testcase_01 | AC | 19 ms
11,600 KB |
testcase_02 | AC | 77 ms
11,652 KB |
testcase_03 | AC | 78 ms
11,608 KB |
testcase_04 | AC | 76 ms
11,596 KB |
testcase_05 | AC | 74 ms
11,520 KB |
testcase_06 | AC | 75 ms
11,628 KB |
testcase_07 | AC | 73 ms
11,584 KB |
testcase_08 | AC | 76 ms
11,636 KB |
testcase_09 | AC | 75 ms
11,628 KB |
testcase_10 | AC | 77 ms
11,596 KB |
testcase_11 | AC | 67 ms
11,632 KB |
testcase_12 | AC | 65 ms
11,536 KB |
testcase_13 | AC | 67 ms
11,548 KB |
testcase_14 | AC | 19 ms
11,472 KB |
testcase_15 | AC | 20 ms
11,612 KB |
testcase_16 | AC | 19 ms
11,672 KB |
testcase_17 | AC | 20 ms
11,668 KB |
testcase_18 | AC | 20 ms
11,544 KB |
testcase_19 | AC | 18 ms
11,800 KB |
ソースコード
#include <stdio.h> #include <bits/stdc++.h> #include <atcoder/all> using namespace atcoder; using mint = modint1000000007; using namespace std; #define rep(i,n) for (int i = 0; i < (n); ++i) #define Inf 1000000001 struct combi{ deque<mint> kaijou; deque<mint> kaijou_; combi(int n){ kaijou.push_back(1); for(int i=1;i<=n;i++){ kaijou.push_back(kaijou[i-1]*i); } mint b=kaijou[n].inv(); kaijou_.push_front(b); for(int i=1;i<=n;i++){ int k=n+1-i; kaijou_.push_front(kaijou_[0]*k); } } mint combination(int n,int r){ if(r>n)return 0; mint a = kaijou[n]*kaijou_[r]; a *= kaijou_[n-r]; return a; } mint junretsu(int a,int b){ mint x = kaijou_[a]*kaijou_[b]; x *= kaijou[a+b]; return x; } mint catalan(int n){ return combination(2*n,n)/(n+1); } }; int main(){ combi C(1000000); int N,M; cin>>N>>M; mint ans = C.junretsu(N,N) * (N*2); rep(_,M){ int t,x,y; scanf("%d %d %d",&t,&x,&y); mint cur = C.junretsu(x,y); if(t==1){ cur *= C.junretsu((N-(x+1)),(N-y)); } else{ cur *= C.junretsu((N-x),(N-(y+1))); } ans -= cur; } cout<<ans.val()<<endl; return 0; }