| # | Time | Username | Problem | Language | Result | Execution time | Memory |
|---|---|---|---|---|---|---|---|
| 669953 | alvingogo | Towns (IOI15_towns) | C++14 | 0 ms | 0 KiB |
This submission is migrated from previous version of oj.uz, which used different machine for grading. This submission may have different result if resubmitted.
#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;
map<pair<long long,long long>,long long> mp;
long long qt=0;
long long f(long long x,long long 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);
}
long long hubDistance(int n,int sub){
qt=0;
mp.clear();
long long a=0;
for(long long i=0;i<n;i++){
if(f(0,i)>f(0,a)){
a=i;
}
}
long long b=a;
for(long long i=0;i<n;i++){
if(f(a,i)>f(a,b)){
b=i;
}
}
assert(a!=b);
map<long long,pair<long long,long long> > m;
long long dd=f(a,0)-((f(a,0)+f(b,0)-f(a,b))/2);
m[dd].fs=max(m[dd].fs,f(a,0)-dd);
for(long long i=0;i<n;i++){
if(i!=a && i!=b){
long long x=min(dd,(f(a,i)-f(0,i)+f(a,0))/2);
m[x].fs=max(m[x].fs,f(i,a)-x);
}
}
vector<long long> g;
long long cz=0;
for(auto &h:m){
g.push_back(h.fs);
h.sc.sc=cz;
cz++;
}
long long p=g.size();
vector<long long> dis(p,0);
for(long long i=0;i<p;i++){
dis[i]=max(dis[i],max(g[i],f(a,b)-g[i]));
for(long long j=i+1;j<p;j++){
long long 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);
}
}
long long y=0;
for(long long i=0;i<p;i++){
if(dis[i]<dis[y]){
y=i;
}
}
vector<long long> dis2(n);
for(long long i=0;i<n;i++){
long long x=min(dd,(f(a,i)-f(0,i)+f(a,0))/2);
dis2[i]=f(a,i)-x+abs(g[y]-x);
//cout << dis2[i] << " ";
}
//cout << '\n';
//assert(qt<=2*n);
//cout << qt << endl;
long long at=1;
vector<long long> vis(n),pp(n);
long long cnt=0,lo=0;
for(long long i=0;i<n;i++){
long long u=0;
for(long long j=0;j<n;j++){
if(f(i,j)!=dis2[i]+dis2[j]){
u++;
}
}
if(u>=n/2){
at=-1;
}
//cout << i << " " << u << "\n";
}
/*
for(long long 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--;
}
}
}
long long cc=0;
for(long long i=0;i<n;i++){
//if(!vis[i]){
long long q=(f(lo,i)!=dis2[i]+dis2[lo]);
if(q){
cc+=pp[i]+1;
}
//}
}
if(cc>n/2){
at=-1;
}
*/
//assert(2*qt<=7*n);
return dis[y]*at;
}