莫比乌斯反演
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∑ninans=∑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 · 171 阅读 · 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 · 125 阅读 · 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 · 155 阅读 · 0 评论