【BZOJ 3569】DZY Loves Chinese II 随机化+线性基

  用到一个结论——[先建树,再给每个非树边一个权值,每个树边的权值为覆盖他的非树边的权值的异或和,然后如果给出的边存在一个非空子集异或和为0则不连通,否则连通](必须保证每条边的出现和消失只能由自己产生,即一个边不能由其他其他边异或得到,这就是我们随机化边权的原因)
  证明:(前置性质:I.只割非树边一定不可以 II.非树边"藏"在树边里 III.非树边是在树上是简单路径 IV.对于一个连通块(只考虑树结构),"藏"在他周围树边里的非树边要么连接他与外界(只在他周围树边里出现一次),要么从他进从他出(在他周围树边里出两次) V.若异或和为0,则所有边出现(非树边藏在树边里)偶数次)I.如果此时图不连通,则删掉的边一定存在一个非空子集异或和为0:我们就关注一个被割出来的块,我们先把他周围的树边割掉,那么现在得到的异或值里只有连接他与外界的非树边,因为他不连通,所以剩下的边都被割,异或进来就好,所以得0 II.如果此时的非空集合异或值为0,则原图一定不连通:先把树边割掉得到了几个连通块(只考虑树结构),并且异或起来,那么异或和里面所有出现次数为奇数的非树边都被割掉,现在能连接这些连通块的只有剩下的出现次数为偶数的非树边,而只用出现次数为偶数的非树边一定不能把这些连通块连城一个整体(你可以先把某个出现次数为偶数的非树边的两端所在连通块连起来,然后就可以得到(在树结构中)在被连接的连通块的中间的连通块中,紧挨着被连接的连通块的连通块,一定无法连接上被连接的联通块,于是证得) 证毕!
  接下来的实现就很容易了吧.
  我一开始用了一种十分神奇的随机化,虽然可以对,但是T掉了(单点1s可过,但是bzoj上0.5s),就是先进行40(30也行)次随机边权不同的最小生成树,得到对于每个边分别在哪几次出现,询问的时候只要不存在[给出所有边都没有选用的生成树]就不连通,否则连通
  (也许用不严格的随机化代替严格的运行是一种方法吧.)

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
char xB[(1<<15)+10],*xS,*xT;
#define gtc (xS==xT&&(xT=(xS=xB)+fread(xB,1,1<<15,stdin),xS==xT)?0:*xS++)
inline void read(int &x){
  register char ch=gtc;
  for(x=0;ch<‘0‘||ch>‘9‘;ch=gtc);
  for(;ch>=‘0‘&&ch<=‘9‘;x=(x<<1)+(x<<3)+ch-‘0‘,ch=gtc);
}
const int N=100010,M=500010,A=30;
int n,m,f[N],size[N],a[M],b[M],key[M],val[N];
inline int find(int x){return x==f[x]?x:f[x]=find(f[x]);}
struct V{int to,next,id;}c[N<<1];
int head[N],t;
inline void add(int x,int y,int z){c[++t].to=y,c[t].next=head[x],head[x]=t,c[t].id=z;}
#define bin(a,b) (((a)>>(b))&1)
struct Gay{
  int gay[A];
  inline void clear(){memset(gay,0,sizeof(gay));}
  inline bool insert(int x){
    for(int i=A-1;i>=0;--i){
      if(!bin(x,i))continue;
      if(gay[i])x^=gay[i];
      else return gay[i]=x,false;
    }return true;
  }
}gay;
inline void dfs(int x,int fa){
  for(int i=head[x];i;i=c[i].next)
    if(c[i].id!=fa)dfs(c[i].to,c[i].id),val[x]^=val[c[i].to];
  key[fa]=val[x];
}
int main(){
  read(n),read(m),srand(45);
  register int i,x;int y,k,temp,T=(1<<10)-1;
  for(i=1;i<=m;++i)read(a[i]),read(b[i]);
  for(i=1;i<=n;++i)f[i]=i,size[i]=1;
  for(i=1,temp=1;i<=m;++i){
    x=find(a[i]),y=find(b[i]);
    if(x==y){key[i]=rand(),val[a[i]]^=key[i],val[b[i]]^=key[i];continue;}
    ++temp,size[x]<size[y]?(x^=y^=x^=y):0,f[y]=x,size[x]+=size[y],add(a[i],b[i],i),add(b[i],a[i],i);
  }
  if(temp!=n){read(T);while(T--)puts("Disconnected");return 0;}
  dfs(1,0),key[0]=0,read(T);
  register bool ans;register int cnt=0;
  while(T--){
    read(k),ans=false,gay.clear();
    while(k--){
      if(ans){read(x);continue;}
      read(x),x^=cnt;
      if(gay.insert(key[x]))ans=true;
    }
    puts(ans?"Disconnected":(++cnt,"Connected"));
  }return 0;
}

