【BZOJ3944】Sum(杜教筛)

题面

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

范围:\(n<2^{31}\)

令$$S(n)=\sum_{i=1}^n\mu(i)$$

随便找个函数\(g\)和\(\mu\)做狄利克雷卷积

\[(g*\mu)(i)=\sum_{d|i}\mu(d)g(\frac{i}{d})
\]

对这个玩意求前缀和

\[\sum_{i=1}^n\sum_{d|i}\mu(d)g(\frac{i}{d})
\]

把\(d\)给提出来

\[\sum_{d=1}^n\sum_{d|i}\mu(d)g(\frac{i}{d})
\]

\[\sum_{d=1}^n\sum_{i=1}^{n/d}\mu(d)g(i)
\]

把\(g(d)\)提出去

\[\sum_{d=1}^ng(d)\sum_{i=1}^{n/d}\mu(i)
\]

把前面设的东西带回去

\[\sum_{d=1}^ng(d)S(\frac{n}{d})
\]

那么,我们发现,我们要求的是\(S(n)\)

\[g(1)S(n)=\sum_{d=1}^ng(d)S(\frac{n}{d})-\sum_{d=2}^ng(d)S(\frac{n}{d})
\]

\[=\sum_{d=1}^n(g*\mu)(d)-\sum_{d=2}^ng(d)S(\frac{n}{d})
\]

我们知道:

\[\sum_{i|n}\mu(i)=[n=1]
\]

令\(g(x)=1\)

前面一堆东西是什么呀?

\[(1*\mu)(i)=\sum_{d|i}\mu(i)=[i=1]
\]

所以,

\[S(n)=1-\sum_{d=2}^nS(\frac{n}{d})
\]

后面的东西很显然可以数论分块

但是,我们不可能算完所有的东西

所以,线性筛预处理\(n^{\frac{2}{3}}\)项的前缀和

对于大于这个范围的值就做数论分块,然后递归处理


再看一下\(\phi\)怎么做

和前面一样,设出前缀和

\[S(n)=\sum_{i=1}^n\phi(i)
\]

找个函数\(g\)来做卷积

\[(g*\phi)(i)=\sum_{d|i}g(d)\phi(\frac{i}{d})
\]

求个前缀和

\[\sum_{i=1}^n(g*\phi)(i)=\sum_{i=1}^n\sum_{d|i}g(d)\phi(\frac{i}{d})
\]

\(g\)提出来

\[\sum_{d=1}^ng(d)\sum_{d|i}\phi(\frac{i}{d})
\]

\[=\sum_{d=1}^ng(d)\sum_{i=1}^{n/d}\phi(i)
\]

\[=\sum_{d=1}^ng(d)S(\frac{n}{d})
\]

根刚才的式子没什么区别

\[S(n)=\sum_{d=1}^ng(d)S(\frac{n}{d})-\sum_{d=2}^ng(d)S(\frac{n}{d})
\]

\[S(n)=\sum_{i=1}^n(g*\phi)(i)-\sum_{d=2}^ng(d)S(\frac{n}{d})
\]

我们又知道一个式子:

\[\sum_{d|i}\phi(d)=i
\]

也就是:$$(1*\phi)(i)=i$$

所以

\[S(n)=\sum_{i=1}^ni-\sum_{d=2}^nS(\frac{n}{d})
\]

\[S(n)=\frac{n(n+1)}{2}-\sum_{d=2}^nS(\frac{n}{d})
\]

和前面一样的,预处理\(n^\frac{2}{3}\)项

然后递归算就可以了


如果紧紧把范围卡在\(n^\frac{2}{3}\)这题会\(TLE\)???

一定是我常数太丑

反正可以多筛点

于是就多筛点

然后就\(AC\)啦???

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define MAX 2500000
#define ll long long
inline ll read()
{
ll x=0,t=1;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=-1,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return x*t;
}
bool zs[MAX+10];
int pri[MAX+10],tot;
ll mu[MAX+10],phi[MAX+10];
map<ll,ll> mm,pp;
void pre()
{
zs[1]=true;mu[1]=phi[1]=1;
for(int i=2;i<=MAX;++i)
{
if(!zs[i])pri[++tot]=i,mu[i]=-1,phi[i]=i-1;
for(int j=1;j<=tot&&i*pri[j]<=MAX;++j)
{
zs[i*pri[j]]=true;
if(i%pri[j])mu[i*pri[j]]=-mu[i],phi[i*pri[j]]=phi[i]*phi[pri[j]];
else{mu[i*pri[j]]=0;phi[i*pri[j]]=phi[i]*pri[j];break;}
}
}
for(int i=1;i<=MAX;++i)phi[i]+=phi[i-1];
for(int i=1;i<=MAX;++i)mu[i]+=mu[i-1];
}
ll Mu(ll x)
{
if(x<=MAX)return mu[x];
if(mm[x])return mm[x];
ll ret=0;
ll i=2,j;
while(i<=x)
{
j=x/(x/i);
ret+=(j-i+1)*Mu(x/i);
i=j+1;
}
return mm[x]=1-ret;
}
ll Phi(ll x)
{
if(x<=MAX)return phi[x];
if(pp[x])return pp[x];
ll ret=0;
ll i=2,j;
while(i<=x)
{
j=x/(x/i);
ret+=(j-i+1)*Phi(x/i);
i=j+1;
}
return pp[x]=x*(x+1)/2-ret;
}
int main()
{
pre();
int T=read();
while(T--)
{
ll n=read();
printf("%lld %lld\n",Phi(n),Mu(n));
}
return 0;
}

