Cod sursa(job #68524)

Utilizator bogdan2412Bogdan-Cristian Tataroiu bogdan2412 Data 28 iunie 2007 12:19:36
Problema Rsir Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.37 kb
#include <cstdio>
#include <algorithm>

using namespace std;

#define MAXMOD 7005

long long N;
int T0, T1, a, b, x, y, z, MOD;

int precalc[MAXMOD];

inline int getNew( int T0, int T1 )
{
	int rez = (precalc[T0] * a) % MOD;
	rez += (precalc[T1] * b) % MOD;
	rez += (T0 * x + T1 * y + z) % MOD;
	return rez % MOD;
}

int main()
{
	freopen("rsir.in", "rt", stdin);
	freopen("rsir.out", "wt", stdout);
	scanf("%d %d %d %d %d %d %d %d %lld", &T0, &T1, &a, &b, &x, &y, &z, &MOD, &N);

	T0 %= MOD; T1 %= MOD;
	a %= MOD; b %= MOD; x %= MOD; y %= MOD; z %= MOD;

	if (N == 0)
	{
		printf("%d\n", T0);
		return 0;
	}
	if (N == 1)
	{
		printf("%d\n", T1);
		return 0;
	}

	for (int i = 0; i < MOD; i++)
		precalc[i] = (i * i) % MOD;

	pair<int, int> cur = make_pair(T0, T1), nxt, goal;
	for (int i = 2; i <= MOD * MOD + 1; cur = nxt, i++)
		nxt = make_pair( cur.second, getNew( cur.first, cur.second ) );
	goal = cur;

	cur = make_pair(T0, T1);
	for (int i = 1; i <= MOD * MOD + 1; cur = nxt, i++)
	{
		if (cur == goal)
		{
			int start = i - 1, len = MOD * MOD - i + 1;

			N -= start; N %= len;

			for (int k = 1; k <= N; cur = nxt, k++)
				nxt = make_pair( cur.second, getNew( cur.first, cur.second ) );

			printf("%d\n", cur.second);
			return 0;
		}
		nxt = make_pair( cur.second, getNew( cur.first, cur.second ) );
		if (i == N)
		{
			printf("%d\n", nxt.second);
			return 0;
		}
	}
		
	return 0;
}