Pagini recente » Cod sursa (job #1973909) | Cod sursa (job #2219068) | Cod sursa (job #902648) | Cod sursa (job #2764881) | Cod sursa (job #917679)
Cod sursa(job #917679)
#include<fstream>
#define MOD 666013
using namespace std;
ifstream f("kfib.in");
ofstream g("kfib.out");
long long a[3][3], p[3][3], aux[3][3];
int n;
void Inmulteste_p_a()
{
int i, j, k;
aux[1][1]=aux[1][2]=aux[2][1]=aux[2][2]=0;
for (i=1; i<3; ++i)
{
for (k=1; k<3; ++k)
for (j=1; j<3; ++j)
aux[i][k]=(aux[i][k]+p[i][j]*a[j][k])%MOD;
}
for (i=1; i<3; ++i)
for (j=1; j<3; ++j) p[i][j]=aux[i][j];
}
void Ridica()
{
int i, j, k;
aux[1][1]=aux[1][2]=aux[2][1]=aux[2][2]=0;
for (i=1; i<3; ++i)
{
for (k=1; k<3; ++k)
for (j=1; j<3; ++j)
aux[i][k]=(aux[i][k]+a[i][j]*a[j][k])%MOD;
}
for (i=1; i<3; ++i)
for (j=1; j<3; ++j) a[i][j]=aux[i][j];
}
int main()
{
f>>n; n-=2;
a[1][2]=a[2][1]=a[2][2]=1;
p[1][1]=p[2][2]=1;
while (n)
{
if (n%2==1)
{
Inmulteste_p_a();
--n;
}
else
{
Ridica();
n/=2;
}
}
g<<(p[1][2]+p[2][2])%MOD<<"\n";
f.close();
g.close();
return 0;
}