Pagini recente » Cod sursa (job #1532018) | Formatare Textile | Cod sursa (job #1856925) | Cod sursa (job #1373560) | Cod sursa (job #2193046)
#include <bits/stdc++.h>
using namespace std;
const int MAXM = 7e3;
typedef pair < int, int > PII;
int amod[MAXM], bmod[MAXM], m;
inline PII nxt(PII x) {
int aux = amod[x.first] + bmod[x.second];
if (aux >= m)
aux -= m;
return {x.second, aux};
}
int main()
{
int t0, t1, a, b, x, y, z;
long long n;
ifstream fin("rsir.in");
fin >> t0 >> t1 >> a >> b >> x >> y >> z >> m >> n;
fin.close();
t0 %= m; t1 %= m; a %= m; b %= m; x %= m; y %= m; z %= m;
for (int i = 0; i < m; ++i) {
amod[i] = (a * i % m * i % m + i * x) % m;
bmod[i] = (b * i % m * i % m + i * y + z) % m;
}
ofstream fout("rsir.out");
if (n == 0) {
fout << t0;
} else if (n == 1) {
fout << t1;
} else {
PII p1 = {t0, t1}, p2 = {t0, t1};
do {
p1 = nxt(p1);
p2 = nxt(nxt(p2));
} while (p1 != p2);
int lcyc = 0;
do {
p1 = nxt(p1);
++lcyc;
} while (p1 != p2);
p1 = p2 = {t0, t1};
for (int i = 0; i < lcyc; ++i)
p2 = nxt(p2);
int lchain = 0;
while (p1 != p2) {
++lchain;
p1 = nxt(p1);
p2 = nxt(p2);
}
if (n > lchain)
n = lchain + (n - lchain) % lcyc;
p1 = {t0, t1};
for (int i = 1; i < n; ++i)
p1 = nxt(p1);
fout << p1.second << '\n';
}
fout.close();
return 0;
}