Pagini recente » Cod sursa (job #2370644) | Cod sursa (job #1717384) | Cod sursa (job #234278) | Cod sursa (job #600327) | Cod sursa (job #3135237)
#include <fstream>
#define MOD 1999999973
using namespace std;
ifstream fin("lgput.in");
ofstream fout("lgput.out");
int n,p;
int sol[3][3],put[3][3];
void pwr(int a[3][3],int b[3][3])
{
long long c[3][3];
for(int i=1;i<=2;i++)
{
for(int j=1;j<=2;j++)
{
c[i][j]=0;
}
}
for(int i = 1; i <=2; i++)
for(int j = 1; j <=2; j++)
{
for(int k = 1; k <=2; k++)
c[i][j] = (c[i][j] + 1LL * a[i][k] * b[k][j])%MOD;
}
for(int i=1;i<=2;i++)
{
for(int j=1;j<=2;j++)
{
a[i][j]=c[i][j];
}
}
}
void solve()
{
sol[1][1]=sol[2][2]=1;
put[1][1]=n;
put[2][2]=1;
while(p)
{
if(p % 2)
{
pwr(sol,put);
}
pwr(put,put);
p = p/ 2;
}
fout<<(sol[1][1]%MOD)<<" ";
}
int main()
{
fin>>n>>p;
solve();
return 0;
}