結果

問題 No.1659 Product of Divisors
ユーザー snrnsidy
提出日時 2021-10-22 22:34:29
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 10 ms / 2,000 ms
コード長 1,495 bytes
コンパイル時間 2,252 ms
コンパイル使用メモリ 201,952 KB
最終ジャッジ日時 2025-01-25 03:58:28
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
 
using namespace std;
 
vector <int> prime;
map <int,int> cnt;
const long long int MOD = 1e9 + 7;

long long int mypow(long long int x,long long int n)
{
	long long int res = 1;
	while(n>0)
	{
		if(n%2==1)
		{
			res = ((res%MOD)*(x%MOD)%MOD);
			res%=MOD;
		}
		x = ((x%MOD)*(x%MOD)%MOD);
		x%=MOD;
		n/=2;
	}
	return res;
}
 
void eratos(void)
{
	bool isprime[1000001];
	memset(isprime,true,sizeof(isprime));
	for(int i=2;i<=1000000;i++)
	{
		if(isprime[i])
		{
			prime.push_back(i);
			for(int j=2*i;j<=1000000;j+=i)
			{
				isprime[j] = false;
			}
		}
	}
}
 
int main(void)
{
	cin.tie(0);
	ios::sync_with_stdio(false);
 
	long long int N,M;
 
	cin >> N >> M;
 
	eratos();

	for(int i=0;i<prime.size();i++)
	{
		if(N==1 || N < prime[i]*prime[i])
		{
			break;
		}
		if(N%prime[i]==0)
		{
			int count = 0;
			while(1)
			{
				if(N%prime[i]!=0)
				{
					break;
				}
				count++;
				N/=prime[i];
			}
			cnt[prime[i]] = count;
		}
	}
 
	if(N!=1)
	{
		cnt[N] = 1;
	}
 
	long long int res = 1;
 
	for(auto it : cnt)
	{
		long long int y = it.second;
		long long int x = M;
		long long int A = x + y;
		long long int B = y;
		long long int val = 1;
		for(int i=0;i<B;i++)
		{
			long long int X = (A-i)%MOD;
			if(X < 0) X += MOD;
			val = ((val%MOD)*(X%MOD)%MOD);
			val%=MOD;
		}
		for(int i=B;i>=1;i--)
		{
			val = ((val%MOD)*(mypow(i,MOD-2)%MOD)%MOD);
			val%=MOD;
		}
		res = ((res%MOD)*(val%MOD)%MOD);
		res%=MOD;
	}
 
	cout << res << '\n';
 
	return 0;	
}
0