結果

問題 No.978 Fibonacci Convolution Easy
ユーザー wait_sushiwait_sushi
提出日時 2020-01-31 22:02:21
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 38 ms / 2,000 ms
コード長 2,644 bytes
コンパイル時間 1,639 ms
コンパイル使用メモリ 170,800 KB
実行使用メモリ 20,188 KB
最終ジャッジ日時 2023-10-19 01:02:35
合計ジャッジ時間 2,802 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 3 ms
4,348 KB
testcase_01 AC 17 ms
11,996 KB
testcase_02 AC 11 ms
8,108 KB
testcase_03 AC 37 ms
20,188 KB
testcase_04 AC 12 ms
8,108 KB
testcase_05 AC 4 ms
4,348 KB
testcase_06 AC 16 ms
11,996 KB
testcase_07 AC 29 ms
20,188 KB
testcase_08 AC 18 ms
11,996 KB
testcase_09 AC 30 ms
20,188 KB
testcase_10 AC 38 ms
20,188 KB
testcase_11 AC 15 ms
11,996 KB
testcase_12 AC 5 ms
4,348 KB
testcase_13 AC 16 ms
11,996 KB
testcase_14 AC 6 ms
5,288 KB
testcase_15 AC 19 ms
11,996 KB
testcase_16 AC 37 ms
20,188 KB
testcase_17 AC 38 ms
20,188 KB
testcase_18 AC 2 ms
4,348 KB
testcase_19 AC 2 ms
4,348 KB
testcase_20 AC 2 ms
4,348 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef long long int ll;
typedef vector<ll> vl;
typedef pair<ll, ll> PP;
#define rep(i, n) for(ll i = 0; i < ll(n); i++)
#define all(v) v.begin(), v.end()
#define inputv(v, n)                                                           \
    vl v;                                                                      \
    rep(i, n) {                                                                \
        ll x;                                                                  \
        cin >> x;                                                              \
        v.push_back(x);                                                        \
    }
bool chmin(ll& a, ll& b) { if (b < a) { a = b; return 1; } return 0; }
bool chmax(ll& a, ll& b) { if (b > a) { a = b; return 1; } return 0; }
const ll INF = 999999999999999;
const ll MOD = 1000000007;
const ll MAX_N = 500010;
ll a, b, c, d, e, f, p, t, x, y, z, q, m, n, r, h, k, w, l, ans;
#include <cstdint>                                                                                                              

template <std::uint_fast64_t Modulus> class modint {
	using u64 = std::uint_fast64_t;

public:
	u64 a;

	constexpr modint(const u64 x = 0) noexcept : a(x% Modulus) {}
	constexpr u64& value() noexcept { return a; }
	constexpr const u64& value() const noexcept { return a; }
	constexpr modint operator+(const modint rhs) const noexcept {
		return modint(*this) += rhs;
	}
	constexpr modint operator-(const modint rhs) const noexcept {
		return modint(*this) -= rhs;
	}
	constexpr modint operator*(const modint rhs) const noexcept {
		return modint(*this) *= rhs;
	}
	constexpr modint operator/(const modint rhs) const noexcept {
		return modint(*this) /= rhs;
	}
	constexpr modint& operator+=(const modint rhs) noexcept {
		a += rhs.a;
		if (a >= Modulus) {
			a -= Modulus;
		}
		return *this;
	}
	constexpr modint& operator-=(const modint rhs) noexcept {
		if (a < rhs.a) {
			a += Modulus;
		}
		a -= rhs.a;
		return *this;
	}
	constexpr modint& operator*=(const modint rhs) noexcept {
		a = a * rhs.a % Modulus;
		return *this;
	}
	constexpr modint& operator/=(modint rhs) noexcept {
		u64 exp = Modulus - 2;
		while (exp) {
			if (exp % 2) {
				*this *= rhs;
			}
			rhs *= rhs;
			exp /= 2;
		}
		return *this;
	}
};

using mint = modint<MOD>;
int main() {
	cin >> n >> p;
	vector<mint> A,B;
	mint ans = 0,wa=0;
	A.push_back(0);
	A.push_back(1);
	rep(i, n) {
		A.push_back(A[i] + A[i + 1] * p);
	}
	rep(i, n) {
		wa += A[i];
	}
	wa = wa * wa;
	wa=wa;
	rep(i, n) {
		wa += A[i]*A[i];
	}
	wa = wa / 2;
	cout << wa.a << endl;
}
0