题目给了一个公式,主要就是让你求这个函数能进去的最大次数 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;
}