CF 615D Multipliers 素因子,费马小定理

题意:给出n的素数表示形式(可以相同的p[i],并且指数都为1):n=p1*p2*...p[m].
m,p[i]<=2e5. 问n的所有因子相乘的结果mod1e9+7 

设y[i]为n的第i个因子,p[i]为n的第i个素数,q[i]为第i个素数的指数.
则 res= y[1]*y[2]....y[k]= (p[1]^q[1] * p[2]^q[2] *....p[n]^q[n]) *(.......) ..(每个因子,素因子分解)

考虑p[i]^q[i] 这个素因子一共出现了多少次，出现次数等于其他素因子任意选择的次数.设其余出现次数为b 
b很大,注意！！!： a^b %p != a^(b%p) %p  

利用Fermat's Little Theorem ： a^(p-1)≡1(mod p) 令b=k(p-1)+r 则 a^b≡a^r (mod p)

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=2e5+5,mod=1e9+7;
ll p[N],pre[N],suf[N],a[N],n;
ll powmod(ll x,ll n)
{
	ll s=1;
	while(n)
	{
		if(n&1)
			s=(s*x)%mod;
		n>>=1;
		x=(x*x)%mod;
	}
	return s;
}
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);
	cin>>n;
	for(int i=1;i<=n;i++)
		cin>>p[i],a[p[i]]++;
	ll s=1,res=1;
	pre[0]=suf[N-1]=1;
	for(int i=1;i<N;i++)
	{
		pre[i]=pre[i-1];
		if(a[i])
			pre[i]=(pre[i]*(a[i]+1))%(mod-1);
	}
	for(int i=N-2;i>=2;i--)
	{
		suf[i]=suf[i+1];
		if(a[i])
			suf[i]=(suf[i+1]*(a[i]+1))%(mod-1);
	}
	for(int i=1;i<N;i++)
	{
		ll b=(pre[i-1]*suf[i+1])%(mod-1);
		ll x=1;
		for(int j=1;j<=a[i];j++)
		{
			x=(x*i)%mod;
			res=(res*powmod(x,b))%mod;
		}
	}
	cout<<res<<'\n';
	return 0;
}