Pagini recente » Cod sursa (job #2023608) | Cod sursa (job #1515302) | Cod sursa (job #1723709) | Cod sursa (job #2277480) | Cod sursa (job #1822366)
#include <iostream>
#include <fstream>
#define MOD 666013
using namespace std;
ifstream f("kfib.in");
ofstream g("kfib.out");
int n;
long long a[4][4],s[4][4];
void inmultire(long long b[4][4],long long c[4][4])
{
int i,j,k;
long long d[4][4];
for(i=1;i<=2;i++)
for(j=1;j<=2;j++)
{
d[i][j]=0LL;
for(k=1;k<=2;k++)
d[i][j]=(0LL+d[i][j]+1LL*b[i][k]*c[k][j])%MOD;
}
for(i=1;i<=2;i++)
for(j=1;j<=2;j++)
b[i][j]=d[i][j];
}
int main()
{
f>>n;
a[1][2]=1;
a[2][1]=1;
a[2][2]=1;
s[1][1]=0;
s[1][2]=1;
n--;
while(n>0)
{
if(n%2!=0)
{
inmultire(s,a);
}
inmultire(a,a);
n=n/2;
}
g<<s[1][2];
return 0;
}