杜教筛

杜教筛

 

1.概述

 

杜教筛是用以解决积性函数前缀和的算法。

在学习了莫比乌斯反演之后,杜教筛的过程就会显得简单而自然。

 

2.基本形式 

 

对于积性函数f,我们定义如下函数:

S(x)=\sum_{i=1}^{x}f(i)

构造积性函数 h,g,使得 f*g=h

显然 :

\sum_{i=1}^{n}h(i)=\sum_{i=1}^{n}(f*g)(i)

进一步转化:

\sum_{i=1}^{n}h(i)=\sum_{d=1}^{n}g(d)\cdot \sum_{i=1}^{\left \lfloor \frac{n}{d} \right \rfloor}f(i)

也就是:

\sum_{i=1}^{n}h(i)=\sum_{d=1}^{n}g(d)\cdot S(\left \lfloor \frac{n}{d} \right \rfloor)

把右式的第一项提出:

\sum_{i=1}^{n}h(i)=g(1)S(n)+\sum_{d=2}^{n}g(d)\cdot S(\left \lfloor \frac{n}{d} \right \rfloor)

进一步整理,得:

g(1)S(n)=\sum_{i=1}^{n}h(i)-\sum_{d=2}^{n}g(d)\cdot S(\left \lfloor \frac{n}{d} \right \rfloor)

显然,如果我们能够快速算出任意的  h(i) ,就能够通过数论分块在O(n^{\frac{2}{3}})的时间内求出S

 

3.一些常见函数的杜教筛形式

    3.1  \mu函数的前缀和

有一个显而易见的结论:\mu*I=\epsilon

设  S(n)=\sum_{i=1}^{n}\mu(i)

由杜教筛易得:

S(n)=1-\sum_{d=2}^{n}S(\left \lfloor \frac{n}{d} \right \rfloor)

    3.2  \varphi函数的前缀和

有一个显而易见的结论:\varphi*I=id

设  S(n)=\sum_{i=1}^{n}\varphi(i)

由杜教筛易得:

S(n)=\frac{n(n-1)}{2}-\sum_{d=2}^{n}S(\left \lfloor \frac{n}{d} \right \rfloor)

4.一些总结

 

杜教筛作为一个数论算法,代码实现不算太难,但有几个小细节需要注意:

  • S(i)  是需要全部记录的,需要用到unordered_map或者 map(多一个log)。
  • 实测unordered_map、map,时间差不多
  • 容易爆long long

平时遇到积性函数求和,并且n高达10^9,10^{10},10^{11},就可以尝试用杜教筛来解。

当然杜教筛有一个弊端:没有一个系统的构造函数的方法。

如果遇到一个难以构造的积性函数,可以尝试使用 min_25筛、洲阁筛,(啥?你不会,我也不会qwq),那自求多福吧(QAQ~)。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/315903.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

ML.NET机器学习、API容器化与Azure DevOps实践(二):案例

在上文中,我简单地介绍了机器学习以及ML.NET的相关知识,从本讲开始,我会基于一个简单的案例:学生成绩预测,来介绍使用ML.NET进行机器学习以及API部署的基本过程。本案例的数据来源为加州大学尔湾分校的机器学习公开样本…

业界萌新对斯坦纳树的小结

业界萌新对斯坦纳树的小结 0.简介 斯坦纳树问题是组合优化问题,与最小生成树相似,是最短网络的一种。最小生成树是在给定的点集和边中寻求最短网络使所有点连通。而最小斯坦纳树允许在给定点外增加额外的点,使生成的最短网络开销最小。 ——…

cf1208E. Let Them Slide

cf1208E. Let Them Slide 题意: 都放在一个长度为W的框里面。有n个序到,第i个序列的长度是1。这些序到并排放在一起,每一个序列都放在一个长度为w的框里 这些序列可以在框里面滑动,但是不能划出框。 对于每一个位置&#xff0…

结合eShopOnWeb全面认识领域模型架构

一.项目分析在上篇中介绍了什么是"干净架构",DDD符合了这种干净架构的特点,重点描述了DDD架构遵循的依赖倒置原则,使软件达到了低藕合。eShopOnWeb项目是学习DDD领域模型架构的一个很好案例,本篇继续分析该项目各层的职…

SOS_dp算法

Codeforces博客 简介: 前置知识:状压dp Sum over Subsets dynamic programming,简称Sos dp,状压dp的一种 用一个列题引出SOS dp: 给你一个由2N2^N2N个整数组成的确定数组A,我们需要计算对于任意的x,F(x)所…

微软开源Bing搜索背后的关键算法

微软今天宣布开源了一项 Bing 搜索背后的关键算法 —— SPTAG,它使 Bing 能够快速将搜索结果返回给用户。仅在几年前,网络搜索很简单,用户输入几个关键词然后浏览结果页面。现如今,这些用户可能会在手机上拍照并将其放入搜索框中&…

Stern-Brocot Tree

