#include "highway.h"
#include <bits/stdc++.h>
#define ll long long
#define pb push_back
#define mp make_pair
using namespace std;
const int maxn = 100100;
int n, m;
vector<int>w;
vector<int>gs[maxn];
vector<pair<int,int> > g[maxn];
vector<int>gt[maxn];
int dist[maxn][2], parent[maxn][2];
bool visited[maxn];
ll dista;
map<pair<int,int>, int> ind;
void bfs(int st, int d) {
memset(visited,false,sizeof(visited));
queue<int>q;
q.push(st);
parent[st][d] = -1;
visited[st] = true;
while(!q.empty()) {
int curr = q.front();
q.pop();
for(auto i:g[curr]) {
if(!visited[i.first]) {
visited[i.first] = true;
dist[i.first][d] = dist[curr][d] + 1;
parent[i.first][d] = i.second;
q.push(i.first);
}
}
}
}
set<pair<int,int> > sta, stb;
void preprocess() {
for(int i=0;i<n;i++) {
if(dist[i][0] == dist[i][1]) continue; // skip this node
if(dist[i][0] < dist[i][1]) {
for(int j:gs[i]) {
if(dist[j][0] < dist[j][1]) {
sta.insert(mp(min(i, j), max(i, j)));
}
}
}
else {
for(int j:gs[i]) {
if(dist[j][0] > dist[j][1]) {
stb.insert(mp(min(i, j), max(i, j)));
}
}
}
}
}
int solve(set<pair<int,int> > edges, int root, int d) {
// build bfs tree
memset(visited,false,sizeof(visited));
memset(parent,-1,sizeof(parent));
for(auto i:edges) {
gt[i.first].pb(i.second);
gt[i.second].pb(i.first);
}
vector<int>v, v2;
queue<int>q;
q.push(root);
visited[root] = true;
parent[root][d] = -1;
while(!q.empty()) {
int curr = q.front();
q.pop();
v.pb(parent[curr][d]);
v2.pb(curr);
for(auto i:gt[curr]) {
if(!visited[i]) {
visited[i] = true;
parent[i][d] = curr;
dist[i][d] = dist[curr][d] + 1;
q.push(i);
}
}
}
if(v.size() == 1) return v[0];
// perform binary search
for(int i=0;i<v.size();i++) {
if(v[i] == -1) continue; // root
int cind = ind[mp(v[i], v2[i])];
w[cind] = 1;
}
int li = 0, ri = v.size() - 1;
while(li < ri) {
int mid = (li + ri) / 2;
// all before mid put them as low cost and other as high cost
for(int i=0;i<=mid;i++) {
if(v[i] == -1) continue; // root
int cind = ind[mp(v[i], v2[i])];
w[cind] = 0;
}
ll cdist = ask(w);
for(int i=0;i<=mid;i++) {
if(v[i] == -1) continue; // root
int cind = ind[mp(v[i], v2[i])];
w[cind] = 1;
}
if(cdist == dista) {
ri = mid;
}
else {
li = mid + 1;
}
}
for(int i=0;i<v.size();i++) {
if(v[i] == -1) continue; // root
int cind = ind[mp(v[i], v2[i])];
w[cind] = 0;
}
return v2[li];
}
void find_pair(int N, std::vector<int> U, std::vector<int> V, int A, int B) {
n = N; m = U.size();
int li=0, ri=m-1;
for(int i=0;i<m;i++) {
ind[mp(U[i], V[i])] = ind[mp(V[i], U[i])] = i;
gs[U[i]].pb(V[i]);
gs[V[i]].pb(U[i]);
g[U[i]].pb(mp(V[i], i));
g[V[i]].pb(mp(U[i], i));
}
for(int i=0;i<m;i++) w.pb(0);
dista = ask(w);
while(li<ri) {
int mid = (li+ri)/2;
for(int i=li;i<=mid;i++) {
w[i] = 1;
}
ll distb = ask(w);
for(int i=li;i<=mid;i++) {
w[i] = 0;
}
if(dista == distb) {
li = mid + 1;
}
else {
ri = mid;
}
}
int x = li;
int u = U[x], v = V[x];
//cout<<u<<" "<<v<<"\n";
bfs(u, 0);
bfs(v, 1);
preprocess();
int xx = solve(sta, u, 0);
int yy = solve(stb, v, 1);
//cout<<"Ans: "<<xx<<" "<<yy<<"\n";
answer(xx, yy);
}
Compilation message
highway.cpp: In function 'int solve(std::set<std::pair<int, int> >, int, int)':
highway.cpp:106:18: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for(int i=0;i<v.size();i++) {
~^~~~~~~~~
highway.cpp:139:18: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for(int i=0;i<v.size();i++) {
~^~~~~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
10 ms |
8184 KB |
Output is correct |
2 |
Correct |
9 ms |
8184 KB |
Output is correct |
3 |
Incorrect |
10 ms |
8232 KB |
Output is incorrect: {s, t} is wrong. |
4 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
13 ms |
8612 KB |
Output is incorrect: {s, t} is wrong. |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
49 ms |
11748 KB |
Output is correct |
2 |
Correct |
90 ms |
15004 KB |
Output is correct |
3 |
Correct |
126 ms |
18456 KB |
Output is correct |
4 |
Correct |
502 ms |
38020 KB |
Output is correct |
5 |
Correct |
430 ms |
38000 KB |
Output is correct |
6 |
Correct |
769 ms |
39160 KB |
Output is correct |
7 |
Correct |
423 ms |
39680 KB |
Output is correct |
8 |
Correct |
633 ms |
38344 KB |
Output is correct |
9 |
Correct |
401 ms |
38052 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
13 ms |
8696 KB |
Output is correct |
2 |
Correct |
77 ms |
11640 KB |
Output is correct |
3 |
Correct |
558 ms |
33684 KB |
Output is correct |
4 |
Correct |
784 ms |
40920 KB |
Output is correct |
5 |
Correct |
702 ms |
41048 KB |
Output is correct |
6 |
Incorrect |
704 ms |
40924 KB |
Output is incorrect: {s, t} is wrong. |
7 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
84 ms |
11844 KB |
Output is correct |
2 |
Correct |
71 ms |
12216 KB |
Output is correct |
3 |
Execution timed out |
2171 ms |
42804 KB |
Time limit exceeded |
4 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
72 ms |
11776 KB |
Output is correct |
2 |
Correct |
114 ms |
12356 KB |
Output is correct |
3 |
Correct |
765 ms |
39680 KB |
Output is correct |
4 |
Incorrect |
1188 ms |
42964 KB |
Output is incorrect: {s, t} is wrong. |
5 |
Halted |
0 ms |
0 KB |
- |