「hdu6608」Fansblog 威尔逊定理

「hdu6608」 Fansblog

今天回顾之前多校联赛的题目：Fansblog ，发现一件有趣的事情，快速积的时间复杂度，比普通*快多了，刚才一直Tle 在用了普通*，卧槽这也能T。

快速积 >> * 速度的原因:

威尔逊定理：

• 威尔逊定理就是对于任意的正质数 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;
    }
}