#pragma GCC optimize("Ofast")
#include <bits/stdc++.h>
using namespace std;
#define int long long
vector <int> edges[407];
int dist[407][407];
int vis[407];
int makod[407];
int zysk[407][407];
queue <int> q;
int d(int u, int v){
return dist[u][v];
}
int32_t main()
{
int n, i, j, a, b, m, c, v, t, l, r, s, k, w, u, req, met;
char zn;
cin >> t;
while(t--){
cin >> n;
for(i=0; i<=n; i++){
for(j=0; j<=n; j++){
zysk[i][j]=0;
dist[i][j]=0;
}
makod[i]=0;
vis[i]=0;
edges[i].clear();
}
for(i=1; i<=n; i++){
for(j=1; j<=n; j++){
cin >> zn;
if(zn-'0'){
edges[i].push_back(j);
edges[j].push_back(i);
}
}
}
// getting distances
l=0;
r=0;
for(i=1; i<=n; i++){
q.push(i);
vis[i]=i;
dist[i][i]=0;
while(!q.empty()){
v=q.front();
q.pop();
for(auto e : edges[v]){
if(vis[e]==i) continue;
q.push(e);
vis[e]=i;
dist[i][e]=dist[i][v]+1;
}
}
for(j=1; j<=n; j++){
makod[i]=max(makod[i], dist[i][j]);
}
r=max(r, makod[i]);
}
r++;
// binsearch
while(l<r){
k=(l+r)/2;
//k=l;
//cout << "########### " << k << '\n';
bool exists=false;
// v is one of the ports
for(v=1; v<=n; v++){
req=0;
for(i=0; i<=n; i++){
for(j=0; j<=makod[i]; j++){
zysk[i][j]=0;
}
}
for(u=1; u<=n; u++){
for(w=u+1; w<=n; w++){
if(d(u,w)<=k) continue;
req++;
if(d(u,v)<=k){
zysk[w][min(k-d(u,v), makod[w])]++;
}
if(d(w,v)<=k){
zysk[u][min(k-d(w,v), makod[u])]++;
}
}
}
for(i=1; i<=n; i++){
/*for(j=0; j<=makod[i]; j++){
cout << zysk[i][j] << ' ';
}*/
for(j=makod[i]; j>=0; j--){
zysk[i][j]+=zysk[i][j+1];
}
//cout << '\n';
}
//cout << '\n';
// looking for second port
for(u=1; u<=n; u++){
met=0;
for(w=1; w<=n; w++){
met+=zysk[w][d(u,w)];
}
if(met>=req){
exists=true;
break;
}
}
if(exists) break;
}
//if(exists) break;
//l++;continue;
if(exists){
r=k;
}
else{
l=k+1;
}
}
cout << l << '\n';
}
return 0;
}
/*
1
10
0111000000
1010100000
1101010000
1010110000
0101011100
0011100010
0000100001
0000100000
0000010001
0000001010
*/
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 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 | #pragma GCC optimize("Ofast") #include <bits/stdc++.h> using namespace std; #define int long long vector <int> edges[407]; int dist[407][407]; int vis[407]; int makod[407]; int zysk[407][407]; queue <int> q; int d(int u, int v){ return dist[u][v]; } int32_t main() { int n, i, j, a, b, m, c, v, t, l, r, s, k, w, u, req, met; char zn; cin >> t; while(t--){ cin >> n; for(i=0; i<=n; i++){ for(j=0; j<=n; j++){ zysk[i][j]=0; dist[i][j]=0; } makod[i]=0; vis[i]=0; edges[i].clear(); } for(i=1; i<=n; i++){ for(j=1; j<=n; j++){ cin >> zn; if(zn-'0'){ edges[i].push_back(j); edges[j].push_back(i); } } } // getting distances l=0; r=0; for(i=1; i<=n; i++){ q.push(i); vis[i]=i; dist[i][i]=0; while(!q.empty()){ v=q.front(); q.pop(); for(auto e : edges[v]){ if(vis[e]==i) continue; q.push(e); vis[e]=i; dist[i][e]=dist[i][v]+1; } } for(j=1; j<=n; j++){ makod[i]=max(makod[i], dist[i][j]); } r=max(r, makod[i]); } r++; // binsearch while(l<r){ k=(l+r)/2; //k=l; //cout << "########### " << k << '\n'; bool exists=false; // v is one of the ports for(v=1; v<=n; v++){ req=0; for(i=0; i<=n; i++){ for(j=0; j<=makod[i]; j++){ zysk[i][j]=0; } } for(u=1; u<=n; u++){ for(w=u+1; w<=n; w++){ if(d(u,w)<=k) continue; req++; if(d(u,v)<=k){ zysk[w][min(k-d(u,v), makod[w])]++; } if(d(w,v)<=k){ zysk[u][min(k-d(w,v), makod[u])]++; } } } for(i=1; i<=n; i++){ /*for(j=0; j<=makod[i]; j++){ cout << zysk[i][j] << ' '; }*/ for(j=makod[i]; j>=0; j--){ zysk[i][j]+=zysk[i][j+1]; } //cout << '\n'; } //cout << '\n'; // looking for second port for(u=1; u<=n; u++){ met=0; for(w=1; w<=n; w++){ met+=zysk[w][d(u,w)]; } if(met>=req){ exists=true; break; } } if(exists) break; } //if(exists) break; //l++;continue; if(exists){ r=k; } else{ l=k+1; } } cout << l << '\n'; } return 0; } /* 1 10 0111000000 1010100000 1101010000 1010110000 0101011100 0011100010 0000100001 0000100000 0000010001 0000001010 */ |
English