【Peaks加强版 BZOJ 3551】你被坑了吗？

这道在没加读入优化时间在20s左右的题终于在大米饼两天的死缠烂打、鬼混、乱整乱撞后艰难地AC了。但惋惜的是，大米饼一号代码其实更加简洁，但至今找不出BUG，我将它放在下面，也许有一天从远方来的另一个大米饼会拯救它。让我们一起念出这道题的关键字：

Kurskal，LCA倍增，Kurskal重构树，dfs序，主席树。

比较新鲜的是重构树，它有一些美妙性质。

红色的是由边化成的节点。用Kurskal建出这棵树，在这一棵树上进行DFS序处理，并且在DFS序上进行主席树（可持久化线段树）处理。有趣的是，这道题读入优化可以加速5s，是不是很美妙！

#include<stdio.h>
#include<algorithm>
#include<cstring>
#define ro(i,a,b) for(int i=a;i>=b;i--)
#define go(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int inf=2147483646;const int logp=18;const int N=1e5+233;
const int M=5e5+233;const int P=N<<1;const int E=P<<1;const int S=P*logp ;
int n,m,q,p,num,cnt,K,d[N],h[N],dad[P],val[P],rt[P],lc[S],rc[S],sum[S];
int si[P],dep[P],dfn[P],pos[P],head[P],nex[E],ver[E],fa[logp+1][P],pr[logp+1],lg[P];
struct edge{int x,y,z;};edge a[M];
bool cmp(edge a, edge b){return a.z<b.z;}
int find(int x){return dad[x]!=x?dad[x]=find(dad[x]):x;}
void ADD(int u, int v){nex[++K]=head[u];head[u]=K;ver[K]=v;}
void dfs(int u)
{
    si[u]=1;dfn[u]=++num;pos[num]=u;
    for(int i=head[u];i;i=nex[i])
    {int v=ver[i];fa[0][v]=u;dep[v]=dep[u]+1;dfs(v);si[u]+=si[v];}
}
void Kruskal()
{
    p=n;sort(a+1,a+1+m,cmp);
    go(i,1,m){int x=find(a[i].x),y=find(a[i].y);
    if(x!=y){dad[++p]=p;val[p]=a[i].z;dad[x]=dad[y]=p;
    ADD(p,x),ADD(p,y);if(p==(n<<1)-1)break;}}
}
void Rmq()
{
    go(k,1,lg[p])go(i,1,p){if(dep[i]<pr[k])continue;
    fa[k][i]=fa[k-1][fa[k-1][i]];}
}
int Add(int p,int l,int r,int x)
{
    int k=++cnt;sum[k]=sum[p]+1;if(l==r)return k;int mid=l+r>>1;
    if(x<=mid)lc[k]=Add(lc[p],l,mid,x),rc[k]=rc[p];
    else rc[k]=Add(rc[p],mid+1,r,x),lc[k]=lc[p];return k;
}
void prepare()
{
    Kruskal();num=0;go(i,1,p)if(!dfn[i])dfs(find(i));pr[0]=1;
    go(i,1,logp){pr[i]=pr[i-1]<<1;if(pr[i]>p)break;lg[pr[i]]=1;}
    go(i,1,p)lg[i]+=lg[i-1];Rmq();go(i,1,num)
    {int x=pos[i];if(x>n)rt[i]=rt[i-1];else rt[i]=Add(rt[i-1],1,n,h[x]);}
}
int Jump(int x, int v){
    int len=lg[dep[x]];ro(i,len,0)if(val[fa[i][x]]<=v)x=fa[i][x];return x;
}
int Query(int a, int b, int l, int r, int k)
{
    if(l==r)return l;int amo=sum[rc[b]]-sum[rc[a]];
    int mi=l+r>>1;if(amo<k)return Query(lc[a],lc[b],l,mi,k-amo);
    return Query(rc[a],rc[b],mi+1,r,k);
}
int Ask(int x, int v, int k)
{
    int anc=Jump(x,v),beg=rt[dfn[anc]-1],end=rt[dfn[anc]+si[anc]-1];
    if(sum[end]-sum[beg]<k)return -1;
    int hi=Query(beg,end,1,n,k);return d[hi];
}
int main()
{
    scanf("%d%d%d",&n,&m,&q);val[0]=inf;
    go(i,1,n)scanf("%d",&h[i]),d[i]=h[i],dad[i]=i;
    sort(d+1,d+1+n);go(i,1,n)h[i]=lower_bound(d+1,d+1+n,h[i])-d;
    go(i,1,m)scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].z);
    prepare();int ans=-1;while(q--)
    {
        int x,v,k;scanf("%d%d%d",&x,&v,&k);
        if(ans != -1)x^=ans,v^=ans,k^=ans;
        ans=Ask(x,v,k);printf("%d
",ans);
    }
}

