Pagini recente » Cod sursa (job #2016884) | Cod sursa (job #1793092) | Cod sursa (job #63946) | Cod sursa (job #93353) | Cod sursa (job #68628)
Cod sursa(job #68628)
#include <cstdio>
#include <cmath>
using namespace std;
long long n;
unsigned t0, t1, a, b, x, y, z, MOD, sqr[10000];
void modif(unsigned &t0, unsigned &t1)
{
unsigned p0 = sqr[t0], p1 = sqr[t1], old = t1;
t1 = (a*p0 + b*p1 + x*t0 + y*t1 + z)%MOD;
t0 = old;
}
void brute(unsigned n)
{
for (unsigned i=0; i<n; i++)
modif(t0, t1);
printf("%u\n", t0);
}
int main()
{
unsigned i;
freopen("rsir.in", "r", stdin);
freopen("rsir.out", "w", stdout);
scanf("%u %u %u %u %u %u %u %u %lld",
&t0, &t1, &a, &b, &x, &y, &z, &MOD, &n);
t0 %= MOD; t1 %= MOD;
unsigned last = MOD*MOD;
for (i=0; i<MOD; i++)
sqr[i] = (i*i) % MOD;
last <?= (1<<19);
if (n <= last){
brute(n);
return 0;
}
for (i=0; i<last; i++)
modif(t0, t1);
n-=last;
unsigned T0 = t0, T1 = t1, lung = 1;
modif(t0, t1);
while (t0!=T0 || t1!=T1) modif(t0, t1), lung++;
brute(n%lung);
return 0;
}