-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathNEWCURR.c
92 lines (92 loc) · 1.31 KB
/
NEWCURR.c
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
89
90
91
92
#include<bits/stdc++.h>
using namespace std;
#define ll long long int
#define MOD 12345678901
ll mull(ll a, ll b)
{
ll res = 0;
while(b)
{
if(b&1) res = (res+a)%MOD;
a = (a+a)%MOD;
b /= 2;
}
return res;
}
ll modexp(ll a, ll b)
{
if(a == 1 || b == 1) return a;
ll res = 1;
while(b)
{
if(b&1) res = mull(res,a);
a = mull(a,a);
b/=2;
}
return res;
}
void mul(ll A[2][2] , ll B[2][2])
{
int i,j,k;
ll C[2][2];
for(i=0;i<2;++i)
{
for(j=0;j<2;++j)
{
C[i][j] = 0;
for(k=0;k<2;++k)
C[i][j] = (C[i][j] + mull(A[i][k],B[k][j]))%MOD;
}
}
for(i=0;i<2;++i)
{
for(j=0;j<2;++j)
A[i][j]=C[i][j];
}
}
void matexp(ll M[2][2], ll b)
{
if(n == 1) return;
ll RES[2][2] = {1,0,0,1};
while(b)
{
if(b&1) mul(RES,M);
mul(M,M);
b/=2;
}
int i,j;
for(i=0;i<2;++i)
{
for(j=0;j<2;++j)
M[i][j] = RES[i][j];
}
}
ll ans(ll n)
{
if(n <= 2) return n;
ll M[2][2] = {1,1,1,0};
matexp(M,n-1);
ll ab = (M[0][0]*2 + M[0][1])%MOD;
if(ab < 0) ab+=MOD;
return ab;
}
int main()
{
/*ll A[2][2]={1,2,0,1};
matexp(A,2);
int i,j;
for(i=0;i<2;++i)
{
for(j=0;j<2;++j)
cout<<A[i][j]<<" ";
cout<<endl;
}*/
int t;
cin>>t;
ll n;
while(t--)
{
scanf("%lld",&n);
printf("%lld\n",ans(n));
}
}