Pagini recente » Cod sursa (job #1138886) | Cod sursa (job #1566169) | Cod sursa (job #2058323) | Cod sursa (job #1640875) | Cod sursa (job #1118348)
#include <fstream>
#include <iostream>
#include <vector>
#include <bitset>
#include <string.h>
#include <algorithm>
#include <iomanip>
#include <math.h>
#include <time.h>
#include <stdlib.h>
#include <set>
#include <map>
#include <string>
#include <queue>
#include <deque>
using namespace std;
const char infile[] = "kfib.in";
const char outfile[] = "kfib.out";
ifstream fin(infile);
ofstream fout(outfile);
const int MAXN = 100005;
const int MOD = 666013;
const int oo = 0x3f3f3f3f;
typedef vector<int> Graph[MAXN];
typedef vector<int> :: iterator It;
const inline int min(const int &a, const int &b) { if( a > b ) return b; return a; }
const inline int max(const int &a, const int &b) { if( a < b ) return b; return a; }
const inline void Get_min(int &a, const int b) { if( a > b ) a = b; }
const inline void Get_max(int &a, const int b) { if( a < b ) a = b; }
int fibo[3][3], Ans[3][3];
inline void Multiply(int A[3][3] , int B[3][3], int C[3][3]) {
for(int i = 0 ; i < 2 ; ++ i)
for(int j = 0 ; j < 2 ; ++ j)
for(int k = 0 ; k < 2 ; ++ k)
C[i][j] = (C[i][j] + 1LL * A[i][k] * B[k][j]) % MOD;
}
void lgPow(int P, int M[3][3]) {
int Aux[3][3], C[3][3];
memcpy(C, fibo, sizeof(fibo));
M[0][0] = M[1][1] = 1;
for( ; P ; P >>= 1) {
if(P & 1) {
memset(Aux, 0, sizeof(Aux));
Multiply(M, C, Aux);
memcpy(M, Aux, sizeof(Aux));
}
memset(Aux, 0, sizeof(Aux));
Multiply(C, C, Aux);
memcpy(C, Aux, sizeof(Aux));
}
}
int main() {
int N;
fin >> N;
fibo[0][0] = 0;
fibo[0][1] = 1;
fibo[1][0] = 1;
fibo[1][1] = 1;
lgPow(N - 1, Ans);
fout << Ans[1][1] << '\n';
fin.close();
fout.close();
return 0;
}