nefu120 梅森素数 Lucas-Lehmer判定法

这道题目，它的题目意思是表达错了，但是看案例就能看明白到底问的什么

下面介绍一下 对梅森数进行判定是否是梅森素数非常有效的方法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");
		}
	}
}