結果
| 問題 | 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;
}
            
            
            
        