[模板][P3377]杜教筛

Description:

求 $ \sum_{i=1}^n \phi(i) ,\sum_{i=1}^n \mu(i)$

Hint:

\(n<=10^{10}?\)

Solution:

考虑积性函数 \(f,g,h?\) 及其前缀和 \(F,G,H?\)

其中 \(h=f*g?\)

首先 \(H(x)=\sum_{n=1}^xh(n)\)

\(=\sum_{n=1}^x \sum_{d|n} f(d) g(\frac{n}{d})\)

枚举倍数转枚举因数

\(=\sum_{k=1}^x \sum_{d=1}^{\lfloor \frac{x}{k} \rfloor} f(d) g(k)?\)

$=\sum_{k=1}^x g(k)\sum_{d=1}^{\lfloor \frac{x}{k} \rfloor} f(d) $

\(=\sum_{k=1}^x g(k) F(\lfloor \frac{x}{k} \rfloor)\)

\(=\sum_{k=1}^{x}g(d)F(\lfloor \frac{x}{k} \rfloor)\)

故 \(g(1)F(n)=H(n)-\sum_{d=2}^{n}g(d)F(\lfloor \frac{n}{d} \rfloor)\)

此式是杜教筛的核心式,适用于非线性求一个积性函数的前缀和

只要能快速求出 \(H(n),\sum g(d)\) 就能在$ O(n^{ \frac{2}{3}} ) $ 求出\(F(n)\)

\(\sum_{i=1}^n \phi(i)=\frac{n(n+1)}{2}-\sum_{i=2}^n \phi(\lfloor \frac{n}{i} \rfloor)\)

\(\sum_{i=1}^n \mu(i)=1-\sum_{i=2}^n \mu(\lfloor \frac{n}{i}\rfloor)\)

先筛出线性数据范围内的,再杜教筛

递归求解即可

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int mxn=5e6+5,inf=2147483647;
int T,n,tot;
int p[mxn],vis[mxn];
ll mu[mxn],ph[mxn];
map<int ,ll > smu,sph;

void init()
{
    vis[1]=mu[1]=ph[1]=1;
    for(int i=2;i<=mxn;++i) {
        if(!vis[i]) mu[i]=-1,ph[i]=i-1,p[++tot]=i;
        for(int j=1;j<=tot&&p[j]*i<=mxn;++j) {
            vis[p[j]*i]=1;
            if(i%p[j]) mu[p[j]*i]=-mu[i],ph[p[j]*i]=ph[p[j]]*ph[i];
            else {ph[p[j]*i]=ph[i]*p[j];break;}
        }
    }
    for(int i=2;i<=mxn;++i) mu[i]+=mu[i-1],ph[i]+=ph[i-1];
}

ll get_mu(int n)
{
    if(n<=mxn) return mu[n];
    if(smu[n]) return smu[n]; ll ans=0;
    for(int l=2,r;r<inf&&l<=n;l=r+1)
        r=n/(n/l),ans+=(r-l+1)*get_mu(n/l);
    return smu[n]=1ll-ans;
}

ll get_ph(int n)
{
    if(n<=mxn) return ph[n];
    if(sph[n]) return sph[n]; ll ans=0;
    for(int l=2,r;r<inf&&l<=n;l=r+1)  //一定是从2开始
        r=n/(n/l),ans+=(r-l+1)*get_ph(n/l);
    return sph[n]=(ull)n*(n+1ll)/2-ans;
}

int main()
{
    cin>>T; init();
    while(T--) {
        scanf("%d",&n);
        printf("%lld %lld\n",get_ph(n),get_mu(n));
    }
    return 0;
}

原文地址:https://www.cnblogs.com/list1/p/10395182.html

时间: 2024-11-09 00:34:12

[模板][P3377]杜教筛的相关文章

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

【模板】杜教筛(Sum)

