/** * code generated by JHelper * More info: https://github.com/AlexeyDmitriev/JHelper * @author */ #include #include #ifndef SOLUTION_COMMON_H #include using namespace std; using ll = long long; using PI = pair; template using V = vector; using VI = V; #define _1 first #define _2 second #ifdef MY_DEBUG # define DEBUG(x) x #else # define DEBUG(x) #endif template inline void debug(T &A) { DEBUG( for (const auto &a : A) { cerr << a << " "; } cerr << '\n'; ) } template inline void debug_with_format(T &A, Func f) { DEBUG( for (const auto &a : A) { cerr << f(a) << " "; } cerr << '\n'; ) } template inline void debug_dim2(T &A) { DEBUG( for (const auto &as : A) { debug(as); } ) } template inline void debug(const char *format, Args const &... args) { DEBUG( fprintf(stderr, format, args ...); cerr << '\n'; ) } template string format(const string &fmt, Args ... args) { size_t len = snprintf(nullptr, 0, fmt.c_str(), args ...); vector buf(len + 1); snprintf(&buf[0], len + 1, fmt.c_str(), args ...); return string(&buf[0], &buf[0] + len); } template string fmtP(pair a) { stringstream ss; ss << "(" << a._1 << "," << a._2 << ")"; return ss.str(); } #define SOLUTION_COMMON_H #endif //SOLUTION_COMMON_H const int MOD = 1000000007; class E { V> pow(V> a, int n, int m) { if (n == 1) return a; V> b = pow(mul(a, a, m), n/2, m); if (n % 2 == 1) b = mul(b, a, m); return b; } V> mul(V> a, V> b, int m) { assert(a[0].size() == b.size()); ll th = (ll)7e18; V> res(a.size(), V(b[0].size())); for (int i = 0; i < a.size(); ++i) { for (int j = 0; j < a[0].size(); ++j) { for (int k = 0; k < a[0].size(); ++k) { res[i][j] += a[i][k] * b[k][j]; if (res[i][j] > th) res[i][j] %= m; } res[i][j] %= m; } } return res; } public: void solve(std::istream& in, std::ostream& out) { int a, b, n; in >> a >> b >> n; V> A = {{a, b}, {1, 0}}; auto B = pow(A, n, MOD); out << B[1][0]; } }; int main() { E solver; std::istream& in(std::cin); std::ostream& out(std::cout); solver.solve(in, out); return 0; }