2019牛客多校第三场 Graph Games

19 篇文章 0 订阅
8 篇文章 0 订阅

在这里插入图片描述
在这里插入图片描述

题意

给出一个无向图, S ( x ) S(x) S(x) 表示点 x x x 所能达到的所有点的集合
有两种操作,一种是修改一个区间的边的状态,即在图上的变为不在,不在的加上去
另一种是询问两个点的 S ( x ) 、 S ( y ) S(x)、S(y) S(x)S(y) 是否相同

题解

利用 H a s h Hash Hash,给每个点附一个随机的64位整数,将所连点的权值异或和作为 S ( x ) S(x) S(x)的值,来判断是否相同

将所有的点分为两种点,重点和轻点
重点是度数大于等于 s q r t ( m ) sqrt(m) sqrt(m) 的点,否则是轻点
建立一个图,特征是,重点连重点,轻点连重点,轻点连轻点,不能重点连轻点
因为,我们这样做就是要让轻点更新重点,重点连轻点的话,复杂度会大大上升

将边的序列分块
对于轻点,修改时,块间打标记,块内暴力更新。查询时,由于轻点连的边少,所以枚举轻点的边,看边所在的块是否有标记,有的话,更新答案。
对于重点,修改时,块间打标记,块内暴力更新。查询时,肯定不能遍历边了,我们要遍历块
考虑每个块对重点的贡献,重点记录下来每个块对它的贡献,遍历块时,加上被标记的块的贡献

代码
#include<bits/stdc++.h>
#define N 100010
#define INF 0x3f3f3f3f
#define eps 1e-10
// #define pi 3.141592653589793
#define mod 998244353
#define P 1000000007
#define LL long long
#define pb push_back
#define fi first
#define se second
#define cl clear
#define si size
#define lb lower_bound
#define ub upper_bound
#define bug(x) cerr<<#x<<"      :   "<<x<<endl
#define mem(x) memset(x,0,sizeof x)
#define sc(x) scanf("%lld",&x)
#define scc(x,y) scanf("%d%d",&x,&y)
#define sccc(x,y,z) scanf("%d%d%d",&x,&y,&z)
using namespace std;

namespace IO{ 
    #define BUF_SIZE 100000 
    bool IOerror=0; 
    inline char nc(){ 
        static char buf[BUF_SIZE],*p1=buf+BUF_SIZE,*pend=buf+BUF_SIZE; 
        if (p1==pend){ 
            p1=buf; pend=buf+fread(buf,1,BUF_SIZE,stdin); 
            if (pend==p1){IOerror=1;return -1;} 
        } 
        return *p1++; 
    } 
    inline bool blank(char ch){return ch==' '||ch=='\n'||ch=='\r'||ch=='\t';} 
    inline void read(int &x){ 
        bool sign=0; char ch=nc(); x=0; 
        for (;blank(ch);ch=nc()); 
        if (IOerror)return; 
        if (ch=='-')sign=1,ch=nc(); 
        for (;ch>='0'&&ch<='9';ch=nc())x=x*10+ch-'0'; 
        if (sign)x=-x; 
    } 
};

using namespace IO;

typedef pair<int,LL> pi;
int n,m,q,block,e[N<<1][2],in[N],fg[N],tag[1010],la[N],cnt=0;
vector<pi> g[N];
LL w[N],ans[N],tp[N];
struct node{
    int x,y,z;
}v[N<<2];
inline void add(int x,int y,int z){
    v[++cnt]=node{y,z,la[x]};
    la[x]=cnt;
}

int main(){
    mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());
    int T;
    read(T);
    while(T--){
        read(n);read(m); cnt=0;
        memset(in,0,sizeof(int)*(n+3));
        memset(la,0,sizeof(int)*(n+3));
        memset(ans,0,sizeof(LL)*(n+3));
        mem(tag);
        for(int i=1;i<=n;i++) w[i]=rng(),g[i].clear();
        for(int i=1,x,y;i<=m;i++) {
            read(x);read(y);
            in[x]++;in[y]++;
            ans[x]^=w[y]; ans[y]^=w[x];
            e[i][0]=x;e[i][1]=y;
        }

        block=sqrt(m); 
        for(int i=1;i<=n;i++) fg[i]=(in[i]>=block);

        int t=0;
        for(int i=0;i<=m;i+=block){
            t++;
            for(int j=(i==0);j<block;j++)if (i+j<=m){
                int x=e[i+j][0],y=e[i+j][1];
                if (fg[x]) tp[x]^=w[y];
                if (fg[y]) tp[y]^=w[x];
                if(fg[x]&&fg[y]){
                    add(x,y,i+j); add(y,x,i+j);
                }
                if (!fg[x]) add(x,y,i+j);
                if (!fg[y]) add(y,x,i+j);
            }
            for(int j=(i==0);j<block;j++)if(i+j<=m){
                int x=e[i+j][0],y=e[i+j][1];
                if (fg[x]&&tp[x])g[x].pb(pi(t,tp[x]));
                if (fg[y]&&tp[y])g[y].pb(pi(t,tp[y]));
                tp[y]=0;tp[x]=0;
            }
        }

        read(q);
        while(q--){
            int op,l,r;
            read(op);read(l);read(r);
            if (op==1){
                int bl=l/block+1,br=r/block+1;
                if (bl==br){
                    for(int i=l;i<=r;i++) {
                    int x=e[i][0],y=e[i][1];
                    ans[x]^=w[y];
                    ans[y]^=w[x];
                    }
                    continue;
                }
                for(int i=bl+1;i<br;i++) tag[i]^=1;
                for(int i=l;i<bl*block;i++) {
                    int x=e[i][0],y=e[i][1];
                    ans[x]^=w[y];
                    ans[y]^=w[x];
                }
                for(int i=(br-1)*block;i<=r;i++){
                    int x=e[i][0],y=e[i][1];
                    ans[x]^=w[y];
                    ans[y]^=w[x];
                }
            }else{
                LL tl=ans[l],tr=ans[r];
                if (fg[l]) {
                    for(auto i:g[l])if (tag[i.fi]) tl^=i.se;
                }else
                    for(int i=la[l];i;i=v[i].z) if(tag[v[i].y/block+1]) tl^=w[v[i].x];
                if (fg[r]) {
                    for(auto i:g[r])if (tag[i.fi]) tr^=i.se;
                }else
                    for(int i=la[r];i;i=v[i].z) if(tag[v[i].y/block+1]) tr^=w[v[i].x];
                if (tl==tr) putchar('1');else putchar('0');
            }
        }
        putchar('\n');
    }
}



  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值