#include using namespace std; //------------------------------------------------------------ // ここからはけんちょんさんのブログからコピペ // https://drken1215.hatenablog.com/entry/2018/06/08/210000 const int MAX = 510000; const int MOD = 1000000007; long long fac[MAX], finv[MAX], inv[MAX]; void COMinit() { fac[0] = fac[1] = 1; finv[0] = finv[1] = 1; inv[1] = 1; for (int i = 2; i < MAX; i++){ fac[i] = fac[i - 1] * i % MOD; inv[i] = MOD - inv[MOD%i] * (MOD / i) % MOD; finv[i] = finv[i - 1] * inv[i] % MOD; } } long long COM(int n, int k){ if (n < k) return 0; if (n < 0 || k < 0) return 0; return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD; } // ここまでコピペ //------------------------------------------------------------ int main() { // 入力 long long a, b; int n, k; cin >> a >> b; cin >> n >> k; // a,bがすでに大きいと大変なのでMODをとる a%=MOD; b%=MOD; // 初期化 COMinit(); // (n,k)を計算 cout << (a*COM(n-1,k-1)+b*COM(n-1,k-2)) % MOD << endl; // n段目の数の2乗の和を計算 long long ans = 0; for(int i = 0; i < n+1; ++i) { long long tmp = (a*COM(n-1,i)+b*COM(n-1,i-1)) % MOD; ans += (tmp*tmp) % MOD; ans %= MOD; } cout << ans << endl; return 0; }