补题:
本题或许是参考 Problem - C - Codeforces
根据题意,f(i)
就是不能被整除的最小的一个质因子。
打表发现,当15个质因子相乘后,长度就大于18。
因此可以知道小于等于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,...k−1)r−lcm(1,...k)r⌋−⌊lcm(1,...k−1)l−lcm(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,...k−1)r−lcm(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
k≥2∑k×(⌊lcm(1,...k−1)r−lcm(1,...k)r⌋)=2×(r−⌊lcm(1,2)r⌋)+3×(⌊lcm(1,2)r−lcm(1,2,3)r⌋)+...=r+k≥1∑⌊lcm(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;
}
总结
赛时对这题容斥没有找到切入点。这个容斥处理极具有思维性,还是需要多做思维题!
参考: