[BZOJ3569] 线性基+随机化

[BZOJ3569]

  • 考虑先构出一棵原图的生成树,对于每条非树边rand()一个权值,每条树边的权值设为覆盖了这条边的非树边边权异或和。
  • 那么每次如果导致图不连通,当且仅当一条树边被删除且覆盖了这条树边的所有边都被删除了,这种情况下这些边的权值异或起来为0。 
  • 线性基搞一搞
  • zbl zbl
  • 前往非洲的单程票拿好

  • 果然非洲人不适合写随机算法 

Code


#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define per(i,a,b) for(int i=(a);i>=(b);i--)
#define ll long long
using namespace std;
const int inf=2333333333;
const int N=1e6+100;
const int lim=30;
struct node{int y,n;}e[N];
int lin[N],flag,len=0,ans=0,x,y,n,m,q;
int a[N],b[N],c[50];
bool v[N];
void read(int x,int y)
{e[++len].y=y,e[len].n=lin[x],lin[x]=len;}
void dfs1(int x){
	v[x]=1;
    for(int i=lin[x];i;i=e[i].n){
        int y=e[i].y;
        if(!v[y]){
            dfs1(y);
        }else{
            b[i]=rand()%inf;
            a[x]^=b[i];
            a[y]^=b[i];
        }
    }
}
void dfs2(int x){
	v[x]=1;
    for(int i=lin[x];i;i=e[i].n){
        int y=e[i].y;
        if(v[y])continue;
        dfs2(y);
        b[i]^=a[y];
        a[x]^=a[y];
    }
}
void insert(int v){
    per(i,lim,0)if((v>>i)&1){
        if(c[i])v^=c[i];
        else{
            per(j,i-1,0)if((v>>j)&1)v^=c[j];
            rep(j,i+1,lim)if((c[j]>>i)&1)c[j]^=v;
            c[i]=v;
            break;
        }
    }
	if(!v)flag=1;
}
int main()
{
    srand(time(0));
    scanf("%d%d",&n,&m);
    rep(i,1,m)scanf("%d%d",&x,&y),read(min(x,y),max(x,y));
    dfs1(1);memset(v,0,sizeof(v));dfs2(1);
    scanf("%d",&q);
    while(q--){
        memset(c,0,sizeof(c));flag=0;
        int k;scanf("%d",&k);
        rep(i,1,k){
            int x;scanf("%d",&x);
            x^=ans;
            insert(b[x]);
        }
        printf("%s\n",flag?"Disconnected":"Connected");
        ans+=!flag;
    }
    return 0;
}

 

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值