题目:点击打开链接
思路:balabala一大堆,其实本意就是考你威尔逊定理,关于威尔逊定理
如果p为素数,则(p-1)! mod p = p-1;
如果p为合数,则(p-1)! mod p = 0;当p等于4时除外,当p=4时,答案为2而不是0;
所以这个题就符合上边所说的定理,唯一的cheak就是p=4时。
解题:那么既然知道了是判断素数的问题,那么就好做了,但是判断素数的方法有这么多,选择哪一个呢,不妨让我们来试试。
1.素数表筛选法:
前边:今天算是见识了,疯狂输出,狂打表,结果疯狂RE,真心不知道为啥RE,最后终于知道为啥了,但是我真的不知道为什么会RE啊,好气啊,今天一直在纠结素数打表的问题,感觉爆炸了。
素数表的模板就在这一起整理了把
最后的结果是180ms左右
代码:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;
const int maxn = 1e6+10;//尽量开1e6,re不要怪我
bool isprime[maxn];
int prime[maxn];
//打印素数表的模板
void build()
{
int i,j;
memset(isprime,true,sizeof(isprime));
isprime[1]=false;
int nprime=0;
for(i=2;i<maxn;i++){ //wtm为啥写sqrt(maxn)就RE
if(isprime[i]){
prime[nprime++]=i;
for(j=i+i;j<maxn;j+=i) isprime[j]=false;//wtm为啥写i*i就RE!!!什么鬼
}
}
}
bool judgeprime(int x)
{
for(int i=0;prime[i]*prime[i]<=x;i++)
if(x%prime[i]==0)
return false;
return true;
}
int main()
{
int T;
build();
scanf("%d",&T);
while(T--){
int x;scanf("%d",&x);
if(judgeprime(x)) printf("%d\n",x-1);
else{
if(x==4) printf("2\n");
else printf("0\n");
}
}
return 0;
}
2.miller_Rabin算法
前边:这个也是很奇怪,感觉算法很快,但是一直T,最后把乘法改进了下,结果速度变得奇快!!!这个乘法的快速计算的方法,没怎么明白。先收藏下来,不过miller_rabin倒是挺快的,但是最后的结果,没比打表快到那里去。速度是280ms
代码:
include<iostream>
#include<cstdio>
#include<ctime>
#include<cstdlib>
using namespace std;
typedef long long ll;
const int Times=3;//MillerRabin检测次数
ll quickmul(ll x,ll y,ll Z)//大数相乘取模,非常神奇的一个快速乘法的计算方法
{
ll tmp=x/(long double)Z*y+1e-3;
return (x*y+Z-tmp*Z)%Z;
}
ll quickmod(ll a,ll b,ll c)
{
ll ans=1;
while(b){
if(b&1) ans=(ans*a)%c;
a=quickmul(a,a,c);
b>>=1;
}
return ans;
}
bool millerRabin(ll x)//只利用费马小定理的测试函数 a^(x-1)%x==1则 x为素数
{
if(x<=1) return false;
if(x==2) return true;
if(!(x&1)) return false;
for(int i=0;i<Times;i++){
ll a=rand()%(x-1)+1;//随机生成一个数 来进行费马小定理检测
if(quickmod(a,x-1,x)!=1) return false;
}
return true;
}
int main()
{
int T;
while(scanf("%d",&T)!=EOF){
while(T--){
ll x;scanf("%lld",&x);
if(millerRabin(x)) printf("%lld\n",x-1);
else{
if(x==4) printf("2\n");
else printf("0\n");
}
}
}
}