Codeforces Round 729（Div.2）C Strange Function 题解

题目大意

让$f(i)$表示最小的正整数$x$，满足$x$不是$i$的一个除数。
计算${\sum }_{i=1}^{n}f(i)$ $mod$ $10^9+7$。

输入描述

第一行包含一个整数$t(1\le t\le 10^4)$，表示测试用例的数量。
每一组测试用例包含一个整数$n(1\le n\le 10^{16})$。

输出描述

对每组测试用例，输出一个整数$ans$，满足$ans={\sum }_{i=1}^{n}f(i)$ $mod$ $10^9+7$。

样例输入

6
1
2
3
4
10
10000000000000000

样例输出

2
5
7
10
26
366580019

思路

先考虑$f(i)=x$的意义。由于$1,2,3,\dots ,x-1$都是$i$的因子，故$1,2,3,\dots ,x-1$的最小公倍数$LCM(1,2,3,\dots ,x-1)$也一定是$i$的因子，而$LCM(1,2,3,\dots ,x)$一定不是$i$的因子。我们得到结论：$i\ge LCM(1,2,3,\dots ,x-1)$。已知$i\le 10^{16}$，因此$LCM(1,2,3,\dots ,x-1)\le i\le 10^{16}$，这样$x$的取值上界得以确定。

在$lcm$函数中，涉及到两个数相乘的操作，我们判断一下会不会超过$long$ $long$。设$a$为累计的$lcm$，$b$为当前要求$lcm$的数，已经限定$a\le 10^{16}$，则$b\le 100$。$100$以内的质数有$25$个，不考虑$2,3,5,7$，有$21$个，只考虑质数则它们的$lcm$也超过了$10^{21}$，因此$b$不可能超过$100$。所以我们只需要枚举不到$100$个数就可以确定$x$的取值上界。实际上取值上界的值为$42$。

由于要计算${\sum }_{i=1}^{n}f(i)$ ，我们单独考虑取值上界下的每一个值$x$，并考虑在累加过程中不同的$x$出现的次数。在$1$~$n$中，能整除$LCM(1,2,3,\dots ,x-1)$的显然有$\lfloor \frac{n}{LCM(1,2,3,\dots ,x-1)}\rfloor$个，能整除$LCM(1,2,3,\dots ,x)$的显然有$\lfloor \frac{n}{LCM(1,2,3,\dots ,x)}\rfloor$个。由于$LCM(1,2,3,\dots ,x-1)$是$LCM(1,2,3,\dots ,x)$的因子，记集合$U$为{$1,2,3,\dots ,n$}，集合$A$为$U$中所有能整除$LCM(1,2,3,\dots ,x-1)$的数，集合$B$为$U$中所有能整除$LCM(1,2,3,\dots ,x)$的数，则$B\subseteq A$。在$Venn$图中我们要求$A\cap {\complement }_{U}B$，发现答案为${\complement }_{A}B$。则${\sum }_{i=1}^{n}f(i)={\sum }_{x=2}^{42}(x\times (\lfloor \frac{n}{LCM(1,2,3,\dots ,x-1)}\rfloor -\lfloor \frac{n}{LCM(1,2,3,\dots ,x)}\rfloor ))$。

考点

最小公倍数的性质
容斥原理的应用

AC代码

#include<iostream>
#define int long long
using namespace std;
const int mod = 1e9 + 7;
const int maxn = 105;
int T, n, ans, LCM[maxn];
inline int gcd(int a, int b) {
	return b == 0 ? a : gcd(b, a%b);
}
int lcm(int a, int b) {
	return a*b / gcd(a, b);
}
void init() {
	ios::sync_with_stdio(0);
	cin.tie(0), cout.tie(0);
	cin >> T;
}
int get_imax(int upp, int i) {
	for (; upp <= 1e16; i++)upp = lcm(upp, i);
	return i;
}
signed main()
{
	init();
	int imax = get_imax(1, 2);
	LCM[1] = 1;
	for (int i = 2; i <= imax; i++) {
		LCM[i] = lcm(LCM[i - 1], i);
	}
	while (T--) {
		cin >> n;
		ans = 0;
		for (int i = 2; n >= LCM[i - 1]; i++) {
			ans = (ans + i*((n / LCM[i - 1] - n / LCM[i]) % mod)) % mod;
		}
		cout << ans << endl;
	}
    return 0;
}

心得

做出这道题目，不仅要推导数学公式，还要考虑模域。由于题目限定了取值范围，因此求最小公倍数时乘积不会越界，不用取模。直到最后计算答案时再取模，这样就避免了除法取模，否则还要再去求逆元，会非常麻烦。