#pragma GCC target("avx2") #pragma GCC optimize("03") #pragma GCC optimize("unroll-loops") #include using namespace std; typedef long double ld; typedef long long ll; typedef unsigned long long ull; #define endl "\n" #define FOR(i,a,b) for(int i=(a);i<=(b);i++) #define rep(i,n) for(int i=0;i<(n);i++) #define PII pair #define PLL pair #define ALL(x) (x).begin(), (x).end() constexpr int INF=1<<30; constexpr ll LINF=1LL<<60; constexpr ll mod=1e9+7; constexpr int NIL = -1; templateinline bool chmax(T &a, const T &b) { if (ainline bool chmin(T &a, const T &b) { if (binline int popcount(T a) {return __builtin_popcount(a);} templateinline T emod(T a, T p) { return (a%p + p) % p;} struct mint { ll x; mint(ll x=0):x(x%mod){} bool operator==(const mint a)const{return x==a.x;} bool operator!=(const mint a)const{return x!=a.x;} bool operator>=(const mint a){return (x >= a.x)? 1: 0;} bool operator<(const mint a){return !(*this>=a);} bool operator>(const mint a){return (x > a.x)? 1:0;} bool operator<=(const mint a){return !(*this>a);} mint& operator+=(const mint a) { if ((x += a.x) >= mod) x -= mod; return *this; } mint& operator-=(const mint a) { if ((x += mod-a.x) >= mod) x -= mod; return *this; } mint& operator*=(const mint a) { (x *= a.x) %= mod; return *this; } mint& operator/=(const mint a) { return (*this) *= a.inv(); } mint operator+(const mint a) const { mint res(*this); return res+=a; } mint operator-(const mint a) const { mint res(*this); return res-=a; } mint operator*(const mint a) const { mint res(*this); return res*=a; } mint operator/(const mint a) const { mint res(*this); return res/=a; } mint pow(ll t) const { if (!t) return 1; mint a = pow(t>>1); a *= a; //2 square if (t&1) a *= *this; return a; } // for prime mod mint inv() const { return pow(mod-2); } };//------------------- int main() { cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(15); ll B,C,D; cin >> B >> C >> D; mint b = mint(B); mint c = mint(C); if(C==1) { auto ans = b*mint(D); cout << ans.x << endl; return 0; } mint dv = b * c / ( c -1); auto s = c.pow(D) * dv; s -= dv; cout << s.x << endl; return 0; }