首页 > 其他分享 >学习笔记_莫比乌斯反演

学习笔记_莫比乌斯反演

时间:2023-08-11 19:56:59浏览次数:46  
标签:lfloor frac 乌斯 sum rfloor mu 反演 莫比

前置知识:整除分块

莫比乌斯函数(\(\mu\))

$\mu(d)=\begin{cases}
1 & (d=1)
\\
(-1)^{k} & \forall C_i=1
\\
0 & \exists C_i>1
\end{cases} $

性质

1.对于任意正整数\(n\),$$\sum_{d|n} \mu(d) = [n=1]$$
[]是一个返回bool型值的判断,当[]内条件成立时返回1,否则返回0.

2.对于任意正整数\(n\),

\[\sum_{d|n}\frac{\mu(d)}{d}=\frac{\phi(n)}{n} \]

代码实现


int num_prime;
//质数个数
bool noprime[N];
//判断每个数是否是质数
int mu[N],prime[N];
//mu[]莫比乌斯函数   prime[]质数

void Wonder_of_U(int n)
{
    mu[1]=1;
    for(int i=2;i<=n;++i)
    {

        if(!noprime[i])
        {
            p[++num_prime]=i;
            mu[i]=-1;
        }
        for(int j=1,k;j<=num_prime && (k=i*prime[j])<=n;++j)
        {
            noprime[k]=true;
            if(!(i % prime[j]))break;
            else mu[k]=-mu[i];
        }
    }
}

莫比乌斯反演

莫比乌斯定理

定义\(F(n)\)和\(f(n)\)是定义在非负整数集合上的两个函数。他们满足:

\[F(n)= \sum_{d|n} f(d) \]

那么一定存在

\[f(n)=\sum_{d|n}\mu(d)F(\lfloor \frac nd \rfloor) \]

此即为莫比乌斯定理.

证明

由\(F(n)\)和\(f(n)\)定义可得

\[F(\lfloor \frac nd \rfloor)=\sum_{i|{\lfloor \frac nd \rfloor}}f(i) \]

那么

\[\sum_{d|n}\mu(d)F(\lfloor \frac nd \rfloor)=\sum_{d|n}\mu(d)\sum_{i|{\lfloor \frac nd \rfloor}}f(i)=\sum_{i|n}f(i) \sum_ {d|{\lfloor \frac ni\rfloor}}\mu(d)=f(n) \]

证毕.


还有个式子.
当\(F(n)\)和\(f(n)\)满足

\[F(n)=\sum_{n|d}f(d) \]

得到

\[f(n)=\sum_{n|d}\mu (\frac dn)F(d) \]

标签:lfloor,frac,乌斯,sum,rfloor,mu,反演,莫比
From: https://www.cnblogs.com/Cayde-6/p/17623827.html

相关文章

  • 单位根反演
    一个非常不阳间的关于单位根的性质:\[\forallk,[n∣k]=\frac{1}{n}\sum_{i=0}^{n−1}\omega^{ik}_{n}\]证明:若\(n∣k\),则有\[Ans=\frac{1}{n}\sum_{i=0}^{n−1}ω^{ink^′}_n=\frac{1}{n}\sum_{i=0}^{n−1}1=1\]否则,等比数列求和有\[Ans=\frac1n\frac{ω^0_n−ω^{nk}_n}......
  • 莫比乌斯反演
    机房卷哥让我写的。大部分是习题总结。不如Wiki两个积性函数\(f,g\),它们的狄利克雷卷积(Dirichlet卷积)为\[h(n)=(f*g)(n)=\sum_{d|n}f(d)g(\frac{n}{d})\]记为\(h=f*g\).狄利克雷卷积满足交换律,结合律,且得到的函数也是积性函数。定义\[\epsilon(n)=\lb......
  • 正泰电力携手图扑:VR 变电站事故追忆反演
    VR(VirtualReality,虚拟现实)技术作为近年来快速发展的一项新技术,具有广泛的应用前景,支持融合人工智能、机器学习、大数据等技术,实现更加智能化、个性化的应用。在电力能源领域,VR技术在高性能计算机和专有设备支持下,已被应用于变电站的培训和安全管理。大力提升了变电站的运行效......
  • 算法学习笔记(24): 狄利克雷卷积和莫比乌斯反演
    狄利克雷卷积和莫比乌斯反演看了《组合数学》,再听了学长讲的……感觉三官被颠覆……目录狄利克雷卷积和莫比乌斯反演狄利克雷卷积特殊的函数函数之间的关系除数函数和幂函数欧拉函数和恒等函数莫比乌斯函数和欧拉函数卷积的逆元莫比乌斯函数与莫比乌斯反演求法数论分块(整除分......
  • ENVI大气校正方法反演Landsat 7地表温度
    本文介绍基于ENVI软件,实现对Landsat7遥感影像加以大气校正方法的地表温度反演操作。目录1图像前期处理与本文理论部分2实际操作2.1植被覆盖度计算2.2地表比辐射率计算2.3相同温度下黑体辐射亮度值计算2.4地表真实温度计算2.5图像导出3专题地图制作4不同地物地表温度对......
  • 莫比乌斯反演学习笔记
    莫比乌斯反演数论函数数论函数是指定义域为正整数的一类函数。基本的数论函数恒等函数\(I(n)=1\)元函数\(e(n)=[n=1]\)单位函数\(id(n)=n\)莫比乌斯函数$$\mu(n)=\begin{cases}0,&n的约数中包含大于1的完全平方数\(-1)^k,&k为x含有的质因子种类数\end{cases}$$欧......
  • 莫比乌斯反演
    函数定义\[\begin{aligned}\\I(n)&=1\\\epsilon(n)&=[n=1]\\id(n)&=n\end{aligned}\]卷积定义\[(f*g)(x)=\sum\limits_{d|n}f(d)g(\dfrac{n}{d})\]有以下性质:\[\begin{aligned}f*g&=g*f\\(f*g)*h&=f*(g*h)\......
  • 拉格朗日反演
    前置引理首先设\(\Bbb{C}[[x]]=\{\sum\limits_{i\ge0}a_ix^i|a_i\in\Bbb{C}\},\Bbb{C}((x))=\{\sum\limits_{i\ge-N}a_ix^i|N\inZ,a_i\in\Bbb{C}\}\)有以下引理:\(f(x)\in\Bbb{C}((x)),[x^{-1}]f'(x)=0\)显然,因为没有一个幂次项求导后是\(-1\)次项。\(f(......
  • 二项式反演
    第一种形式:\[f(n)=\sum\limits_{i=0}^n\dbinom{n}i{g(i)}\Leftrightarrowg(n)=\sum\limits_{i=0}^n(-1)^{n+i}\dbinomnif(i)\]证明:\[\begin{aligned}f(n)&=\sum\limits_{i=0}^n\dbinom{n}i{g(i)}=\sum\limits_{i=0}^n\dbinomni\sum\limits_{j=0......
  • 单位根反演
    命题如下:$$\forallk\inZ,[n|k]=\frac{1}{n}\sum\limits_{i=0}{n-1}\omega_n$$证明:设$[n|k]=1$,则根据单位根性质,我们可以得到:$$\sum\limits_{i=0}{n-1}\omega_n=n$$设$[n|k]=0$,则:$$\sum\limits_{i=0}{n-1}\omega_n=\frac{\omega_n{nk}-1}{\omega_n-1}=0$$由此可知......