Pagini recente » Cod sursa (job #572614) | Cod sursa (job #2171997) | Cod sursa (job #1323692) | Cod sursa (job #861462) | Cod sursa (job #68798)
Cod sursa(job #68798)
#include <cstdio>
#include <cmath>
using namespace std;
long long n;
unsigned t0, t1, a, b, x, y, z, MOD, v1[8000], v2[8000];
void modif(unsigned &t0, unsigned &t1)
{
unsigned old = t1;
t1 = (v1[t0] + v2[t1] + z);
if (t1>=MOD) t1-=MOD;
if (t1>=MOD) t1-=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++){
v1[i] = (i*i)%MOD;
v1[i] = (a*v1[i] + x*i)%MOD;
v2[i] = (i*i)%MOD;
v2[i] = (b*v2[i] + y*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;
}