#line 1 "main.cpp" #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include /* template start */ using i64 = std::int_fast64_t; using u64 = std::uint_fast64_t; #define rep(i, a, b) for (i64 i = (a); (i) < (b); (i)++) #define all(i) i.begin(), i.end() #ifdef LOCAL #define debug(...) \ std::cerr << "LINE: " << __LINE__ << " [" << #__VA_ARGS__ << "]:", \ debug_out(__VA_ARGS__) #else #define debug(...) #endif void debug_out() { std::cerr << std::endl; } template void debug_out(Head h, Tail... t) { std::cerr << " " << h; if (sizeof...(t) > 0) std::cout << " :"; debug_out(t...); } template std::ostream& operator<<(std::ostream& os, std::pair pa) { return os << pa.first << " " << pa.second; } template std::ostream& operator<<(std::ostream& os, std::vector vec) { for (std::size_t i = 0; i < vec.size(); i++) os << vec[i] << (i + 1 == vec.size() ? "" : " "); return os; } template inline bool chmax(T1& a, T2 b) { return a < b && (a = b, true); } template inline bool chmin(T1& a, T2 b) { return a > b && (a = b, true); } template constexpr Num mypow(Num a, u64 b, Num id = 1) { if (b == 0) return id; Num x = id; while (b > 0) { if (b & 1) x *= a; a *= a; b >>= 1; } return x; } template std::vector> enumerate(const std::vector& data) { std::vector> ret; for (std::size_t index = 0; index < data.size(); index++) ret.emplace_back(index, data[index]); return ret; } /* template end */ #line 2 "/home/imperi/cp-library/lib/utility/modint.hpp" /* author:noshi91 reference:https://noshi91.hatenablog.com/entry/2019/03/31/174006 https://github.com/noshi91/Library/blob/master/other/modint.cpp 1つめのmodintにoperator==,!=を追加したものです thx @noshi91!! */ #line 14 "/home/imperi/cp-library/lib/utility/modint.hpp" template 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 bool operator==(const modint rhs) const noexcept { return a == rhs.a; } constexpr bool operator!=(const modint rhs) const noexcept { return !(*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) 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; } }; #line 100 "main.cpp" constexpr u64 MOD=1000000007; using mint=modint; #line 2 "/home/imperi/cp-library/lib/math/Matrix.hpp" #line 6 "/home/imperi/cp-library/lib/math/Matrix.hpp" template class Matrix { public: using value_t = Field; using size_t = std::size_t; static constexpr size_t row = ROW; static constexpr size_t column = COLUMN; std::array, ROW> a; explicit Matrix(const value_t init = value_t()) { std::fill(a[0].begin(), a[ROW - 1].end(), init); } Matrix(const std::array, ROW>& a_) : a(a_) {} constexpr value_t& operator()(size_t i, size_t j) { return a[i][j]; } constexpr bool operator==(const Matrix& rhs) { for (size_t i = 0; i < ROW; i++) { for (size_t j = 0; j < COLUMN; j++) { if (a[i][j] != rhs.a[i][j]) return false; } } return true; } constexpr bool operator!=(const Matrix& rhs) { return !(*this == rhs); } constexpr Matrix operator+(const Matrix rhs) { return Matrix(*this) += rhs; } constexpr Matrix operator-(const Matrix rhs) { return Matrix(*this) -= rhs; } template constexpr Matrix operator*( const Matrix rhs) { Matrix tmp; for (size_t r = 0; r < ROW; r++) { for (size_t c = 0; c < NEWCOLUMN; c++) { for (size_t i = 0; i < COLUMN; i++) { tmp.a[r][c] += a[r][i] * rhs.a[i][c]; } } } return tmp; } constexpr Matrix& operator+=(const Matrix& rhs) { for (size_t i = 0; i < ROW; i++) { for (size_t j = 0; j < COLUMN; j++) { a[i][j] += rhs.a[i][j]; } } return *this; } constexpr Matrix& operator-=(const Matrix& rhs) { for (size_t i = 0; i < ROW; i++) { for (size_t j = 0; j < COLUMN; j++) { a[i][j] -= rhs.a[i][j]; } } } constexpr Matrix& operator*=(const Matrix& rhs) { std::array, ROW> tmp; std::fill(tmp[0].begin(), tmp[ROW - 1].end(), 0); for (size_t r = 0; r < ROW; r++) { for (size_t c = 0; c < COLUMN; c++) { for (size_t i = 0; i < COLUMN; i++) { tmp[r][c] += a[r][i] * rhs.a[i][c]; } } } a = std::move(tmp); return *this; } constexpr static Matrix id() { static_assert(ROW == COLUMN, "ROW must be equal to COLUMN"); Matrix tmp; for (size_t i = 0; i < ROW; i++) tmp.a[i][i] = 1; return tmp; } void print(){ for(size_t r=0;r; int main(){ std::cin.tie(nullptr); std::ios::sync_with_stdio(false); i64 a,b,n; std::cin>>a>>b>>n; mat ret = mat::id(); mat x({{ {a,b},{1,a}}}); u64 bit = n; while(bit>0){ if(bit&1)ret *= x; x*=x; bit>>=1; } std::cout<<(ret(0,0)*2).value()<<"\n"; return 0; }