Codeforces 1299D - Around the World(线性基+图论+dp)
一道线性基的综合题 %%%%%%
首先注意到“非简单路径”“异或和”等字眼,可以本能地想到线性基。根据线性基那一套理论,一个值 \(v\) 可以表示为某个 \(1\to 1\) 的非简单回路上边权的异或和当前节点它可以表示为 \(1\) 所在连通块的若干个 \(\ge 1\) 简单环上权值的异或和。其次我们还可以注意到本题至于很小,最高不过 \(2^5-1=31\),而稍微打个表即可发现大小为 \(5\) 的本质不同的线性基个数只有 \(374\) 个,因此可以暴搜求出所有本质不同的线性基,具体来说我们每次枚举 \(i=1,2,3,\cdots,31\),如果当前线性基里能够插入 \(i\) 就插入 \(i\) 并继续 DFS 下去,复杂度 \(374\times 31\times 5\),不会出问题。
然后考虑以线性基为状态的 DP,我们假设去掉 \(1\) 号点之后原图有 \(m\) 个连通块,注意到我们只能删除与 \(1\) 相连的边,故每个连通块内部环的情况是不受影响的,因此我们可以先按照 P4151 [WC2011]最大XOR和路径 的套路求出每个连通块中所有环的异或和组成的线性基,当然如果这个线性基是线性有关集那说明如果这个连通块如果与 \(1\) 相连就会存在不合法的路径,也就是说这个连通块必须与 \(1\) 断开,这个维护一个数组 \(ok_i\) 表示编号为 \(i\) 的连通块是否合法即可。不难发现每个连通块只可能有两种可能,一是连通块内部存在一个与 \(1\) 相连的点,二是连通块内部存在两个与 \(1\) 相连的点,不可能出现三个及以上与 \(1\) 相连的点,否则就会出现长度 \(>3\) 的经过 \(1\) 的环了,记 \(dp_{i,j}\) 为考虑前 \(i\) 个连通块,当前线性基为 \(j\)(我们将大小为 \(5\) 的本质不同的线性基编号 \(1,2,3,\cdots,374\))的方案数,考虑按照上面的分析过程分情况讨论:
- 若 \(ok_i=0\),那么 \(1\) 与当前连通块之间的边必须断开,即 \(dp_{i,j}=dp_{i-1,j}\)
- 若 \(ok_i=1\),继续分两种情况:
- 连通块内部存在一个与 \(1\) 相连的点,记 \(b_i\) 为连通块 \(i\) 的线性基,那么如果断开与 \(1\) 相连的边线性基不会发生变化,即 \(dp_{i,j}\leftarrow dp_{i-1,j}\),否则相当于将线性基 \(j\) 与 \(b_i\) 进行合并,即 \(dp_{i,j\cup b_i}\leftarrow dp_{i-1,j}\)
- 连通块内部存在两个与 \(1\) 相连的点,如果断开两条与 \(1\) 相连的边,线性基不会发生变化,即 \(dp_{i,j}\leftarrow dp_{i-1,j}\),如果断开一条与 \(1\) 相连的边,相当于将线性基 \(j\) 与 \(b_i\) 进行合并,即 \(dp_{i,j\cup b_i}\leftarrow 2·dp_{i-1,j}\),如果与 \(1\) 相连的边都不断开,那么这部分会多出一个环没有加入线性基,这个环就是经过 \(1\) 的三元环,记 \(w\) 为该三元环的权值异或和,那么 \(dp_{i,j\cup b_i\cup w}\leftarrow dp_{i-1,j}\)
预处理线性基的合并关系即可在常数时间内实现转移,时间复杂度 \(155\times 374+374^2+374n\)
const int MAXN=1e5;
const int MAXB=380;
const int MOD=1e9+7;
int n,m;
struct lbase{
int a[5];
lbase(){memset(a,0,sizeof(a));}
bool insert(int x){
for(int i=4;~i;i--) if(x>>i&1){
if(a[i]) x^=a[i];
else{
a[i]=x;
for(int j=0;j<i;j++) if(a[i]>>j&1) a[i]^=a[j];
for(int j=i+1;j<=4;j++) if(a[j]>>i&1) a[j]^=a[i];
return 1;
}
} return 0;
}
int hash(){return (a[4]<<10)|(a[3]<<6)|(a[2]<<3)|(a[1]<<1)|a[0];}
} b[MAXB+5],c[MAXN+5];
int rid[MAXN+5],num=0;
void dfsfind(lbase cur){
int hs=cur.hash();if(rid[hs]) return;
else rid[hs]=++num,b[num]=cur;
for(int i=1;i<=31;i++){
lbase nxt=cur;
if(nxt.insert(i)) dfsfind(nxt);
}
}
int trs[MAXB+5][MAXB+5];
void init_trs(){
for(int i=1;i<=num;i++) for(int j=1;j<=num;j++){
lbase tmp=b[i];bool ok=1;
for(int k=0;k<=4;k++) if(b[j].a[k]) ok&=tmp.insert(b[j].a[k]);
if(ok) trs[i][j]=rid[tmp.hash()];
}
}
int hd[MAXN+5],to[MAXN*2+5],val[MAXN*2+5],nxt[MAXN*2+5],ec=0;
void adde(int u,int v,int w){to[++ec]=v;val[ec]=w;nxt[ec]=hd[u];hd[u]=ec;}
int bel[MAXN+5],dis[MAXN+5],ok[MAXN+5],cmp=0,con[MAXN+5],fst[MAXN+5];
bool is[MAXN+5];int dfn[MAXN+5],tim=0;
void dfs(int x,int f){
bel[x]=cmp;dfn[x]=++tim;
for(int e=hd[x];e;e=nxt[e]){
int y=to[e],z=val[e];if(y==1) continue;
if(!bel[y]) dis[y]=dis[x]^z,dfs(y,x);
else if(y!=f&&dfn[x]>dfn[y]) ok[cmp]&=c[cmp].insert(dis[x]^dis[y]^z);
}
}
int dp[MAXN+5][MAXB+5];
int main(){
scanf("%d%d",&n,&m);dfsfind(*new(lbase));init_trs();
for(int i=1,u,v,w;i<=m;i++){
scanf("%d%d%d",&u,&v,&w);
adde(u,v,w);adde(v,u,w);
} //printf("%d\n",num);
for(int e=hd[1];e;e=nxt[e]){
int y=to[e],z=val[e];
if(!bel[y]){
cmp++;ok[cmp]=1;con[cmp]=z;
fst[cmp]=y;dfs(y,-1);
} else {
for(int ee=hd[y];ee;ee=nxt[ee]){
int v=to[ee],w=val[ee];
if(v==fst[bel[y]]){
is[bel[y]]=1;con[bel[y]]^=w^z;
break;
}
}
}
}
dp[0][rid[0]]=1;
for(int i=1;i<=cmp;i++){
for(int j=1;j<=num;j++) dp[i][j]=dp[i-1][j];
if(!ok[i]) continue;
if(!is[i]){
int id=rid[c[i].hash()];
for(int j=1;j<=num;j++) if(trs[j][id]){
dp[i][trs[j][id]]=(dp[i][trs[j][id]]+dp[i-1][j])%MOD;
}
} else {
int id1=rid[c[i].hash()];
bool ok=c[i].insert(con[i]);
int id2=rid[c[i].hash()];
for(int j=1;j<=num;j++){
if(trs[j][id1]) dp[i][trs[j][id1]]=(dp[i][trs[j][id1]]+2*dp[i-1][j]%MOD)%MOD;
if(ok&&trs[j][id2]) dp[i][trs[j][id2]]=(dp[i][trs[j][id2]]+dp[i-1][j]%MOD)%MOD;
}
}
} int ans=0;
for(int i=1;i<=num;i++) ans=(ans+dp[cmp][i])%MOD;
printf("%d\n",ans);
return 0;
}
Codeforces 1299D - Around the World(线性基+图论+dp)的更多相关文章
- 洛谷P4151 最大XOR和路径 [WC2011] 线性基+图论
正解:线性基+图论 解题报告: 传送门 首先可以思考一下有意义的路径会是什么样子,,,那就一定是一条链+一些环 挺显然的因为一条路径原路返回有没有意义辣?所以一定是走一条链+一些环(当然也可以麻油环, ...
- 【做题】CF388D. Fox and Perfect Sets——线性基&数位dp
原文链接https://www.cnblogs.com/cly-none/p/9711279.html 题意:求有多少个非空集合\(S \subset N\)满足,\(\forall a,b \in ...
- Codeforces 938G 线段树分治 线性基 可撤销并查集
Codeforces 938G Shortest Path Queries 一张连通图,三种操作 1.给x和y之间加上边权为d的边,保证不会产生重边 2.删除x和y之间的边,保证此边之前存在 3.询问 ...
- 【题解】 Codeforces 662A Gambling Nim (线性基)
662A,戳我戳我 Solution: 我们先取\(ans=a[1] \bigoplus a[2] \bigoplus ... \bigoplus a[n]\),然后我们定义\(c[i]=a[i] \ ...
- codeforces 1100F Ivan and Burgers 线性基 离线
题目传送门 题意: 给出 n 个数,q次区间查询,每次查询,让你选择任意个下标为 [ l , r ] 区间内的任意数,使这些数异或起来最大,输出最大值. 思路:离线加线性基. 线性基学习博客1 线性基 ...
- [bzoj 2115]线性基+图论
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=2115 给定一个带权无向图,要找出从1到n路径权值异或和最大的那一条的路径异或和. 考虑1到 ...
- BZOJ CF388D. Fox and Perfect Sets [线性基 数位DP]
CF388D. Fox and Perfect Sets 题意:求最大元素\(le n\)的线性空间的个数 给神题跪了 orz 容易想到 每个线性基对应唯一的线性空间,我们可以统计满足条件的对应空间不 ...
- bzoj 3811: 玛里苟斯【线性基+期望dp】
这个输出可是有点恶心啊--WA*inf,最后抄了别人的输出方法orz 还有注意会爆long long,要开unsigned long long 对于k==1,单独考虑每一位i,如果这一位为1则有0.5 ...
- codeforces 388D Fox and Perfect Sets(线性基+数位dp)
#include<bits/stdc++.h> using namespace std; #define fi first #define se second #define mp mak ...
随机推荐
- WeakMap与Map,使用WeakMap实现深拷贝循环引用问题
1.Map可以使用任意类型的key值,不限字符串,对象等. 2.WeakMap只能使用对象作为key值,是弱引用,当从WeakMap中移除时,会自动垃圾回收 3.Object只能用基本类型作为key值 ...
- Java:并发笔记-07
Java:并发笔记-07 说明:这是看了 bilibili 上 黑马程序员 的课程 java并发编程 后做的笔记 6. 共享模型之不可变 本章内容 不可变类的使用 不可变类设计 无状态类设计 6.1 ...
- Java:Set接口小记
Java:Set接口小记 对 Java 中的 Set接口 与 其实现类,做一个微不足道的小小小小记 概述 public interface Set<E> extends Collectio ...
- [技术博客] 通过ItemTouchHelper实现侧滑删除功能
通过ItemTouchHelper实现侧滑删除功能 一.效果 二.具体实现 demo中演示的这种左滑删除的效果在手机APP中比较常用,安卓也为我们提供了专门的辅助类ItemTouchHelper来帮助 ...
- Spring Cloud Gateway + Jwt + Oauth2 实现网关的鉴权操作
Spring Cloud Gateway + Jwt + Oauth2 实现网关的鉴权操作 一.背景 二.需求 三.前置条件 四.项目结构 五.网关层代码的编写 1.引入jar包 2.自定义授权管理器 ...
- 【http】https加速优化
目录 前言 HTTPS 的连接很慢 https 步骤简要划分 握手耗时 证书验证 CRL OCSP 硬件优化 软件优化 软件升级 协议优化 证书优化 会话复用 会话票证 预共享密钥 前言 主要记录 h ...
- 【数据结构&算法】05-线性表之数组
目录 前言 线性结构与非线性结构 数组 数组的两个限制 数组的随机访问特性 数组的操作 插入操作 删除操作 数组越界 容器 数组下标 前言 本笔记主要记录数组的一些基础特性及操作. 顺便解答下为什么大 ...
- 【python】以souhu邮箱为例学习DDT数据驱动测试
前言 DDT(Data-Driven Tests)是针对 unittest 单元测试框架设计的扩展库.允许使用不同的测试数据来运行一个测试用例,并将其展示为多个测试用例.通俗理解为相同的测试脚本使用不 ...
- Screenshot 库和Collections 库
一.screenShot 是 robot framework的标准类库,用于截取当前窗口,需要手动加载. 示例: 运行结果: 二.Collections 库 Collections 库同样为 Robo ...
- redis-sentinel "DENIED Redis is running in protected mode"
protected-mode no in sentinel.conf https://github.com/antirez/redis/issues/3106