#include <bits/stdc++.h>
using namespace std;
bool check(int d, int n, vector <vector<int>> &dist) {
int numFarPairs = 0;
for (int i = 0; i < n; i++) for (int j = 0; j < n; j++) {
numFarPairs += dist[i][j] > d;
}
vector <vector<int>> coveredFarPairs(n, vector <int> (n, 0));
vector <int> atDist(d + 1, 0);
for (int jTel = 0; jTel < n; jTel++) for (int i = 0; i < n; i++) {
for (int x = 0; x <= d; x++) {
atDist[x] = 0;
}
for (int j = 0; j < n; j++) if (dist[i][j] > d && dist[jTel][j] <= d) {
atDist[dist[jTel][j]]++;
}
for (int x = 1; x <= d; x++) {
atDist[x] += atDist[x - 1];
}
if (atDist[d] == 0) {
continue;
}
int iTel = 0;
while (iTel + 4 <= n) {
int dLeft1 = d - dist[i][iTel];
int dLeft2 = d - dist[i][iTel + 1];
int dLeft3 = d - dist[i][iTel + 2];
int dLeft4 = d - dist[i][iTel + 3];
if (dLeft1 >= 0) coveredFarPairs[jTel][iTel] += atDist[dLeft1];
if (dLeft2 >= 0) coveredFarPairs[jTel][iTel + 1] += atDist[dLeft2];
if (dLeft3 >= 0) coveredFarPairs[jTel][iTel + 2] += atDist[dLeft3];
if (dLeft4 >= 0) coveredFarPairs[jTel][iTel + 3] += atDist[dLeft4];
iTel += 4;
}
while (iTel < n) {
int dLeft = d - dist[i][iTel];
if (dLeft >= 0) coveredFarPairs[jTel][iTel] += atDist[dLeft];
iTel++;
}
}
for (int iTel = 0; iTel < n; iTel++) for (int jTel = iTel + 1; jTel < n; jTel++) {
int covered = coveredFarPairs[iTel][jTel] + coveredFarPairs[jTel][iTel];
if (covered == numFarPairs) {
return true;
}
}
return false;
}
int solve(int n, vector <vector<int>> &g) {
vector <vector<int>> dist(n, vector <int> (n, n));
for (int i = 0; i < n; i++) {
dist[i][i] = 0;
}
for (int i = 0; i < n; i++) for (int j = 0; j < n; j++) if (g[i][j]) {
dist[i][j] = 1;
}
for (int j = 0; j < n; j++) for (int i = 0; i < n; i++) for (int k = 0; k < n; k++) {
dist[i][k] = min(dist[i][k], dist[i][j] + dist[j][k]);
}
int maxDist = 0;
for (int i = 0; i < n; i++) for (int j = 0; j < n; j++) {
maxDist = max(maxDist, dist[i][j]);
}
int l = 0, r = maxDist;
while (r - l > 1) {
int mid = (l + r) / 2;
if (check(mid, n, dist)) {
r = mid;
} else {
l = mid;
}
}
return r;
}
int main() {
ios_base::sync_with_stdio(false);
int t;
cin >> t;
while (t--) {
int n;
cin >> n;
vector <vector<int>> g(n, vector <int> (n));
for (int i = 0; i < n; i++) {
string s;
cin >> s;
for (int j = 0; j < n; j++) {
g[i][j] = s[j] == '1';
}
}
cout << solve(n, g) << '\n';
}
return 0;
}
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 | #include <bits/stdc++.h> using namespace std; bool check(int d, int n, vector <vector<int>> &dist) { int numFarPairs = 0; for (int i = 0; i < n; i++) for (int j = 0; j < n; j++) { numFarPairs += dist[i][j] > d; } vector <vector<int>> coveredFarPairs(n, vector <int> (n, 0)); vector <int> atDist(d + 1, 0); for (int jTel = 0; jTel < n; jTel++) for (int i = 0; i < n; i++) { for (int x = 0; x <= d; x++) { atDist[x] = 0; } for (int j = 0; j < n; j++) if (dist[i][j] > d && dist[jTel][j] <= d) { atDist[dist[jTel][j]]++; } for (int x = 1; x <= d; x++) { atDist[x] += atDist[x - 1]; } if (atDist[d] == 0) { continue; } int iTel = 0; while (iTel + 4 <= n) { int dLeft1 = d - dist[i][iTel]; int dLeft2 = d - dist[i][iTel + 1]; int dLeft3 = d - dist[i][iTel + 2]; int dLeft4 = d - dist[i][iTel + 3]; if (dLeft1 >= 0) coveredFarPairs[jTel][iTel] += atDist[dLeft1]; if (dLeft2 >= 0) coveredFarPairs[jTel][iTel + 1] += atDist[dLeft2]; if (dLeft3 >= 0) coveredFarPairs[jTel][iTel + 2] += atDist[dLeft3]; if (dLeft4 >= 0) coveredFarPairs[jTel][iTel + 3] += atDist[dLeft4]; iTel += 4; } while (iTel < n) { int dLeft = d - dist[i][iTel]; if (dLeft >= 0) coveredFarPairs[jTel][iTel] += atDist[dLeft]; iTel++; } } for (int iTel = 0; iTel < n; iTel++) for (int jTel = iTel + 1; jTel < n; jTel++) { int covered = coveredFarPairs[iTel][jTel] + coveredFarPairs[jTel][iTel]; if (covered == numFarPairs) { return true; } } return false; } int solve(int n, vector <vector<int>> &g) { vector <vector<int>> dist(n, vector <int> (n, n)); for (int i = 0; i < n; i++) { dist[i][i] = 0; } for (int i = 0; i < n; i++) for (int j = 0; j < n; j++) if (g[i][j]) { dist[i][j] = 1; } for (int j = 0; j < n; j++) for (int i = 0; i < n; i++) for (int k = 0; k < n; k++) { dist[i][k] = min(dist[i][k], dist[i][j] + dist[j][k]); } int maxDist = 0; for (int i = 0; i < n; i++) for (int j = 0; j < n; j++) { maxDist = max(maxDist, dist[i][j]); } int l = 0, r = maxDist; while (r - l > 1) { int mid = (l + r) / 2; if (check(mid, n, dist)) { r = mid; } else { l = mid; } } return r; } int main() { ios_base::sync_with_stdio(false); int t; cin >> t; while (t--) { int n; cin >> n; vector <vector<int>> g(n, vector <int> (n)); for (int i = 0; i < n; i++) { string s; cin >> s; for (int j = 0; j < n; j++) { g[i][j] = s[j] == '1'; } } cout << solve(n, g) << '\n'; } return 0; } |
English