結果
| 問題 |
No.995 タピオカオイシクナーレ
|
| コンテスト | |
| ユーザー |
Shibuyap
|
| 提出日時 | 2020-02-21 22:17:04 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 43 ms / 2,000 ms |
| コード長 | 1,899 bytes |
| コンパイル時間 | 1,632 ms |
| コンパイル使用メモリ | 167,200 KB |
| 実行使用メモリ | 5,248 KB |
| 最終ジャッジ日時 | 2024-10-08 22:49:02 |
| 合計ジャッジ時間 | 3,157 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 23 |
ソースコード
#include <bits/stdc++.h>
#define rep(i,n) for(int i = 0; i < (n); ++i)
#define srep(i,s,t) for (int i = s; i < t; ++i)
#define drep(i,n) for(int i = (n)-1; i >= 0; --i)
using namespace std;
typedef long long int ll;
typedef pair<int,int> P;
#define yn {puts("Yes");}else{puts("No");}
#define MAX_N 200005
// ax + by = gcd(a, b) となるような (x, y) を求める
// 多くの場合 a と b は互いに素として ax + by = 1 となる (x, y) を求める
long long extGCD(long long a, long long b, long long &x, long long &y) {
if (b == 0) {
x = 1;
y = 0;
return a;
}
long long d = extGCD(b, a%b, y, x); // 再帰的に解く
y -= a / b * x;
return d;
}
// 負の数にも対応した mod (a = -11 とかでも OK)
inline long long mod(long long a, long long m) {
return (a % m + m) % m;
}
// 逆元計算 (ここでは a と m が互いに素であることが必要)
long long modinv(long long a, long long m) {
long long x, y;
extGCD(a, m, x, y);
return mod(x, m); // 気持ち的には x % m だが、x が負かもしれないので
}
ll mod_pow(ll x, ll n, ll mod){ // x ^ n % mod
ll res = 1;
while(n > 0){
if (n & 1) res = res * x % mod;
x = x * x % mod;
n >>= 1;
}
return res;
}
int main() {
ll n, m, k, p, q;
cin >> n >> m >> k >> p >> q;
ll MOD = 1e9+7;
ll b[n];
ll v = 0;
ll w = 0;
rep(i,n){
cin >> b[i];
if(i < m){
v += b[i];
}else{
w += b[i];
}
}
v %= MOD;
w %= MOD;
ll qn = mod_pow(q,k,MOD);
ll q2p = q - 2*p + MOD * 123456;
q2p %= MOD;
ll q2pn = mod_pow(q2p,k,MOD);
ll ans = ((v * q2pn % MOD + v * qn % MOD + w * qn % MOD - w * q2pn % MOD + MOD * 100) % MOD);
ans *= modinv(2*qn, MOD);
ans %= MOD;
cout << ans << endl;
return 0;
}
Shibuyap