題意:給你一個素數(shù)Q,讓你找到一個最大的P<Q, 求P! mod Q。
思路:威爾遜定理就是對于一個素數(shù)P, 滿足(P - 1) ! ≡ -1 (mod P) , 其實也就是(P - 1) ! mod P == P - 1, 所以我們就能很容易的解決這道題了,我是打了1e7的素數(shù)表,然后相乘的時候用了快速乘防止爆long long,另外不要忘了除法取模要用逆元。
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 1e7 + 10;
int prime[2000000], a[maxn];
int num = 0;
void initprime() //線性篩
{ memset(a, 0, sizeof(a));
memset(prime, 0, sizeof(prime));
for(int i = 2; i < maxn; ++i)
{
if(!a[i]) prime[num++] = i;
for(int j = 0; j < num && i * prime[j] < maxn; ++j)
{
a[i * prime[j]] = 1;
if(!(i % prime[j])) break;
}
}
}
ll mul(ll a, ll b, ll c)//快速乘 防止ll 相乘會炸ll
{
ll res = 0;
while(b)
{
if(b & 1)
res = (res + a) % c;
b >>= 1;
a = (a << 1) % c;
}
return res % c;
}
ll powermod(ll a, ll b, ll c)
{
a = a % c;
ll res = 1;
while(b)
{
if(b & 1)
res = mul(res, a, c);
b >>= 1;
a = mul(a, a, c);
}
return res % c;
}
bool check(ll x) //檢查x是否為素數(shù)
{
for(int i = 0; i < num && prime[i] < x; ++i)
{
if(x % prime[i] == 0)
return false;
}
return true;
}
int main()
{
initprime();
int t;
scanf("%d", &t);
while(t--)
{
ll Q, P;
scanf("%lld", &Q);
bool flag = false;
ll ans = Q - 1;
for(P = Q - 1; P > 1; --P)
{
if(check(P))
{
flag = true;
break;
}
}
if(flag)
{
for(ll i = P + 1; i <= Q - 1; ++i)
{
ans = mul(ans, powermod(i, Q - 2, Q), Q);
}
}
printf("%lld\n", ans);
}
return 0;
}
本站文章版權歸原作者及原出處所有 。內容為作者個人觀點, 并不代表本站贊同其觀點和對其真實性負責,本站只提供參考并不構成任何投資及應用建議。本站是一個個人學習交流的平臺,網(wǎng)站上部分文章為轉載,并不用于任何商業(yè)目的,我們已經(jīng)盡可能的對作者和來源進行了通告,但是能力有限或疏忽,造成漏登,請及時聯(lián)系我們,我們將根據(jù)著作權人的要求,立即更正或者刪除有關內容。本站擁有對此聲明的最終解釋權。