洛谷 P3935 Calculating

虽然对这道题没有什么帮助，但是还是记一下：约数个数也是可以线性筛的

http://www.cnblogs.com/xzz_233/p/8365414.html

测正确性题目：https://www.luogu.org/problemnew/show/P1403

这个好像叫d函数
看$d=(a_1+1)(a_2+1)\cdots(a_k+1)$
然而还不行，你还要记这个数的$a_1$（定义在上面）记为f
首先，如果p是质数，那么d(p)=2,f(p)=1
然后，将合数n分解成n=px（p是n最小的质因子），
若$p\nmid x$则d(n)=2d(x),f(n)=1（d乘2相当于是要不要新选p）
否则$f(n)=f(x)+1$,$d(n)=d(x)*\frac{f(n)+1}{f(x)+1}$

---

https://www.luogu.org/problemnew/show/P3935

题目给的f(x)就是x的约数个数。。。

那么，$\sum_{i=1}^n(\sum_{d|n}1)=\sum_{i=1}^n({\lfloor}{\frac{n}{i}}{\rfloor})$

数论分块即可。。。

 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #include<vector>
 using namespace std;
 #define fi first
 #define se second
 #define mp make_pair
 #define pb push_back
 typedef long long ll;
 typedef unsigned long long ull;
 typedef pair<int,int> pii;
 #define md 998244353
 ll x,y,ans;
 int main()
 {
     ll i,j;
     scanf("%lld%lld",&x,&y);
     for(i=;i<=y;i=j+)
     {
         j=min(y,y/(y/i));
         ans+=(y/i)*(j-i+);
     }
     x--;
     for(i=;i<=x;i=j+)
     {
         j=min(x,x/(x/i));
         ans-=(x/i)*(j-i+);
     }
     printf("%lld",ans%md);
     return ;
 }