【BZOJ 3545】【ONTAK 2010】Peaks & 【BZOJ 3551】【ONTAK 2010】Peaks加强版 Kruskal重构树

sunshine的A题我竟然调了一周！！！

把循环dfs改成一个dfs就可以，，，我也不知道为什么这样就不会RE，但它却是A了，，，

这周我一直在调这个题，总结一下智障错误：

1.倍增的范围设成了n而不是n*2-1，，，

2.重构树的顶点是n*2-1，而我一开始设成了n，，，

3.define里的for3和for4的i--打成i++，，，，，，，，，，，，

4.dfs爆栈了，找CA爷问的编译命令里手动扩栈，真是愚蠢的问题，，，，

比赛时绝不会有太多时间，在这么犯逗就得滚粗了QAQ

3545：

#include<cstdio> 
#include<cstring> 
#include<algorithm> 
#define for1(i,a,n) for(int i=(a);i<=(n);i++) 
#define for2(i,a,n) for(int i=(a);i<(n);i++) 
#define for3(i,a,n) for(int i=(a);i>=(n);i--) 
#define for4(i,a,n) for(int i=(a);i>(n);i--) 
#define CC(i,a) memset(i,a,sizeof(i)); 
#define read(x) x=getint() 
using namespace std; 
inline const int getint(){char c=getchar();int k=1,r=0;for(;c<'0'||c>'9';c=getchar())if(c=='-')k=-1;for(;c>='0'&&c<='9';c=getchar())r=r*10+c-'0';return k*r;} 
inline const int max(const int &a,const int &b){return a>b?a:b;} 
inline const int min(const int &a,const int &b){return a<b?a:b;} 
const int N=1E5+10; 
const int M=5*1E5+10; 
struct node{int x,y,z;}E[M]; 
struct snnn{int l,r,s;}T[N*30]; 
int n,m,fa[N<<1],num[N<<1],HH[N],id[N],H[N],cnt,f[N<<1][20]; 
int root[N],L[N<<1],R[N<<1],lch[N<<1],rch[N<<1],ST[N]; 
inline bool cmp(node X,node Y){return X.z<Y.z;} 
inline bool cmp2(int X,int Y){return HH[X]<HH[Y];} 
inline void init(){CC(fa,0);CC(num,0);CC(HH,0);CC(id,0);CC(H,0);CC(f,0);CC(root,0);CC(L,0);CC(R,0);CC(lch,0);CC(rch,0);CC(ST,0);} 
inline int find(int x){ 
    if (x==fa[x]) return x; 
    else{fa[x]=find(fa[x]);return fa[x];} 
} 
inline void LCA(int dd){ 
    for1(j,1,19) 
        for1(i,1,dd) 
            f[i][j]=f[f[i][j-1]][j-1]; 
} 
inline void dfs(int x){ 
    if (lch[x]&&rch[x]){ 
        L[x]=cnt+1; 
        dfs(lch[x]); 
        dfs(rch[x]); 
        R[x]=cnt; 
    }else{ 
        cnt++; ST[cnt]=x; 
        L[x]=cnt; R[x]=cnt; 
    } 
} 
inline void update(int l,int r,int &pos,int key){ 
    T[++cnt]=T[pos]; pos=cnt; T[pos].s++; 
    if (l==r) return; 
    int mid=(l+r)>>1; 
    if (key<=mid) update(l,mid,T[pos].l,key); else update(mid+1,r,T[pos].r,key); 
} 
inline int LCA_find(int x,int y){ 
    for3(i,19,0) 
     if ((f[x][i]!=0)&&(num[f[x][i]]<=y)) 
      x=f[x][i]; 
    return x; 
} 
inline int query(int LL,int RR,int key){ 
    int mid,s,x=root[LL],y=root[RR],l=1,r=n; 
    while (l<r){ 
        mid=(l+r)>>1; s=T[T[y].l].s-T[T[x].l].s; 
        if (key<=s) r=mid,x=T[x].l,y=T[y].l; 
        else l=mid+1,key-=s,x=T[x].r,y=T[y].r; 
    }return l; 
} 
int main(){
    init(); 
    read(n); read(m); int Q; read(Q); 
    for1(i,1,n) read(HH[i]),id[i]=i; 
    sort(id+1,id+n+1,cmp2); 
    for1(i,1,n) H[id[i]]=i; 
    for1(i,1,m){read(E[i].x);read(E[i].y);read(E[i].z);} 
    sort(E+1,E+m+1,cmp); 
    cnt=n+1; 
    for2(i,1,n<<1) fa[i]=i; 
    for1(i,1,m){ 
        int fx=find(E[i].x),fy=find(E[i].y); 
        if (fx!=fy){ 
            fa[fx]=cnt; fa[fy]=cnt; 
            f[fx][0]=cnt; f[fy][0]=cnt; 
            lch[cnt]=fx; rch[cnt]=fy; 
            num[cnt]=E[i].z; 
            cnt++; 
            if (cnt>((n<<1)-1)) break; 
        } 
    } 
    int dd=cnt-1; 
    LCA(dd); 
    cnt=0; 
    //for2(i,1,n<<1) if (L[i]==0) dfs(find(i)); 
    dfs(dd);
    cnt=0; 
    for1(i,1,n){ 
        root[i]=root[i-1]; 
        update(1,n,root[i],H[ST[i]]); 
    } 
    int a,b,c,la=0; 
    for1(i,1,Q){ 
        read(a); read(b); read(c); 
        //a^=la; b^=la; c^=la; 
        int rt=LCA_find(a,b),nn; 
        if (R[rt]-L[rt]+1<c) nn=-1; 
        else nn=HH[id[query(L[rt]-1,R[rt],R[rt]-L[rt]+2-c)]],la=nn; 
        printf("%d
",nn); 
    }
    return 0; 
}

3551：

#include<cstdio> 
#include<cstring> 
#include<algorithm> 
#define for1(i,a,n) for(int i=(a);i<=(n);i++) 
#define for2(i,a,n) for(int i=(a);i<(n);i++) 
#define for3(i,a,n) for(int i=(a);i>=(n);i--) 
#define for4(i,a,n) for(int i=(a);i>(n);i--) 
#define CC(i,a) memset(i,a,sizeof(i)); 
#define read(x) x=getint() 
using namespace std; 
inline const int getint(){char c=getchar();int k=1,r=0;for(;c<'0'||c>'9';c=getchar())if(c=='-')k=-1;for(;c>='0'&&c<='9';c=getchar())r=r*10+c-'0';return k*r;} 
inline const int max(const int &a,const int &b){return a>b?a:b;} 
inline const int min(const int &a,const int &b){return a<b?a:b;} 
const int N=1E5+10; 
const int M=5*1E5+10; 
struct node{int x,y,z;}E[M]; 
struct snnn{int l,r,s;}T[N*30]; 
int n,m,fa[N<<1],num[N<<1],HH[N],id[N],H[N],cnt,f[N<<1][20]; 
int root[N],L[N<<1],R[N<<1],lch[N<<1],rch[N<<1],ST[N]; 
inline bool cmp(node X,node Y){return X.z<Y.z;} 
inline bool cmp2(int X,int Y){return HH[X]<HH[Y];} 
inline void init(){CC(fa,0);CC(num,0);CC(HH,0);CC(id,0);CC(H,0);CC(f,0);CC(root,0);CC(L,0);CC(R,0);CC(lch,0);CC(rch,0);CC(ST,0);} 
inline int find(int x){ 
    if (x==fa[x]) return x; 
    else{fa[x]=find(fa[x]);return fa[x];} 
} 
inline void LCA(int dd){ 
    for1(j,1,19) 
        for1(i,1,dd) 
            f[i][j]=f[f[i][j-1]][j-1]; 
} 
inline void dfs(int x){ 
    if (lch[x]&&rch[x]){ 
        L[x]=cnt+1; 
        dfs(lch[x]); 
        dfs(rch[x]); 
        R[x]=cnt; 
    }else{ 
        cnt++; ST[cnt]=x; 
        L[x]=cnt; R[x]=cnt; 
    } 
} 
inline void update(int l,int r,int &pos,int key){ 
    T[++cnt]=T[pos]; pos=cnt; T[pos].s++; 
    if (l==r) return; 
    int mid=(l+r)>>1; 
    if (key<=mid) update(l,mid,T[pos].l,key); else update(mid+1,r,T[pos].r,key); 
} 
inline int LCA_find(int x,int y){ 
    for3(i,19,0) 
     if ((f[x][i]!=0)&&(num[f[x][i]]<=y)) 
      x=f[x][i]; 
    return x; 
} 
inline int query(int LL,int RR,int key){ 
    int mid,s,x=root[LL],y=root[RR],l=1,r=n; 
    while (l<r){ 
        mid=(l+r)>>1; s=T[T[y].l].s-T[T[x].l].s; 
        if (key<=s) r=mid,x=T[x].l,y=T[y].l; 
        else l=mid+1,key-=s,x=T[x].r,y=T[y].r; 
    }return l; 
} 
int main(){
    init(); 
    read(n); read(m); int Q; read(Q); 
    for1(i,1,n) read(HH[i]),id[i]=i; 
    sort(id+1,id+n+1,cmp2); 
    for1(i,1,n) H[id[i]]=i; 
    for1(i,1,m){read(E[i].x);read(E[i].y);read(E[i].z);} 
    sort(E+1,E+m+1,cmp); 
    cnt=n+1; 
    for2(i,1,n<<1) fa[i]=i; 
    for1(i,1,m){ 
        int fx=find(E[i].x),fy=find(E[i].y); 
        if (fx!=fy){ 
            fa[fx]=cnt; fa[fy]=cnt; 
            f[fx][0]=cnt; f[fy][0]=cnt; 
            lch[cnt]=fx; rch[cnt]=fy; 
            num[cnt]=E[i].z; 
            cnt++; 
            if (cnt>((n<<1)-1)) break; 
        } 
    } 
    int dd=cnt-1; 
    LCA(dd); 
    cnt=0; 
    //for2(i,1,n<<1) if (L[i]==0) dfs(find(i)); 
    dfs(dd);
    cnt=0; 
    for1(i,1,n){ 
        root[i]=root[i-1]; 
        update(1,n,root[i],H[ST[i]]); 
    } 
    int a,b,c,la=0; 
    for1(i,1,Q){ 
        read(a); read(b); read(c); 
        a^=la; b^=la; c^=la; 
        int rt=LCA_find(a,b),nn; 
        if (R[rt]-L[rt]+1<c) nn=-1,la=0;
        else nn=HH[id[query(L[rt]-1,R[rt],R[rt]-L[rt]+2-c)]],la=nn; 
        printf("%d
",nn); 
    }
    return 0; 
}

 然后就可以了O(∩_∩)O~~