Pagini recente » Cod sursa (job #296229) | Cod sursa (job #1648427) | Cod sursa (job #1452167) | Cod sursa (job #1976158) | Cod sursa (job #1922662)
#include <cstdio>
#include <vector>
using namespace std;
const int N=2,mod=666013;
struct matrice
{
vector<vector<int> > v;
matrice()
{
v.resize(N,vector<int>(N,0));
}
matrice operator *(const matrice &aux) const
{
matrice ret;
for(int i=0;i<N;i++)
for(int j=0;j<N;j++)
for(int k=0;k<N;k++)
ret.v[i][k]=(ret.v[i][k]+1LL*v[i][j]*aux.v[j][k])%mod;
return ret;
}
void unitate()
{
for(int i=0;i<N;i++) v[i][i]=1;
}
};
matrice rid_put(matrice a,int n)
{
matrice p;
p.unitate();
for(;n>=1;n>>=1)
{
if(n&1) p=p*a;
a=a*a;
}
return p;
}
int main()
{
freopen("kfib.in", "r", stdin);
freopen("kfib.out", "w", stdout);
int n;
scanf("%d",&n);
matrice a;
a.v[0][0]=0;
a.v[0][1]=a.v[1][0]=a.v[1][1]=1;
matrice sol=rid_put(a,n);
printf("%d",sol.v[1][0]);
return 0;
}