正解

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
typedef long long LL;
char xB[(1<<15)+10],*xS=xB,*xT=xB;
#define gtc (xS==xT&&(xT=(xS=xB)+fread(xB,1,1<<15,stdin),xS==xT)?0:*xS++)
inline void read(int &x){
    register char ch=gtc;
    for(x=0;ch<‘0‘||ch>‘9‘;ch=gtc);
    for(;ch>=‘0‘&&ch<=‘9‘;x=(x<<1)+(x<<3)+ch-‘0‘,ch=gtc);
}
const int N=100010;
const int M=500010;
int n,m,f[N],size[N],a[M],b[M],sort[M],key[M],buc[M],full=(1<<16)-1;
LL pos[M],Full;
inline int find(int x){
  return x==f[x]?x:f[x]=find(f[x]);
}
int main(){
  read(n),read(m);
  register int i;int k,x,y,t,id;LL temp;
  for(i=1;i<=m;++i)read(a[i]),read(b[i]);
  Full=(1LL<<40LL)-1;
  for(k=1;k<=40;++k){
    temp=1LL<<(k-1LL);
    srand(666233+k);
    for(i=0;i<=full;++i)buc[i]=0;
    for(i=1;i<=m;++i)key[i]=rand()&full,++buc[key[i]];
    for(i=1;i<=full;++i)buc[i]+=buc[i-1];
    if(k&1)for(i=1;i<=m;++i)sort[buc[key[i]]--]=i;
    else for(i=m;i>0;--i)sort[buc[key[i]]--]=i;
    for(i=1;i<=n;++i)f[i]=i,size[i]=1;
    for(i=1,id=sort[1],t=1;i<=m&&t<n;++i,id=sort[i]){
      x=find(a[id]),y=find(b[id]);
      if(x==y)continue;
      ++t,pos[id]|=temp;
      if(size[x]<size[y])std::swap(x,y);
      f[y]=x,size[x]+=size[y];
    }
    if(k==1&&t!=n){
      read(t);
      while(t--)puts("Disconnected");
      return 0;
    }
  }
  read(t);LL ans;y=0;
  while(t--){
    ans=0,read(k);
    while(k--)read(x),x^=y,ans|=pos[x];
    puts(ans==Full?"Disconnected":(++y,"Connected"));
  }return 0;
}

一开始的TLE解法

原文地址:https://www.cnblogs.com/TSHugh/p/8379531.html

时间: 2024-11-07 00:18:46

【BZOJ 3569】DZY Loves Chinese II 随机化+线性基的相关文章

[BZOJ3569]DZY Loves Chinese II(随机化+线性基)

3569: DZY Loves Chinese II Time Limit: 5 Sec  Memory Limit: 64 MBSubmit: 1515  Solved: 569[Submit][Status][Discuss] Description 神校XJ之学霸兮,Dzy皇考曰JC. 摄提贞于孟陬兮,惟庚寅Dzy以降. 纷Dzy既有此内美兮,又重之以修能. 遂降临于OI界,欲以神力而凌♂辱众生. 今Dzy有一魞歄图,其上有N座祭坛,又有M条膴蠁边. 时而Dzy狂WA而怒发冲冠,神力外溢,

BZOJ 3563 DZY Loves Chinese / BZOJ 3569 DZY Loves Chinese II 随机化+高斯消元解异或方程组

题目大意:给出一个无向图,问删掉k条边的时候,图是否联通. 思路:虽然我把这两个题放在了一起,但是其实这两个题可以用完全不同的两个解法来解决. 第一个题其实是DZY出错了...把每次的边数也异或了,那就直接用这个性质一个一个往后推就行了..最后一个暴力求一下.. 第二个题才是本意啊. 听到做法的时候我惊呆了.. 首先是将整个图中拆出一个树,那么所有边就分为树边和非树边.将所有非树边都加一个随机权值.树边的权值是所有能够覆盖它的非树边的权值的异或和. 把整个图拆开的充要条件是拆掉一条树边,同时将所

