杭电多校第八场 HDU6858- Discovery of Cycles(LCT+快读)

题意:

给定一个图有m条边,边的序号为1-m,q个强制在线询问(假的强制在线,不过在线能做就是了),每次问第l到第r条边组成的子图是否存在环。

思路:

首先每次询问对应的图不一样,要建立一个新图的时间复杂度必然是不够的,考虑先进行预处理。思维一下,一个区间左端点固定了,右端点越大,图上有环的概率越大。实际上,对于一个固定的左端点,如果找出一个最小右端点使子图上有环,那么右端点再往右移也必然有环。因此只要对所有左端点预处理出一个最小的右端点即可O(1)判断给定区间是否可以构成有环的子图。

如何预处理?思维一下想到了滑动区间,左端点l固定,右端点r一直向右移动,将边加入子图,直到图上可以构成环为止,此时得到pre[l]=r,先不将边r加入图中,将左端点右移一位,即将图上一条边删去,再向右移动右端点,直至图上有环……如此往复,只要进行O(n)次加边和删边的操作,就可以得到所有左端点对应的最小右端点。如何进行删边加边和判断图上有没有环?用LCT即可,若新加入的一条边对应的两个节点的根相同,则加入这条边后图上就会有环(其实这里就是将LCT当成可以删边的并查集进行使用)

这题读入数据量非常大,有将近1e7个int,可以使用超级快读模板,使用之后快1000ms

