結果
| 問題 |
No.1102 Remnants
|
| コンテスト | |
| ユーザー |
🍮かんプリン
|
| 提出日時 | 2020-07-03 22:24:21 |
| 言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 2,766 bytes |
| コンパイル時間 | 2,074 ms |
| コンパイル使用メモリ | 159,524 KB |
| 実行使用メモリ | 6,948 KB |
| 最終ジャッジ日時 | 2024-09-17 03:57:02 |
| 合計ジャッジ時間 | 9,096 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | WA * 3 |
| other | WA * 25 |
ソースコード
/**
* @FileName a.cpp
* @Author kanpurin
* @Created 2020.07.03 22:24:14
**/
#include "bits/stdc++.h"
using namespace std;
typedef long long ll;
int MOD = 1e9 + 7;
struct mint {
private:
long long x;
public:
mint(long long x = 0) : x((MOD + x) % MOD) {}
mint(std::string& s) {
long long z = 0;
for (int i = 0; i < s.size(); i++) {
z *= 10;
z += s[i] - '0';
z %= MOD;
}
this->x = z;
}
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) {
long long n = MOD - 2;
mint u = 1, b = a;
while (n > 0) {
if (n & 1) {
u *= b;
}
b *= b;
n >>= 1;
}
return *this *= u;
}
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 operator/(const mint& a) const {
mint res(*this);
return res /= a;
}
friend std::ostream& operator<<(std::ostream& os, const mint& n) {
return os << n.x;
}
bool operator==(const mint& a) const {
return this->x == a.x;
}
};
long long powMod(long long k, long long n, long long mod) {
long long x = 1;
while (n > 0) {
if (n & 1) {
x = x * k % mod;
}
k = k * k % mod;
n >>= 1;
}
return x;
}
long long comb(ll n, ll r, ll mod) {
long long ret = 1;
while (true) {
if (r == 0) break;
long long N = n % mod;
long long R = r % mod;
if (N < R) return 0;
for (int i = 0; i < R; i++) {
ret = ret * (N - i) % mod;
}
long long imul = 1;
for (int i = 0; i < R; i++) {
imul = imul * (i + 1) % mod;
}
ret = ret * powMod(imul, mod - 2, mod) % mod;
n /= mod;
r /= mod;
}
return ret;
}
int main() {
int n, k;
cin >> n >> k;
mint ans = 0;
mint c1 = comb(n + k - 1, n - 1, MOD), c2 = 1;
for (int i = 0; i < n; i++) {
int a;
cin >> a;
ans += c1 * c2 * a;
cout << ans << endl;
c1 *= n - i - 1;
c1 /= n - i + k - 1;
c2 *= k + i + 1;
c2 /= i + 1;
}
cout << ans << endl;
return 0;
}
🍮かんプリン