「hdu6608」 Fansblog
今天回顾之前多校联赛的题目:Fansblog ,发现一件有趣的事情,快速积
的时间复杂度,比普通*
快多了,刚才一直Tle
在用了普通*
,卧槽这也能T,当然也有乘法爆long long
的问题。
快速积 >> * 速度的原因:
威尔逊定理:
- 威尔逊定理就是对于任意的正质数 K 有:
((k−1)!)%k=k−1
题解:逆元+威尔逊定理
- 首先我们用普通的素数筛,筛到小于p的最大素数
- 然后利于威尔逊定理推得:
((p−1)!)%p=q!×(q+1)×(q+2)×...×(p−1)%p=(p−1)
推 |
导 |
得 V
q!≡ (1/((q+1)*(q+2)..... *(p-2)))(mod p)
- 然后利用逆元,求: q! (mod p)
AC代码:
#include<bits/stdc++.h>
using namespace std;
#define rep(i,a,b) for(long long i=(a);i<=(b);i++)
#define ll long long
ll ksc(ll a,ll b,ll mod)
{
ll res=0;
while(b)
{
if(b&1)
res=res+a%mod;
a=a+a%mod;
b>>=1;
}
return res%mod;
}
ll ksm(ll a,ll b,ll mod)
{
ll res=1;
while(b)
{
if(b&1)
res=ksc(res,a,mod);
a=ksc(a,a,mod);
b>>=1;
}
return res%mod;
}
bool shai(ll x)
{
for(ll i=2;i<=sqrt(x);i++)
{
if(x%i==0)
return false;
}
return true;
}
int main()
{
ios::sync_with_stdio(false);
cin.tie(0);
ll t;
cin>>t;
while(t--)
{
ll p;
cin>>p;
ll q=p;
while(q--)
{
if(shai(q))
break;
}
ll res=1;
rep(i,q+1,p-2)
{
res=ksc(res,ksm(i,p-2,p),p);
}
cout<<res<<endl;
}
}