結果
問題 |
No.1127 変形パスカルの三角形
|
ユーザー |
|
提出日時 | 2020-08-20 17:10:57 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 39 ms / 1,500 ms |
コード長 | 1,168 bytes |
コンパイル時間 | 2,233 ms |
コンパイル使用メモリ | 192,780 KB |
最終ジャッジ日時 | 2025-01-13 04:26:42 |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 30 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define REP(i,n) for(int i=0; i<(int)(n); i++) #define FOR(i,b,e) for (int i=(int)(b); i<(int)(e); i++) #define ALL(x) (x).begin(), (x).end() const double PI = acos(-1); const long long MOD = 1e9 + 7; long long f[200000], r[200000]; long long modpow(long long x, long long p, long long mod) { long long ret = 1; while (p) { if (p & 1) ret = ret * x % mod; x = x * x % mod; p >>= 1; } return ret; } long long comb(int n, int k) { if (0 <= k && k <= n) return f[n] * r[k] % MOD * r[n-k] % MOD; return 0; } int main() { ios_base::sync_with_stdio(0); cin.tie(0); r[0] = f[0] = 1; for (int i = 1; i < 200000; i++) { f[i] = i * f[i-1] % MOD; r[i] = modpow(f[i], MOD - 2, MOD); } int n, k; long long a, b; cin >> a >> b >> n >> k; a %= MOD; b %= MOD; long long val = (a * comb(n-1, k-1) + b * comb(n-1, k-2)) % MOD; long long sum = 0; for (int k = 1; k <= n+1; k++) { long long tmp = (a * comb(n-1, k-1) + b * comb(n-1, k-2)) % MOD; sum += tmp * tmp % MOD; sum %= MOD; } cout << val << endl; cout << sum << endl; return 0; }