本文搬运自本人初中博客园博客,若图片加载不出来,可到原文查看:https://www.cnblogs.com/zhangtingxi/p/15780513.html

题目链接

题目

定义f(n)=ndn,dndf(n) = |n - \sum_{d|n,d\not =n}d|

每次给出A,BA,B,求i=ABf(i)\sum_{i=A}^B f(i)

对于100%100\%的数据,A,B107A,B\le 10^7

思路

看到 ABA,B 的范围,想着能不能快速求 f(n)f(n)

难点就在于求 dn,dnd\sum_{d|n,d\not =n}d,就是求一个数除自己的因子的和。

可以用类似筛的思想,一个 O(nlogn)O(n\log n) 解决。

总结

模拟赛时没打出这题,归根到底还是胆小了。

看到题目的数据范围,即使下面有个3秒,第一个反应还是 O(n)O(n),完全没打 O(nlogn)O(n\log n),很崩溃。

最后就只能乱打了个 O(nn)O(n\sqrt n) 交上去。

从中我也得到一个启示:有些时候 O(nlogn)O(n\log n) 能过 1e71e7 的数据。

就像NOIP2021T1一样。

orz HJH O(n)O(n) 做法。

Code

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
#include<bits/stdc++.h>
using namespace std;
#define int long long
inline int read(){int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;
ch=getchar();}while(ch>='0'&&ch<='9'){x=(x<<1)+
(x<<3)+(ch^48);ch=getchar();}return x*f;}
#define N 10000010
//#define M
//#define mo
int n, m, i, j, k;
int a, b, f[N];

signed main()
{
// freopen("tiaoshi.in","r",stdin);
// freopen("tiaoshi.out","w",stdout);
a=read(); b=read();
for(i=1; i<=b; ++i)
for(j=i+i; j<=b; j+=i) f[j]+=i;
for(i=a; i<=b; ++i) m+=abs(i-f[i]);
printf("%lld", m);
return 0;
}