Stern-Brocot Tree 0.简介 Stern-Brocot Tree,俗称SB树(滑稽)。它能够表示出所有的最简分数,如下图。 1.一些规律 显然,对于两个相邻的最简分数 可以得到另一个最简分数 这样就可以在Stern-Brocot Tree上表示出所有…

FWT(快速沃尔什变换)

文章目录引入:or卷积and卷积xor卷积IFWT模板:例题:引入: FFT/NTT是用来解决∑ijkA[i]B[j]\sum_{ijk}A[i]B[j]∑ijk​A[i]B[j]的式子 而FWT是用来解决Ci∑j⊕kiAjBkC_i\sum_{j⊕ki}A_jB_kCi​∑j⊕ki​Aj​Bk​ ​ FWT是一种用于处…

教你自制.NET Core Global Tools

点击上方蓝字关注“汪宇杰博客”命令行是程序员装逼利器,.NET Core也可以写命令行程序,但是如何分发给其他程序员使用,一直是个问题。现在,有了.NET Core Global Tools,可以很方便的解决分发问题,我们来看看…

三点间LCA

三点间LCA 1.直接上题——jzoj5883. 【NOIP2018模拟A组9.25】到不了 Dscription wy 和 wjk 是好朋友。 今天他们在一起聊天,突然聊到了以前一起唱过的《到不了》。 “说到到不了,我给你讲一个故事吧。” “嗯?” “从前,神和凡人…

微软拥抱开源,Win10为啥要引入真Linux4.X内核?

来源 | 异步 | 文末赠书2019 年微软 Build 开发者大会在雷德蒙德召开。继将 Bash shell、原生 OpenSSH、WSL 引入 Windows,以及在微软商店提供 Ubuntu、SUSE Linux 和 Fedora 等发行版,微软又宣布了一个重大的决定 —— 将完整的 Linux 内核引入 Windows…

F.孤独(牛客小白月赛39)

F.孤独(牛客小白月赛39) 题意: 给定一棵树,寻找一个路径,将断掉所有与这个路径上的点相连的边,使得剩下的最大连通块的大小最小 题解: 这题有点印象,感觉做过,至少这…

分布式 - 分布式系统的特点

20世纪60年代,IBM研发了System 360架构大型机,与同时期的波音707、福特汽车誉为商业三大成就,凭借其卓越的性能和良好的稳定性,开启了大型机的时代,诞生了非常多的集中式系统,采用单机架构,有非…

[WC2011][BZOJ2115] Xor

BZOJ2115 Xor 题目描述: 题目大意: 给定一张 n 个点 m 条边的无向带权连通图,求一条从点 1 到点 n 的路径,使得经过的边权异或和最大。 路径可以经过重复点和重复边,当一条边被重复经过时也会相应地被 xor 多次。 s…

.NET Core 3.0 可回收程序集加载上下文

.NET诞生以来,程序集的动态加载和卸载都是一个Hack的技术,之前的NetFx都是使用AppDomain的方式去加载程序集,然而AppDomain并没有提供直接卸载一个程序集的API,而是要卸载整个AppDomain才能卸载包含在其中的所有程序集。然而卸载整…

ADPC2-G 希望

希望 题意: 有A,B两棵树,对于一个1到n的全排列a[i],让树A中的点i和树B的节点a[i]连一条边,希望指数:两棵树和新加入的边构成的图中,环长为m的环的个数。数组a[]可以任意交换位置,且任意&#…

.Net Core中依赖注入服务使用总结

一、依赖注入引入依赖注入的目的是为了解耦和。说白了就是面向接口编程,通过调用接口的方法,而不直接实例化对象去调用。这样做的好处就是如果添加了另一个种实现类,不需要修改之前代码,只需要修改注入的地方将实现类替换。上面的…

[BZOJ3944] Sum

[BZOJ3944] Sum 题目描述&#xff1a; solution 裸的杜教筛。 唯一的坑点在于卡常。 似乎ans1和ans2都杜教筛超时了。 然而用杜教筛求出ans2&#xff0c;并用求出ans1不超时&#xff1f;&#xff01;&#xff1f;。 ​ ​ ​ #include<bits/stdc.h> using namespac…

SuperSocket 2.0 Preview1 发布,.NET Socket服务器框架

今天&#xff0c;SuperSocket的作者发布了2.0版本的第一个预览版。SuperSocket 2.0 是一个经过全新设计的&#xff0c;第一个完全基于.NET Core的版本。作者正在积极尝试提供更简单易用的API的同时&#xff0c;尽量保证与老版本相似的原汁原味的开发体验。新的版本中亦删除了一…

[LG P2519][BZOJ2298][HAOI2011]problem a

[LG P2519][BZOJ2298][HAOI2011]problem a 题目描述 一次考试共有n个人参加 第i个人说&#xff1a;"有ai个人分数比我高&#xff0c;bi个人分数比我低。" 问最少有几个人没有说真话(可能有相同的分数) 输入输出格式 输入格式&#xff1a; 第一行一个整数n&…