欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页

2020牛客暑期多校训练营(第四场)B-Basic Gcd Problem

程序员文章站 2022-04-01 16:11:11
...

传送门

Description

2020牛客暑期多校训练营(第四场)B-Basic Gcd Problem

  • fci(ni)mod(109+7)f_{c_i}(n_i)mod(10^9+7)

题解

  • n=i=1npibin=\prod_{i = 1}^{n}p_i^{b_i}

  • sumpow=i=1nbisumpow=\sum_{i = 1}^{n}b_i

  • fci(ni)=nisumpowf_{c_i}(n_i)=n_i^{sumpow}

实现

#include <bits/stdc++.h>
#define ll long long
const int N=1e6+10,mod=1e9+7;
using namespace std;
ll T,n,c;
bool np[N];
ll tot,p[N],mn[N];
ll ksm(ll a,ll b){ll r=1;while(b){if(b&1)r=r*a%mod;a=a*a%mod,b>>=1;}return r;}
void sieve(){
	np[0]=np[1]=mn[1]=1;
	for(int i=2;i<N;i++){
		if(!np[i]) p[++tot]=i,mn[i]=i;
		for(int j=1;j<=tot&&i*p[j]<N;j++){
			np[i*p[j]]=1,mn[i*p[j]]=p[j];
			if(i%p[j]==0) break;
		}
	}
}int main(){
	sieve(),scanf("%lld",&T);
	while(T--){
		scanf("%lld%lld",&n,&c);
		ll p=0;
		while(n^1) n/=mn[n],++p;
		printf("%lld\n",ksm(c,p));
	}
}