結果
| 問題 |
No.3226 2×2行列累乗
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2025-08-08 23:07:45 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 2 ms / 2,000 ms |
| コード長 | 1,245 bytes |
| コンパイル時間 | 1,627 ms |
| コンパイル使用メモリ | 195,624 KB |
| 実行使用メモリ | 7,716 KB |
| 最終ジャッジ日時 | 2025-08-08 23:07:49 |
| 合計ジャッジ時間 | 2,856 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 27 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const ll mod = 1e9 + 7;
const int N = 200005;
const int INF = 0x3f3f3f3f;
struct Mat {
ll a, b, c, d;
};
ll cal(ll x, ll k) {
return (x % k + k) % k;
}
Mat mul(Mat A, Mat B, ll K) {
Mat C;
C.a = (ll)(((__int128)A.a * B.a + (__int128)A.b * B.c) % K);
C.b = (ll)(((__int128)A.a * B.b + (__int128)A.b * B.d) % K);
C.c = (ll)(((__int128)A.c * B.a + (__int128)A.d * B.c) % K);
C.d = (ll)(((__int128)A.c * B.b + (__int128)A.d * B.d) % K);
return C;
}
Mat powmod(Mat x, unsigned long long y, ll K) {
Mat res = {1 % K, 0, 0, 1 % K};
while (y) {
if (y & 1) res = mul(res, x, K);
x = mul(x, x, K);
y >>= 1;
}
return res;
}
int main() {
ll A, B, C, D;
ll S, T;
unsigned long long n;
ll K;
cin >> A >> B >> C >> D;
cin >> S >> T;
cin >> n >> K;
A = cal(A, K);
B = cal(B, K);
C = cal(C, K);
D = cal(D, K);
S = cal(S, K);
T = cal(T, K);
Mat M = {A, B, C, D};
Mat r = powmod({A, B, C, D}, n, K);
cout << cal((ll)(((__int128)r.a * S + (__int128)r.b * T) % K), K) << " " << cal((ll)(((__int128)r.c * S + (__int128)r.d * T) % K), K) << "\n";
return 0;
}