-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathRABBIT1.cpp
63 lines (63 loc) · 953 Bytes
/
RABBIT1.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
#include<bits/stdc++.h>
using namespace std;
#define ll long long int
//#define MOD 1048576
int MOD;
void multiply(ll A[2][2], ll B[2][2])
{
ll RES[2][2];
int i,j,k;
for(i=0;i<2;++i)
{
for(j=0;j<2;++j)
{
RES[i][j]=0;
for(k=0;k<2;++k)
{
RES[i][j] = (RES[i][j]%MOD + (A[i][k] * B[k][j])%MOD)%MOD;
}
}
}
for(i=0;i<2;++i)
{
for(j=0;j<2;++j)
A[i][j] = RES[i][j];
}
}
void matexp(ll A[2][2], ll p)
{
if(p==1) return;
ll RES[2][2]={1,0,0,1};
while(p)
{
if(p&1) multiply(RES,A);
multiply(A,A);
p/=2;
}
for(int i=0;i<2;++i)
{
for(int j=0;j<2;++j)
A[i][j]=RES[i][j]%MOD;
}
}
ll ans(ll n)
{
if(n<=2) return n;
ll M[2][2]={1,1,1,0};
matexp(M,n-2);
return ((2*M[0][0])%MOD + M[0][1])%MOD;
}
int main()
{
int t;
cin>>t;
while(t--)
{
int n,m;
cin>>n>>m;
MOD = 1<<m;
//cout<<MOD<<endl;
cout<<ans(n)%MOD<<endl;
}
return 0;
}