結果
問題 | No.794 チーム戦 (2) |
ユーザー | やむなく |
提出日時 | 2019-02-22 22:11:36 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 100 ms / 1,500 ms |
コード長 | 1,386 bytes |
コンパイル時間 | 1,229 ms |
コンパイル使用メモリ | 165,240 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-11-25 08:49:41 |
合計ジャッジ時間 | 3,758 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 32 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define rep(i,n) for(int i=0;i<(int)(n);i++) #define repe(i,n) rep(i,(n)+1) #define per(i,n) for(int i=(int)(n)-1;i>=0;i--) #define pere(i,n) rep(i,(n)+1) #define all(x) (x).begin(),(x).end() #define SP <<" "<< #define MOD 1000000007 #define IINF 1000000000 #define LINF 1000000000000000000 typedef long long LL; typedef long double LD; long long modpow(long long a, long long n, long long mod) { long long res = 1; while (n > 0) { if (n & 1) res = res * a % mod; a = a * a % mod; n >>= 1; } return res; } // a^{-1} mod を計算する long long modinv(long long a, long long mod) { return modpow(a, mod - 2, mod); } int main(){ LL n,k; cin >> n >> k; vector<int> a(n); vector<LL> fact(n+1); fact[0]=1; for(int i=1;i<=n;i++){ fact[i]=fact[i-1]*i%MOD; } for(int i=0;i<n;i++) cin >> a[i]; sort(all(a)); LL ans=1; for(int i=n-1;i>=0;i--){ LL idx=distance(a.begin(),upper_bound(all(a),k-a[i])); if(idx-(n-1-i)<=0){ cout << 0 << endl; return 0; } if(idx>i){ LL nokori=n-2*(n-1-i); LL tmp=nokori/2; while(nokori){ ans=(ans*(nokori*(nokori-1)/2%MOD))%MOD; nokori-=2; } ans=ans*modinv(fact[tmp],MOD)%MOD; break; } ans=(ans*(idx-(n-1-i)))%MOD; } cout << ans << endl; return 0; }