2020牛客多校第4场 B

题目给了一个公式，主要就是让你求这个函数能进去的最大次数 cnt。然后求c*c的(cnt-1)次方,即是c的cnt次方

想要函数进去的最大，必然是让进去后得到的x最大，对于gcd能得到最大值就一定会是x的最大因数

列如 x = 24 . 则 24 -> 12 -> 6 -> 3 -> 1 进去了4次。 可以知道，对于x和它的最大因数（y）相除，可以得到它的最小质因数，同理，对于y的最大因数(z)相除，也可以得到y的最小质因数这也是x的最小质因数。递归下去，发现就是求出尽可能小的质因子，而对于任何一个大于1的数，有唯一分解定律。
故即是求出x用唯一分解定律后，所有质因子次数的个数

可以先用欧拉筛把质数全部筛出来。（正解）
下面我是直接循环。
x为long long 会超时。因为这样写的话复杂度上界是 1e9（虽然到不了）
代码如下：

#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#include<bits/stdc++.h>
//#define int long long
#define sc scanf
#define pf printf
using namespace std;
typedef pair<int,int> pii;
typedef long long ll;
  
const int INF = 0x3f3f3f3f;
const double eps = 1e-5;
const int mod = 1e9+7;
const int N = 200010;
 
ll qmi(ll a,ll b){
    ll res = 1;
    while(b){
        if(b&1) res = res*a%mod;
        a = a*a%mod;
        b >>= 1;
    }
    return res;
}
 
signed main(){
//  IOS;
    #ifdef ddgo
        freopen("C:/Users/asus/Desktop/ddgoin.txt","r",stdin);
    #endif
     
    int tt; sc("%d",&tt);
    while(tt --){
        int x,c; sc("%d%d",&x,&c);
        int cnt = 0;
        for(int i=2;i*i<=x;i++)
            while(x%i == 0) x/=i,cnt++;
        if(x > 1) cnt ++;
        pf("%lld\n",qmi(c,cnt)%mod);
    }
     
    return 0;
}

欧拉筛，很多地方要处理一下：

#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define CLOCK printf("Time used:%.2lf\n",(double)clock()/CLOCKS_PER_SEC);
#include<bits/stdc++.h>
//#define int long long
#define sc scanf
#define pf printf
using namespace std;
typedef pair<int,int> pii;
typedef long long ll;
 
const int INF = 0x3f3f3f3f;
const double eps = 1e-5;
const int mod = 1e9+7;
const int N = 1000010;

int is_prime[N],primes[N],k;

void ols(int n) 
{
    is_prime[1] = 1;
    for (int i = 2; i <= n; i++) {
        if (!is_prime[i]) primes[k++] = i;
        for (int j = 0; j<k&&primes[j] <= n / i; j++) {
            is_prime[primes[j] * i] = 1;
            if (i % primes[j] == 0) break;
        }
    }
}

ll qmi(ll a,ll b){
	ll res = 1;
	while(b){
		if(b&1) res = res*a%mod;
		a = a*a%mod;
		b >>= 1;
	}
	return res;
}

signed main(){
//	IOS;
	#ifdef ddgo
		freopen("C:/Users/asus/Desktop/ddgoin.txt","r",stdin);
	#endif
	ols(1000000);
	int tt; sc("%d",&tt);
	while(tt --){
		int x,c; sc("%d%d",&x,&c);
		int cnt = 0;
		if(!is_prime[x]){
			cnt = 1;
		}else{
			for(int i=0;i<k&&primes[i] <= x&&is_prime[x];i++)
			while(x%primes[i] == 0) x/=primes[i],cnt++;
			if(x > 1) cnt ++;
		}
		pf("%lld\n",qmi(c,cnt)%mod);
	}
    return 0;
}