結果

問題 No.1596 Distance Sum in 2D Plane
コンテスト
ユーザー irona
提出日時 2021-07-10 09:44:06
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 455 ms / 2,000 ms
コード長 1,908 bytes
コンパイル時間 4,025 ms
コンパイル使用メモリ 230,548 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-02 01:46:50
合計ジャッジ時間 10,623 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
#include <atcoder/all>
#define ll long long int
#define vec vector<ll>
#define mat vector<vector<ll>>
#define pll pair<ll,ll>
#define count __builtin_popcountll

using namespace std;
using namespace atcoder;
using mint = modint1000000007;

const ll mod=1000000007;//998244353;
const ll inf=1000000000000000000;
ll dx4[4]={1,0,-1,0};
ll dy4[4]={0,-1,0,1};
ll dx8[8]={1,0,-1,1,-1,1,0,-1};
ll dy8[8]={1,1,1,0,0,-1,-1,-1};

ll powmod(ll a,ll n,ll m) {
    if(n == 0)
        return 1;

    if(n % 2 ==0){
        ll r = powmod(a,n/2,m);
        return r*r % m;
    }

    return a*powmod(a,n-1,m)%m;
}

ll factmod(ll a,ll b,ll m) {
  if(a<b)return 0;
  ll s=1;
  while(a>=b){
    s = s*a % m;
    a -= 1;
  }
  return s;
}

ll modinv(ll a, ll m) {
  ll b = m, u = 1, v = 0;
  while (b) {
    ll t = a / b;
    a -= t * b; swap(a, b);
    u -= t * v; swap(u, v);
  }
  u %= m;
  if (u < 0) u += m;
  return u;
}

ll combmod(ll a,ll b,ll m){
  if(b==0)return 1;
  ll c=factmod(a,a-b+1,m);
  ll d=modinv(factmod(b,1,m),m);
  ll e=(c*d)%m;
  return e;
}

void modsum(ll &a,ll b,ll m){
  a+=b;
  a%=m;
}

void modsub(ll &a,ll b,ll m){
  a-=b;
  a+=m;
  a%=m;
}

void modmul(ll &a,ll b,ll m){
  a*=b;
  a%=m;
}

void modquo(ll &a,ll b,ll m){
  a*=modinv(b,m);
  a%=m;
}

int main(){
  cout << fixed << setprecision(15);

  ll n,m;
  cin >> n >> m;

  mint ans=0;
  ans=combmod(2*n,n,mod);
  ans*=(2*n);

  vector<ll> f(2*n+1);
  f[0]=1;
  for(ll i=1;i<2*n+1;i++){
    f[i]=f[i-1]*i%mod;
  }


  for(ll i=0;i<m;i++){
    ll t,x,y;
    cin >> t >> x >> y;
    mint ret=1;
    ret*=f[x+y];
    ret*=modinv(f[x],mod);
    ret*=modinv(f[y],mod);
    ret*=f[2*n-x-y-1];    
    if(t==1){
      ret*=modinv(f[n-x-1],mod);
      ret*=modinv(f[n-y],mod);
    }
    else{
      ret*=modinv(f[n-x],mod);
      ret*=modinv(f[n-y-1],mod);
    }
    ans-=ret;
  }

  cout << ans.val() << endl;

  return 0;
}
0