#include using namespace std; #define inf INT_MAX #define INF LLONG_MAX #define ll long long #define ull unsigned long long #define M (int)(1e9+7) #define P pair #define FOR(i,m,n) for(int i=(int)m;i<(int)n;i++) #define RFOR(i,m,n) for(int i=(int)m;i>=(int)n;i--) #define rep(i,n) FOR(i,0,n) #define rrep(i,n) RFOR(i,n,0) #define all(a) a.begin(),a.end() const int vx[4] = {0,1,0,-1}; const int vy[4] = {1,0,-1,0}; #define F first #define S second #define PB push_back #define EB emplace_back #define int ll #define vi vector #define IP pair #define PI pair #define PP pair #define Yes(f){cout<<(f?"Yes":"No")<= 0 ? y % M : (M - (-y) % M) % M) {} Mint &operator+=(const Mint &p) { if((x += p.x) >= M) x -= M; return *this; } Mint &operator-=(const Mint &p) { if((x += M - p.x) >= M) x -= M; return *this; } Mint &operator*=(const Mint &p) { x = (int) (1LL * x * p.x % M); return *this; } Mint &operator/=(const Mint &p) { *this *= p.inverse(); return *this; } Mint operator-() const { return Mint(-x); } Mint operator+(const Mint &p) const { return Mint(*this) += p; } Mint operator-(const Mint &p) const { return Mint(*this) -= p; } Mint operator*(const Mint &p) const { return Mint(*this) *= p; } Mint operator/(const Mint &p) const { return Mint(*this) /= p; } bool operator==(const Mint &p) const { return x == p.x; } bool operator!=(const Mint &p) const { return x != p.x; } Mint inverse() const { int a = x, b = M, u = 1, v = 0, t; while(b > 0) { t = a / b; swap(a -= t * b, b); swap(u -= t * v, v); } return Mint(u); } Mint power(int n) const { Mint ret(1), mul(x); while(n > 0) { if(n & 1) ret *= mul; mul *= mul; n >>= 1; } return ret; } friend ostream &operator<<(ostream &os, const Mint &p) { return os << p.x; } friend istream &operator>>(istream &is, Mint &a) { int t; is >> t; a = Mint(t); return (is); } }; using mint = Mint; template using Matrix = vector>; template Matrix mat_add(Matrix A, Matrix B) { int m = A.size(); int n = B[0].size(); Matrix C(m, vector(n)); for (int i = 0; i < m; i++) { for (int j = 0; j < n; j++) { C[i][j] = A[i][j] + B[i][j]; } } return C; } template Matrix mat_sub(Matrix A, Matrix B) { int m = A.size(); int n = B[0].size(); Matrix C(m, vector(n)); for (int i = 0; i < m; i++) { for (int j = 0; j < n; j++) { C[i][j] = A[i][j] - B[i][j]; } } return C; } template Matrix mat_mul(Matrix A, Matrix B) { int m = A.size(); int l = B.size(); int n = B[0].size(); Matrix C(m, vector(n)); for (int i = 0; i < m; i++) { for (int j = 0; j < n; j++) { for (int k = 0; k < l; k++) { C[i][j] += A[i][k] * B[k][j]; } } } return C; } template Matrix mat_pow(Matrix A, ll p) { Matrix ret(A.size(), vector(A.size(), 0)); for (int i = 0; i < A.size(); i++) ret[i][i] = 1; while (p > 0) { if (p & 1) ret = mat_mul(ret, A); A = mat_mul(A, A); p >>= 1; } return ret; } signed main(){ cin.tie(0); ios::sync_with_stdio(false); cout<>a>>b>>n; if(n==0){ cout<<0< A(2,vector(2)); A[0][0]=a; A[0][1]=b; A[1][0]=1; A[1][1]=0; Matrix B(2,vector(1)); B[0][0]=1; B[0][1]=0; auto ans=mat_mul(mat_pow(A,n-1),B); cout<