数论专题(1)数论函数,整数分块
admin
2024-05-01 13:32:12
0次

从今天起,我们将要开始数论的学习,是不是感觉很难?难的话就听我讲吧,讲了后就不难了(bushi)

数论函数定义 (数论函数)

数论函数的定义:在全体正整数(或者整数)上定义的函数称作数论函数。 积性的定义:若 gcd(a,b)=1,则f(ab)=f(a)f(b)。 举个栗子: 欧拉函数\varphi (x) :                                  ​​​​​​​        1(x)=1,Id(x)=x,I(x)=[x=1]。 接下来,我们再举一个函数,狄利克雷卷积: 数论函数f(n)和g(n)的狄利克雷卷积h(n)定义为:         ​​​​​​​        ​​​​​​​        ​​​​​​​                        h(n)=\sum_{d|n}^{}f(d)g(\frac{n}{d}) 记作h=f*g; 定理:两个积性函数的狄利克雷卷积也是积性函数。 那我们是如何证明的呢,那可就说来话长了: h(xy)=\sum_{d_{1}|x}^{}f(d_{1})g(\frac{x}{d_{1}})=\sum_{d_{1}|x}^{}\sum_{d_{2}|y}^{}f(d_{1}d_{2})g(\frac{xy}{d_{1}d_{2}})=\sum_{d_{1}|x}^{}f(d_{1})g(\frac{x}{d_{1}})\sum_{d_{2}|y}^{}f(d_{2})g(\frac{y}{d_{2}})=h(x)h(y) 看到这么长一串的函数,我的心就发抖……… 定理:两个积性函数的对应位置相乘也是积性函数。 这个就很显然了,证明我就不给出来了,你们自己想着吧。 所以我们可以定义更多的数论函数,以及用卷积描述它们之间的关系: 例如:         ​​​​​​​        ​​​​​​​        ​​​​​​​        ​​​​​​​        ​​​​​​​        d(x)=\sum_{d|n}1,即d=1*1         ​​​​​​​        ​​​​​​​        ​​​​​​​        ​​​​​​​        ​​​​​​​        n=\sum_{d|n}\varphi (d),即Id=\varphi *1 接下来,我们来研究一下,卷积的性质: 定理(卷积运算律): 交换律:设有两个积性函数f,g则f*g=g*f 结合律:设有两个积性函数f,g,h则(f*g)*h=f*(h*g) 这不就是乘法结合律吗?小学生都会好不好…… 交换律的证明很显然。 结合律的证明可以把式子改写成矩阵形式,然后用矩阵的结合律 来证明。应该大力推式子也可以。 单位元的定理: I是单位元:对任意的数论函数f,f*I=I*f=f; 证明过程如下: ​​​​​​​        ​​​​​​​        ​​​​​​​        ​​​​​​​        ​​​​​​​        ​​​​​​​        f(n)=\sum_{d|n}^{})f(d)I(\frac{n}{d}) 接下来又来到了下一个函数,狄利克雷逆: 若f*g=I,则数论函数f,g互为彼此的狄利克雷逆。

整数分块

接下来,我们进入下一个章节的学习:整数分块

定理:对于任意的i, \left \lfloor \frac{n}{i} \right \rfloor 只有O(\sqrt{n})种取值

证明过程如下:

当1\leq i\leq \sqrt{n}时,\frac{n}{i}只有​​​​​​​​​​​​​​O(\sqrt{n})种取值 

当i> \sqrt{n}​​​​​​​时,\frac{n}{i}\leq \sqrt{n},也只有O(\sqrt{n})种取值

 综上,只有O(\sqrt{n})种取值

下面举个栗子:

计算下取整分式的和式,计算\sum_{i}^{n}=1\left \lfloor \frac{n}{i} \right \rfloor

由于\left \lfloor \frac{n}{i} \right \rfloor只有O(\sqrt{n})种取值 ,并且,使得 \left \lfloor \frac{n}{i} \right \rfloor取相同的取值的i也是一段一段的,所以我们只需要一段一段地计算即可。

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

哎哟,终于讲完了,累死我了,咱么今天就讲到这里,下篇博文我们会讲莫比乌斯反演,记得来收看哦!! 

 

 

 

 

相关内容

热门资讯

两只大牛股,紧急发声:不属实!... 5天3板超声电子:“超声电子高频板通过英伟达认证”情况不属实,目前无产品供货给英伟达9月13日晚间,...
科创板创新药企BD授权交易总额... 2026年半年度财报季已收官,在创新药企业以业绩扭亏和BD兑现吸引市场目光的同时,以CXO为代表的上...
蔚来李斌1-8月交付增近58%... 9月9日消息,上海蔚来汽车有限公司副总裁品牌与传播资深专家马麟在微博透露,9月8日下午,蔚来创始人、...
京东物流部署过万前置仓应对手机... IT之家 9 月 13 日消息,据京东物流黑板报,9 月,多个品牌旗舰手机集中上新。针对消费者希望尽...
特朗普让泽连斯基勿袭击俄柴油燃... 9月13日,当地时间9月13日,美国总统特朗普要求乌克兰总统泽连斯基停止袭击俄罗斯柴油燃料设施,称相...
下周,年内最便宜新股启动打新 ... 2026.09.13本文字数:1467,阅读时长大约2分钟作者 |第一财经 安卓封图 |AI生成下周...
北京:激活地理信息要素赋能城市... 9月13日,从北京市规划和自然资源委员会获悉,《北京市测绘地理信息条例》已由北京市第十六届人民代表大...
接下来该如何逃顶?建议看看这篇... 最近这段时间,很多朋友的心里其实都是悬着的。行情走到这个位置,大家见面聊得最多的就是,接下来到底该怎...
武汉:到2035年武汉区域科技... 9月13日,武汉市委十四届十三次全会召开,审议通过了《关于打造具有全球竞争力的开放创新生态 高水平建...
卡塔尔股市收盘上涨0.25% 9月13日,卡塔尔证券交易所主要股指今日收盘上涨0.25%,报9848.95点,上涨24.46点。全...