题意

∑ k = 1 n ∑ i ∣ k ∑ j ∣ i λ ( i ) λ ( j ) \sum_{k = 1}^{n}\sum_{i \mid k} \sum_{j \mid i} \lambda(i) \lambda(j) k=1∑n​i∣k∑​j∣i∑​λ(i)λ(j)

对 998244353 998244353 998244353 取模

其中 λ ( x ) = ( − 1 ) ∑ i e i , x = ∏ i p i e i \lambda(x) = (-1)^{\sum\limits_{i}e_i},x=\prod\limits_{i}p_i^{e_i} λ(x)=(−1)i∑​ei​,x=i∏​piei​​

分析:

λ ( x ) \lambda(x) λ(x) 为刘维尔函数,可以打表发现 ∑ d ∣ n λ ( d ) = [ n = a 2 , a ∈ N + ] \sum_{d \mid n}\lambda(d) =[n = a^2,a \in N^+] d∣n∑​λ(d)=[n=a2,a∈N+]

也就是 n n n 是否为完全平方数

把式子中的 λ ( i ) \lambda(i) λ(i) 提到前面

∑ k = 1 n ∑ i ∣ k λ ( i ) ∑ j ∣ i λ ( j ) \sum_{k = 1}^{n}\sum_{i \mid k} \lambda(i)\sum_{j \mid i} \lambda(j) k=1∑n​i∣k∑​λ(i)j∣i∑​λ(j)

那么就变为

∑ k = 1 n ∑ i ∣ k λ ( i ) [ i = a 2 , a ∈ N + ] \sum_{k = 1}^{n}\sum_{i \mid k} \lambda(i)[i= a^2,a \in N^+] k=1∑n​i∣k∑​λ(i)[i=a2,a∈N+]

那么完全平方数的刘维尔函数为 1 1 1,再设 f ( x ) = [ i = a 2 , a ∈ N + ] f(x)=[i= a^2,a \in N^+] f(x)=[i=a2,a∈N+] 得

∑ i = 1 n ∑ d ∣ i f ( d ) \sum_{i = 1}^{n}\sum_{d \mid i}f(d) i=1∑n​d∣i∑​f(d)

交换求和次序

∑ d = 1 n f ( d ) ⌊ n d ⌋ \sum_{d = 1}^{n}f(d) \lfloor\frac{n}{d}\rfloor d=1∑n​f(d)⌊dn​⌋

这样直接枚举平方数即可,时间复杂度 O ( n ) O(\sqrt{n}) O(n ​)

代码:

#include <bits/stdc++.h>
#define int long long
using namespace std;
const int mod = 998244353;
int n, res;
signed main() {cin >> n;for (int i = 1; i * i <= n; i ++) {res = (res + n / (i * i)) % mod;}cout << res << endl;
}

[2021CCPC威海热身赛] Number Theory (打表 推公式)相关推荐

  1. hdu 4215 Number Theory?(打表)

    题目: E(N) = |{i | gcd(N, i) = 1, 1 <= i <= N}| F(N) = |{i | N % i = 0, 1 <= i <= N}| 求有多少 ...

  2. 2016级算法第二次上机-F.ModricWang's Number Theory II

    891 ModricWang's Number Theory II 思路 使得序列的最大公约数不为1,就是大于等于2,就是找到一个大于等于2的数,它能够整除序列中的所有数. 考虑使得一个数d整除数组中 ...

  3. concrete maths ch4 number theory

    ch4 number theory 数论研究正数的性质 1.整除 gcd lcm 扩展欧几里得. 整除求和\(\sum_{n|m}\)的几个公式.ch2的知识会很有用. 2.质数 Fundamenta ...

  4. Number Theory Problem(The 2016 ACM-ICPC Asia China-Final Contest 找规律)

    题目: Mr. Panda is one of the top specialists on number theory all over the world. Now Mr. Panda is in ...

  5. NUMTRYE - Number Theory (Easy)

    NUMTRYE - Number Theory (Easy) Hard 版本就是用 pollard_rho 分解质因子. f(n)=∏(pi2ei+1+1)f(n) = \prod(p_i ^{2e_ ...

  6. Codeforces Round #628 (Div. 2) E. Ehab‘s REAL Number Theory Problem 巧妙的质因子建图

    传送门 文章目录 题意: 思路: 题意: 给你nnn个数,每个数的因子个数不超过777个,选出最少的数使其乘积为平方数. n≤1e5n\le 1e5n≤1e5 思路: 由于因子不超过777个,所以由约 ...

  7. Segment 2:Introduction Number Theory——Fermat and Euler【费马定理和欧拉定理】

    Segment 2:Introduction Number Theory--Fermat and Euler[费马定理和欧拉定理] 这是整个数论简介的内容,下面是这个的主目录,其链接为:https:/ ...

  8. Elliptic Curves Number Theory And Cryptography——pairing learning

    <Elliptic Curves Number Theory And Cryptography 2n>中Example 11.5, magma脚本: clear; q:=7; Fq:=GF ...

  9. CodeForces - 466C Number of Ways(推公式/dp)

    题目链接:点击查看 题目大意:给出一个长度为 n 的数列,现在要求出满足条件的 ( i , j ) 的匹配数量,满足: 题目分析:训练时推的公式,简单说一下吧,维护前缀和 sum,则确定两个断点 ( ...

最新文章

  1. gis 数据框裁剪_BIM+GIS的八大挑战!大挑战,见未来
  2. 火币矿池参选EOS超级节点, 不只为分一杯羹
  3. 多分辨率适配常用目录
  4. TiDB 源码阅读系列文章(十五)Sort Merge Join
  5. 《卓有成效的程序员》----读书笔记二
  6. 什么是透传通道?(透明传输通道,就是当中继使)
  7. python僵尸进程和孤儿进程_python中多进程应用及僵尸进程、孤儿进程
  8. logo自动旋转 html,HTML/CSS3 Logo选择器
  9. python mysqldb 查询不到最新记录_python – MySQLdb是否缓存SELECT结果?
  10. 《linux c编程指南》学习手记2
  11. GoLang笔记—容器篇
  12. Unity资源加载闪退问题深度分析
  13. 各大云服务厂商 轻量应用服务器 性能评测对比,阿里云、腾讯云、华为云、Ucloud
  14. 刚上大一,C 语言压根听不懂怎么办?
  15. 数据类型,栈内存、堆内存
  16. 2023五一旅游必备物品清单!快记到手机待办APP里
  17. 2022计算机软件毕业设计大全
  18. html 设置请求头,http请求头和响应头设置
  19. 限时免费 | 学习成本超低的Python数据分析手册
  20. 英文Essay写作诀窍讲解

热门文章

  1. 申请idc许可证的条件、流程、需要的材料
  2. 数字冰雹 数字孪生城市智能运营中心(IOC)可视化决策系统
  3. antd+react ---手机号校验(流程中使用)
  4. Java物联网开发(一) —— MQTT协议
  5. 毕业设计-基于微信小程序的智能膳食分析系统
  6. 谷粒商城六商品服务三级分类
  7. 记一次win to go不能启动的问题
  8. 2021-12-11 WPF上位机 110-欧姆龙串口转网口的封装
  9. 条形码(barcode)code128生成代码
  10. 步进电机噪声和振动过大