QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#286490 | #5439. Meet in the Middle | ushg8877 | WA | 120ms | 37624kb | C++14 | 4.8kb | 2023-12-17 22:46:03 | 2023-12-17 22:46:03 |
Judging History
answer
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define MP make_pair
mt19937 rnd(time(0));
const int MAXN=1e5+5;
int n,q;
struct edge{
int v,w;
edge(int _v=0,int _w=0){v=_v,w=_w;}
};
struct query{
int v,id;
query(int _v=0,int _i=0){v=_v,id=_i;}
};
struct node{
ll w,c,fr;
node(ll _w=-1e18,ll _c=-1,ll _f=0){w=_w,c=_c,fr=_f;}
};
vector<edge> edg1[MAXN],edg2[MAXN];
vector<int> vir[MAXN];
vector<query> que[MAXN];
ll dep2[MAXN],dep1[MAXN],in[MAXN],tot;
ll ex[MAXN];// 该点属性
node b1[MAXN],b2[MAXN];// 维护最优点
int siz[MAXN],maxx[MAXN],col[MAXN],son,tsiz;
ll ans[MAXN<<3];
bool vis[MAXN];
struct LCA_tech{
int a[20][MAXN<<1],lg[MAXN<<1];
int chk(int x,int y){
return dep2[x]<dep2[y]?x:y;
}
void build(){
for(int i=2;i<(MAXN<<1);i++) lg[i]=lg[i>>1]+1;
for(int i=1;i<20;i++)
for(int j=1;j+(1<<i)-1<=tot;j++)
a[i][j]=chk(a[i-1][j],a[i-1][j+(1<<i-1)]);
}
int lca(int x,int y){
x=in[x],y=in[y];
if(x>y) swap(x,y);
int k=lg[y-x+1];
return chk(a[k][x],a[k][y-(1<<k)+1]);
}
}T2;
void dfs(int u,int fa){
T2.a[0][++tot]=u;in[u]=tot;
for(auto e:edg2[u]) if(e.v!=fa){
dep2[e.v]=dep2[u]+e.w;
dfs(e.v,u);
T2.a[0][++tot]=u;
}
}// init for T2
void get_size(int u,int fa){
maxx[u]=siz[u]=1;
for(auto e:edg1[u]) if(e.v!=fa&&!vis[e.v]){
get_size(e.v,u);
maxx[u]=max(maxx[u],siz[e.v]);
siz[u]+=siz[e.v];
}
maxx[u]=max(maxx[u],tsiz-siz[u]);
if(maxx[u]<maxx[son]) son=u;
}
void build_vir(vector<int> a){
static int st[MAXN],top;
st[top=1]=1;
sort(a.begin(),a.end(),[&](int x,int y){return in[x]<in[y];});
a.erase(unique(a.begin(),a.end()),a.end());
auto link=[&](int u,int v){
vir[u].push_back(v);
vir[v].push_back(u);
};
for(int u:a) if(u!=1){
int p=T2.lca(u,st[top]);
while(st[top]!=p){
if(dep2[st[top-1]]>=dep2[p]) link(st[top],st[top-1]),top--;
else{
link(st[top],p);
st[top]=p;
}
}
st[++top]=u;
}
while(top>1) link(st[top-1],st[top]),top--;
}
void dfz(int r){
vis[r]=true;
vector<int> nd;
auto dfs=[&](auto self,int u,int fa)->void{
nd.push_back(u);ex[u]=dep1[u];
for(query q:que[u]) nd.push_back(q.v);
for(edge e:edg1[u]) if(e.v!=fa&&!vis[e.v]){
dep1[e.v]=dep1[u]+e.w;
col[e.v]=(u==r?e.v:col[u]);
self(self,e.v,u);
}
};dep1[r]=0;dfs(dfs,r,0);col[r]=r;
// cout<<"BUILD VIR: ";
// for(int j:nd) cout<<j<<' ';cout<<endl;
build_vir(nd);
// cout<<"VIRTURE TREE "<<r<<endl;
// for(int i=1;i<=n;i++){
// for(int j:vir[i]) if(i<j) cout<<"EDGE "<<i<<' '<<j<<endl;
// }
// cout<<"VALUE: "<<ex[1]<<' '<<ex[2]<<' '<<dep1[1]<<' '<<dep1[2]<<endl;
auto push=[&](int u,node b){
if(b1[u].fr==b.fr||b2[u].fr==b.fr) return;
if(b1[u].w<b.w){
if(b1[u].c!=b.c) b2[u]=b1[u];
b1[u]=b;
}else if(b1[u].c!=b.c&&b2[u].w<b.w) b2[u]=b;
};
auto dfs1=[&](auto self,int u,int fa)->void{
b1[u]=node(ex[u],col[u],u);b2[u]=node();ex[u]=-1e18;
for(int v:vir[u])if(v!=fa){
self(self,v,u);
push(u,node(b1[v].w+dep2[v]-dep2[u],b1[v].c,b1[v].fr));
push(u,node(b2[v].w+dep2[v]-dep2[u],b2[v].c,b2[v].fr));
}
};dfs1(dfs1,1,0);// 自下而上
// cout<<"FIRST"<<endl;
// cout<<b1[1].w<<' '<<b1[1].c<<endl;
// cout<<b1[2].w<<' '<<b1[2].c<<' '<<dep1[1]<<endl;
// cout<<col[1]<<' '<<col[2]<<' '<<dep1[r]<<endl;
auto dfs2=[&](auto self,int u,int fa)->void{
for(int v:vir[u])if(v!=fa){
push(v,node(b1[u].w+dep2[v]-dep2[u],b1[u].c,b1[u].fr));
push(v,node(b2[u].w+dep2[v]-dep2[u],b2[u].c,b2[u].fr));
self(self,v,u);
}
vir[u].clear();
};dfs2(dfs2,1,0);// 自上而下
// cout<<b1[1].w<<' '<<b1[1].c<<endl;
// cout<<b1[2].w<<' '<<b1[2].c<<' '<<dep1[1]<<endl;
// cout<<col[1]<<' '<<col[2]<<' '<<dep1[r]<<endl;
auto dfs3=[&](auto self,int u,int fa)->void{
for(query q:que[u]){
// cout<<"GO! "<<u<<' '<<q.v<<endl;
if(u==r||b1[q.v].c!=col[u]) ans[q.id]=max(ans[q.id],dep1[u]+b1[q.v].w);
if(u==r||b2[q.v].c!=col[u]) ans[q.id]=max(ans[q.id],dep1[u]+b2[q.v].w);
}
for(edge e:edg1[u]) if(e.v!=fa&&!vis[e.v]) self(self,e.v,u);
};dfs3(dfs3,r,0);
for(edge e:edg1[r]) if(!vis[e.v]){
tsiz=siz[e.v];son=0;
get_size(e.v,0);get_size(son,0);
dfz(son);
}
}
int main(){
ios::sync_with_stdio(false);
memset(ex,-0x3f,sizeof(ex));
cin>>n>>q;
for(int i=1;i<n;i++){
int u,v,w;cin>>u>>v>>w;
edg1[u].push_back(edge(v,w));
edg1[v].push_back(edge(u,w));
}
for(int i=1;i<n;i++){
int u,v,w;cin>>u>>v>>w;
edg2[u].push_back(edge(v,w));
edg2[v].push_back(edge(u,w));
}
for(int i=1;i<=q;i++){
int u,v;cin>>u>>v;
que[u].push_back(query(v,i));
}
dfs(1,0);T2.build();
// cout<<"DEEP: "<<dep2[1]<<' '<<dep2[2]<<endl;
maxx[0]=1e9;tsiz=n;
get_size(1,0);get_size(son,0);dfz(son);
for(int i=1;i<=q;i++) cout<<ans[i]<<endl;
}
/*
2 4
1 2 1
1 2 3
1 1
1 2
2 1
2 2
*/
詳細信息
Test #1:
score: 100
Accepted
time: 0ms
memory: 27048kb
input:
3 4 1 2 1 2 3 2 1 2 2 2 3 1 1 1 1 2 2 1 2 2
output:
6 4 5 3
result:
ok 4 number(s): "6 4 5 3"
Test #2:
score: 0
Accepted
time: 2ms
memory: 26524kb
input:
2 1 1 2 1 1 2 1 1 1
output:
2
result:
ok 1 number(s): "2"
Test #3:
score: 0
Accepted
time: 2ms
memory: 26968kb
input:
2 1 1 2 1 1 2 1 1 2
output:
1
result:
ok 1 number(s): "1"
Test #4:
score: -100
Wrong Answer
time: 120ms
memory: 37624kb
input:
10000 50000 8101 5996 108427744 5996 7605 870838849 5996 5599 603303696 8101 3006 339377417 8101 6463 442516687 6463 5560 109174079 5560 4063 127596224 3006 1682 947915262 5996 1986 130416311 6463 5455 279771516 6463 2634 516688796 4063 3034 217886863 7605 5092 742375061 5599 2266 193804402 5092 140...
output:
647838384844 362326016102 387098802897 637066393138 456072330399 645842915960 641537859258 400174767683 459407279331 573316121269 674370549986 434572336505 435282975629 620172136910 553499885160 526743824759 471341504556 689550247042 549161302822 498449365451 396706495314 393003299495 701575393972 6...
result:
wrong answer 2nd numbers differ - expected: '626539793176', found: '362326016102'