目录

前言:发现最近都没怎么写博客,,,赶紧发篇以前记的笔记凑凑数

卷积

卷积定义:

  如果有数论函数\(f, g\), 那么它们卷积的第\(n\)项为\((f * g) (n)\),设这个卷出来的新函数为h,那么有
  \[h(n) = \sum_{k | n} f(k) g(n / k) = \sum_{ij = n}f(i) g(j)\]
  

一些性质:

  1,积性函数的卷积还是积性函数
  证明: 现有\(f, g\)为积性函数,且\(gcd(p, q) == 1\),求证\(h(p) \cdot h(q) = h(qp).\)
  \[h(p) \cdot h(q) = \sum_{ab = p}f(a) g(b) \cdot \sum_{cd = q}f(c) g(d)\]
  \[= \sum_{a} f(a) g(\frac{p}{a}) \sum_{c} f(c) g(\frac{q}{c})\]
  \[= \sum_{a} \sum_{c} f(a) f(c) g(\frac{p}{a}) g(\frac{q}{c})\]
  \[= \sum_{a} \sum_{c} f(ac) g(\frac{pq}{ac})\]
  因为\(ac \cdot \frac{pq}{ac} = pq\),所以
  \[\sum_{a} \sum_{c} f(ac)g(\frac{pq}{ac}) = \sum_{ac | pq} f(ac)g(\frac{pq}{ac})\]
  \[= \sum_{x|pq}f(x)g(\frac{pq}{x}) = h(pq)\]
  

常见狄利克雷卷积:

\(id = \phi * 1\)
\(d = 1 * 1\)
\(\sigma = id * 1\)
\(e = 1 * \mu\)
\(\phi = id * \mu\)
  

有趣的事情:

  1,FFT(多项式乘法)就是在求卷积。
    可以把\(f,g\)分别看做两个多项式,\(f(i)\)表示\(x^i\)的系数,\(g(i)\)同理。
    那么卷出来的\(h(i)\)就等于两式相乘后\(x^{i}\)的系数,因为卷积定义里面\(h(k)\)就要加上满足\(ij = k\)的\(f(i) g(j)\).

杜教筛

目的:在低于线性的时间内求积性函数前缀和

算法:

  设\(S(n) = \sum_{i = 1}^{n}f_{i}\).
  设有一个积性函数\(g(i)\).则:
  \[\sum_{i = 1}^{n}(g * f)(i) = \sum_{i = 1}^{n} \sum_{d | i} g(d)f(\frac{i}{d})\]
  因为当\(d | i\)时会统计到\(g(d)f(\frac{i}{d})\),因此直接枚举\(d\)和\(i = \frac{i}{d}\),\(i\)的最大值为\(\frac{n}{d}\),因为要满足\(i \cdot d \le n\).所以:
  \[原式= \sum_{d = 1}^{n}g(d) \sum_{i = 1}^{\frac{n}{d}}f(i) = \sum_{d = 1}^{n}g(d)S(\frac{n}{d})\]
  所以:
  \[\sum_{i = 1}^{n}(g * f)(i) = \sum_{d = 1}^{n}g(d)S(\frac{n}{d})\]
  
  因为我们有:
  \[g(1)S(n) = \sum_{i = 1}^{n}g(i)S(\frac{n}{i}) - \sum_{i = 2}^{n}g(i)S(\frac{n}{i})\]
  带入上式得:
  \[g(1)S(n) = \sum_{i = 1}^{n}(g * f)(i) - \sum_{i = 2}^{n}g(i)S(\frac{n}{i})\]
  因此我们只需要找到一个合适的\(g\),带入上式,使得\(\sum_{i = 1}^{n}(g * f)(i)\)可以快速计算,那么就只需要对\(\sum_{i = 2}^{n}g(i)S(\frac{n}{i})\)进行求解就可以得到\(S(n)\),而这个式子中的\(S(\frac{n}{i})\)是可以整数分块求的。
  于是就可以做到\(O(n^{\frac{3}{4}})\)求解了。
  同时为了降低复杂度,可以先预处理一部分\(S\)值。通过证明可得:设我们当前预处理了前\(k\)个S,那么当\(k\)取到\(n^{\frac{2}{3}}\)时可以使得整体复杂度降为\(n^{\frac{2}{3}}\).(后面部分记忆化求解)
  
  关于记忆化:
  1,可以使用map/hash解决。
  2,或者观察到对于同一个\(S(n)\),我们所有要求的\(S\)值都是形如\(S(\frac{n}{x})\)的,并且我们知道\(\lfloor{ \frac{\lfloor{\frac{n}{x}}\rfloor}{y} }\rfloor = \lfloor {\frac{n}{xy}} \rfloor\),因此对于同一个\(n\),我们求的所有\(S(x)\)都可以表示为\(S(\frac{n}{x})\),因此设\(now = \frac{n}{x}\)。那么:
    1,对于\(now <= n^{\frac{2}{3}}\),我们已经预处理出来了。
    2,对于\(now > n^{\frac{2}{3}}\),因为\(now > n^{\frac{2}{3}}\),所以\(x = \frac{n}{now} <= n ^ {\frac{1}{3}}\),因此我们可以令\(sum[x] = S(\frac{n}{x})\)。即若我们有\(S(now)\) 其中$now > n ^ {\frac{2}{3}} $,那么我们可以把这个值存入 \(sum[\frac{n}{now}]\).
    注意:但这样的话,虽然减少一个log,但若有多个不同n值,则每次处理一个新的n值时就需要清空S数组。
    因为虽然对于同一个n值,我们要求的所有S的下标都可以表示为 \(\frac{n}{x}\),但不是任意一个数都可以被表示为 \(\frac{n}{x}\) 的,
    如果这个时候的另一个\(n\)恰好不能被表示,那么答案就会出错

