結果
問題 |
No.1084 積の積
|
ユーザー |
|
提出日時 | 2020-12-11 19:03:48 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 31 ms / 2,000 ms |
コード長 | 3,710 bytes |
コンパイル時間 | 1,342 ms |
コンパイル使用メモリ | 131,236 KB |
最終ジャッジ日時 | 2025-01-16 22:19:54 |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 5 |
other | AC * 27 |
ソースコード
#include <iostream> #include <vector> #include <algorithm> #include <cmath> #include <queue> #include <string> #include <map> #include <set> #include <stack> #include <tuple> #include <deque> #include <array> #include <numeric> #include <bitset> #include <iomanip> #include <cassert> #include <chrono> #include <random> #include <limits> #include <iterator> #include <functional> #include <sstream> #include <fstream> #include <complex> #include <cstring> #include <unordered_map> using namespace std; using ll = long long; using P = pair<int, int>; constexpr int INF = 1001001001; constexpr int mod = 1000000007; // constexpr int mod = 998244353; template<class T> inline bool chmax(T& x, T y){ if(x < y){ x = y; return true; } return false; } template<class T> inline bool chmin(T& x, T y){ if(x > y){ x = y; return true; } return false; } struct mint { int x; mint() : x(0) {} mint(int64_t y) : x(y >= 0 ? y % mod : (mod - (-y) % mod) % mod) {} mint& operator+=(const mint& p){ if((x += p.x) >= mod) x -= mod; return *this; } mint& operator-=(const mint& p){ if((x -= p.x) < 0) x += mod; return *this; } mint& operator*=(const mint& p){ x = (int)(1LL * x * p.x % mod); return *this; } mint& operator/=(const mint& p){ *this *= p.inverse(); return *this; } mint operator-() const { return mint(-x); } mint operator+(const mint& p) const { return mint(*this) += p; } mint operator-(const mint& p) const { return mint(*this) -= p; } mint operator*(const mint& p) const { return mint(*this) *= p; } mint operator/(const mint& p) const { return mint(*this) /= p; } bool operator==(const mint& p) const { return x == p.x; } bool operator!=(const mint& p) const { return x != p.x; } mint pow(int64_t n) const { mint res = 1, mul = x; while(n > 0){ if(n & 1) res *= mul; mul *= mul; n >>= 1; } return res; } mint inverse() const { return pow(mod - 2); } friend ostream& operator<<(ostream& os, const mint& p){ return os << p.x; } friend istream& operator>>(istream& is, mint& p){ int64_t val; is >> val; p = mint(val); return is; } }; int main(){ ios::sync_with_stdio(false); cin.tie(nullptr); int N; cin >> N; vector<int> A(N); for(int i = 0; i < N; ++i) cin >> A[i]; constexpr int bound = 1e+9; vector<ll> R(N); int l = 0, r = 0; ll prod = 1; while(r < N){ while(l < r && prod * A[r] >= bound){ R[l] = r; prod /= A[l++]; } prod *= A[r++]; } while(l < r) R[l++] = r; vector<ll> coeff_1(N + 1), coeff_0(N + 1); // (R - L) - (i - L) = -i + R for(int i = 0; i < N; ++i){ coeff_1[i] += -1, coeff_1[R[i]] -= -1; coeff_0[i] += R[i], coeff_0[R[i]] -= R[i]; } for(int i = 0; i < N; ++i){ coeff_1[i + 1] += coeff_1[i]; coeff_0[i + 1] += coeff_0[i]; } vector<ll> W(N + 1); for(int i = 0; i < N; ++i){ W[i] = coeff_1[i] * i + coeff_0[i]; } // for(int i = 0; i < N; ++i) cout << R[i] << ' '; // cout << endl; // for(int i = 0; i < N; ++i) cout << coeff_1[i] << ' '; // cout << endl; // for(int i = 0; i < N; ++i) cout << coeff_0[i] << ' '; // cout << endl; // for(int i = 0; i < N; ++i) cout << W[i] << ' '; // cout << endl; mint ans = 1; for(int i = 0; i < N; ++i){ ans *= mint(A[i]).pow(W[i]); } cout << ans << '\n'; return 0; }