結果

問題 No.856 増える演算
ユーザー kakira9618
提出日時 2019-03-28 12:19:10
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 5,219 bytes
コンパイル時間 1,937 ms
コンパイル使用メモリ 179,604 KB
実行使用メモリ 30,276 KB
最終ジャッジ日時 2024-10-09 16:03:07
合計ジャッジ時間 29,715 ms
ジャッジサーバーID
(参考情報)
judge3 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 53 WA * 27
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

// O((N + max(A_i))log max(A_i))
// A_i > 50000
#include <bits/stdc++.h>
using namespace std;
#define all(x) (x).begin(),(x).end()
#define rep(i, n) for (int i = 0; i < (n); i++)
#define chmin(x, y) (x) = min((x), (y))
#define chmax(x, y) (x) = max((x), (y))
#define endl "\n"
typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
typedef complex<double> C;
template <typename T> ostream &operator<<(ostream &os, const vector<T> &vec) {os << "["; for (const auto &v : vec) {os << v << ","; } os << "]";
    return os;}
template <typename T, typename U> ostream &operator<<(ostream &os, const pair<T, U> &p) {os << "(" << p.first << ", " << p.second << ")"; return os;}
const int mod = 1e9 + 7;
const int MAX_A = 100000;
ll mod_pow(ll x, ll n, ll mod) {
ll res = 1;
while(n > 0) {
if(n & 1) res = res * x % mod;
x = x * x % mod;
n >>= 1;
}
return res;
}
ll inv(ll x) {
return mod_pow(x, mod - 2, mod);
}
ll f(ll x, ll y) {
return (((x + y) % mod) * mod_pow(x, y, mod)) % mod;
}
// (x + y)x^y < (s + t)s^t
// log(x + y) + ylog(x) < log(s + t) + tlog(s)
double g(ll x, ll y) {
return log(x + y) + y * log(x);
}
void fft(vector<C> &f, int n, int dir) {
if (n == 1) return;
vector<C> f0, f1;
for(int i = 0; i < n; i++) {
if (i % 2 == 0) {
f0.push_back(f[i]);
} else {
f1.push_back(f[i]);
}
}
fft(f0, n / 2, dir);
fft(f1, n / 2, dir);
C zeta = polar(1.0, 2 * M_PI * dir / n);
C now = C(1.0, 0.0);
for(int i = 0; i < n; i++) {
f[i] = f0[i % (n / 2)] + now * f1[i % (n / 2)];
now *= zeta;
}
}
vector<C> comb(vector<C> A, vector<C> B) {
int size = A.size() + B.size();
int n = 1;
while(n < size) n <<= 1;
A.resize(n);
B.resize(n);
fft(A, n, 1);
fft(B, n, 1);
vector<C> ret;
for(int i = 0; i < n; i++) {
ret.push_back(A[i] * B[i]);
}
fft(ret, n, -1);
for(int i = 0; i < n; i++) {
ret[i] /= n;
}
return ret;
}
vector<ll> getPairCnt(vector<ll> &K) {
vector<C> A, B;
for (int i = 0; i < K.size(); i++) {
A.push_back(C(K[i], 0));
B.push_back(C(K[i], 0));
}
vector<ll> ret;
vector<C> ret_ = comb(A, B);
for (int i = 0; i < ret_.size(); i++) {
ret.push_back((ll)(real(ret_[i]) + 0.1));
}
for (int i = 0; i < ret.size(); i++) {
if (i % 2 == 0) ret[i] -= (i / 2 < K.size() ? K[i / 2] : 0);
ret[i] >>= 1;
}
return ret;
}
void solve() {
int N;
cin >> N;
assert(N >= 1 && N <= MAX_A);
vector<ll> A(N), K(MAX_A + 1);
for (int i = 0; i < N; i++) {
cin >> A[i];
assert(A[i] >= 1 && A[i] <= MAX_A);
A[i] %= mod;
K[A[i]]++;
}
// 1. miO(N)
// iji<jA[j]
// iO(1)O(N)
vector<ll> mi_table(N + 1, 1LL << 60);
for (int i = N - 1; i >= 0; i--) {
mi_table[i] = min(mi_table[i + 1], A[i]);
}
double mi = g(A[0], mi_table[1]);
ll mi_f = f(A[0], mi_table[1]);
for (int i = 1; i < N; i++) {
if (mi > g(A[i], mi_table[i + 1])) {
mi = g(A[i], mi_table[i + 1]);
mi_f = f(A[i], mi_table[i + 1]);
}
}
// 2. \prod_i \prod_j f(A[i], A[j]) O(N)
// 2a.
// i A[i]^{\sum_{i<j}{A[j]}}
// \sum_{i<j}{A[j]}
vector<ll> acc_back(N + 1);
for (int i = N - 1; i >= 0; i--) {
acc_back[i] = acc_back[i + 1] + A[i];
acc_back[i] %= mod;
}
ll prod = 1;
for (int i = 0; i < N; i++) {
prod *= mod_pow(A[i], acc_back[i + 1], mod);
prod %= mod;
}
// 2b.
// A[i] + A[j] = k (i < j)(i, j)
// A[i]k
// vector<ll> pair_cnt(2001);
// for (int k = 2; k <= 2000; k++) {
// ll c = 0;
// for(int a = 0; a <= k; a++) {
// int b = k - a;
// c += K[a] * K[b];
// c %= mod;
// }
// if (k % 2 == 0) c -= K[k / 2];
// pair_cnt[k] += c / 2;
// }
// 2b-2.
// FFT
// k
vector<ll> pair_cnt = getPairCnt(K);
for (int k = 0; k < pair_cnt.size(); k++) {
prod *= mod_pow(k, pair_cnt[k], mod);
prod %= mod;
}
cout << (prod * inv(mi_f)) % mod << endl;
}
int main() {
#ifdef LOCAL_ENV
cin.exceptions(ios::failbit);
#endif
cin.tie(0);
ios::sync_with_stdio(false);
cout.setf(ios::fixed);
cout.precision(16);
solve();
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0