【ACM数论】和式变换技术,也许是最好的讲解之一
在做数论题时,往往需要进行和式变换,然后变换成我们可以处理的和式,再针对和式做筛法、整除分块等操作。
本文将介绍一些常见的和式变换技术。
以下出现的概念大部分为个人总结,未必是学术界/竞赛界的统一说法,有不严谨的地方请谅解。
🎈 作者:Eriktse
🎈 简介:19岁,211计算机在读,现役ACM银牌选手🏆力争以通俗易懂的方式讲解算法!❤️欢迎关注我,一起交流C++/Python算法。(优质好文持续更新中……)🚀
🎈 原文链接(阅读原文获得更好阅读体验):https://www.eriktse.com/algorithm/1101.html
和式的基本形式
和式一般有两种:区间枚举型和整除枚举型。
区间枚举型
我们的前缀和就是一个典型的区间枚举型和式。
假设我们有一个定义域为
求和符号中,如果没有特殊说明,一般枚举的都是整数,且步长为1。
整除枚举型
约数个数是一个典型的整除枚举型和式,我们可以容易的写出它的表达式:
其中
表示 可以整除 ,即 是 的因子。
约数之和也是一个整除枚举型和式,表达式如下:
和式的基本性质
可拆分性质
第一种拆分如下:
这是显然的,但是基本上用不着。
第二种拆分如下:
这也是显然的。
常数可提取
当我们的和式里面乘上了一个常数
整除枚举型变换为区间枚举型(重要)
就比如上面那个约数之和的函数:
我们知道
我们称枚举的东西为“指标”,例如上面和式中
d|n
中的d
,i=1
中的i
。
指标变换(重要)
给定一个整数
现在令
所以我们可以不一个个枚举
然后变换枚举范围和布尔函数,注意这里
现在我们可以发现后面这个布尔函数就变成了一个常见的积性函数
交换求和次序(重要)
上式进行莫比乌斯反演后可以得到如下的和式(如果不懂莫比乌斯反演可以暂时先不管,之后再学),设
我们可以发现
接下来我们进行一次枚举类型的转换:
接下来我们将
转换为整除分块形式(十分重要)
上式转换完成后,我们可以发现后面两坨是可以进行整除分块的。
怎么理解呢?这个式子表达的就是当
那么和式就可转换为:
例题
luogu P2257 YY的GCD:https://www.luogu.com.cn/problem/P2257
阅读题意我们可以知道题目所求为,不妨设
接下来开始变换:
莫比乌斯反演:
注意这里
后面两坨可以进行整除分块,同时换一下
令
再交换求和次序:
现在发现
我们设一个函数:
那么
做完了。
Code:
#include <bits/stdc++.h>
#define int long long
using namespace std;
const int N = 1e7 + 9;
int sum[N], mu[N];
void init(int n = N - 2)
{
bitset<N> vis;
vector<int> prim;
vis[1] = true;
mu[1] = 1;
for(int i = 2;i <= n; ++ i)
{
if(!vis[i])prim.push_back(i), mu[i] = -1;
for(int j = 0;j < prim.size() and i * prim[j] <= n; ++ j)
{
vis[i * prim[j]] = true;
if(i % prim[j] == 0)break;//此时i * prim[j]含有平方因子
mu[i * prim[j]] = -mu[i];//此时i * prim[j]的本质不同质因子+1,或已经含有平方因子
}
}
for(int i = 0;i < prim.size(); ++ i)
{
for(int j = 1; prim[i] * j <= n; ++ j)
{
sum[prim[i] * j] += mu[j];
}
}
for(int i = 1;i <= n; ++ i)sum[i] += sum[i - 1];
}
void solve()
{
int n, m;scanf("%lld %lld", &n, &m);
if(n > m)swap(n, m);
int ans = 0;
for(int l = 1, r;l <= n; l = r + 1)
{
r = min(n / (n / l), m / (m / l));
ans += (sum[r] - sum[l - 1]) * (n / l) * (m / l);
}
printf("%lld\n", ans);
}
signed main()
{
init();
int _;scanf("%lld", &_);
while(_ --)solve();
return 0;
}
结束
🎈 本文由eriktse原创,创作不易,如果对您有帮助,欢迎小伙伴们点赞👍、收藏⭐、留言💬
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 阿里最新开源QwQ-32B,效果媲美deepseek-r1满血版,部署成本又又又降低了!
· 单线程的Redis速度为什么快?
· SQL Server 2025 AI相关能力初探
· AI编程工具终极对决:字节Trae VS Cursor,谁才是开发者新宠?
· 展开说说关于C#中ORM框架的用法!