【BZOJ3944】Sum(杜教筛)的更多相关文章

  1. [BZOJ3944]Sum(杜教筛)

    3944: Sum Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 6201  Solved: 1606[Submit][Status][Discuss ...

  2. [bzoj3944] sum [杜教筛模板]

    题面: 传送门 就是让你求$ \varphi\left(i\right) $以及$ \mu\left(i\right) $的前缀和 思路: 就是杜教筛的模板 我们把套路公式拿出来: $ g\left( ...

  3. bzoj3944: Sum 杜教筛板子题

    板子题(卡常) 也可能是用map太慢了 /************************************************************** Problem: 3944 Us ...

  4. 3944: Sum[杜教筛]

    3944: Sum Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 3471  Solved: 946[Submit][Status][Discuss] ...

  5. 【Bzoj3944】杜教筛模板(狄利克雷卷积搞杜教筛)

    题目链接 哇杜教筛超炫的 有没有见过$O(n^\frac{2}{3})$求欧拉函数前缀和的算法?没有吧?蛤蛤蛤 首先我们来看狄利克雷卷积是什么 首先我们把定义域是整数,陪域是复数的函数叫做数论函数. ...

  6. 洛谷P4213 Sum(杜教筛)

    题目描述 给定一个正整数N(N\le2^{31}-1)N(N≤231−1) 求ans_1=\sum_{i=1}^n\phi(i),ans_2=\sum_{i=1}^n \mu(i)ans1​=∑i=1 ...

  7. bzoj 3944 Sum —— 杜教筛

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3944 杜教筛入门题! 看博客:https://www.cnblogs.com/zjp-sha ...

  8. BZOJ 3944: Sum [杜教筛]

    3944: Sum 贴模板 总结见学习笔记(现在还没写23333) #include <iostream> #include <cstdio> #include <cst ...

  9. BZOJ3944: Sum(杜教筛模板)

    BZOJ3944: Sum(杜教筛模板) 题面描述 传送门 题目分析 求\(\sum_{i=1}^{n}\mu(i)\)和\(\sum_{i=1}^{n}\varphi(i)\) 数据范围线性不可做. ...

随机推荐

  1. ASP.NET Core Razor 页面使用指南

    ASP.NET Core Razor 页面作为 ASP.NET Core 2.0的一部分发布,它是基于页面的全新的Web开发框架.如果您想学习如何使用 ASP.NET Core Razor 页面,可以 ...

  2. WPF ----在UserControl的xaml里绑定依赖属性

    场景:在定义wpf 用户控件的时候,希望使用时设置自定义的属性来改变用户控件里的状态或内容等. 下面直接上实例代码: 用户控件的后台代码,定义依赖属性 public partial class MyU ...

  3. 在Swift项目中使用OC,在OC项目中使用Swift

    几天前,我开始新的App的开发了.终于有机会把swift用在实战中了,也学到了之前纯学语法时没有机会获得的知识. 这篇博文中,我就如何使用swift.OC混编做一个介绍. OC中使用Swift 首先, ...

  4. 洛谷P2832 行路难 分析+题解代码【玄学最短路】

    洛谷P2832 行路难 分析+题解代码[玄学最短路] 题目背景: 小X来到了山区,领略山林之乐.在他乐以忘忧之时,他突然发现,开学迫在眉睫 题目描述: 山区有n座山.山之间有m条羊肠小道,每条连接两座 ...

  5. 华为云照片的爬虫程序更新(python3.6)

    一.背景: 每年终都有一个习惯,就是整理资料进行归档,结果发现手机照片全备份在华为云里,在官网上找了一圈,没找到官方的pc工具用来同步照片. 于是找出上次写的程序,看看能不能爬到数据,然而……果然不好 ...

  6. vue学习问题总结(一)

    使用comopontent组件报错错误信息:vue.js:491 [Vue warn]: Unknown custom element: <todo-item> - did you reg ...

  7. Java long类型和Long类型的那些事

    还记得最近做了一个项目使用的是Long类型作为主键Id坑死人了,对于我们来说Long类型一样是一个包装类型,类似String类型,使用==符号进行比较的时候有时候会出现问题,建议适应equal()方法 ...

  8. mysql2 - 基础

    一.SQL 练习 在java1701下,创建Stuednt 表,并插入以下数据: 增加创建时间字段,如下: 更改所有时间为当前时间: 二.数据库基础知识 1. 关系模型 1.1 表 table.列 c ...

  9. Redis限制在规定时间范围内登陆错误次数限制

    在博客里之前有过一篇文章是 <PHP结合Redis来限制用户或者IP某个时间段内访问的次数>,这篇文章的思路也是一样的.看下代码吧 //登录错误次数校验 $key = "logi ...

  10. numpy模块中的sum(axis)方法

    1.sum函数声明 sum(a, axis=None, dtype=None, out=None, keepdims=<class 'numpy._globals._NoValue'>) ...