結果
問題 | No.856 増える演算 |
ユーザー |
![]() |
提出日時 | 2019-03-17 19:57:28 |
言語 | C++11 (gcc 13.3.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 2,416 bytes |
コンパイル時間 | 1,300 ms |
コンパイル使用メモリ | 162,456 KB |
実行使用メモリ | 12,068 KB |
最終ジャッジ日時 | 2024-10-09 16:00:12 |
合計ジャッジ時間 | 9,347 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 TLE * 1 -- * 49 |
ソースコード
// TLE解,O(N^2)// 「最も小さいもの」はO(N)で求められる.#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;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;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 solve() {int N;cin >> N;vector<ll> A(N);for (int i = 0; i < N; i++) {cin >> A[i];A[i] %= mod;}// 考察1. miをO(N)で求める// iを固定すると,jとしては(i<jの中で)もちろんA[j]が小さいのを選ぶべき.// i以降の要素で最も小さい要素をO(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]);}}ll prod = 1;for (int i = 0; i < N; i++) {for (int j = i + 1; j < N; j++) {prod *= f(A[i], A[j]);prod %= mod;}}cout << (prod * inv(mi_f)) % mod << endl;}int main() {#ifdef LOCAL_ENVcin.exceptions(ios::failbit);#endifcin.tie(0);ios::sync_with_stdio(false);cout.setf(ios::fixed);cout.precision(16);solve();}