結果
| 問題 |
No.1102 Remnants
|
| コンテスト | |
| ユーザー |
chocono2230
|
| 提出日時 | 2020-07-04 11:26:33 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 129 ms / 2,000 ms |
| コード長 | 2,118 bytes |
| コンパイル時間 | 1,709 ms |
| コンパイル使用メモリ | 169,368 KB |
| 実行使用メモリ | 6,528 KB |
| 最終ジャッジ日時 | 2024-09-18 20:02:22 |
| 合計ジャッジ時間 | 4,376 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 25 |
ソースコード
#include <bits/stdc++.h>
#define rep(i,n) for(int i = 0; i < (int)(n); i++)
#define rrep(ri,n) for(int ri = (int)(n-1); ri >= 0; ri--)
#define rep2(i,x,n) for(int i = (int)(x); i < (int)(n); i++)
#define rrep2(ri,x,n) for(int ri = (int)(n-1); ri >= (int)(x); ri--)
#define repit(itr,x) for(auto itr = x.begin(); itr != x.end(); itr++)
#define rrepit(ritr,x) for(auto ritr = x.rbegin(); ritr != x.rend(); ritr++)
#define ALL(x) x.begin(), x.end()
using ll = long long;
using namespace std;
// auto mod int
// https://youtu.be/L8grWxBlIZ4?t=9858
// https://youtu.be/ERZuLAxZffQ?t=4807 : optimize
// https://youtu.be/8uowVvQ_-Mo?t=1329 : division
const int mod = 1000000007;
struct mint {
ll x; // typedef long long ll;
mint(ll x=0):x((x%mod+mod)%mod){}
mint operator-() const { return mint(-x);}
mint& operator+=(const mint a) {
if ((x += a.x) >= mod) x -= mod;
return *this;
}
mint& operator-=(const mint a) {
if ((x += mod-a.x) >= mod) x -= mod;
return *this;
}
mint& operator*=(const mint a) {
(x *= a.x) %= mod;
return *this;
}
mint operator+(const mint a) const {
mint res(*this);
return res+=a;
}
mint operator-(const mint a) const {
mint res(*this);
return res-=a;
}
mint operator*(const mint a) const {
mint res(*this);
return res*=a;
}
mint pow(ll t) const {
if (!t) return 1;
mint a = pow(t>>1);
a *= a;
if (t&1) a *= *this;
return a;
}
// for prime mod
mint inv() const {
return pow(mod-2);
}
mint& operator/=(const mint a) {
return (*this) *= a.inv();
}
mint operator/(const mint a) const {
mint res(*this);
return res/=a;
}
};
void fc(vector<mint> &comb, int n, int k){
comb.at(0) = 1;
rep2(i, 1, n){
comb.at(i) = comb.at(i-1) * (k+i) / i;
}
}
int main(){
int n, k;
cin >> n >> k;
vector<mint> a(n);
rep(i, n){
int in;
cin >> in;
a.at(i) = in;
}
vector<mint> comb(n);
fc(comb, n, k);
mint ans = 0;
rep(i, n){
mint add = a.at(i);
add *= comb.at(i) * comb.at(n-1-i);
ans += add;
}
cout << ans.x << endl;
return 0;
}
chocono2230