bzoj 3569: DZY Loves Chinese II

链接 3569: DZY Loves Chinese II 题目大意:给出一张\(n\)个点\(m\)条边的无向图,进行\(q\)次询问,问删掉某\(k\)条边后图是否联通,强制在线. \(N≤100000 \ M≤500000\ Q≤50000\ 1≤K≤15\) 先考虑一下离线怎么做: \(cdq\)分治. 首先把所有没有影响的边都建出来 分治过程: 1.把左边没有右边有的边建出来 2.分治左边 3.把并查集恢复至初始的样子 4.把右边没有左边有的边建出来 5.分治右边 虽然每次\(cdq\

BZOJ 3569 DZY Loves Chinese II 高斯消元

题目大意:给定一个[魞歄连通图],多次询问当图中某k条边消失时这个图是否联通 强制在线 我们找到这个图的任意一棵生成树 然后对于每条非树边将其的权值赋为一个随机数 对于每条树边 我们将这条树边的权值设为所有覆盖这条树边的边权的异或和 那么图不连通当且仅当删除一条树边和覆盖这条树边的所有边集 而由于刚才的处理一条树边和覆盖这条边的所有边集的异或和为零 于是问题转化成了对于给定的k条边是否存在一个边权的异或和为零的子集 果断高斯消元 由于使用了随机化所以碰撞率极低 好方法学习了...构思真是巧妙 记

bzoj 3569 DZY Loves Chinese II 随机算法 树上倍增

题意:给你一个n个点m条边的图,有若干组询问,每次询问会选择图中的一些边删除,删除之后问此图是否联通?询问之间相互独立.此题强制在线. 思路:首先对于这张图随便求一颗生成树,对于每一条非树边,随机一个权值.树边的权值为所有覆盖这条树边的非树边的权值异或和.覆盖是指这条边是个返祖边,并且一端在父节点方向,一端在子节点方向.这样,我们选出若干条边,看一下他们异或起来是不是0,如果是0,那么相当于把一条树边和它的所有子节点方向的返祖边全部断开,那么图就不连通了. 代码: #include <bits/

【BZOJ3563/3569】DZY Loves Chinese II 线性基神题

[BZOJ3563/3569]DZY Loves Chinese II Description 神校XJ之学霸兮,Dzy皇考曰JC. 摄提贞于孟陬兮,惟庚寅Dzy以降. 纷Dzy既有此内美兮,又重之以修能. 遂降临于OI界,欲以神力而凌♂辱众生. 今Dzy有一魞歄图,其上有N座祭坛,又有M条膴蠁边. 时而Dzy狂WA而怒发冲冠,神力外溢,遂有K条膴蠁边灰飞烟灭. 而后俟其日A50题则又令其复原.(可视为立即复原) 然若有祭坛无法相互到达,Dzy之神力便会大减,于是欲知其是否连通. Input 第

BZOJ 3563 DZY Loves Chinese 并查集

题目大意:给定一个无向联通图,q次询问当图中某k条边消失时图是否联通 强制在线 逗比题233 不明白什么意思的去看DZY Loves Chinese II的红字就明白这题为何逗比了0.0 #include <cstdio> #include <cstring> #include <iostream> #include <algorithm> #define M 100100 using namespace std; struct edges{ int x,y

[BZOJ3563&amp;3569]DZY Loves Chinese

bzoj 加强版 sol 其实前一题还有一种解法的,具体方法请见bzoj讨论版. 以下是正解(?) 建一棵生成树. 考虑什么时候图会不连通:当且仅当存在一条树边被删除,同时所有覆盖了他的非树边也被删除. 所以考虑给每条非树边随机一个权值,然后把每条树边的权值设为所以覆盖了这条边的非树边的边权异或和. 这样一来在要删除的边集中,若存在一个异或和为零的子集则说明图会不连通. 线性基维护一下.复杂度\(O(n+Qk*30)\) code 我WA了半天结果是随机种子出锅了... #include<cst

BZOJ 3563 DZY Loves Chinese

题解:为每条非树边赋一个权值 每条树边的权值为覆盖他的非树边权值异或和 如果边集的子集线性相关,相当于把树边和非树边拦腰砍断,则不连通 用线性基判断线性相关 问题:为什么srand(time(0))会RE? #include<iostream> #include<cstdio> #include<cstring> #include<ctime> #include<cstdlib> using namespace std; const int ma