这道题目,它的题目意思是表达错了,但是看案例就能看明白到底问的什么
下面介绍一下 对梅森数进行判定是否是梅森素数非常有效的方法Lucas-Lehmer判定法
设p是素数,第p个梅森数为Mp=2^p-1,r1=4,对于k>=2;利用
rk%Mp=r^2(k-1)-2,0<=rk<Mp,注明 这里的 k是下标,不好打出来,博客还要留着自己看呢,以后不要混淆了
这样可以得到rk的序列,则有Mp是素数,当且仅当r(p-1)%Mp=0,这里的p也是下标
#include<iostream>
#include<cstdio>
#include<list>
#include<algorithm>
#include<cstring>
#include<string>
#include<queue>
#include<stack>
#include<map>
#include<vector>
#include<cmath>
#include<memory.h>
#include<set>
#define ll long long
#define LL __int64
#define eps 1e-8
#define e 2.718281828
//const ll INF=9999999999999;
#define M 400000100
#define inf 0xfffffff
using namespace std;
//vector<pair<int,int> > G;
//typedef pair<int,int> P;
//vector<pair<int,int>> ::iterator iter;
//
//map<ll,int>mp;
//map<ll,int>::iterator p;
//
//vector<int>G[30012];
ll detal(ll x,ll y,ll m)//x*y%m
{
ll temp=0;
while(y>0)//把乘法换成加法,直接用乘法平方会溢出
{
if(y&1)//y为奇数
temp=(temp+x)%m;
y>>=1;
x=(x<<1)%m;
}
return temp;
}
int main(void)
{
int t;
int p;
cin>>t;
while(t--)
{
cin>>p;
ll m=1;
ll a[112];
a[1]=4;
ll temp;
m<<=p;
m-=1;
for(int i=2;i<=p-1;i++)
{
temp=detal(a[i-1],a[i-1],m);
a[i]=(temp-2)%m;
}
if(p==2)
puts("yes");
else
{
if(a[p-1]==0)
puts("yes");
else
puts("no");
}
}
}