結果
問題 | No.995 タピオカオイシクナーレ |
ユーザー |
![]() |
提出日時 | 2020-02-21 21:47:30 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 242 ms / 2,000 ms |
コード長 | 1,205 bytes |
コンパイル時間 | 1,594 ms |
コンパイル使用メモリ | 169,576 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-10-08 22:39:56 |
合計ジャッジ時間 | 4,876 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 23 |
ソースコード
#include <bits/stdc++.h>using namespace std;typedef long long ll;typedef pair<ll, ll> P;#define DUMP(x) cout << #x << " = " << (x) << endl;#define FOR(i, m, n) for (ll i = m; i < n; i++)#define IFOR(i, m, n) for (ll i = n - 1; i >= m; i--)#define REP(i, n) FOR(i, 0, n)#define IREP(i, n) IFOR(i, 0, n)#define FOREACH(x, a) for (auto&(x) : (a))#define ALL(v) (v).begin(), (v).end()#define SZ(x) ll(x.size())const int MOD = 1000000007;ll modpow(ll a, ll b, ll m) {if (b == 0) return 1;ll ret;if (b % 2 == 0) {ll t = modpow(a, b / 2, m) % m;ret = t * t % m;} else {ll t = modpow(a, (b - 1) / 2, m) % m;ret = a * t % m * t % m;}ret %= m;return ret;}ll inv(ll a) { return modpow(a, MOD - 2, MOD); }int main() {ll N, M, K, p, q;cin >> N >> M >> K >> p >> q;M--;vector<ll> b(N);REP(i, N) { cin >> b[i]; }ll r = p * inv(q) % MOD;vector<ll> pr(N);REP(i, N) {ll t = 2 * ((MOD - r) % MOD) + 1 % MOD;ll p0 = i <= M ? 1 : 0;pr[i] = (modpow(t, K, MOD) * (p0 + ((MOD - inv(2)) % MOD)) % MOD) % MOD +inv(2) % MOD;}ll ans = 0;REP(i, N) { (ans += b[i] * pr[i] % MOD) %= MOD; }cout << ans << endl;}