2020牛客多校第四场 B

2020牛客多校第四场 B

解法：
如果要最优，必须要使gcd的次数最多。
已知x=p1^k1p2^k2…
我们可以让y=gcd(i,x)=p1^k1-1p2^k2…
以此类推，最多可以执行(k1+k2+…kn)次。

#include <bits/stdc++.h>
 
using namespace std;
typedef long long ll;
const int N = 1e6;
bool vis[N+10];
int p[N+10],a[N+10];
void init(int n){
   int m=0;
   for(int i=2;i<=n;i++){
        if(!vis[i])p[++m]=i,a[i]=1;
        for(int j=1;j<=m&&p[j]*i<=n;j++){
            vis[p[j]*i]=true;
            a[p[j]*i]=a[i]+1;
            if(i%p[j]==0)break;
        }
   }
}
 
const ll mod = 1e9 +7;
 
ll ksm(ll a, ll p)
{
    ll res = 1;
    while(p)
    {
        if (p & 1) res = res * a % mod;
        a = a * a % mod;
        p >>= 1;
    }
    return res;
}
 
int main()
{
    init(1e6);
    ll t, n, c;
    scanf("%lld", &t);
    while(t--)
    {
        scanf("%lld%lld", &n, &c);
        printf("%lld\n", ksm(c, a[n]));
    }
}