結果
| 問題 |
No.978 Fibonacci Convolution Easy
|
| ユーザー |
kibuna
|
| 提出日時 | 2020-01-31 21:59:52 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 61 ms / 2,000 ms |
| コード長 | 2,221 bytes |
| コンパイル時間 | 1,585 ms |
| コンパイル使用メモリ | 170,172 KB |
| 実行使用メモリ | 65,728 KB |
| 最終ジャッジ日時 | 2024-09-18 20:58:00 |
| 合計ジャッジ時間 | 2,707 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 21 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
using lint = long long;
const lint inf = 1LL << 60;
const lint mod = 1000000007;
struct mint {
lint v;
lint _mod;
mint() : v(0) {}
mint(signed v, lint _mod = mod) : v(v), _mod(_mod) {}
mint(lint t, lint _mod = mod) : _mod(_mod) {
v = t % _mod;
if (v < 0)
v += _mod;
}
mint pow(lint k) {
mint res(1), tmp(v);
while (k) {
if (k & 1)
res *= tmp;
tmp *= tmp;
k >>= 1;
}
return res;
}
static mint add_identity() { return mint(0); }
static mint mul_identity() { return mint(1); }
mint inv() { return pow(_mod - 2); }
mint &operator+=(mint a) {
v += a.v;
if (v >= _mod)
v -= _mod;
return *this;
}
mint &operator-=(mint a) {
v += _mod - a.v;
if (v >= _mod)
v -= _mod;
return *this;
}
mint &operator*=(mint a) {
v = v * a.v % _mod;
return *this;
}
mint &operator/=(mint a) { return (*this) *= a.inv(); }
mint operator+(mint a) const { return mint(v) += a; };
mint operator-(mint a) const { return mint(v) -= a; };
mint operator*(mint a) const { return mint(v) *= a; };
mint operator/(mint a) const { return mint(v) /= a; };
mint operator-() const { return v ? mint(_mod - v) : mint(v); }
bool operator==(const mint a) const { return v == a.v; }
bool operator!=(const mint a) const { return v != a.v; }
bool operator<(const mint a) const { return v < a.v; }
};
ostream &operator<<(ostream &os, mint m) { return os << m.v; }
int main() {
cin.tie(nullptr);
ios::sync_with_stdio(false);
lint n, p;
cin >> n >> p;
if (n == 1) {
cout << 0 << "\n";
return 0;
}
vector<mint> fib(n), cum(n);
fib[0] = 0;
fib[1] = 1;
for (int i = 2; i < n; ++i) {
fib[i] = fib[i - 1] * p + fib[i - 2];
}
cum[0] = 0;
for (int i = 1; i < n; ++i) {
cum[i] = fib[i] + cum[i - 1];
}
mint ret = 0;
for (int i = 0; i < n; ++i) {
ret += fib[i] * cum[i];
}
cout << ret << "\n";
return 0;
}
kibuna