#include <bits/stdc++.h>
const int maxn=3e5+5;
using namespace std;
namespace Fast_IO{ //orz laofu
    const int MAXL((1 << 18) + 1);int iof, iotp;
    char ioif[MAXL], *ioiS, *ioiT, ioof[MAXL],*iooS=ioof,*iooT=ioof+MAXL-1,ioc,iost[55];
    char Getchar(){
        if (ioiS == ioiT){
            ioiS=ioif;ioiT=ioiS+fread(ioif,1,MAXL,stdin);return (ioiS == ioiT ? EOF : *ioiS++);
        }else return (*ioiS++);
    }
    void Write(){fwrite(ioof,1,iooS-ioof,stdout);iooS=ioof;}
    void Putchar(char x){*iooS++ = x;if (iooS == iooT)Write();}
    inline int read(){
        int x=0;for(iof=1,ioc=Getchar();(ioc<'0'||ioc>'9')&&ioc!=EOF;)iof=ioc=='-'?-1:1,ioc=Getchar();
        if(ioc==EOF)Write(),exit(0);
        for(x=0;ioc<='9'&&ioc>='0';ioc=Getchar())x=(x<<3)+(x<<1)+(ioc^48);return x*iof;
    }
    inline long long read_ll(){
        long long x=0;for(iof=1,ioc=Getchar();(ioc<'0'||ioc>'9')&&ioc!=EOF;)iof=ioc=='-'?-1:1,ioc=Getchar();
        if(ioc==EOF)Write(),exit(0);
        for(x=0;ioc<='9'&&ioc>='0';ioc=Getchar())x=(x<<3)+(x<<1)+(ioc^48);return x*iof;
    }
    void Getstr(char *s, int &l){
        for(ioc=Getchar();ioc==' '||ioc=='\n'||ioc=='\t';)ioc=Getchar();
        if(ioc==EOF)Write(),exit(0);
        for(l=0;!(ioc==' '||ioc=='\n'||ioc=='\t'||ioc==EOF);ioc=Getchar())s[l++]=ioc;s[l] = 0;
    }
    template <class Int>void Print(Int x, char ch = '\0'){
        if(!x)Putchar('0');if(x<0)Putchar('-'),x=-x;while(x)iost[++iotp]=x%10+'0',x/=10;
        while(iotp)Putchar(iost[iotp--]);if (ch)Putchar(ch);
    }
    void Putstr(const char *s){for(int i=0,n=strlen(s);i<n;++i)Putchar(s[i]);}
} // namespace Fast_IO
using namespace Fast_IO;
namespace LCT {
    struct node {
        int f, son[2] ;
        bool fz ;
    } tr[maxn] ;
    inline void init(int n){
        for(int i=0;i<=n;i++)
            tr[i].f=tr[i].son[0]=tr[i].son[1]=tr[i].fz=0;
    }
    inline void reverse ( int x ) {
        tr[x].fz = 0 ;
        swap ( tr[x].son[0], tr[x].son[1] ) ;
        int lc = tr[x].son[0] ;
        int rc = tr[x].son[1] ;
        tr[lc].fz ^= 1, tr[rc].fz ^= 1 ;
    }
    inline void rotate ( int x, int w ) {
        int f = tr[x].f, ff = tr[f].f, r, R ;
        r = tr[x].son[w] ;
        R = f ;
        tr[R].son[1-w] = r ;
        if ( r )
            tr[r].f = R ;
        r = x ;
        R = ff ;
        if ( tr[R].son[0] == f )
            tr[R].son[0] = r ;
        else if ( tr[R].son[1] == f )
            tr[R].son[1] = r ;
        tr[r].f = R ;
        r = f ;
        R = x ;
        tr[R].son[w] = r ;
        tr[r].f = R ;
    }
    int tmp[maxn] ;
    inline void splay ( int x, int rt ) {
        int s = 0, i = x ;
        while ( tr[i].f && (tr[tr[i].f].son[0]==i||tr[tr[i].f].son[1]==i) ) {
            tmp[++s] = i ;
            i = tr[i].f ;
        }
        tmp[++s] = i ;
        while ( s ) {
            i = tmp[s] ;
            s -- ;
            if ( tr[i].fz )
                reverse(i) ;
        }
        while ( tr[x].f!=rt && (tr[tr[x].f].son[0]==x||tr[tr[x].f].son[1]==x) ) {
            int f = tr[x].f, ff = tr[f].f ;
            if ( ff == rt || (tr[ff].son[0]!=f&&tr[ff].son[1]!=f) ) {
                if ( x == tr[f].son[0] )
                    rotate(x,1) ;
                else
                    rotate(x,0) ;
            } else {
                if ( tr[ff].son[0] == f && tr[f].son[0] == x )
                    rotate(f,1), rotate(x,1) ;
                else if ( tr[ff].son[1] == f && tr[f].son[1] == x )
                    rotate(f,0), rotate(x,0) ;
                else if ( tr[ff].son[0] == f && tr[f].son[1] == x )
                    rotate(x,0), rotate(x,1) ;
                else if ( tr[ff].son[1] == f && tr[f].son[0] == x )
                    rotate(x,1), rotate(x,0) ;
            }
        }
    }
    inline void access ( int x ) {
        int y = 0 ;
        while ( x ) {
            splay ( x, 0 ) ;
            tr[x].son[1] = y ;
            if ( y != 0 )
                tr[y].f = x ;
            y = x ;
            x = tr[x].f ;
        }
    }
    inline void makeroot ( int x ) {
        access(x) ;
        splay(x,0) ;
        tr[x].fz ^= 1 ;
    }
    inline void link ( int x, int y ) {
        makeroot(x) ;
        tr[x].f = y ;
        access(x) ;
    }
    inline void cut ( int x, int y ) {
        makeroot(x) ;
        access(y) ;
        splay(y,0) ;
        tr[tr[y].son[0]].f = 0 ;
        tr[y].son[0] = 0 ;
    }
    int find_root ( int x ) {
        access(x) ;
        splay(x,0) ;
        while ( tr[x].son[0] != 0 )
            x = tr[x].son[0] ;
        return x ;
    }
}
using namespace LCT;
int u[maxn],v[maxn];
int pre[maxn];
int main () {
    int T;
    T=read();
    while(T--){
        int n,m,q;
        n=read();
        m=read();
        q=read();
        init(m);
        memset(pre,0,sizeof(int)*(m+1));
        for(int i=1;i<=m;i++){
            u[i]=read();
            v[i]=read();
        }
        int p=1;
        for(int i=1;i<=m;i++){
            if(i>1){
                cut(u[i-1],v[i-1]);
            }
            for(p;p<=m;p++){
                if(find_root(u[p])!=find_root(v[p])){
                    link(u[p],v[p]);
                }
                else{
                    pre[i]=p;
                    break;
                }
            }
        }
        int lastans=0;
        for(int i=1;i<=q;i++){
            int l1,r1;
            l1=read();
            r1=read();
            int k1=(l1^lastans)%m+1;
            int k2=(r1^lastans)%m+1;
            l1=min(k1,k2);
            r1=max(k1,k2);
            if(pre[l1] && r1>=pre[l1]){
                lastans=1;
                puts("Yes");
            }
            else{
                lastans=0;
                puts("No");
            }
        }
    }
}

参考与引用: https://jkchen.blog.csdn.net/article/details/107991684

posted @ 2020-08-14 13:34  UCPRER  阅读(288)  评论(0编辑  收藏  举报