博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
【文文殿下】【BZOJ4804】欧拉心算
阅读量:6992 次
发布时间:2019-06-27

本文共 739 字,大约阅读时间需要 2 分钟。

题解

显然有 \(ans=\sum _{i=1} ^{n} \lfloor \frac{n}{i} \rfloor \sum _{d|i} \mu(d) \phi (\frac{i}{d})\)

前半部分就是个整除分块,后半部分直接让脑子受到了冲击。

但是,我们知道,两个积性函数的\(\text{Direchlet}\)卷积还是积性函数,我们考虑构造线性筛。

我们看到其中有\(\text{Mobius}\)函数,鲁迅曾经说过“一看到\(\mu\),就想到积性函数,就想到唯一分解,就想到\(>2\)的幂的贡献不计,就想到线性筛,\(\text{OIer}\)只有在这个方面思想如此跃进。”

我们尝试对一个数\(T\)因数分解。

\(h(x)=\sum _{d|x} \mu(d) \phi (\frac{i}{d})\)对于其中一个素因子\(p\)和他的出现次数\(q\)\(x=p^q\)我们发现

\(h(x)=1,q=0\)

\(h(x)=p-2,q=1\)

\(h(x)=p^{q-2}(p-1)^2,q>=2\)

显然$h(T)=\prod _{i=1}^{t} h(p_i^{q_i}) $我们根据这个构造线性筛

#include
typedef long long ll;const int maxn = 1e7+10;ll h[maxn];int prime[maxn],cnt,n;bool vis[maxn];inline void prelude() { h[1]=1; for(int i = 2;i

转载于:https://www.cnblogs.com/Syameimaru/p/10233101.html

你可能感兴趣的文章
ORA-12514 TNS 监听程序当前无法识别连接描述符中请求服务 的解决方法
查看>>
PPP中的PAP和CHAP的区别
查看>>
基于CentOS5.5的SVN服务器搭建
查看>>
maven使用笔记
查看>>
JBoss配置使项目能在局域网其他机子上访问项目
查看>>
VIO概述 On-Manifold Preintegration for Real-Time Visual--Inertial Odometry
查看>>
CocoaPods升级安装三方库报错
查看>>
SpringBoot整合RabbitMQ实现微服务间的异步消息沟通
查看>>
pku1338 Ugly Numbers
查看>>
程序算法与人生选择 分类: 转载收藏 2013...
查看>>
牛客网校招全国统一模拟笔试(三月场)- Java方向
查看>>
Apache主站点配置
查看>>
[转]蓝牙开发
查看>>
C语言程序举例
查看>>
$.param()的实例应用
查看>>
web安全:xss && csrf
查看>>
数据保存(永久保存)方式
查看>>
POJ 3320 尺取法(基础题)
查看>>
如何使表格中的文字不换行?多出的字用“..."代替
查看>>
c# 进程间通信
查看>>