結果
問題 | No.1773 Love Triangle |
ユーザー | hitonanode |
提出日時 | 2021-10-06 23:47:34 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 254 ms / 2,000 ms |
コード長 | 1,460 bytes |
コンパイル時間 | 1,457 ms |
コンパイル使用メモリ | 134,816 KB |
実行使用メモリ | 6,948 KB |
最終ジャッジ日時 | 2024-07-03 20:57:27 |
合計ジャッジ時間 | 16,259 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 90 |
ソースコード
#include <chrono> #include <iostream> #include <random> #include <tuple> #include <vector> using namespace std; #include <atcoder/modint> using mint = atcoder::modint1000000007; mt19937 mt(530629); uniform_int_distribution<int> rndgen(0, mint::mod() - 1); template <class T> int rank_of_matrix(std::vector<std::vector<T>> &M) { const int N = M.size(); int rank = 0; for (int i = 0; i < N; ++i) { int ti = i; while (ti < N and M[ti][i] == 0) ti++; if (ti == N) continue; ++rank; M[i].swap(M[ti]); T inv = T(1) / M[i][i]; for (int j = i + 1; j < N; ++j) M[i][j] *= inv; for (int h = 0; h < N; ++h) { if (i == h) continue; const T c = -M[h][i]; for (int j = i + 1; j < N; ++j) M[h][j] += M[i][j] * c; } } return rank; } int solve(int N, const vector<tuple<int, int, int>> &uvws) { vector M(N, vector<mint>(N)); for (auto [u, v, w] : uvws) { mint x = rndgen(mt); M[u][v] += x; M[v][w] += x; M[w][u] += x; M[v][u] -= x; M[w][v] -= x; M[u][w] -= x; } return rank_of_matrix(M) / 2; } int main() { int N, M; cin >> N >> M; vector<tuple<int, int, int>> uvws; while (M--) { int u, v, w; cin >> u >> v >> w; u--, v--, w--; uvws.emplace_back(u, v, w); } cout << max(solve(N, uvws), solve(N, uvws)) << '\n'; }