整除分块

整除分块

问题提出

最近在学习莫比乌斯反演,发现基本上所有有关莫比乌斯反演的题目,都会涉及一个小知识点整除分块

所以在学习毒瘤莫比乌斯之前学习一下整除分块是很有必要的

然后我们就先来学习这个问题

我们来讨论一个有趣的问题

\(\sum_{i=1}^N \lfloor \frac Ni \rfloor ,N \leq 10^{12}\)

问题解决

显然不能直接做

经过一番冷静推理暴力打表,我们发现以下性质:

·1.\(\large \lfloor \frac Ni \rfloor\)向下取整最多只有\(2\sqrt{N}\)种取值

证明:对于\(i\le \sqrt{N}\),只有\(\sqrt{N}\)种,对于\(i>\sqrt{N},\large{\frac Ni}<\sqrt{N}\),也只有\(\sqrt {N}\)种取值,共计\(2 \sqrt{N}\)种取值

·2.设\(\large \lfloor \frac N{i'} \rfloor\)\(\large \lfloor \frac N{i} \rfloor\)相等,则\(i'\)的最大值\(max(i')= \large \left \lfloor \frac N{\left \lfloor \frac Ni \right \rfloor } \right \rfloor\)

证明:设\(\large{ \lfloor \frac Ni \rfloor}=k\),于是可以写成\(ki+p=N,1\le p<i\)的形式,若有\[\large{\lfloor \frac N{i+d} \rfloor}=k\],于是\[k(i+d)+p'=N\] ,即有\(p'=p-kd\),则\(能得到的最大值为d能得到的最大值为 \large \lfloor \frac pk \rfloor\)

于是有

\(\begin{aligned}i'&=i+d_{max} \\ &=i+\lfloor \frac pk \rfloor \\&=i+\left \lfloor \frac {N \;mod\; i}{\lfloor \frac Ni \rfloor} \right \rfloor \\ &=i+\left \lfloor \frac {N-\lfloor \frac Ni\rfloor i}{\lfloor \frac Ni \rfloor} \right \rfloor \\ &=\left \lfloor i + \frac {N-\lfloor \frac Ni\rfloor i}{\lfloor \frac Ni \rfloor} \right \rfloor \\ &=\left \lfloor \frac{\lfloor \frac Ni \rfloor i}{\lfloor \frac Ni \rfloor} + \frac {N-\lfloor \frac Ni\rfloor i}{\lfloor \frac Ni \rfloor} \right \rfloor \\ &=\left \lfloor \frac N{\lfloor \frac Ni \rfloor} \right \rfloor \quad \quad\Box\end{aligned}\)

然后,设两个指针L,R,L的初始值为1,每次令\(\large R=\left \lfloor \frac N{\lfloor \frac NL \rfloor} \right \rfloor\)\(\large (R-L+1)\cdot \lfloor \frac NL \rfloor\)由于\(\large \lfloor \frac NL \rfloor\)只有\(2 \sqrt{N}\)种取值,且单调递减,这最多只有\(2 \sqrt{N}\)个取值不同的,时间复杂度\(O(\sqrt{N})\).

for(int l=1,r;l<=n;l=r+1)
{
    r=n/(n/l);
    ans+=(r-l+1)*(n/l);
}

与其他函数的联系

  • 有时候,可能推出来的式子不一定就是一个很裸的整除分块,可能会与某些积性函数相乘,如:μ,φ...... 这时候,我们就需要对这些函数统计一个前缀和。因为,每当我们使用整除分块跳过一个区间的时候,其所对应的函数值也跳过了一个区间。所以此时,就需要乘上那一个区间的函数值。
  • (当然,如果当出题人想要考考你的数论能力的话,这时就不是统计前缀和这么简单了。可能\(O(N)\)线筛都会TLE,那么我们就需要杜教筛了)
  • (PS:关于杜教筛的内容,我没过多久就会写,那时候再填坑吧!)

转载于:https://www.cnblogs.com/My-snowing/p/10390340.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
SQLAlchemy 是一个 SQL 工具包和对象关系映射(ORM)库,用于 Python 编程语言。它提供了一个高级的 SQL 工具和对象关系映射工具,允许开发者以 Python 类和对象的形式操作数据库,而无需编写大量的 SQL 语句。SQLAlchemy 建立在 DBAPI 之上,支持多种数据库后端,如 SQLite, MySQL, PostgreSQL 等。 SQLAlchemy 的核心功能: 对象关系映射(ORM): SQLAlchemy 允许开发者使用 Python 类来表示数据库表,使用类的实例表示表中的行。 开发者可以定义类之间的关系(如一对多、多对多),SQLAlchemy 会自动处理这些关系在数据库中的映射。 通过 ORM,开发者可以像操作 Python 对象一样操作数据库,这大大简化了数据库操作的复杂性。 表达式语言: SQLAlchemy 提供了一个丰富的 SQL 表达式语言,允许开发者以 Python 表达式的方式编写复杂的 SQL 查询。 表达式语言提供了对 SQL 语句的灵活控制,同时保持了代码的可读性和可维护性。 数据库引擎和连接池: SQLAlchemy 支持多种数据库后端,并且为每种后端提供了对应的数据库引擎。 它还提供了连接池管理功能,以优化数据库连接的创建、使用和释放。 会话管理: SQLAlchemy 使用会话(Session)来管理对象的持久化状态。 会话提供了一个工作单元(unit of work)和身份映射(identity map)的概念,使得对象的状态管理和查询更加高效。 事件系统: SQLAlchemy 提供了一个事件系统,允许开发者在 ORM 的各个生命周期阶段插入自定义的钩子函数。 这使得开发者可以在对象加载、修改、删除等操作时执行额外的逻辑。
SQLAlchemy 是一个 SQL 工具包和对象关系映射(ORM)库,用于 Python 编程语言。它提供了一个高级的 SQL 工具和对象关系映射工具,允许开发者以 Python 类和对象的形式操作数据库,而无需编写大量的 SQL 语句。SQLAlchemy 建立在 DBAPI 之上,支持多种数据库后端,如 SQLite, MySQL, PostgreSQL 等。 SQLAlchemy 的核心功能: 对象关系映射(ORM): SQLAlchemy 允许开发者使用 Python 类来表示数据库表,使用类的实例表示表中的行。 开发者可以定义类之间的关系(如一对多、多对多),SQLAlchemy 会自动处理这些关系在数据库中的映射。 通过 ORM,开发者可以像操作 Python 对象一样操作数据库,这大大简化了数据库操作的复杂性。 表达式语言: SQLAlchemy 提供了一个丰富的 SQL 表达式语言,允许开发者以 Python 表达式的方式编写复杂的 SQL 查询。 表达式语言提供了对 SQL 语句的灵活控制,同时保持了代码的可读性和可维护性。 数据库引擎和连接池: SQLAlchemy 支持多种数据库后端,并且为每种后端提供了对应的数据库引擎。 它还提供了连接池管理功能,以优化数据库连接的创建、使用和释放。 会话管理: SQLAlchemy 使用会话(Session)来管理对象的持久化状态。 会话提供了一个工作单元(unit of work)和身份映射(identity map)的概念,使得对象的状态管理和查询更加高效。 事件系统: SQLAlchemy 提供了一个事件系统,允许开发者在 ORM 的各个生命周期阶段插入自定义的钩子函数。 这使得开发者可以在对象加载、修改、删除等操作时执行额外的逻辑。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值