传送门 Description 给定一个正整数\(N(N\le2^{31}-1)\) 求 \[ans1=\sum_{i=1}^n \varphi(i)\] \[ans_2=\sum_{i=1}^n \mu(i)\] Solution 总算是写了一个不会\(TLE\)的杜教筛,不想用\(map\),因此上了一个很丑的\(Hash\)-- Code #include<bits/stdc++.h> #define ll long long #define max(a,b) ((a)>(b)?(

P4213【模板】杜教筛(Sum)

思路:杜教筛 提交:\(2\)次 错因:\(\varphi(i)\)的前缀和用\(int\)存的 题解: 对于一类筛积性函数前缀和的问题,杜教筛可以以低于线性的时间复杂度来解决问题. 先要构造\(h=f*g\),并且\(h\)的前缀和易求,\(g\)的区间和易求. 具体地: \[\sum_{i=1}^{n}h(i)=\sum_{i=1}^{n}\sum_{d|i}g(d)\cdot f(\frac{i}{d})\] \[\sum_{i=1}^{n}h(i)=\sum_{d=1}^{n}g(d)\

P4213 【模板】杜教筛(Sum) min_25筛

\(\color{#0066ff}{ 题目描述 }\) 给定一个正整数\(N(N\le2^{31}-1)\) 求 \(ans_1=\sum_{i=1}^n\varphi(i)\) \(ans_2=\sum_{i=1}^n \mu(i)\) \(\color{#0066ff}{输入格式}\) 一共T+1行 第1行为数据组数T(T<=10) 第2~T+1行每行一个非负整数N,代表一组询问 \(\color{#0066ff}{输出格式}\) 一共T行,每行两个用空格分隔的数ans1,ans2 \(\c

P4213 【模板】杜教筛

题面:https://www.luogu.org/problem/P4213 #include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> #include<cmath> #include<algorithm> #include<ctime> #include<queue> #include<map> using na

P4213【模板】杜教筛

题面 https://www.luogu.org/problem/P4213 题解 #include<cstdio> #include<iostream> #include<cstring> #pragma GCC optimize(2) #define ri register int #define maxn 2147483647 #define maxm 3000000 #define maxk 2000 #define LL long long using nam

CCPC 2019 网络赛 HDU huntian oy (杜教筛)

1005 huntian oy (HDU 6706) 题意: 令,有T次询问,求 f(n, a, b). 其中 T = 10^4,1 <= n,a,b <= 1e9,保证每次 a,b互质. 思路: 首先我们需要知道 公式: gcd(a^n - b^n, a^m - b^m) = a^(gcd(m,n)) - b^(gcd(m,n)) 由a,b互质,原式即为 f(n, a, b) = ∑∑ (i-j)*[(i,j)=1] = ∑ (i*∑ [(i, j)=1] ) - ∑∑ j*[(i, j)=

【数论】狄利克雷卷积及其快速计算方法及杜教筛

目录(假的 狄利克雷卷积基础知识 数论函数 狄利克雷卷积定义 狄利克雷卷积性质 常用卷积 卷积计算方法 最暴力的暴力 稍好的暴力 优美的暴力 莫比乌斯反演(待填坑) 杜教筛 经典杜教筛 第二种杜教筛 第三种杜教筛 背景 本人即将去CTS&APIO2019,由于一些特殊原因,发现自己数论突然变得很菜. 就决定在去的前一天,翻出来以前的数论学习资料看一看.翻到了czgj的校内狄利克雷卷积课件,发现其中提到了的任意数列\(f(n)\)和\(g(n)\)的狄利克雷卷积\((f*g)(n)\)(从1到n,

luoguP4213 [模板]杜教筛

https://www.luogu.org/problemnew/show/P4213 同 bzoj3944 考虑用杜教筛求出莫比乌斯函数前缀和,第二问随便过,第一问用莫比乌斯反演来做,中间的整除分块里的莫比乌斯前缀和刚好用第二问来做 杜教筛的时候先线性筛出前 N 个数的莫比乌斯函数前缀和,其余的用 map 记忆化搜索,实测 N 取 3670000 最佳(其实我只测了3次) #include <bits/stdc++.h> using namespace std; typedef unsign