結果

問題 No.3225 2×2行列相似判定 〜easy〜
ユーザー tnakao0123
提出日時 2025-08-10 17:21:26
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 233 ms / 2,000 ms
コード長 2,567 bytes
コンパイル時間 330 ms
コンパイル使用メモリ 46,112 KB
実行使用メモリ 7,716 KB
最終ジャッジ日時 2025-08-10 17:21:32
合計ジャッジ時間 5,548 ms
ジャッジサーバーID
(参考情報)
judge3 / judge4
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 33
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In static member function ‘static Mat2<T> Mat2<T>::read() [with T = MI<67>]’:
main.cpp:64:10: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   64 |     scanf("%d%d%d%d", &w, &x, &y, &z);
      |     ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~

ソースコード

diff #

/* -*- coding: utf-8 -*-
 *
 * 3225.cc:  No.3225 2テ・陦悟・逶ク莨シ蛻、螳・縲彳asy縲・- yukicoder
 */

#include<cstdio>
#include<algorithm>

using namespace std;

/* constant */

const int MOD = 67;

/* typedef */

template<const int MOD>
struct MI {
  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); }
};

using mi = MI<MOD>;

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);
  }

  Mat2 operator^(int 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() {
  mat2 ma = mat2::read();
  mat2 mb = mat2::read();

  for (int pa = 0; pa < MOD; pa++)
    for (int pb = 0; pb < MOD; pb++)
      for (int pc = 0; pc < MOD; pc++)
	for (int pd = 0; pd < MOD; pd++) {
	  mat2 mp(pa, pb, pc, pd);
	  if (mp.det() != 0 && mp * ma == mb * mp) {
	    puts("Yes");
	    return 0;
	  }
	}
	  
  puts("No");
  return 0;
}
0