快速求整个数组最大公约数

试题链接：2021年牛客寒假集训营第四场–J题

如果我们要求一堆数的最大公约数，我们只要对每一个质数求出在这些数的最大次数中的最小值，最后把这些幂次累乘起来，就是这些数的最大公约数。

#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=1e4+5;
const int mod=1e9+7;
int x[N],p[N];
int main(){
	int n;
	scanf("%d",&n);
	for(int i=0;i<n;i++) scanf("%d",&x[i]);
	for(int i=0;i<n;i++) scanf("%d",&p[i]);
	ll res=1;
	for(int i=2;i<=10000;i++){//枚举因子 
		int cnt=0x3f3f3f3f;//每个数、这个因子的最小次幂 
		for(int j=0;j<n;j++){
			int s=0;
			while(x[j]%i==0){
				x[j]/=i;
				s++;
			}
			cnt=min(s*p[j],cnt);
		}
		for(int k=1;k<=cnt;k++){
			res*=i;
			res%=mod;
		}
	}
	printf("%d\n",res);
}