Codeforces Round #589 (Div. 2) C. Primes and Multiplication （数学）

题目链接：https://codeforces.com/contest/1228/problem/C

反思：做这个题的时候没想出思路，补题的时候还把质因数分解那里打错了，找了好久错。

思路：先对x进行质因数分解，对于每个质数，从f(1,p) …f(n,p),先算出能整除p的数字个数，这里肯定还包含能被p^k（k >= 2）整除的，他们的p肯定是要算在答案里的，所以先除掉是可以的，令n = n / p，之前的被整除的数变成了 p ^ （k - 1），就这样每次除p，最后可算出一个总次数，用快速幂即可。

代码：

#include <bits/stdc++.h>

using namespace std;

const int maxn = 1e5 + 5;
const int mod = 1e9 + 7;
vector <int> v;

long long qpow(long long x , long long y ) {
	long long cnt = 1 , t = x;
	while(y) {
		if(y & 1)cnt = t * cnt % mod;
		t = t * t % mod;
		y /= 2;
	}
	return cnt;
}

int main() {
	long long x , n;
	ios::sync_with_stdio(0);
	cin >> x >> n;
	int y = sqrt(x);
	long long t = x;
	for(int i = 2 ; i <= y + 1 ; i++) {
		if(t && t % i == 0) {
			while(t % i == 0) {
				t /= i;
			}
			v.push_back(i);	
		}
	}
	if(t > 1)v.push_back(t);
	long long ans = 1;
	for(int i = 0 ; i < v.size() ; i++) {
		t = n;
		long long cnt = 0;
		while(t) {
			cnt += t / v[i];
			t /= v[i];
		}
		ans = ( ans % mod )* ( qpow(v[i] , cnt) % mod) % mod;
	}
	cout << ans % mod << "\n";
	return 0;
}