結果

問題 No.1127 変形パスカルの三角形
ユーザー forest3
提出日時 2020-07-26 21:00:11
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 43 ms / 1,500 ms
コード長 1,087 bytes
コンパイル時間 1,575 ms
コンパイル使用メモリ 169,488 KB
実行使用メモリ 9,472 KB
最終ジャッジ日時 2024-06-28 19:22:48
合計ジャッジ時間 3,493 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

int main()
{
	long long a, b;
	int n, k;
	cin >> a >> b >> n >> k;

	const long long MOD = 1000000000 + 7;
	a %= MOD;
	b %= MOD;

	auto power = [&]( long long x, long long nn ) -> long long
	{
		long long res = 1;
		while( nn > 0 ) {
			if( nn % 2 ) res = res * x % MOD;
			x = x * x % MOD;
			nn /= 2;
		}
		return res;
	};

	vector<long long> fac( n );
	vector<long long> ifac( n );
	fac[0] = 1;
	ifac[0] = power( 1, MOD - 2 );
	for( int i = 1; i < n; i++ ) {
		fac[i] = fac[i - 1] * i % MOD;
		ifac[i] = power( fac[i], MOD - 2 );
	}

	vector<long long> na( n + 1 );
	na[0] = 1;
	na[n] = 0;
	vector<long long> nb( n + 1 );
	nb[0] = 0;
	nb[n] = 1;
	for( int i = 1; i < n; i++ ) {
		int nn = n - 1;
		na[i] = fac[nn] * ifac[i] % MOD * ifac[nn - i] % MOD;
		nb[n - i] = na[i];
	}
	long long ans = (na[k - 1] * a % MOD + nb[k - 1] * b % MOD) % MOD;

	cout << ans << endl;

	ans = 0;
	for( int i = 0; i < n + 1; i++ ) {
		long long x = (na[i] * a % MOD + nb[i] * b % MOD) % MOD;
		ans += x * x % MOD;
		ans %= MOD;
	}

	cout << ans << endl;
}
0