标签:
0.前置知识
一些函数
- \(1(n)=1\)
- \(id(n)=n\)
- \(\sigma(n)\)为\(n\)的约数和
狄利克雷卷积
定义两个数论函数运算\(*\),
若\(h=f*g\),则
\[
h(n)=\sum_{d|n}f(d)g(\frac nd)
\]
它满足一些性质:
\(f*g=g*f\)
\(f*(g*h)=(f*g)*h\)
因为\(\sum_{(ij)k=n}(f(i)g(j))h(k)=\sum_{i(jk)=n}f(i)(g(j)h(k))\)
\((f+g)*h=f*h+g*h\)
\((xf)*g=x(f*g)\)
设单位元\(\epsilon(n) = [n = 1]\),则\(\epsilon*f=f\)
对于每个\(f(1)\neq 0\)的函数\(f\),存在逆元\(g\)使得\(f*g=\epsilon\)
如何求一个函数的逆元呢?
令\(g(n)\)满足以下式子:
\[
g(n)=\frac 1{f(1)}\left(\epsilon(n)-\sum_{i|n,i\neq1}f(i)g(\frac ni)\right)
\]
这样的话:
\[
\sum_{i|n}f(i)g(\frac ni) \\
=f(1)g(n) + \sum_{i|n, i\neq 1} f(i)g(\frac ni) \\
=\epsilon(n) \\
\therefore f*g=\epsilon
\]
1.莫比乌斯反演
定义一个函数\(\mu\)使得\(\mu*1=\epsilon\)
这样的话,如果\(g*1=f\),则\(f*\mu=g\)
即:如果\(f(n)=\sum_{d|n}g(d)\),则\(g(n)=\sum_{d|n}\mu(d)f(\frac nd)\)
是不是很简单
2.题目
若\(n<m\)求:
#1
\[ \sum_{i=1}^n\sum_{j=1}^m[gcd(i,j)=1] \]
设\(f=\epsilon\),\(\because \epsilon=1*\mu,\;\therefore g=\mu\)
于是原式变为:
\[
\sum_{i=1}^n\sum_{j=1}^m\sum_{d|i,d|j}\mu(d) \\
=\sum_{d=1}^n \mu(d)\sum_{i=1}^n\sum_{j=1}^m[d|i][d|j] \\
=\sum_{d=1}^n \mu(d)\left[\frac nd\right]\left[\frac md\right]
\]
预处理\(\mu\),数论分块就可以\(O(\sqrt n)\)求了
#2
\[ \sum_{i=1}^n\sum_{j=1}^m gcd(i,j) \]
设\(f=id\),\(\because id=1*\varphi,\;\therefore g=\varphi\)
于是原式变为:
\[
\sum_{d=1}^n\varphi(d)\left[\frac nd\right]\left[\frac md\right]
\]
#3
\[ \sum_{i=1}^n\sum_{j=1}^m\sigma(gcd(i,j)) \]
设\(f=\sigma\),\(\because\sigma=1*(\mu*\sigma),\;\therefore g=\mu*\sigma\)
于是原式变为:
\[
\sum_{d=1}^ng(d)\left[\frac nd\right]\left[\frac md\right]
\]
因为积性函数可以线性筛,所以也可以预处理
#4
\[ \sum_{i=1}^n\sum_{j=1}^mf(gcd(i,j)) \]
设\(g=\mu*f\),如果我们能预处理出\(f\),那么可以这样求\(g\)
(当然是蒯的啦)
void get_g_1(int N, const int *f, int *g)
{
for (int i = 1; i <= N; i++) g[i] = 0;
for (int i = 1; i <= N; i++)
for (int j = 1; i * j <= N; j++)
g[i * j] = (g[i * j] + mu[i] * f[j]) % mod;
} // 依照定义,O(nlogn)
void get_g_2(int N, const int *f, int *g)
{
for (int i = 1; i <= N; i++) g[i] = f[i];
for (int i = 1; i <= N; i++)
for (int j = 2; i * j <= N; j++)
g[i * j] = (g[i * j] - g[i]) % mod;
} // 类似求狄利克雷卷积逆的方式,不需要线性筛 mu ,O(nlogn)
void get_g_3(int N, const int *f, int *g)
{
for (int i = 1; i <= N; i++) g[i] = f[i];
for (int i = 0; i < prime_count; i++)
for (int j = N / prime[i]; j >= 1; j--)
g[j * prime[i]] = (g[j * prime[i]] - g[j]) % mod;
} // Magic! O(nloglogn)
于是原式变为:
\[
\sum_{d=1}^ng(d)\left[\frac nd\right]\left[\frac md\right]
\]
标签: 来源: https://www.cnblogs.com/cj-xxz/p/10182999.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。