結果
問題 | No.2615 ペアの作り方 |
ユーザー |
|
提出日時 | 2024-01-26 21:27:13 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 64 ms / 2,000 ms |
コード長 | 3,442 bytes |
コンパイル時間 | 2,639 ms |
コンパイル使用メモリ | 221,196 KB |
最終ジャッジ日時 | 2025-02-18 22:55:55 |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 21 |
ソースコード
#pragma GCC optimize("Ofast")#include <bits/stdc++.h>using namespace std;typedef long long int ll;typedef unsigned long long int ull;mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());ll myRand(ll B) {return (ull)rng() % B;}inline double time() {return static_cast<long double>(chrono::duration_cast<chrono::nanoseconds>(chrono::steady_clock::now().time_since_epoch()).count()) * 1e-9;}template <int mod>struct static_modint {using mint = static_modint;int x;static_modint() : x(0) {}static_modint(int64_t y) : x(y >= 0 ? y % mod : (mod - (-y) % mod) % mod) {}mint& operator+=(const mint& rhs) {if ((x += rhs.x) >= mod) x -= mod;return *this;}mint& operator-=(const mint& rhs) {if ((x += mod - rhs.x) >= mod) x -= mod;return *this;}mint& operator*=(const mint& rhs) {x = (int) (1LL * x * rhs.x % mod);return *this;}mint& operator/=(const mint& rhs) { return *this = *this * rhs.inv(); }mint pow(long long n) const {mint _x = *this, r = 1;while (n) {if (n & 1) r *= _x;_x *= _x;n >>= 1;}return r;}mint inv() const { return pow(mod - 2); }mint operator+() const { return *this; }mint operator-() const { return mint() - *this; }friend mint operator+(const mint& lhs, const mint& rhs) {return mint(lhs) += rhs;}friend mint operator-(const mint& lhs, const mint& rhs) {return mint(lhs) -= rhs;}friend mint operator*(const mint& lhs, const mint& rhs) {return mint(lhs) *= rhs;}friend mint operator/(const mint& lhs, const mint& rhs) {return mint(lhs) /= rhs;}friend bool operator==(const mint& lhs, const mint& rhs) {return lhs.x == rhs.x;}friend bool operator!=(const mint& lhs, const mint& rhs) {return lhs.x != rhs.x;}friend ostream &operator<<(ostream &os, const mint &p) {return os << p.x;}friend istream &operator>>(istream &is, mint &a) {int64_t t; is >> t;a = static_modint<mod>(t);return (is);}};const unsigned int mod = 998244353;using modint = static_modint<mod>;modint mod_pow(ll n, ll x) { return modint(n).pow(x); }modint mod_pow(modint n, ll x) { return n.pow(x); }template <typename T>struct Comination {vector<T> p, invp;Comination(int sz) : p(sz+1), invp(sz+1) {p[0] = 1;for (int i = 1; i <= sz; ++i) {p[i] = p[i-1] * i;}invp[sz] = p[sz].inv();for (int i = sz-1; i >= 0; --i) {invp[i] = invp[i+1] * (i+1);}}T comb(int n, int r) {if (r < 0 or n < r) return 0;return p[n]*invp[n-r]*invp[r];}T big_comb(T n, int r) {T res = invp[r];for (int i = 0; i < r; ++i) {res *= (n-i);}return res;}};using Comb = Comination<modint>;Comb p(1<<20);int main(){cin.tie(nullptr);ios::sync_with_stdio(false);int n; cin >> n;vector<pair<int,int>> v;for (int i = 0; i < 2*n; ++i) {int a; cin >> a;v.push_back({a, i/n});}sort(v.begin(), v.end());int cnt = 0;for (int i = 0; i < n; ++i) {cnt += v[i].second;}cout << p.p[cnt]*p.p[n-cnt] << endl;}