結果

問題 No.1596 Distance Sum in 2D Plane
ユーザー kanzakihikarukanzakihikaru
提出日時 2021-07-09 22:59:06
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 86 ms / 2,000 ms
コード長 2,144 bytes
コンパイル時間 1,591 ms
コンパイル使用メモリ 167,808 KB
実行使用メモリ 12,560 KB
最終ジャッジ日時 2023-09-14 10:32:05
合計ジャッジ時間 4,210 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 16 ms
12,452 KB
testcase_01 AC 15 ms
12,172 KB
testcase_02 AC 85 ms
12,312 KB
testcase_03 AC 85 ms
12,232 KB
testcase_04 AC 86 ms
12,308 KB
testcase_05 AC 82 ms
12,560 KB
testcase_06 AC 82 ms
12,444 KB
testcase_07 AC 84 ms
12,316 KB
testcase_08 AC 82 ms
12,444 KB
testcase_09 AC 85 ms
12,340 KB
testcase_10 AC 83 ms
12,408 KB
testcase_11 AC 64 ms
12,304 KB
testcase_12 AC 65 ms
12,344 KB
testcase_13 AC 63 ms
12,244 KB
testcase_14 AC 17 ms
12,548 KB
testcase_15 AC 17 ms
12,416 KB
testcase_16 AC 17 ms
12,508 KB
testcase_17 AC 16 ms
12,456 KB
testcase_18 AC 16 ms
12,312 KB
testcase_19 AC 16 ms
12,240 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <string>

#define rep(i,n) for(ll i=0;i < (n);i++)
#define rep2(i, s, n) for (ll i = (s); i < (n); i++)
#define rrep(i,n) for(ll i=n-1;i>=0;i--)
#define ff first
#define ss second
#define pb push_back
#define ALL(a) (a).begin(),(a).end()

typedef long long ll;

const ll INF=1000000000000000000ll;
const ll MOD=1000000007ll;
const int MAX=5100000;

using namespace std;

// aよりもbが大きいならばaをbで更新する
// (更新されたならばtrueを返す)
template <typename T>
bool chmax(T &a, const T& b) {
  if (a < b) {
    a = b;  // aをbで更新
    return true;
  }
  return false;
}
// aよりもbが小さいならばaをbで更新する
// (更新されたならばtrueを返す)
template <typename T>
bool chmin(T &a, const T& b) {
  if (a > b) {
    a = b;  // aをbで更新
    return true;
  }
  return false;
}

vector<long long> fact, fact_inv, inv;
/*  init_nCk :二項係数のための前処理
    計算量:O(n)
*/
void init_nCk(int SIZE) {
    fact.resize(SIZE + 5);
    fact_inv.resize(SIZE + 5);
    inv.resize(SIZE + 5);
    fact[0] = fact[1] = 1;
    fact_inv[0] = fact_inv[1] = 1;
    inv[1] = 1;
    for (int i = 2; i < SIZE + 5; i++) {
        fact[i] = fact[i - 1] * i % MOD;
        inv[i] = MOD - inv[MOD % i] * (MOD / i) % MOD;
        fact_inv[i] = fact_inv[i - 1] * inv[i] % MOD;
    }
}
/*  nCk :MODでの二項係数を求める(前処理 int_nCk が必要)
    計算量:O(1)
*/
long long nCk(int n, int k) {
    if(n < k){return 0;};
    if(n < 0 || k < 0){return 0;};
    return fact[n] * (fact_inv[k] * fact_inv[n - k] % MOD) % MOD;
}

int main()
{
	cin.tie(0);	ios::sync_with_stdio(false);
  init_nCk(400000);
  ll N, M; cin >> N >> M;

  ll ans = nCk(N*2, N);
  ans *= N; ans %= MOD;
  ans *= 2; ans %= MOD;

  ll t, x, y, tmp;
  rep(i, M){
    cin >> t >> x >> y;
    if(t == 1){
      tmp = nCk(x+y, x) * nCk(2*N-x-y-1, N-y); tmp %= MOD;
      ans -= tmp; if(ans < 0){ans += MOD;}
    }
    else{
      tmp = nCk(x+y, x) * nCk(2*N-x-y-1, N-x); tmp %= MOD;
      ans -= tmp; if(ans < 0){ans += MOD;}
    }
  }

  cout << ans << endl;

  return 0;
}
0