下面是暂时RE的代码：

#include<stdio.h>
#include<algorithm>
#define go(i,a,b) for(int i=a;i<=b;i++)
#define ro(i,a,b) for(int i=a;i>=b;i--)
#define rule 1000000000
using namespace std;const int N=100003;struct TT{int siz,l,r;}t[N*20];
struct G{int u,v,w;bool operator<(const G&X)const{return w<X.w;};}g[N*20];
int n,m,q,fa[N],ch[N][2],SZ,sz,val[N],root,v,x,k,ans=0,T[N];
int up[N][20],dis[N][20],start[N],end[N],dfs_clock,DFS[N],ROOT[N];
int find(int x){return fa[x]==x?x:fa[x]=find(fa[x]);}
int get_root(int x,int val){ro(j,19,0)if(up[x][j]&&dis[x][j]<=val)x=up[x][j];return x;} 
void dfs(int u){if(!u)return;start[u]=dfs_clock;
    if(ch[u][1]+ch[u][0]<1)DFS[++dfs_clock]=val[u];
    go(i,0,1)dfs(ch[u][i]);end[u]=dfs_clock;
}
void insert(int& x,int Val,int l,int r){
    t[++sz]=t[x];x=sz;++t[x].siz;if(l==r)return;
    int mid=(l+r)/2;if(Val<=mid)insert(t[x].l,Val,l,mid);
    else insert(t[x].r,Val,mid+1,r);
}
void Que(int i,int j,int k,int l,int r){
    if( l==r ){ans=l;return;}
    int tmp=t[t[j].l].siz-t[t[i].l].siz,mid=(r+l)/2;
    if(k<=tmp)Que(t[i].l,t[j].l,k,l,mid);else Que(t[i].r,t[j].r,k-tmp,mid+1,r);
}
int main(){
    freopen("in.in","r",stdin);
    scanf("%d%d%d",&n,&m,&q);SZ=n;
    go(i,1,n)scanf("%d",&val[i]),fa[i]=i;
    go(i,1,m)scanf("%d%d%d",&g[i].u,&g[i].v,&g[i].w);
    sort(g+1,g+m+1);int tot=n;go(i,1,m)
    {
        int Fa=find(g[i].u),Fb=find(g[i].v);
        if(Fa==Fb)continue;if(!--tot)break;
        int u=++SZ;fa[u]=fa[Fa]=fa[Fb]=u;
        ch[u][0]=Fa,ch[u][1]=Fb;
        up[Fa][0]=up[Fb][0]=u;
        dis[Fa][0]=dis[Fb][0]=g[i].w;
    }
    
    dfs(SZ);go(i,1,n)ROOT[i]=ROOT[i-1],insert(ROOT[i],DFS[i],0,rule);  
    go(j,1,19)go(i,1,SZ)up[i][j]=up[up[i][j-1]][j-1],  
    dis[i][j]=max(dis[i][j-1],dis[up[i][j-1]][j-1]); 
     
    go(i,1,q)  
    {  
        scanf("%d%d%d",&v,&x,&k);x^=ans;v^=ans;k^=ans;
        root=get_root(x,v);
        int gap=end[root]-start[root];
        if(gap<k){puts("-1");ans=0;continue;}  
        Que(ROOT[start[root]],ROOT[end[root]],gap-k+1,0,rule);  
        if(!ans){puts("-1");ans=0;continue;}printf("%d
",ans);  
    } 
    return 0;
}//Paul_Guderian

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

你没有达到顶尖！但记得给我买杯柠檬汁。