結果
| 問題 |
No.2255 Determinant Sum
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2023-04-06 12:21:13 |
| 言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 33 ms / 2,000 ms |
| コード長 | 2,026 bytes |
| コンパイル時間 | 3,848 ms |
| コンパイル使用メモリ | 258,404 KB |
| 実行使用メモリ | 6,824 KB |
| 最終ジャッジ日時 | 2024-10-02 11:55:25 |
| 合計ジャッジ時間 | 4,006 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 23 |
ソースコード
#include <bits/stdc++.h>
#include <atcoder/modint>
using namespace std;
using namespace atcoder;
using mint = modint;
int main() {
ios_base::sync_with_stdio(false);
cin.tie(nullptr);
int t; cin >> t;
while (t--) {
int n, p; cin >> n >> p;
mint::set_mod(p);
vector a(n, vector<mint>(n));
vector<int> row(n), col(n);
for (int i = 0; i < n; i++) {
for (int j = 0; j < n; j++) {
int u; cin >> u;
row[i] += (u == -1);
col[j] += (u == -1);
a[i][j] = u;
}
}
int mir = *min_element(row.begin(), row.end()),
mxr = *max_element(row.begin(), row.end()),
mic = *min_element(col.begin(), col.end()),
mxc = *max_element(col.begin(), col.end());
if (mxr > 1 || mxc > 1) {
cout << "0\n";
} else {
for (int j = n - 1; j >= 0; j--) {
if (col[j]) {
for (int i = 0; i < n; i++) {
a[i].erase(a[i].begin() + j);
}
}
}
for (int i = n - 1; i >= 0; i--) {
if (row[i]) {
a.erase(a.begin() + i);
}
}
mint ans = (p == 2);
n = a.size();
for (int i = 0; i < n; i++) {
for (int j = i; j < n; j++) {
if (a[j][i] != 0) {
swap(a[i], a[j]);
break;
}
}
if ((ans *= a[i][i]) == 0) {
break;
}
for (int j = i + 1; j < n; j++) {
mint d = a[j][i] / a[i][i];
for (int k = i; k < n; k++) {
a[j][k] -= a[i][k] * d;
}
}
}
// cerr << n << '\n';
cout << ans.val() << '\n';
}
}
}