| # | 제출 시각 | 아이디 | 문제 | 언어 | 결과 | 실행 시간 | 메모리 |
|---|---|---|---|---|---|---|---|
| 670062 | alvingogo | 도시들 (IOI15_towns) | C++14 | 0 ms | 0 KiB |
이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include <bits/stdc++.h>
//#include <towns.h>
#pragma GCC optimize("Ofast")
#define AquA cin.tie(0);ios_base::sync_with_stdio(0);
#define fs first
#define sc second
#define p_q priority_queue
using namespace std;
int nn;
vector<vector<int> > dis;
int getDistance(int x,int y){
assert(x>=0 && x<nn);
assert(y>=0 && y<nn);
return dis[x][y];
}
map<pair<int,int>,int> mp;
int qt=0;
int f(int x,int y){
if(x==y){
return 0;
}
if(x>y){
swap(x,y);
}
if(mp.find({x,y})!=mp.end()){
return mp[{x,y}];
}
qt++;
return mp[{x,y}]=getDistance(x,y);
}
int hubDistance(int n,int sub){
qt=0;
mp.clear();
int a=0;
for(int i=0;i<n;i++){
if(f(0,i)>f(0,a)){
a=i;
}
}
int b=a;
for(int i=0;i<n;i++){
if(f(a,i)>f(a,b)){
b=i;
}
}
assert(a!=b);
map<int,pair<int,vector<int> > > m;
int dd=(f(0,a)-f(0,b)+f(a,b))/2;
for(int i=0;i<n;i++){
if(i!=a){
int x=min(dd,(f(a,i)-f(0,i)+f(a,0))/2);
m[x].fs=max(m[x].fs,f(i,a)-x);
m[x].sc.push_back(i);
}
}
vector<int> g;
int cz=0;
for(auto &h:m){
g.push_back(h.fs);
cz++;
}
int p=g.size();
vector<int> dis(p,0);
for(int i=0;i<p;i++){
dis[i]=max(dis[i],max(g[i],f(a,b)-g[i]));
dis[i]=max(dis[i],m[g[i]].fs);
for(int j=i+1;j<p;j++){
int u=g[j]-g[i];
dis[i]=max(dis[i],m[g[j]].fs+u);
dis[j]=max(dis[j],m[g[i]].fs+u);
}
}
int ans=2e9;
for(int i=0;i<p;i++){
ans=min(ans,dis[i]);
}
return ans;
int at=-1;
auto cal=[&](int y){
vector<int> dis2(n);
for(int i=0;i<n;i++){
int x=min(dd,(f(a,i)-f(0,i)+f(a,0))/2);
dis2[i]=f(a,i)-x+abs(g[y]-x);
}
vector<int> vis(n),pp(n);
int cnt=0,lo=0;
for(int i=0;i<n;i++){
if(cnt==0){
lo=i;
cnt++;
}
else{
if(f(lo,i)!=dis2[i]+dis2[lo]){
cnt++;
vis[i]=1;
pp[lo]++;
}
else{
cnt--;
}
}
}
int cc=0;
for(int i=0;i<n;i++){
if(!vis[i]){
int q=(f(lo,i)!=dis2[i]+dis2[lo]);
if(q){
cc+=pp[i]+1;
}
}
}
if(cc<=n/2){
at=1;
}
};
int nw=1;
int rr=0;
for(int i=0;i<p;i++){
int c=m[g[i]].sc.size();
int u=n-c-nw;
if(nw<=n/2 && u<=n/2 && dis[i]==ans){
rr++;
cal(i);
}
nw+=c;
}
assert(rr<=1);
//assert(2*qt<=7*n);
return ans*at;
}
int main(){
cin >> nn;
dis.resize(nn,vector<int>(nn));
for(int i=0;i<nn;i++){
for(int j=0;j<nn;j++){
cin >> dis[i][j];
}
}
cout << hubDistance(nn,1) << "\n";
return 0;
}