VJudge - Product and GCD(数论)

传送门：https://vjudge.net/problem/AtCoder-caddi2018_a

题意：

给定N，P表示： ${\prod }_{i=1}^N{a}_i=P$

求 m a x { g c d i = 1 N   a i } max\{gcd_{i=1}^{N}\ a_i\} max{gcdi=1N​ ai​}

思路:

数论
要求P的所有可能的n个乘积的因子的最大公因数

显然，当$n=1$的时候，答案就是P
当$n!=1$的时候：
要求最大公因数，我们可以把所有对这个最大公因数存在贡献的质因子$p_i$求出来
假设有m个质因子$p_i$,每个质因子对答案贡献了$k_i$次，那么最大公因数答案就是 $ans={\prod }_{i=1}^m(p_i{)}^{k_i}$

然后，我们再来分析一下$p_i$和$k_i$如何求出
对于给定的P,先对它质因数分解，
假设$c_i$表示P的一个质因数，$a_i$表示$c_i$对P贡献了多少次，可以得到的是:
$P={c_1}^{a_1}\ast {c_2}^{a_2}\ast {c_3}^{a_3}\ast ......$
要使P的某个质因数对答案ans有贡献，那么该质因数$c_i$对P的贡献次数$a_i$应满足:$a_i>=n$

为什么呢？
因为只有$a_i>=n$的情况下，才有可能分给P的N个因数，对N个因素的最大公因数才有贡献
如果$a_i<n$，那么必然不会对答案有贡献
分析到这里，$k_i$也就很容易得出，$k_i=a_i/n$,也就是$c_i$平均下来后对N个因数的贡献次数

---

内心OS:
老师说这道题草稿纸上画画就会了，然而我拿到这道题的却无从下手，毫无思路
终究还是太菜了，什么时候才能变成真正的大佬/(ㄒoㄒ)/~~
勤能补拙、勤能补拙。。。。。

若有不足有误之处希望大佬们不吝赐教

AC代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;

template<class T> T read() {
	T f=1,k=0;
	char c=getchar();
	while(c<'0' || c>'9') {
		if(c=='-') f=-1;
		c=getchar();
	}
	while(c>='0' && c<='9') {
		k=k*10+c-'0';
		c=getchar();
	}
	return f*k;
}

int main() {
	ll n=read<ll>(),p=read<ll>();
	if(n==1) {
		cout<<p<<endl;
		return 0;
	}
	ll ans=1;
	for(ll i=2; i*i<=p; i++) {
		ll cnt=0;
		if(p%i==0) {
			while(p%i==0)  cnt++,p/=i;
			ll k=cnt/n;//均分质因数
			while(k--) ans*=i;
		}
	}
	cout<<ans<<endl;
	return 0;
}