社区讨论
97pts代码求调
CF487ETourists参与者 8已保存回复 9
讨论操作
快速查看讨论及其快照的属性,并进行相关操作。
- 当前回复
- 9 条
- 当前快照
- 1 份
- 快照标识符
- @locvw5pc
- 此快照首次捕获于
- 2023/10/30 20:36 2 年前
- 此快照最后确认于
- 2023/11/02 10:38 2 年前
#37WA了
我的思路是圆方树+树链剖分+线段树,各位dalao拜托了
CPP我的思路是圆方树+树链剖分+线段树,各位dalao拜托了
#include <bits/stdc++.h>
#define int long long
#define MAXN 200010
using namespace std;
const int INF = 0x3f3f3f3f3f3f3f3f;
int n,m,q;
multiset<int>tran[200010];
struct Edge{
int nxt,to;
}e[MAXN<<1],edge[MAXN<<1];
int h[MAXN],head[MAXN];
int tt,tot;
int d[MAXN],val[MAXN];
inline void add_edge(int u,int v){
edge[++tot].to=v;
edge[tot].nxt=head[u];
head[u]=tot;
}
inline void add(int u,int v){
e[++tt].to=v;
e[tt].nxt=h[u];
h[u]=tt;
}
int low[MAXN],dfn[MAXN],num;
int stk[MAXN],tp,col;
void Tarjan(int x){
low[x]=dfn[x]=++num;
stk[++tp]=x;
for(int i=h[x];i;i=e[i].nxt){
int v=e[i].to;
if(!dfn[v]){
Tarjan(v);
low[x]=min(low[v],low[x]);
if(low[v]==dfn[x]){
col++;
while(stk[tp]!=x){
add_edge(n+col,stk[tp]);
add_edge(stk[tp],n+col);
tp--;
}
add_edge(n+col,stk[tp]);
add_edge(stk[tp],n+col);
}
}
else low[x]=min(low[x],dfn[v]);
}
}
int fa[MAXN],siz[MAXN],dep[MAXN],son[MAXN];
void dfs_1(int now,int fath){
fa[now]=fath;
siz[now]=1;
son[now]=0;
dep[now]=dep[fath]+1;
for(int i=head[now];i;i=edge[i].nxt){
int v=edge[i].to;
if(v==fath) continue;
dfs_1(v,now);
siz[now]+=siz[v];
if(siz[son[now]]<siz[v]) son[now]=v;
}
}
int top[MAXN],per[MAXN],cnt;
void dfs_2(int now,int topx){
top[now]=topx;
per[now]=++cnt;
if(son[now]) dfs_2(son[now],topx);
else return;
for(int i=head[now];i;i=edge[i].nxt){
int v=edge[i].to;
if(v==fa[now]||v==son[now]) continue;
dfs_2(v,v);
}
}
void prework(){
dfs_1(1,0);
dfs_2(1,1);
}
struct SG{
int l,r;
int minn;
}t[MAXN<<2];
void build(int p,int l,int r){
t[p].l=l;t[p].r=r;
if(l==r){
t[p].minn=val[l];
return;
}
int mid=(l+r)>>1;
build(p*2,l,mid);
build(p*2+1,mid+1,r);
t[p].minn=min(t[p*2].minn,t[p*2+1].minn);
}
void change(int p,int pos,int w){
if(t[p].l==t[p].r){
t[p].minn=w;
return;
}
int mid=(t[p].l+t[p].r)>>1;
if(mid>=pos) change(p*2,pos,w);
else change(p*2+1,pos,w);
t[p].minn=min(t[p*2].minn,t[p*2+1].minn);
}
int query(int p,int l,int r){
if(l<=t[p].l&&t[p].r<=r){
return t[p].minn;
}
int mid=(t[p].l+t[p].r)>>1;
int ans=INF;
if(mid>=l) ans=min(ans,query(p*2,l,r));
if(mid<r) ans=min(ans,query(p*2+1,l,r));
return ans;
}
inline int LCA(int x,int y){
int ans=INF;
while(top[x]!=top[y]){
if(dep[top[x]]<dep[top[y]]) swap(x,y);
ans=min(ans,query(1,per[top[x]],per[x]));
x=fa[top[x]];
}
if(dep[y]<dep[x]) swap(x,y);
ans=min(ans,query(1,per[x],per[y]));
// if(x>n) ans=min(ans,query(1,per[fa[x]],per[fa[x]]));
if(x>n) ans=min(ans,val[per[fa[x]]]);
return ans;
}
signed main(){
scanf("%lld%lld%lld",&n,&m,&q);
for(int i=1;i<=n;i++){
scanf("%lld",&d[i]);
}
for(int i=1;i<=m;i++){
int u,v;
scanf("%lld%lld",&u,&v);
add(u,v);
add(v,u);
}
for(int i=1;i<=n;i++){
if(!per[i]) Tarjan(i);
}
prework();
for(int i=1;i<=n;i++){
val[per[i]]=d[i];
if(i!=1) tran[fa[i]-n].insert(d[i]);
}
for(int i=n+1;i<=n+col;i++){
if(tran[i-n].size())val[per[i]]=(*tran[i-n].begin());
else val[per[i]]=INF;
}
build(1,1,n+col);
for(int i=1;i<=q;i++){
char c[5];
scanf("%s",c);
if(c[0]=='A'){
int a,b;
scanf("%lld%lld",&a,&b);
printf("%lld\n",LCA(a,b));
}
else if(c[0]=='C'){
int a,w;
scanf("%lld%lld",&a,&w);
change(1,per[a],w);
if(a==1){
val[per[a]]=w;
continue;
}
tran[fa[a]-n].erase(tran[fa[a]-n].find(val[per[a]]));
tran[fa[a]-n].insert(w);
int minn=(*tran[fa[a]-n].begin());
if(minn==val[per[fa[a]]]){
val[per[a]]=w;
continue;
}
change(1,per[fa[a]],minn);
val[per[fa[a]]]=minn;
val[per[a]]=w;
}
}
return 0;
}
回复
共 9 条回复,欢迎继续交流。
正在加载回复...