結果
| 問題 |
No.3226 2×2行列累乗
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2025-08-13 23:06:41 |
| 言語 | C# (.NET 8.0.404) |
| 結果 |
AC
|
| 実行時間 | 64 ms / 2,000 ms |
| コード長 | 1,981 bytes |
| コンパイル時間 | 8,826 ms |
| コンパイル使用メモリ | 168,380 KB |
| 実行使用メモリ | 187,272 KB |
| 最終ジャッジ日時 | 2025-08-13 23:06:59 |
| 合計ジャッジ時間 | 12,510 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 27 |
コンパイルメッセージ
復元対象のプロジェクトを決定しています... /home/judge/data/code/main.csproj を復元しました (121 ミリ秒)。 main -> /home/judge/data/code/bin/Release/net8.0/main.dll main -> /home/judge/data/code/bin/Release/net8.0/publish/
ソースコード
using System;
using static System.Console;
using System.Linq;
using System.Collections.Generic;
class Program
{
static int NN => int.Parse(ReadLine());
static long[] NList => ReadLine().Split().Select(long.Parse).ToArray();
static long[][] NArr(long n) => Enumerable.Repeat(0, (int)n).Select(_ => NList).ToArray();
public static void Main()
{
Solve();
}
static void Solve()
{
var m = NArr(2);
var c = NList;
var s = new long[2][] { new long[1] { c[0] }, new long[1] { c[1] } };
c = NList;
var (n, r) = (c[0], (int)c[1]);
for (var i = 0; i < 2; ++i) for (var j = 0; j < 2; ++j) m[i][j] = (m[i][j] % r + r) % r;
for (var i = 0; i < 2; ++i) s[i][0] = (s[i][0] % r + r) % r;
var ans = Matrix.Mul(Matrix.Pow(m, n, r), s, r);
WriteLine($"{ans[0][0]} {ans[1][0]}");
}
class Matrix
{
// 行列の累乗
public static long[][] Pow(long[][] m, long k, int mod)
{
var multi = m;
var r = new long[m.Length][];
for (var i = 0; i < m.Length; ++i)
{
r[i] = new long[m.Length];
r[i][i] = 1;
}
while (k > 0)
{
if ((k & 1) == 1) r = Mul(r, multi, mod);
multi = Mul(multi, multi, mod);
k >>= 1;
}
return r;
}
// 行列の積
public static long[][] Mul(long[][] x, long[][] y, int mod)
{
var r = new long[x.Length][];
for (var i = 0; i < x.Length; ++i) r[i] = new long[y[0].Length];
for (var i = 0; i < x.Length; ++i) for (var k = 0; k < x[0].Length; ++k)
{
for (var j = 0; j < y[0].Length; ++j) r[i][j] = (r[i][j] + x[i][k] * y[k][j]) % mod;
}
return r;
}
}
}