| # | 제출 시각 | 아이디 | 문제 | 언어 | 결과 | 실행 시간 | 메모리 |
|---|---|---|---|---|---|---|---|
| 123243 | imyujin | 고속도로 설계 (CEOI12_highway) | C++14 | 0 ms | 0 KiB |
이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include "highway.h"
#include<vector>
#include<queue>
#include<bitset>
using namespace std;
#define MAXN 90005
typedef pair<int, int> pii;
int N, M, A, B;
int d, se;
vector<pii> ed[MAXN];
int da[MAXN], db[MAXN];
bitset<MAXN> chk;
int num[MAXN];
void bfs(int x, int dx[]) {
queue<int> q;
for(int i = 0; i < N; i++) dx[i] = -1;
dx[x] = 0;
q.push(x);
while(!q.empty()) {
int t = q.front();
q.pop();
for(auto a : ed[t]) if(dx[a.first] == -1) {
dx[a.first] = dx[t] + 1;
q.push(a.first);
}
}
}
int f(int x, int dx[], int dy[]) {
queue<int> q;
vector<int> con;
for(int i = 0; i < N; i++) chk[i] = false;
q.push(x);
chk[x] = true;
while(!q.empty()) {
int t = q.front();
q.pop();
con.push_back(t);
for(auto a : ed[t]) if(dx[a.first] < dy[a.first] && !chk[a.first]) {
chk[a.first] = true;
q.push(a.first);
}
}
for(int i = 0; i < N; i++) num[i] = -1;
for(int i = 0; i < con.size(); i++) num[con[i]] = i;
/*
for(auto a : con) printf("%d ", a);
printf("\n");
*/
vector<int> w(M, 0);
int s = 0, e = con.size() - 1;
int m = s + e >> 1;
for(int i = m + 1; i <= e; i++) for(auto a : ed[con[i]]) w[a.second] = 1;
while(s < e) {
//printf("s = %d, e = %d\n", s, e);
if(ask(w) == d * A) {
e = m;
m = s + e >> 1;
for(int i = m + 1; i <= e; i++) for(auto a : ed[con[i]]) w[a.second] = 1;
}
else {
s = m + 1;
m = s + e >> 1;
for(int i = s + 1; i <= m; i++) for(auto a : ed[con[i]]) if(num[a.first] < i) w[a.second] = 0;
}
}
return con[s];
}
void find_pair(int n, vector<int> U, vector<int> V, int a, int b) {
N = n;
M = U.size();
A = a;
B = b;
vector<int> w(M, 0);
for(int i = 0; i < M; i++) {
ed[U[i]].push_back(make_pair(V[i], i));
ed[V[i]].push_back(make_pair(U[i], i));
}
d = ask(w) / A;
int s = 0, e = M - 1, m = N - 1 >> 1;
for(int i = m + 1; i <= e; i++) w[i] = 1;
while(s < e) {
if(ask(w) != A * d) {
s = m + 1;
m = s + e >> 1;
for(int i = s; i <= m; i++) w[i] = 0;
}
else {
e = m;
m = s + e >> 1;
for(int i = m + 1; i <= e; i++) w[i] = 1;
}
}
se = s;
int alpha = U[se], beta = V[se];
//printf("se=%d(%d, %d)\n", se, alpha, beta);
bfs(alpha, da);
bfs(beta, db);
//for(int i = 0; i < N; i++) printf("[%d %d]\n", da[i], db[i]);
answer(f(alpha, da, db), f(beta, db, da));
}