Pagini recente » Cod sursa (job #2673561) | Cod sursa (job #2530926) | Cod sursa (job #291020) | Cod sursa (job #1688299) | Cod sursa (job #1796975)
#include <bits/stdc++.h>
#define for0(i, n) for(int i = 0; i < n; i++)
#define for1(i, n) for(int i = 1; i <= n; i++)
#define pb push_back
#define mp make_pair
#define all(v) v.begin(), v.end()
#define V vector<int>
#define VP vector<pair<int, int> >
#define clr(A, x) memset(A, x, sizeof(A))
#define cpy(A, B) memcpy(A, B, sizeof(B))
#define g(s) getline(cin, s) ///ai grija la fin/cin ///
#define FASTIO ios_base::sync_with_stdio(0)
const long long INFLL = (1LL<<62);
const int INFINT = 2000000000;
typedef long long ll;
typedef unsigned long long ull;
using namespace std;
/*template <typename T>
string to_string(const T& n){
ostringstream os;
os << n;
return os.str();
}*/
/*void invers_modular(int a, int b, int &d, int &x, int &y)
{
if(!b)
{
d=a;
x=1;
y=0;
return ;
}
int x0, y0;
invers_modular(b, a%b, d, x0, y0);
x=y0;
y=x0-a/b*y0;
}*/ // daca x<0 se aduna cu mod pana e mai mare, x fiind rezultatul
/*ull putere(ull baza, ull exponent, ull MOD)
{
if(exponent == 0) return 1;
if(exponent % 2 == 0) return putere((baza * baza) % MOD, exponent / 2, MOD) % MOD;
return ((baza % MOD) * (putere(baza, exponent - 1, MOD) % MOD) % MOD);
}*/
ifstream fin("kfib.in"); /// modifica cu numele corespunzator
ofstream fout("kfib.out"); /// modifica cu numele corespunzator
const int MOD = 666013;
ll fibo[4], aux[4];
void putere(ll fibo[], int n)
{
if(n == 1)
fibo[1] = fibo[2] = fibo[3] = 1;
else
{
if(n % 2)
{
putere(fibo, n - 1);
aux[0] = fibo[1];
aux[1] = (fibo[0] + fibo[1]) % MOD;
aux[2] = fibo[3];
aux[3] = (fibo[2] + fibo[3]) % MOD;
for0(i, 4)
fibo[i] = aux[i];
}
else
{
putere(fibo, n / 2);
aux[0] = (fibo[0] * fibo[0] + fibo[1] * fibo[2]) % MOD;
aux[1] = (fibo[0] * fibo[1] + fibo[1] * fibo[3]) % MOD;
aux[2] = (fibo[0] * fibo[2] + fibo[2] * fibo[3]) % MOD;
aux[3] = (fibo[1] * fibo[2] + fibo[3] * fibo[3]) % MOD;
for0(i, 4)
fibo[i] = aux[i];
}
}
}
int main()
{
int k;
fin >> k;
putere(fibo, k + 1);
fout << fibo[0];
return 0;
}