卷积 & 杜教筛的更多相关文章

  1. 中国剩余定理 & 欧拉函数 & 莫比乌斯反演 & 狄利克雷卷积 & 杜教筛

    ssplaysecond的博客(请使用VPN访问): 中国剩余定理: https://ssplaysecond.blogspot.jp/2017/04/blog-post_6.html 欧拉函数: h ...

  2. 洛谷P3768 简单的数学题(莫比乌斯反演+狄利克雷卷积+杜教筛)

    传送门 不会…… 两篇加在一起都看不懂…… https://www.cnblogs.com/cellular-automaton/p/8241128.html https://www.luogu.or ...

  3. [基本操作] Mobius 反演, Dirichlet 卷积和杜教筛

    Dirichlet 卷积是两个定义域在正整数上的函数的如下运算,符号为 $*$ $(f * g)(n) = \sum_{d|n}f(d)g(\frac{n}{d})$ 如果不强调 $n$ 可简写为 $ ...

  4. 【Bzoj3944】杜教筛模板(狄利克雷卷积搞杜教筛)

    题目链接 哇杜教筛超炫的 有没有见过$O(n^\frac{2}{3})$求欧拉函数前缀和的算法?没有吧?蛤蛤蛤 首先我们来看狄利克雷卷积是什么 首先我们把定义域是整数,陪域是复数的函数叫做数论函数. ...

  5. 我也不知道什么是"莫比乌斯反演"和"杜教筛"

    我也不知道什么是"莫比乌斯反演"和"杜教筛" Part0 最近一直在搞这些东西 做了将近超过20道题目吧 也算是有感而发 写点东西记录一下自己的感受 如果您真的 ...

  6. 【Luogu3768】简单的数学题(莫比乌斯反演,杜教筛)

    [Luogu3768]简单的数学题(莫比乌斯反演,杜教筛) 题面 洛谷 \[求\sum_{i=1}^n\sum_{j=1}^nijgcd(i,j)\] $ n<=10^9$ 题解 很明显的把\( ...

  7. 【BZOJ3944】Sum(杜教筛)

    [BZOJ3944]Sum(杜教筛) 题面 求\[\sum_{i=1}^n\mu(i)和\sum_{i=1}^n\phi(i)\] 范围:\(n<2^{31}\) 令\[S(n)=\sum_{i ...

  8. 杜教筛:Bzoj3944: sum

    题意 求\(\sum_{i=1}^{n}\varphi(i)和\sum_{i=1}^{n}\mu(i)\) \(n <= 2^{31}-1\) 不会做啊... 只会线性筛,显然不能线性筛 这个时 ...

  9. 杜教筛进阶+洲阁筛讲解+SPOJ divcnt3

    Part 1:杜教筛进阶在了解了杜教筛基本应用,如$\sum_{i=1}^n\varphi(i)$的求法后,我们看一些杜教筛较难的应用.求$\sum_{i=1}^n\varphi(i)*i$考虑把它与 ...

随机推荐

  1. Docker入门篇(一)之docker基础

    1.Docker 架构 http://blog.csdn.net/u012562943/article/category/6048991/1Docker 使用客户端-服务器 (C/S) 架构模式,使用 ...

  2. SpringBoot-08:SpringBoot采用json的方式实现前后台通用的配置文件

    ------------吾亦无他,唯手熟尔,谦卑若愚,好学若饥------------- 一.需求: 本篇博客是最近笔者做的一个项目,已经上线但是还在不断开发,有些页面上的配置,测试服务器和正式服务器 ...

  3. STM32运行FreeRTOS出现prvTaskExitError错误死机

    文件port.c prvTaskExitError();任务退出错误,一个可能在任务里面写了return,另一个可能任务切换退出问题,入栈和出栈的时候出了问题. static void prvTask ...

  4. Azkaban 工作流调度器

    Azkaban 工作流调度器 1 概述 1.1 为什么需要工作流调度系统 a)一个完整的数据分析系统通常都是由大量任务单元组成,shell脚本程序,java程序,mapreduce程序.hive脚本等 ...

  5. ConfigurationProperties cannot be resolved to a type

    pom.xml 中报错之前: <parent> <groupId>org.springframework.boot</groupId> <artifactId ...

  6. 第四篇 前端学习之JQuery基础

    一 jQuery是什么? jQuery就是一个JavaScript的库. <1> jQuery由美国人John Resig创建,至今已吸引了来自世界各地的众多 javascript高手加入 ...

  7. TPO 02 - Early Cinema

    TPO 02 - Early Cinema NOTE: 主要意思(大概就是主谓宾)用粗体标出:重要的其它用斜体: []中的是大致意思,可能与原文有关也可能无关,但不会离题 目的为训练句子/段落总结能力 ...

  8. IDE看代码,挺好

    初学编程的时候总是收到各种警告:“刚学习编程千万不要用IDE,否则会有xxxxxx的后果”.现在工作后发现使用IDE可以方便编写和查看代码,对于较大的项目来说有很多代码,代码之间的关系也比较复杂,ID ...

  9. CentOS7.x安装Docker1.11.1

    原文发表于cu:2016-05-30 本文属于重发,当前Docker已经分为EE与CE版本,CE版本是17.06.0-ce,最新的CE版本安装方式已略有不同:但可以指定安装版本,如1.11.1,1.1 ...

  10. Python 内置函数介绍

    作者博文地址:http://www.cnblogs.com/spiritman/ Python Built-in Functions