Cod sursa(job #1458482)

Utilizator DuxarFII-Stefan-Negrus Duxar Data 7 iulie 2015 16:57:33
Problema Al k-lea termen Fibonacci Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 2.65 kb
#ifdef ONLINE_JUDGE
	#include <bits/stdc++.h>
#else
	#include <algorithm>
	#include <bitset>
	#include <cassert>
	#include <cmath>
	#include <cstdio>
	#include <cstdlib>
	#include <cstring>
	#include <iostream>
	#include <map>
	#include <set>
	#include <stack>
	#include <string>
	#include <utility>
	#include <vector>
	#include <queue>
	#include <climits>
#endif

using namespace std;

	// lambda : [] (int a, int b) -> bool { body return }
	// string r_str = R"(raw string)"

#define mp make_pair
#define mt make_tuple
#define eb emplace_back
#define pb push_back
#define LL long long
#define ULL unsigned long long
#define BASE 73
#define NMAX 10000
#define NMAX2 20001
#define MOD1 1000000007
#define MOD2 666013
#define ALL(V) (V).begin(), (V).end()
#define ALLR(V) (V).rbegin(), (V).rend()
#define CRLINE Duxar(__LINE__);
#define SHOWME(x) cerr << __LINE__ << ": " << #x << " = " << (x) << endl;
#define ENTER putchar('\n');

int dx4[] = {-1, 0, 1, 0};
int dy4[] = {0, 1, 0, -1};

int dx8[] = {-1, -1, 0, 1, 1, 1, 0, -1};
int dy8[] = {0, 1, 1, 1, 0, -1, -1, -1};

void Duxar(int _this_line) {
#ifndef ONLINE_JUDGE
	printf("\n . . . . . . . . . . . . . Passed line - %d\n", _this_line);
#endif
}

template <class T>
void ReadNo(T &_value) {
	T _sign = 1;
	char ch;
	_value = 0;
	while(!isdigit(ch = getchar())) {
		(ch == '-') && (_sign = -1);
	}
	do {
		_value = _value * 10 + (ch - '0');
	} while(isdigit(ch = getchar()));
	_value *= _sign;
}

int K;
vector <vector <int> > mat;

vector <vector <int> > MultiplyMatrix(vector <vector <int> > mat1, vector <vector <int> > mat2) {
	vector <vector <int> > ans(mat1);
	
	ans[0][0] = ((LL) mat1[0][0] * mat2[0][0] + (LL) mat1[0][1] * mat2[1][0]) % MOD2;
	ans[0][1] = ((LL) mat1[0][0] * mat2[0][1] + (LL) mat1[0][1] * mat2[1][1]) % MOD2;
	ans[1][0] = ((LL) mat1[1][0] * mat2[0][0] + (LL) mat1[1][1] * mat2[1][0]) % MOD2;
	ans[1][1] = ((LL) mat1[1][0] * mat2[0][1] + (LL) mat1[1][1] * mat2[1][1]) % MOD2;
	
	return ans;
}

vector <vector <int> > RiseMatrix(vector <vector <int> > mat, int k) {
	if (k == 1) {
		return mat;
	}
	
	vector <vector <int> > ans;
	ans = RiseMatrix(mat, k / 2);
	ans = MultiplyMatrix(ans, ans);
	if (k & 1) {
		ans = MultiplyMatrix(ans, mat);
	}
	
	return ans;
}

int main(){
	string file_input = "kfib";
#ifdef INFOARENA
	freopen((file_input + ".in").c_str(), "r", stdin);
	freopen((file_input + ".out").c_str(), "w", stdout);
#else
#ifndef ONLINE_JUDGE
	freopen("input.cpp", "r", stdin);
#endif
#endif
	
	int i;
	
	ReadNo(K);
	if (K < 2) {
		printf("%d\n", K);
		return 0;
	}
	
	mat.resize(2, vector <int> (2, 1));
	mat[0][0] = 0;

	vector <vector <int> > mat2 = RiseMatrix(mat, K - 1);
	printf("%d\n", mat2[1][1]);
	
	return 0;
}