Pagini recente » Cod sursa (job #1009272) | Cod sursa (job #1181074) | Cod sursa (job #1679873)
#include <iostream>
#include <fstream>
#define MOD 666013
using namespace std;
ifstream f("kfib.in");
ofstream g("kfib.out");
void product(long long a[][3],long long b[][3])
{
long long c[3][3]; c[0][0]=0;c[0][1]=0;c[1][0]=0;c[1][1]=0;
int i,j,k;
for(i=0;i<=1;i++)for(j=0;j<=1;j++)for(k=0;k<=1;k++) c[i][j]=(c[i][j]+a[i][k]*b[k][j])%MOD;
a[0][0]=c[0][0];a[0][1]=c[0][1];a[1][0]=c[1][0];a[1][1]=c[1][1];
}
void put(long long a[][3],long long sol[][3], long long p)
{
long long x;for(x=0;(1<<x)<=p;x++){if((1<<x)&p)product(sol,a);product(a,a);}
}
void solve()
{
long long a[3][3],sol[3][3],n;
a[1][0]=a[0][1]=a[0][0]=1;a[1][1]=0;
sol[1][1]=sol[0][0]=1;sol[1][0]=sol[0][1]=0;
f>>n;put(a,sol,n-1);
g<<sol[0][0];
}
int main()
{
solve();
return 0;
}