P4213 【模板】杜教筛(Sum)
\(\color{#0066ff}{题 目 描 述}\)
给定一个正整数\(N(N\le2^{31}-1)\)
求
\(\begin{aligned} ans_1=\sum_{i=1}^n\varphi(i) \end{aligned}\)
\(\begin{aligned} ans_2=\sum_{i=1}^n \mu(i) \end{aligned}\)
\(\color{#0066ff}{输 入 格 式}\)
一共T+1行
第1行为数据组数T(T<=10)
第2~T+1行每行一个非负整数N,代表一组询问
\(\color{ #0066ff }{ 输 出 格 式 }\)
一共T行,每行两个用空格分隔的数ans1,ans2
\(\color{#0066ff}{输入样例}\)
6
1
2
8
13
30
2333
\(\color{#0066ff}{ 输 出 样 例}\)
1 1
2 0
22 -2
58 -3
278 -3
1655470 2
\(\color{#0066ff}{数 据 范 围 与 提 示}\)
\(N \leq 2^{31}\)
\(\color{#0066ff}{题 解}\)
前置知识1 : 狄利克雷卷积
对于任意函数f,g,有\(\begin{aligned} h(i) = \sum_{d|i}f(d)*g(\frac{n}{d})\end{aligned}\)
h即为f和g的卷积
常用函数
1、\(i(n) = 1\)
2、\(id(n) = n\)
3、\(e(n)=\left\{\begin{aligned}1\ \ \ n = 1 \\ 0 \ \ \ n \neq 1\end{aligned}\right.\)
4、欧拉函数\(\varphi(n)\)
5、懵逼钨丝函数\(\mu(n)=\left\{\begin{aligned}1\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ n = 1 \\ (-1)^k \ \ \ n由k个不同质数相乘得到\\ 0\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 其它情况\end{aligned}\right.\)
6、\(\sigma(n)=n的约数和\)
7、\(d(n)=n的约数个数\)
常用卷积
1、\(i*\mu = e\)
2、\(e*a=a\)
3、\(\mu * id= \varphi\)
4、\(i*id=\sigma\)
5、\(i*i=d\)
6、\(i*\varphi=id\)
杜教筛
已知\(f(i)\)
用来求\(\begin{aligned}\sum_{i = 1}^n f(i)\end{aligned},n\leq 2^{31}\)
定义\(h(i)=(f*g)(i)=\begin{aligned}\sum_{d|i}f(d)*g(\frac{i}{d})\end{aligned}\)
\(\displaystyle\sum_{i=1}^nh(i)\)
用定义展开
\(=\displaystyle\sum_{i=1}^n\sum_{d|i}g(d)f\left(\frac i d\right)\)
d的范围也是【1.n】的,所以改成枚举d,找它的倍数,这个式子是在求和,找全了就行
\(=\displaystyle \sum_{d=1}^ng(d)\sum_{d|i}f\left(\frac i d \right)\)
把后面变一下
\(=\displaystyle \sum_{d=1}^ng(d)\sum_{i=1}^{\left\lfloor\frac n d \right \rfloor}f( i)\)
然后
\(=\displaystyle \sum_{i=1}^ng(i)S\left(\left\lfloor\frac n i\right\rfloor\right)\)
所以
\(\displaystyle \sum_{i=1}^nh(i)=\sum_{i=1}^ng(i)S\left(\left\lfloor\frac n i\right\rfloor\right)\)
有一个好像没用的式子
\(\displaystyle g(1)S(n)=\sum_{i=1}^ng(i)S\left(\left\lfloor\frac n i\right\rfloor\right)-\sum_{i=2}^ng(i)S\left(\left\lfloor\frac n i\right\rfloor\right)\)
上式把后面移项就成恒等式了
我们把右面第一项用刚刚的结论换走
\(\displaystyle g(1)S(n)=\sum_{i=1}^nh(i)-\sum_{i=2}^ng(i)S\left(\left\lfloor\frac n i\right\rfloor\right)\)
这。。是个递归式
就没了
对于S的递归,用数列分块
一般的h和g都很好求(构造)
对于本题来说
\(i*\varphi=id\)
所以对于\(\varphi\)
\(\displaystyle S(n)=\frac{n*(n+1)}{2}-\sum_{i=2}^nS\left(\left\lfloor\frac n i\right\rfloor\right)\)
刚刚有\(i*\mu=e\)
所以
\(\displaystyle S(n)=1-\sum_{i=2}^nS\left(\left\lfloor\frac n i\right\rfloor\right)\)
没了。。。
把前\(4*10^6\)的东西线性筛一下
最后的复杂度\(O(n^{\frac{2}{3}})\)不会证
#include <bits/stdc++.h>
typedef long long LL;
const int maxn = 4e6;
const int maxx = 4e6 + 10;
int in() {
char ch; int x = 0, f = 1;
while(!isdigit(ch = getchar()))(ch == '-') && (f = -f);
while(isdigit(ch)) x = (x << 1) + (x << 3) + (ch ^ 48), ch = getchar();
return x * f;
}
bool vis[maxx];
LL phi[maxx];
int mu[maxx], pri[maxx], tot;
std::map<int, LL> P;
std::map<int, int> M;
void predoit() {
phi[1] = mu[1] = 1LL;
for(int i = 2; i <= maxn; i++) {
if(!vis[i]) {
pri[++tot] = i;
phi[i] = i - 1;
mu[i] = -1;
}
for(int j = 1; j <= tot && i * pri[j] <= maxn; j++) {
vis[i * pri[j]] = true;
if(i % pri[j] == 0) {
phi[i * pri[j]] = phi[i] * pri[j];
mu[i * pri[j]] = 0;
break;
}
else {
phi[i * pri[j]] = phi[i] * (pri[j] - 1);
mu[i * pri[j]] = -mu[i];
}
}
}
for(int i = 2; i <= maxn; i++) {
phi[i] += phi[i - 1];
mu[i] += mu[i - 1];
}
}
LL workphi(int now)
{
if(now <= maxn) return phi[now];
if(P.count(now)) return P[now];
LL ans = now * (now + 1LL) / 2;
for(int i = 2, lst; i <= now; i = lst + 1) {
lst = now / (now / i);
ans -= 1LL * (lst - i + 1LL) * workphi(now / i);
}
return P[now] = ans;
}
int workmu(int now)
{
if(now <= maxn) return mu[now];
if(M.count(now)) return M[now];
int ans = 1;
for(int i = 2, lst; i <= now; i = lst + 1) {
lst = now / (now / i);
ans -= workmu(now / i) * (lst - i + 1);
}
return M[now] = ans;
}
int main() {
predoit();
for(int T = in(); T --> 0;) {
int n = in();
printf("%lld %d\n", workphi(n), workmu(n));
}
return 0;
}
P4213 【模板】杜教筛(Sum)的更多相关文章
- p4213 【模板】杜教筛(Sum)
传送门 分析 我们知道 $\varphi * 1 = id$ $\mu * 1 = e$ 杜教筛即可 代码 #include<iostream> #include<cstdio> ...
- [模板] 杜教筛 && bzoj3944-Sum
杜教筛 浅谈一类积性函数的前缀和 - skywalkert's space - CSDN博客 杜教筛可以在\(O(n^{\frac 23})\)的时间复杂度内利用卷积求出一些积性函数的前缀和. 算法 ...
- luoguP4213 [模板]杜教筛
https://www.luogu.org/problemnew/show/P4213 同 bzoj3944 考虑用杜教筛求出莫比乌斯函数前缀和,第二问随便过,第一问用莫比乌斯反演来做,中间的整除分块 ...
- 洛谷P4213(杜教筛)
#include <bits/stdc++.h> using namespace std; typedef long long LL; const int maxn = 3e6 + 3; ...
- LG4213 【模板】杜教筛(Sum)和 BZOJ4916 神犇和蒟蒻
P4213 [模板]杜教筛(Sum) 题目描述 给定一个正整数$N(N\le2^{31}-1)$ 求 $$ans_1=\sum_{i=1}^n\varphi(i)$$ $$ans_2=\sum_{i= ...
- 51NOD 1222 最小公倍数计数 [莫比乌斯反演 杜教筛]
1222 最小公倍数计数 题意:求有多少数对\((a,b):a<b\)满足\(lcm(a,b) \in [1, n]\) \(n \le 10^{11}\) 卡内存! 枚举\(gcd, \fra ...
- [洛谷P4213]【模板】杜教筛(Sum)
题目大意:给你$n$,求:$$\sum\limits_{i=1}^n\varphi(i),\sum\limits_{i=1}^n\mu(i)$$最多$10$组数据,$n\leqslant2^{31}- ...
- P4213【模板】杜教筛(Sum)
思路:杜教筛 提交:\(2\)次 错因:\(\varphi(i)\)的前缀和用\(int\)存的 题解: 对于一类筛积性函数前缀和的问题,杜教筛可以以低于线性的时间复杂度来解决问题. 先要构造\(h= ...
- BZOJ3944: Sum(杜教筛模板)
BZOJ3944: Sum(杜教筛模板) 题面描述 传送门 题目分析 求\(\sum_{i=1}^{n}\mu(i)\)和\(\sum_{i=1}^{n}\varphi(i)\) 数据范围线性不可做. ...
随机推荐
- JavaScript的运算符
一.什么是表达式??? 是ECMScript中的一个短语,解释器可以通过计算把它转成一个值,最简单的表达式是字面量或者变量名,单一字面量和组合字面量统称为表达式. 二.一元运算符 1.delete 运 ...
- NLTK词性标注解释
1. CC Coordinating conjunction 连接词2. CD Cardinal number 基数词3. DT Determin ...
- Struts2框架02 消息传递
1 消息传递 浏览器和服务器之间的数据传递 2 服务器项浏览器发送数据 2.1 在控制器类中增加一个属性,该属性的值就是服务器需要向浏览器发送的数据 2.2 为该属性增加 get 方法 2.3 在处理 ...
- CloudNote
目录 406错误 spring+mybatis建立连接失败 生成文档的乱码问题 服务器接收ajax发送的请求参数时出现乱码 获取不到cookie数据 1 利用Jackson返回json字符串时出现错误 ...
- C++ 成员函数前和函数后加const修饰符区别
博客转载自: https://www.iteblog.com/archives/214.html 分析以下一段程序,阐述成员函数后缀const 和 成员函数前const 的作用 #include< ...
- Linux cpulimit命令
一.简介 http://www.3lian.com/edu/2015/06-12/221261.html 二.安装 http://toutiao.com/a6184908783490073090/ 三 ...
- Android之对话框Dialog
首先是确认对话框 //确认对话框 private void showLog1() { AlertDialog.Builder dialog = new AlertDialog.Builder(this ...
- Python程序设计2——列表和元组
数据结构:更好的说法是从数据角度来说,结构化数据,就是说数据并不是随便摆放的,而是有一定结构的,这种特别的结构会带来某些算法上的性能优势,比如排序.查找等. 在Python中,最基本的数据结构是序列( ...
- Umbraco中获取UmbracoContext
在Umbraco项目中,获取当前的UmbracoContext几乎是都需要用到的,我们一般通过一个静态方法来获取,代码如下 public class ContextHelpers { public s ...
- Cookie seesion 赋值
后台赋值: Session["uid"] = userid;//后台 后台取值: Label1.Text = Session["userid"].ToStrin ...