-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathwitmath.cpp
88 lines (88 loc) · 1.24 KB
/
witmath.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
#include <bits/stdc++.h>
using namespace std;
#define ll long long int
int A[]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83,89,97};
bool isprime(ll);
ll pwer(ll a,ll b,ll mod);
ll expo(ll a, ll b,ll mod);
int main()
{
int T;
ll N,i,j;
cin >> T;
while(T--)
{
cin >> N;
for(i=N;i>=2;i--){
if(isprime(i)){
cout << i << endl;
break;
}
}
}
return 0;
}
bool isprime(ll n)
{
if(n==2 || n==3)
return true;
if(n%2==0)
return false;
ll d = n-1;
for(int i=0;i<25;i++)
{
if(A[i]==n)
return true;
else if(n%A[i]==0)
return false;
}
int freq=0;
while(d%2==0){
d=d/2;
freq++;
}
for(int i=0;i<3;i++)
{
int flag=0;
ll a = 1+(rand()%(n-1));
ll pwr = expo(a,d,n);
if(pwr == 1 || pwr==n-1)
continue;
for(int j=1;j<freq;j++)
{
pwr = pwer(pwr,pwr,n);
if(pwr==1)
return false;
if(pwr==n-1){
flag = 1;
break;
}
}
if(flag==1)
continue;
return false;
}
return true;
}
ll pwer(ll a,ll b,ll mod){
a%=mod;b%=mod;
long double res=a;
res*=b;
ll c=(ll)(res/mod);
a*=b;
a-=c*mod;
a%=mod;
if(a<0)a+=mod;
return a;
}
ll expo(ll a, ll b,ll mod)
{
ll result=1;
while(b)
{
if(b&1)result=pwer(result,a,mod);
a=pwer(a,a,mod);
b>>=1;
}
return result;
}