当前位置: 首页 > news >正文

容斥 C. Strange Function改编题

补题:

题目详情 - 9.段坤爱取模%%% - SUSTOJ

本题或许是参考 Problem - C - Codeforces

根据题意,f(i)就是不能被整除的最小的一个质因子。

打表发现,当15个质因子相乘后,长度就大于18。

image-20231119231908296

因此可以知道小于等于1e16内的正整数x,f(x)一定是前20个质因子之一,且合数一定不行。

前20个质因子:2 3 5 7 11 13 17 19 23 29 31 37 41 43 。在第15个质因子相乘时就大于1e16,因此max(f(i))是小于40的。

现在就是要求:第1个质因子在[l, r]的个数 乘上 第1个质因子,加上 第2个质因子在[l, r]的个数 乘上 第2个质因子个数, … 。需要保证i在[l, r]内仅且一次被计算。考虑容斥。

对于f(i) = k来说,i可以被1 ... k - 1任何一个整除,但是不能被k整除。

f()值为k的个数有:
⌊ r l c m ( 1 , . . . k − 1 ) − r l c m ( 1 , . . . k ) ⌋ − ⌊ l l c m ( 1 , . . . k − 1 ) − l l c m ( 1 , . . . k ) ⌋ \lfloor { \frac {r} {lcm(1, ... k-1)} } - { \frac {r} {lcm(1, ... k)} } \rfloor - \lfloor { \frac {l} {lcm(1, ... k-1)} } - { \frac {l} {lcm(1, ... k)} } \rfloor lcm(1,...k1)rlcm(1,...k)rlcm(1,...k1)llcm(1,...k)l
如果将[l, r]分成两部分[1, l][1, r]来处理的话。

以r为例,f()为k的个数:
⌊ r l c m ( 1 , . . . k − 1 ) − r l c m ( 1 , . . . k ) ⌋ \lfloor { \frac {r} {lcm(1, ... k-1)} } - { \frac {r} {lcm(1, ... k)} } \rfloor lcm(1,...k1)rlcm(1,...k)r
那么,[1, r]sum就是:
∑ k ≥ 2 k × ( ⌊ r l c m ( 1 , . . . k − 1 ) − r l c m ( 1 , . . . k ) ⌋ ) = 2 × ( r − ⌊ r l c m ( 1 , 2 ) ⌋ ) + 3 × ( ⌊ r l c m ( 1 , 2 ) − r l c m ( 1 , 2 , 3 ) ⌋ ) + . . . = r + ∑ k ≥ 1 ⌊ r l c m ( 1 , . . . , k ) ⌋ \sum_{k \ge 2}k \times(\lfloor { \frac {r} {lcm(1, ... k-1)} } - { \frac {r} {lcm(1, ... k)} } \rfloor) \\ = 2 \times (r - \lfloor \frac r {lcm(1,2)} \rfloor) + 3 \times (\lfloor { \frac {r} {lcm(1,2)} } - { \frac {r} {lcm(1,2,3)} } \rfloor) + ... \\ = r + \sum_{k \ge 1} \lfloor \frac r {lcm(1, ..., k)} \rfloor k2k×(⌊lcm(1,...k1)rlcm(1,...k)r⌋)=2×(rlcm(1,2)r⌋)+3×(⌊lcm(1,2)rlcm(1,2,3)r⌋)+...=r+k1lcm(1,...,k)r
同理,可得[1, l]内的sum,两者相减即为答案。对于cf上的C,[1, r]即可。

时间复杂度:前面稠密,后面稀疏,最大为O(41 * 2)

const int mod = 998244353;
// https://codeforces.com/contest/1542/problem/C
void solve() {int l,r; cin>>l>>r;int sum = 0;int v = 1;auto lcm = [&](int a, int b) {return a * b / __gcd(a,b);};for(int i = 1; v <= r; v = lcm(i, v), ++i) {sum = (sum + r / v) % mod;}v = 1;for(int i = 1; v <= l - 1; v = lcm(i,v), ++i) {sum = (sum - (l - 1) / v + mod) % mod;}cout<<sum<<endl;
}

总结

赛时对这题容斥没有找到切入点。这个容斥处理极具有思维性,还是需要多做思维题!

参考:

[Codeforces] number theory (R1600) Part.11 - 知乎 (zhihu.com)

Submission #121200660 - Codeforces

http://www.lryc.cn/news/240571.html

相关文章:

  • C++笔记
  • python-opencv 培训课程笔记(1)
  • 【C++初阶】STL详解(七)Stack与Queue的模拟实现
  • 校园报修抢修小程序系统开发 物业小区报修预约上门维修工单系统
  • 【Android】Hilt比Android好在哪里
  • 计算方法 期末总结
  • 【面试】jvm中堆是分配对象存储的唯一选择吗
  • 音视频同步笔记 - 以音频时间为基
  • JavaScript 原始数据类型和对应的对象类型(内置对象)之间的关系
  • 报错For debugging consider passing CUDA_LAUNCH_BLOCKING=1.
  • whisper使用方法
  • 通过easyexcel实现数据导入功能
  • Springboot_文件下载功能(前端后端)
  • Vue框架学习笔记——v-bind数据单向绑定和v-model数据双向绑定
  • 将对象转成URL参数
  • 【论文阅读】MAG:一种用于航天器遥测数据中有效异常检测的新方法
  • 超级武器!深入LoadRunner性能测试流程及极速分析结果!
  • 解决requests库进行爬虫ip请求时遇到的错误的方法
  • 大语言模型领域的重要术语解释
  • 无需公网IP,使用内网穿透实现公网访问本地OpenWRT管理界面
  • 利用 docker 实现JMeter分布式压测
  • Sentinel 分布式系统
  • Cookie、Session、CBV加装饰器的三种方法
  • 【高级渗透篇】网络安全面试
  • 【Flink】Process Function
  • 【源码分析】zeebe actor模型源码解读
  • python3实现类似expect shell的交互式与SFTP的脚本
  • java游戏制作-飞翔的鸟游戏
  • NodeMCU ESP8266构建Web Server网页端控制设备
  • 搭建区块链