結果

問題 No.1158 GCD Products easy
ユーザー Kome_soudou
提出日時 2022-11-09 14:16:59
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 971 bytes
コンパイル時間 1,928 ms
コンパイル使用メモリ 169,256 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-22 20:17:00
合計ジャッジ時間 3,078 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

template <typename T> T modpow(T a, T b, T p)
{
	T ans = 1;
	while(b > 0)
	{
		if(b & 1) ans = ans * a % p; //bの最下位bitが1ならa^(2^i)をかける
		a = a * a % p;
		b >>= 1; //bを1bit右にシフト
	}
	return ans;
}

int main()
{
	int64_t a, b, n;
	cin >> a >> b >> n;
	const int64_t mod = 1000000007;
	
	vector<int64_t> v(11);
	for(int64_t i = a; i <= b; i++) v[i] = 1;
	
	vector<int64_t> dp(11);
	for(int64_t i = 10; i >= 6; i--) if(v[i] == 1) dp[i] = 1;
	
	dp[5] = (modpow(v[5] + v[10], n, mod) - dp[10] + mod) % mod;
	dp[4] = (modpow(v[4] + v[8], n, mod) - dp[8] + mod) % mod;
	dp[3] = (modpow(v[3] + v[6] + v[9], n, mod) - dp[9] - dp[6] +  2 * mod) % mod;
	dp[2] = (modpow(v[2] + v[4] + v[6] + v[8] + v[10], n, mod) - dp[10] - dp[6] - dp[8] - dp[4] + 3 * mod) % mod;
	
	int64_t ans = 1;
	for(int64_t i = 10; i >= 2; i--)
	{
		ans *= modpow(i, dp[i], mod);
		ans %= mod;
	}
	cout << ans << endl;
	return 0;
}
0