链接:https://codeforces.com/contest/1228/problem/C
来源:Codeforces Round #589 (Div. 2)
题意:给出一些定义,对于一个 x ,prime(x) ={p1, p2…, pn},g(x, p) 表示 x 最大可以整除 p 的多少次幂,这个 pk 就是 g(x, p) ,f(x, n) = g(n, p1) × g(n, p2) ×···× g(n,pn),prime(x) ={p1, p2…, pn},p1, p2…, pn表示 x 的质因子。最后求出 f(x, 1) × f(x, 2) × ··· × f(x, n)。
思路: f(x, 1) × f(x, 2) × ··· × f(x, n) = g(1, p1) × g(1, p2) × ··· ×g(1, pn) × g(2, p1) × g(2, p2) × ··· ×g(1, pn) × ··· × g(n, p1) × g(n, p2) × ··· ×g(n, pn) ,从这个式子可以看出对于p1的贡献,在[1, n] 中有 p1,p1 × 某个数(这些数可以贡献一个p1),对于 p12 的贡献就是p12 × 某个数(这些数可以贡献出两个p1)···综上所述,我们反过来想这么一个过程其实就是在寻找对于[1, n]之间的数字能够贡献多少p1,p2···,pn。对于 n 来说,n / p1,在这些数之间有 n / p1个数可以整除 p1,故有n / p1 个数可以贡献出一个 p1,n / p12,在这些数之间有 n / p12个数可以整除 p1,故有n / p12 个数可以在之前的基础上贡献出一个 p1,一次类推。对于p,他出现的次数就是 cnt = n / p + n / p2+ n / p3+ ···· + n / pn,这个值就是 pcnt。这个问题就这样转换成求快速幂的过程。对于此题找出所有的 p ,计算出次数 快速幂累乘即可解决。
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod = 1e9 + 7;
const int Max_n = 1e6 + 10;
int prime[Max_n];
ll qpow(ll a, ll b) {
ll ans = 1, res = a;
while(b) {
if(b & 1) ans = ans * res % mod;
res = res * res % mod;
b >>= 1;
}
return ans % mod;
}
int main() {
ll x, n, inx = 0;
scanf("%lld%lld", &x, &n);
for(int i = 2; i * i <= x; i++) {
if(!(x % i)) {
prime[inx++] = i;
while(!(x % i)) x /= i;
}
}
if(x > 1) prime[inx++] = x;
ll ans = 1;
for(int i = 0; i < inx; i++) {
ll t = n, cnt = 0;
while(t) {//求次数
t /= prime[i];
cnt += t;
}
ans = (ans * qpow(prime[i], cnt)) % mod;
}
printf("%lld\n", ans);
return 0;
}