結果

問題 No.3226 2×2行列累乗
コンテスト
ユーザー ID 21712
提出日時 2026-05-02 14:51:31
言語 Go
(1.26.1)
コンパイル:
env GOCACHE=/tmp go build _filename_
実行:
./Main
結果
AC  
実行時間 1 ms / 2,000 ms
コード長 1,269 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 14,056 ms
コンパイル使用メモリ 281,032 KB
実行使用メモリ 6,400 KB
最終ジャッジ日時 2026-05-02 14:51:47
合計ジャッジ時間 14,829 ms
ジャッジサーバーID
（参考情報） 		judge3_1 / judge2_1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 27
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

package main

import . "fmt"
import . "math/big"

type Modulo struct {
	k *Int
}

// Goの剰余演算子 % は負を割る場合…
// big.IntにはDiv/ModとQuo/Remの2つあり、演算子 % は後者Rem… 
// WHY!?
func (m Modulo) Modi(a int) int {
	return int(new(Int).Mod(NewInt(int64(a)), m.k).Int64())
}

type Mat struct {
	a,b int
	c,d int
}

type Vec struct {
	s int
	t int
}

func (m1 *Mat) MulMat(m2 *Mat, k Modulo) *Mat {
	return &Mat{
		a: k.Modi(k.Modi(m1.a*m2.a) + k.Modi(m1.b*m2.c)),
		b: k.Modi(k.Modi(m1.a*m2.b) + k.Modi(m1.b*m2.d)),
		c: k.Modi(k.Modi(m1.c*m2.a) + k.Modi(m1.d*m2.c)),
		d: k.Modi(k.Modi(m1.c*m2.b) + k.Modi(m1.d*m2.d)),
	}
}

func (m *Mat) MulVec(v *Vec, k Modulo) *Vec {
	return &Vec{
		s: k.Modi(k.Modi(m.a*v.s)+k.Modi(m.b*v.t)),
		t: k.Modi(k.Modi(m.c*v.s)+k.Modi(m.d*v.t)),
	}
}

func (m *Mat) Pow(n int, k Modulo) *Mat{
	if n == 0 {
		return &Mat{
			1, 0,
			0, 1,
		}
	} else if n == 1 {
		return m
	} else if n % 2 == 0 {
		dm := m.MulMat(m, k)
		return dm.Pow(n/2, k)
	} else {
		dm := m.MulMat(m, k)
		return dm.Pow(n/2, k).MulMat(m, k)
	}
}

func main() {
	var m Mat
	var v Vec
	var n int
	var k Int
	Scan(&m.a,&m.b,&m.c,&m.d,&v.s,&v.t,&n,&k)
	ans := m.Pow(n, Modulo{&k}).MulVec(&v, Modulo{&k})
	Println(ans.s, ans.t)
}
0