/* -*- coding: utf-8 -*- * * 3226.cc: No.3226 2テ・陦悟・邏ッ荵・- yukicoder */ #include #include using namespace std; /* constant */ /* typedef */ using ll = long long; struct MI { static int MOD; int v; MI(): v() {} MI(int _v): v(_v % MOD) { if (v < 0) v += MOD; } MI(long long _v): v(_v % MOD) { if (v < 0) v += MOD; } explicit operator int() const { return v; } MI operator+(const MI m) const { return MI(v + m.v); } MI operator-(const MI m) const { return MI(v + MOD - m.v); } MI operator-() const { return MI(MOD - v); } MI operator*(const MI m) const { return MI((long long)v * m.v); } MI &operator+=(const MI m) { return (*this = *this + m); } MI &operator-=(const MI m) { return (*this = *this - m); } MI &operator*=(const MI m) { return (*this = *this * m); } bool operator==(const MI m) const { return v == m.v; } bool operator!=(const MI m) const { return v != m.v; } MI pow(int n) const { // a^n % MOD MI pm = 1, a = *this; while (n > 0) { if (n & 1) pm *= a; a *= a; n >>= 1; } return pm; } MI inv() const { return pow(MOD - 2); } MI operator/(const MI m) const { return *this * m.inv(); } MI &operator/=(const MI m) { return (*this = *this / m); } }; int MI::MOD = 1; using mi = MI; template struct Vec2 { T a, b; Vec2(): a(), b() {}; Vec2(T _a, T _b): a(_a), b(_b) {}; static Vec2 read() { int x, y; scanf("%d%d", &x, &y); return Vec2(x, y); } }; using vec2 = Vec2; template struct Mat2 { T a, b, c, d; Mat2(): a(), b(), c(), d() {} Mat2(T _a, T _b, T _c, T _d): a(_a), b(_b), c(_c), d(_d) {} static Mat2 unit() { return Mat2(1, 0, 0, 1); } static Mat2 read() { int w, x, y, z; scanf("%d%d%d%d", &w, &x, &y, &z); return Mat2(w, x, y, z); } bool operator==(const Mat2 &e) const { return a == e.a && b == e.b && c == e.c && d == e.d; } Mat2 operator*(const Mat2 &e) const { return Mat2(a * e.a + b * e.c, a * e.b + b * e.d, c * e.a + d * e.c, c * e.b + d * e.d); } Vec2 operator*(const Vec2 &v) const { return Vec2(a * v.a + b * v.b, c * v.a + d * v.b); } Mat2 operator^(long long b) const { Mat2 p = Mat2::unit(); Mat2 a = *this; while (b > 0) { if (b & 1) p = p * a; a = a * a; b >>= 1; } return p; } T det() const { return a * d - b * c; } }; using mat2 = Mat2; /* global variables */ /* subroutines */ /* main */ int main() { int a, b, c, d, s, t, k; ll n; scanf("%d%d%d%d%d%d%lld%d", &a, &b, &c, &d, &s, &t, &n, &k); mi::MOD = k; mat2 ma(a, b, c, d); vec2 va(s, t); vec2 vb = (ma ^ n) * va; printf("%d %d\n", (int)vb.a, (int)vb.b); return 0; }