結果
問題 | No.794 チーム戦 (2) |
ユーザー |
![]() |
提出日時 | 2023-04-09 18:17:08 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 63 ms / 1,500 ms |
コード長 | 2,284 bytes |
コンパイル時間 | 3,900 ms |
コンパイル使用メモリ | 258,204 KB |
最終ジャッジ日時 | 2025-02-12 04:17:33 |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 32 |
ソースコード
#include <bits/stdc++.h>#include <atcoder/all>using namespace std;using namespace atcoder;using mint = modint1000000007;const int mod = 1000000007;//using mint = modint998244353;//const int mod = 998244353;//const int INF = 1e9;//const long long LINF = 1e18;#define rep(i, n) for (int i = 0; i < (n); ++i)#define rep2(i,l,r)for(int i=(l);i<(r);++i)#define rrep(i, n) for (int i = (n-1); i >= 0; --i)#define rrep2(i,l,r)for(int i=(r-1);i>=(l);--i)#define all(x) (x).begin(),(x).end()#define allR(x) (x).rbegin(),(x).rend()#define endl "\n"#define P pair<int,int>template<typename A, typename B> inline bool chmax(A & a, const B & b) { if (a < b) { a = b; return true; } return false; }template<typename A, typename B> inline bool chmin(A & a, const B & b) { if (a > b) { a = b; return true; } return false; }// combination mod prime// https://www.youtube.com/watch?v=8uowVvQ_-Mo&feature=youtu.be&t=1619struct combination {vector<mint> fact, ifact;combination(int n) :fact(n + 1), ifact(n + 1) {assert(n < mod);fact[0] = 1;for (int i = 1; i <= n; ++i) fact[i] = fact[i - 1] * i;ifact[n] = fact[n].inv();for (int i = n; i >= 1; --i) ifact[i - 1] = ifact[i] * i;}mint operator()(int n, int k) { return com(n, k); }mint com(int n, int k) {//負の二項係数を考慮する場合にコメントアウトを外す//if (n < 0) return com(-n, k) * (k % 2 ? -1 : 1);if (k < 0 || k > n) return 0;return fact[n] * ifact[k] * ifact[n - k];}mint inv(int n, int k) {//if (n < 0) return inv(-n, k) * (k % 2 ? -1 : 1);if (k < 0 || k > n) return 0;return ifact[n] * fact[k] * fact[n - k];}mint p(int n, int k) { return fact[n] * ifact[n - k]; }}c(2000006);int main() {ios::sync_with_stdio(false);cin.tie(nullptr);int n, k; cin >> n >> k;vector<int>v0, v1;int x = k / 2 + 1;rep(i, n) {int a; cin >> a;if (a >= x)v0.push_back(a);else v1.push_back(a);}sort(allR(v0));sort(all(v1));mint ans = 1;rep(i, v0.size()) {int x = k - v0[i];int count = upper_bound(all(v1), x) - v1.begin();ans *= count - i;}int diff = v1.size() - v0.size();if (diff >= 0) {ans *= c.fact[diff];ans /= pow_mod(2, diff / 2, mod);ans *= c.ifact[diff / 2];}cout << ans.val() << endl;return 0;}