【bzoj4028】【HEOI2015】公约数数列【分块暴力】

公约数数列的分块解决方案
这篇博客探讨了如何解决公约数数列问题,采用了分块策略,并详细解释了利用GCD和XOR以及哈希表进行高效查询的方法。虽然暴力解决复杂度较高,但通过优化可以达到O(n^(1/2) * logn)的时间复杂度。

传送门http://www.lydsy.com/JudgeOnline/problem.php?id=4028
这题十分神奇……
一开始我考虑线段树,后来又考虑分块。。
但是我死在了这么一个问题上:
知道每一块的GCD和XOR,那怎么查询?
相当于gcd(之前的GCD,这一块某处的前缀GCD)*(之前的XOR^这一块某处的XOR)=x。。。
然后我就爆炸了- -
据zzh和tdl等大神的指点,我终于明白了怎么回事。。
他们也是分块,每一块记录GCD和XOR,同时(注意!)开一个哈希表(map也可)记录这一块里的每个前缀XOR值(若有重复取位置靠前的)。
因为整个数列的前缀GCD往后是单调不上升的,且最多有logN种取值(每次变化至少处以2),所以我们可以采取如下策略:
从左至右扫描每个块,要是当前块和前一个块的GCD不同,则暴力扫一遍;
否则可知这一个块里的每个前缀GCD都和前一个块的GCD相等,做如下变换:
XOR*GCD=x,XOR=到当前块之前的异或和(tmpXOR)^当前块里某个前缀异或和xor,GCD=到当前块之前的gcd(tmpGCD)。
tmpXOR ^ xor=xtmpGCD
tmpXOR ^ xtmpGCD=xor
只需在哈希表内查询有没有这个xor即可。
因为最多 logn 次暴力,每次暴力 O(n) ,所以询问复杂度为
O(nlogn) ,修改复杂度为 O(n)

不知道为什么我的程序总是要跑两秒半才出解啊啊啊啊啊啊!!!!
在没解决这个问题之前先不放代码了哈……

### NOIP2015 运输计划 BZOJ4326 题解分析 #### 问题背景 该问题是经典的图论优化问题之一,主要考察树结构上的路径操作以及高效的数据处理能力。题目要求在一个由 $n$ 个节点组成的无向连通树中找到最优的一条边将其改造为虫洞(通过此边不需要耗费时间),从而使得给定的 $m$ 条运输路径中的最长耗时最小化。 --- #### 解决方案概述 解决这一问题的核心在于利用 **二分答案** 和 **树上差分技术** 的组合来实现高效的计算过程。以下是具体的技术细节: 1. **二分答案**: 设当前目标是最小化的最大路径长度为 $T_{\text{max}}$。我们可以通过二分的方式逐步逼近最终的结果。每次尝试验证是否存在一种方式将某条边改为虫洞后使所有路径的最大值不超过当前设定的目标值 $mid$[^1]。 2. **路径标记与统计**: 使用树上差分的思想对每一条路径进行标记并快速统计受影响的情况。假设两点之间的最近公共祖先 (Lowest Common Ancestor, LCA) 是 $r = \text{lca}(u_i, v_i)$,则可以在三个位置分别施加影响:增加 $(u_i + 1), (v_i + 1)$ 同时减少 $(r - 2)$。这种操作能够有效覆盖整条路径的影响范围,并便于后续统一查询和判断[^1]。 3. **数据结构支持**: 结合线段树或者 BIT (Binary Indexed Tree),可以进一步加速区间修改和单点查询的操作效率。这些工具帮助我们在复杂度范围内完成大量路径的同时更新和检索需求[^2]。 4. **实际编码技巧**: 实现过程中需要注意一些边界条件和技术要点: - 正确维护 DFS 序列以便映射原树节点到连续编号序列; - 准备好辅助函数用于快速定位 LCA 节点及其对应关系; - 编码阶段应特别留意变量初始化顺序及循环终止逻辑以防潜在错误发生。 下面给出一段基于上述原理的具体 Python 实现代码作为参考: ```python from collections import defaultdict, deque class Solution: def __init__(self, n, edges): self.n = n self.graph = defaultdict(list) for u, v, w in edges: self.graph[u].append((v, w)) self.graph[v].append((u, w)) def preprocess(self): """Preprocess the tree to get dfs order and lca.""" pass def binary_search_answer(self, paths): low, high = 0, int(1e9) best_possible_time = high while low <= high: mid = (low + high) // 2 if self.check(mid, paths): # Check feasibility with current 'mid' best_possible_time = min(best_possible_time, mid) high = mid - 1 else: low = mid + 1 return best_possible_time def check(self, limit, paths): diff_array = [0]*(self.n+1) for path_start, path_end in paths: r = self.lca(path_start, path_end) # Apply difference on nodes based on their relationship. diff_array[path_start] += 1 diff_array[path_end] += 1 diff_array[r] -= 2 suffix_sum = [sum(diff_array[:i]) for i in range(len(diff_array)+1)] # Verify whether any edge can be modified within given constraints. possible_to_reduce_max = False for node in range(1, self.n+1): parent_node = self.parent[node] if suffix_sum[node]-suffix_sum[parent_node]>limit: continue elif not possible_to_reduce_max: possible_to_reduce_max=True return possible_to_reduce_max # Example usage of class methods would follow here... ``` --- #### 总结说明 综上所述,本题的关键突破点在于如何巧妙运用二分策略缩小搜索空间,再辅以恰当的树形结构遍历技术和差分手段提升整体性能表现。这种方法不仅适用于此类特定场景下的最优化求解任务,在更广泛的动态规划领域也有着广泛的应用前景[^3]。 ---
评论 7
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值