P4213 【模板】杜教筛

2023-07-23 23:03:12

题面:https://www.luogu.org/problem/P4213

#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;
const int mod=1000000007,N=5000005;
long long phi[N];
int t,prime[N],cnt,mobius[N];
bool vis[N];
unordered_mapans_mobius;
unordered_map ans_phi;
inline int read(){int x=0,f=1;char c=getchar();while(c<'0'||c>'9'){if(c=='-'){f=-1;c=getchar();}}while(c>='0'&&c<='9'){x=x*10+c-48;c=getchar();}return x*f;
}
inline long long get_phi(long long x){if(x<=N){return phi[x];}if(ans_phi[x]){return ans_phi[x];}long long ans=((1+x)*x)/2;for(register int l=2,r;l<=x;l=r+1){r=x/(x/l);ans-=(r-l+1)*get_phi(x/l);}return ans_phi[x]=ans;
}
inline int get_mobius(int x){if(x<=N){return mobius[x];}if(ans_mobius[x]){return ans_mobius[x];}int ans=1;for(register int l=2,r;l<=x;l=r+1){r=x/(x/l);ans-=(r-l+1)*get_mobius(x/l);}return ans_mobius[x]=ans;
}
int main(){mobius[1]=phi[1]=1;for(register int i=2;i<=N;++i){if(!vis[i]){prime[++cnt]=i;mobius[i]=-1;phi[i]=i-1;}for(register int j=1;j<=cnt&&prime[j]*i<=N;++j){vis[prime[j]*i]=1;if(i%prime[j]==0){phi[i*prime[j]]=phi[i]*prime[j];break;}phi[i*prime[j]]=phi[prime[j]]*phi[i];mobius[i*prime[j]]=-mobius[i];}}for(register int i=1;i<=N;++i){mobius[i]+=mobius[i-1];phi[i]+=phi[i-1];}int x;t=read();while(t--){x=read();printf("%lld %d\n",get_phi(x),get_mobius(x));}return 0;
}

转载于:https://www.cnblogs.com/ukcxrtjr/p/11311269.html

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！

标签：技术

上一篇 > P2324 骑士精神
下一篇 > MM的Windows 7 登录密码忘记之后

Duilib中list控件支持ctrl和shif多行选中的实现

[ICML2015]Batch Normalization:Accelerating Deep Network Training by Reducing Internal Covariate Shif

win10系统微软输入法于eclipse ctrl+shif+f冲突间接处理办法

Codeforces Round #259 (Div. 2) B. Little Pony and Sort by Shif

读LDD3，内存映射与DMA--PAGE_SHIF…

VMware虚拟机安装XP【要先分区，再设置BOOT 启动CD，shif+上移】

更换iBus五笔的左与右Shif

sublime ctrl+shif+f 没用解决办法

idea 对 ctrl + z 的撤销是 ctrl + shif + z

计算机最早的设计师应用于,计算机应用基础选择题doc.doc

win10自带截图神器：Win+Shift+S

Python基础之文件目录操作

python简述目录_Python基础之文件目录操作(示例代码)

tp5 如何做数据采集

任务2-7(服务器字体+阿里巴巴矢量库)

html标签（1)：h1~h6,p,br,pre,hr

TI 电量计介绍与芯片选型指南

几款TI电源芯片简介

TI DSP芯片C2000系列读取FLASH数据

德州仪器(Ti)平台嵌入式开发基础

TI三相电机智能栅极驱动芯片特点分类

省选模拟（12.08） T3 圈圈圈圈圈圈圈圈

Hadoop生态圈技术栈（上）

大数据开发基础入门与项目实战（三）Hadoop核心及生态圈技术栈之6.Impala交互式查询

小猿圈之Linux下Mysql 操作命令

大数据Hadoop生态圈常用面试题

大数据开发基础入门与项目实战（三）Hadoop核心及生态圈技术栈之4.Hive DDL、DQL和数据操作

备战Noip2018模拟赛11（B组）T3 Monogatari 物语

【智能优化算法-圆圈搜索算法】基于圆圈搜索算法Circle Search Algorithm求解单目标优化问题附matlab代码

NYOJ 78 圈水池

递归问题跑道汽车绕圈问题 Python实现

Hadoop生态圈（三）：MapReduce

P4213 【模板】杜教筛

相关文章