【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[(<<)+],*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=;ch<''||ch>'';ch=gtc);
for(;ch>=''&&ch<='';x=(x<<)+(x<<)+ch-'',ch=gtc);
}
const int N=,M=,A=;
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<<];
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,,sizeof(gay));}
inline bool insert(int x){
for(int i=A-;i>=;--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();
register int i,x;int y,k,temp,T=(<<)-;
for(i=;i<=m;++i)read(a[i]),read(b[i]);
for(i=;i<=n;++i)f[i]=i,size[i]=;
for(i=,temp=;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):,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 ;}
dfs(,),key[]=,read(T);
register bool ans;register int cnt=;
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 ;
}
正解
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
typedef long long LL;
char xB[(<<)+],*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=;ch<''||ch>'';ch=gtc);
for(;ch>=''&&ch<='';x=(x<<)+(x<<)+ch-'',ch=gtc);
}
const int N=;
const int M=;
int n,m,f[N],size[N],a[M],b[M],sort[M],key[M],buc[M],full=(<<)-;
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=;i<=m;++i)read(a[i]),read(b[i]);
Full=(1LL<<40LL)-;
for(k=;k<=;++k){
temp=1LL<<(k-1LL);
srand(+k);
for(i=;i<=full;++i)buc[i]=;
for(i=;i<=m;++i)key[i]=rand()&full,++buc[key[i]];
for(i=;i<=full;++i)buc[i]+=buc[i-];
if(k&)for(i=;i<=m;++i)sort[buc[key[i]]--]=i;
else for(i=m;i>;--i)sort[buc[key[i]]--]=i;
for(i=;i<=n;++i)f[i]=i,size[i]=;
for(i=,id=sort[],t=;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==&&t!=n){
read(t);
while(t--)puts("Disconnected");
return ;
}
}
read(t);LL ans;y=;
while(t--){
ans=,read(k);
while(k--)read(x),x^=y,ans|=pos[x];
puts(ans==Full?"Disconnected":(++y,"Connected"));
}return ;
}
一开始的TLE解法
【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][S ...
- BZOJ 3569 DZY Loves Chinese II 树上差分+线性基
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=3569 Description 神校XJ之学霸兮,Dzy皇考曰JC. 摄提贞于孟陬兮,惟庚寅 ...
- BZOJ 3569 DZY Loves Chinese II ——线性基
[题目分析] 腊鸡题目卡题面. 大概的意思就是给一张无向图,每次删掉其中一些边,问是否联通. 首先想到的是Bitset,可以做到n^2/64.显然过不了. 然而这是lyd在给我们讲线性基的时候的一道题 ...
- BZOJ 3569: DZY Loves Chinese II(线性基)
传送门 解题思路 首先构造出一个生成树,考虑不连接的情况.假设连通两点的非树边和树边都断掉后不连通,那么可以给所有的非树边随机一个互不相同的值,然后树边的权值为过他两端点的非树边权值的异或和,这个可以 ...
- BZOJ 3569: DZY Loves Chinese II [高斯消元XOR 神题]
http://www.lydsy.com/JudgeOnline/problem.php?id=3569 题意:多次询问一个无向连通图当图中某k条边消失时这个图是否联通 强制在线 太神啦啦啦啦啦啦啦啦 ...
- BZOJ 3569 DZY Loves Chinese II
Description 神校XJ之学霸兮,Dzy皇考曰JC. 摄提贞于孟陬兮,惟庚寅Dzy以降. 纷Dzy既有此内美兮,又重之以修能. 遂降临于OI界,欲以神力而凌♂辱众生. 今Dzy有一魞歄图,其上 ...
- bzoj 3569 DZY Loves Chinese II 随机算法 树上倍增
题意:给你一个n个点m条边的图,有若干组询问,每次询问会选择图中的一些边删除,删除之后问此图是否联通?询问之间相互独立.此题强制在线. 思路:首先对于这张图随便求一颗生成树,对于每一条非树边,随机一个 ...
- BZOJ 3563 DZY Loves Chinese
Description 神校XJ之学霸兮,Dzy皇考曰JC. 摄提贞于孟陬兮,惟庚寅Dzy以降. 纷Dzy既有此内美兮,又重之以修能. 遂降临于OI界,欲以神力而凌♂辱众生. 今Dzy有一魞歄图,其上 ...
- 【BZOJ3563/3569】DZY Loves Chinese II 线性基神题
[BZOJ3563/3569]DZY Loves Chinese II Description 神校XJ之学霸兮,Dzy皇考曰JC. 摄提贞于孟陬兮,惟庚寅Dzy以降. 纷Dzy既有此内美兮,又重之以 ...
随机推荐
- PHP使用Redis消息队列
1.redis安装 参考:菜鸟教程http://www.runoob.com/redis/redis-install.html 2.安装php的redis扩展 1)phpinfo()查看php版本信息 ...
- 栈和队列ADT -数据结构(C语言实现)
数据结构与算法分析 栈模型 限制插入和删除只能在表的末端的表 表的末端叫做栈顶(top) 支持Push进栈和Pop入栈操作 //LIFO后进先出表 栈的实现 链表实现 类型声明 struct Node ...
- Python3 Tkinter-Text
1.创建 from tkinter import * root=Tk() t=Text(root) t.pack() root.mainloop() 2.添加文本 from tkinter impor ...
- VBA基础之Excel 工作薄(Book)的操作(三)
三. Excel 工作薄(Book)的操作1. Excel 创建工作薄(Book) Sub addWorkbook() Workbooks.Add End Sub 2. Excel 打开工作薄(Boo ...
- Python中的list
list的创建 1 字面量 >>>L = [1, 2, 3] [1, 2, 3] 2 通过iterable可迭代对象,比如str对象,range对象,map对象 >>&g ...
- DAY4敏捷冲刺
站立式会议 工作安排 (1)服务器配置 已完成对微信小程序登录凭证储存至云端数据库,计划使用微信接口返回的session_id进行转化返回本地,以保持登录态. (2)数据库配置 单词学习记录+用户信息 ...
- 第三章——供机器读取的数据(CSV与JSON)
本书使用的文件.代码:https://github.com/huangtao36/data_wrangling 机器可读(machine readable)文件格式: 1.逗号分隔值(Comma-Se ...
- Debian 7 amd64问题
Debian 7 发布了有1段时间,最近才在自己的电脑硬盘安装,用户体验还算可以.在安装Debian的过程中,有问题还是要记录一下的. 注意:遇到的问题跟硬件体系相关,可能在个别电脑没法重现. 1.默 ...
- React Components Template
React Components Template "use strict"; /** * * @author xgqfrms * @license MIT * @copyrigh ...
- [转]MATLAB cell数据类型
细胞型数据类型(cell)使不同类型和不同维数的数组可以共存,细胞型数组实际上可以认为是一种以任意形式的数组为分量的多维数组. 1.细胞型数据的定义 1)直接赋值定义:细胞型变量在定义时需要使用大括号 ...