結果

問題 No.3226 2×2行列累乗
ユーザー tnakao0123
提出日時 2025-08-10 17:30:20
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1 ms / 2,000 ms
コード長 2,788 bytes
コンパイル時間 440 ms
コンパイル使用メモリ 45,376 KB
実行使用メモリ 7,720 KB
最終ジャッジ日時 2025-08-10 17:30:22
合計ジャッジ時間 1,675 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 27
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:122:8: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
  122 |   scanf("%d%d%d%d%d%d%lld%d", &a, &b, &c, &d, &s, &t, &n, &k);
      |   ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

ソースコード

diff # 

/* -*- coding: utf-8 -*-
 *
 * 3226.cc:  No.3226 2ﾃ・陦悟・邏ｯ荵・- yukicoder
 */

#include<cstdio>
#include<algorithm>

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 <typename T>
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<mi>;
  
template <typename T>
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<T> operator*(const Vec2<T> &v) const {
    return Vec2<T>(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<mi>;

/* 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;
}
0