莫比乌斯反演
pubgoso
这个作者很懒,什么都没留下…
展开
-
牛客挑战赛46 E.反演
设d(n)=∑i=1n[nmod i=0],s(n)=∑i=1nnid(n)=\sum\limits_{i=1}^n[n\mod i=0],s(n)=\sum\limits_{i=1}^n\frac{n}{i}d(n)=i=1∑n[nmodi=0],s(n)=i=1∑nin ans=∑i=1n∑j∣md(ij)=∑i=1n∑j∣m∑x∣i∑y∣j[(x,y)=1]ans=\sum\limits_{i=1}^n \sum\limits_{j|m}d(ij)=\sum\limits_{i=1}^n \s原创 2020-12-15 19:41:00 · 182 阅读 · 2 评论 -
hdu5528 Count a*b
题目链接 #pragma GCC optimize(2) #pragma GCC optimize(3) #include <bits/stdc++.h> using namespace std; typedef unsigned long long LL; const int N = 2e5 + 10; #define fi first #define se second #define pb push_back #define wzh(x) cerr<<#x<<'原创 2020-12-03 10:22:16 · 159 阅读 · 0 评论 -
牛客练习赛62 D.brz的函数
∑i=1n∑j=1nμ(ij)=∑i=1nμ(i)∑j=1nμ(j)[gcd(i,j)=1]=∑i=nnμ(i)∑j=1nμ(j)∑d∣gcd(i,j)μ(d)更改枚举顺序得:∑i=nnμ(i)∑d∣iμ(d)∑j=1⌊nd⌋μ(jd)再次更改枚举顺序得:∑i=1nμ(i)∑j=1⌊ni⌋μ(ij)∑j=1⌊ni⌋μ(ij)设S(n,m)=∑j=1mμ(jn)代入原式得:∑i=1nμ(i)S2(i,ni)通过预处理之后可以O(n)计算每组的答案。但是题目显然是要处理出所有答案。我们观察发现计算答案一定是这样原创 2020-11-06 22:54:02 · 169 阅读 · 0 评论