Pagini recente » Cod sursa (job #747692) | Cod sursa (job #856108) | Cod sursa (job #3122723) | Cod sursa (job #2882706) | Cod sursa (job #3121980)
#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
#define ar array
#define vt vector
#define pq priority_queue
#define pu push
#define pub push_back
#define em emplace
#define emb emplace_back
#define mt make_tuple
#define all(x) x.begin(), x.end()
#define allr(x) x.rbegin(), x.rend()
#define allp(x, l, r) x.begin() + l, x.begin() + r
#define len(x) (int)x.size()
#define uniq(x) unique(all(x)), x.end()
using namespace std;
using namespace __gnu_pbds;
using ll = long long;
using ld = long double;
using ull = unsigned long long;
template <class T, size_t N>
void re(array <T, N>& x);
template <class T>
void re(vt <T>& x);
template <class T>
void re(T& x) {
cin >> x;
}
template <class T, class... M>
void re(T& x, M&... args) {
re(x); re(args...);
}
template <class T>
void re(vt <T>& x) {
for(auto& it : x)
re(it);
}
template <class T, size_t N>
void re(array <T, N>& x) {
for(auto& it : x)
re(it);
}
template <class T, size_t N>
void wr(array <T, N> x);
template <class T>
void wr(vt <T> x);
template <class T>
void wr(T x) {
cout << x;
}
template <class T, class ...M> void wr(T x, M... args) {
wr(x), wr(args...);
}
template <class T>
void wr(vt <T> x) {
for(auto it : x)
wr(it, ' ');
}
template <class T, size_t N>
void wr(array <T, N> x) {
for(auto it : x)
wr(it, ' ');
}
template<class T, class... M>
auto mvt(size_t n, M&&... args) {
if constexpr(sizeof...(args) == 1)
return vector<T>(n, args...);
else
return vector(n, mvt<T>(args...));
}
void set_fixed(int p = 0) {
cout << fixed << setprecision(p);
}
void set_scientific() {
cout << scientific;
}
inline void Open(const string Name) {
#ifndef ONLINE_JUDGE
(void)!freopen((Name + ".in").c_str(), "r", stdin);
(void)!freopen((Name + ".out").c_str(), "w", stdout);
#endif
}
const int mod = 666013;
struct matrix {
int mat[2][2];
matrix () {
memset(mat, 0, sizeof(mat));
}
void init() {
for (int i = 0; i < 2; ++i)
mat[i][i] = 1;
}
matrix operator * (const matrix& b) const {
matrix res;
for (int i = 0; i < 2; ++i)
for (int j = 0; j < 2; ++j)
for (int k = 0; k < 2; ++k) {
int tot = (1LL * mat[i][k] * b.mat[k][j]) % mod;
res.mat[i][j] = (res.mat[i][j] + tot) % mod;
}
return res;
}
};
void solve() {
int n; re(n);
matrix res, a;
res.init();
a.mat[0][1] = a.mat[1][0] = a.mat[1][1] = 1;
while (n > 0) {
if (n & 1) {
res = res * a;
}
a = a * a;
n >>= 1;
}
wr(res.mat[0][1]);
}
int main() {
ios_base::sync_with_stdio(false);
cin.tie(nullptr);
Open("kfib");
int t = 1;
for(;t;t--) {
solve();
}
return 0;
}