結果
問題 | No.1596 Distance Sum in 2D Plane |
ユーザー |
![]() |
提出日時 | 2021-07-09 21:53:34 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 197 ms / 2,000 ms |
コード長 | 1,501 bytes |
コンパイル時間 | 2,112 ms |
コンパイル使用メモリ | 194,280 KB |
最終ジャッジ日時 | 2025-01-22 21:36:02 |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
#include <bits/stdc++.h> #define rep(i, n) for (int i = 0; i < (n); ++i) #define srep(i, s, t) for (int i = s; i < t; ++i) #define drep(i, n) for (int i = (n)-1; i >= 0; --i) using namespace std; typedef long long int ll; typedef pair<int, int> P; #define yn \ { puts("YES"); } \ else { \ puts("NO"); \ } const int MAX = 510000; const int MOD = 1000000007; long long fac[MAX], finv[MAX], inv[MAX]; // テーブルを作る前処理 void COMinit() { fac[0] = fac[1] = 1; finv[0] = finv[1] = 1; inv[1] = 1; for (int i = 2; i < MAX; i++) { fac[i] = fac[i - 1] * i % MOD; inv[i] = MOD - inv[MOD % i] * (MOD / i) % MOD; finv[i] = finv[i - 1] * inv[i] % MOD; } } // 二項係数計算 long long COM(int n, int k) { if (n < k) return 0; if (n < 0 || k < 0) return 0; return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD; } long long FINV(int n) { if (n < 0) return 0; return finv[n]; } int main() { COMinit(); ll n, m; cin >> n >> m; int t[m], x[m], y[m]; rep(i, m) cin >> t[i] >> x[i] >> y[i]; ll ans = COM(n * 2, n) * n * 2 % MOD; rep(i, m) { if (t[i] == 1) { ans += MOD - COM(x[i] + y[i], x[i]) * COM(n - x[i] - 1 + n - y[i], n - y[i]) % MOD; } else { ans += MOD - COM(x[i] + y[i], x[i]) * COM(n - y[i] - 1 + n - x[i], n - x[i]) % MOD; } ans %= MOD; } cout << ans << endl; return 0; }