Pagini recente » Cod sursa (job #658242) | Cod sursa (job #1153127) | Cod sursa (job #1356601) | Cod sursa (job #773529) | Cod sursa (job #2237245)
#include<fstream>
#include<bitset>
#define DN 7000
#define x first
#define y second
using namespace std;
ifstream fin("rsir.in");
ofstream fout("rsir.out");
int r1,r2,l,t0,t1,a,b,c,d,e,M,t2,t3,nr,l0,l1;
long long n;
bitset<DN> r[DN];
int main()
{
fin>>t0>>t1>>a>>b>>c>>d>>e>>M>>n;
l0=t0;
l1=t1;
if(n==0)
{
fout<<t0;
return 0;
}
while(1)
{
nr++;
r[t0][t1]=1;
if(nr==n)
{
fout<<t1;
return 0;
}
t2=(1LL*a*((1LL*t0*t0)%M)+1LL*b*((1LL*t1*t1)%M)+1LL*c*t0+1LL*d*t1+e)%M;
t0=t1;
t1=t2;
if(r[t0][t1]==1)
{
r1=t0;
r2=t1;
break;
}
}
t0=l0;
t1=l1;
n++;
l=nr+1;
while(1)
{
n--;
l--;
if(t0==r1&&t1==r2)
break;
t2=(1LL*a*((1LL*t0*t0)%M)+1LL*b*((1LL*t1*t1)%M)+1LL*c*t0+1LL*d*t1+e)%M;
t0=t1;
t1=t2;
}
n%=l;
if(n==0)
n=l;
nr=0;
while(1)
{
nr++;
if(nr==n)
{
fout<<t1;
return 0;
}
t2=(1LL*a*((1LL*t0*t0)%M)+1LL*b*((1LL*t1*t1)%M)+1LL*c*t0+1LL*d*t1+e)%M;
t0=t1;
t1=t2;
}
}