結果
問題 | No.1596 Distance Sum in 2D Plane |
ユーザー |
|
提出日時 | 2021-07-09 21:36:56 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 202 ms / 2,000 ms |
コード長 | 1,556 bytes |
コンパイル時間 | 954 ms |
コンパイル使用メモリ | 110,628 KB |
実行使用メモリ | 15,744 KB |
最終ジャッジ日時 | 2024-07-01 15:28:13 |
合計ジャッジ時間 | 5,584 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
#include "iostream" #include "climits" #include "list" #include "queue" #include "stack" #include "set" #include "functional" #include "algorithm" #include "string" #include "map" #include "unordered_map" #include "unordered_set" #include "iomanip" #include "cmath" #include "random" #include "bitset" #include "cstdio" #include "numeric" #include "cassert" #include "ctime" using namespace std; constexpr long long int MOD = 1000000007; //constexpr int MOD = 1000000007; //constexpr int MOD = 998244353; //constexpr long long int MOD = 998244353; constexpr double EPS = 1e-8; //int N, M, K, T, H, W, L, R; long long int N, M, K, T, H, W, L, R; long long by[777777]; long long inv[777777]; long long int power(long long int x, long long int n, long long int M) { long long int ret = 1; long long int by = x; while (n) { if (n & 1) { ret *= by; ret %= M; } by *= by; by %= M; n >>= 1; } return ret; } long long ncr(long long n, long long r) { if (n < r || r < 0)return 0; if (n < 0)return 0; return by[n] * inv[r] % MOD * inv[n - r] % MOD; } int main() { ios::sync_with_stdio(false); cin.tie(0); cin >> N >> M; by[0] = inv[0] = 1; for (int i = 1; i < 777777; i++) { by[i] = by[i - 1] * i % MOD; inv[i] = power(by[i], MOD - 2, MOD); } long long ans = N * 2 * ncr(N * 2, N) % MOD; for (int i = 0; i < M; i++) { int a, b, c; cin >> a >> b >> c; long long add = ncr(b + c, b); if (a == 1)b++; else c++; add *= ncr(N * 2 - (b + c), N - b); ans += MOD - add % MOD; ans %= MOD; } cout << ans << endl; }