数学知识——数论分块

对于一个式子：
$$\sum _{i=1}^n\lfloor \frac{n}{i}\rfloor$$
我们容易发现在某段区间$[l,r]$, $\lfloor \frac{n}{i}\rfloor$的值是一样的。那么我们如何获取这个区间呢。
我们假设区间左端点为$l$，同时令$k=\lfloor \frac{n}{l}\rfloor$，$R\ge l$，那么有$\frac{n}{R}\ge k$

要使得$R$最大：也就是区间右端点$r$应该为：$\lfloor \frac{n}{k}\rfloor$

最终结果应该为：$r=\lfloor \frac{n}{\lfloor \frac{n}{l}\rfloor }\rfloor$

话不多说直接看例题：

一：题目链接：The Fool

题意：

Sol:

其实打个表能发现，这个式子有规律 这是规律的AC代码

来看看分块怎么解决，我们发现$\left[\frac{n}{i}\right]$就是$\lfloor \frac{n}{i}\rfloor$
然后我们就可以利用分块算出这个值。
时间复杂度是：$O(\sqrt{N})$

Code:

	int T;
    cin >> T;
    int tt  = 1;
    while(T -- )
    {
        int n;
        cin >> n;
        int res = 0;
        for(int l = 1, r; l <= n; l = r+ 1)
        {
            if(n / l) r = min(n / (n / l), n); // 保险起见（n/l)不能为0
            else r = n;
            res += (n / l) * (r - l + 1);
        }
        if(res & 1) cout << "Case " << tt ++ << ": odd\n";
        else cout << "Case " << tt ++ << ": even\n";
    }

二.题目链接：Ice Rain

题意：

给出$n$和$k$求出下式的值:
$$\sum _{i=1}^n(k\mathrm{m}\mathrm{o}\mathrm{d}i)$$

Sol:

我们知道：$$k\mathrm{m}\mathrm{o}\mathrm{d}i=k-\lfloor \frac{k}{i}\rfloor \ast i$$
那么问题转化为：$$n\ast k-\sum _{i=1}^n(i\ast \lfloor \frac{k}{i}\rfloor )$$
对于后面这个式子：$$\sum _{i=1}^n(i\ast \lfloor \frac{k}{i}\rfloor )$$
我们通过分块能发现$\lfloor \frac{k}{i}\rfloor$在$[l,,r]$区间上值相同，那么这个区间的值就是：$(l+r)\ast (r-l+1)\ast \lfloor \frac{k}{i}\rfloor /2$
那么问题就迎刃而解了~~

Code：

int n, k;
    while(cin >> n >> k)
    {
        int res = n * k;
        for(int l = 1, r; l <= n; l = r + 1)
        {
            if(k / l ) r = min(k / (k / l), n);
            else r = n;
            res -= (l + r) * (r - l + 1) / 2 * (k / l);
        }
        cout << res << endl;
    }

三. 这是链接. 模积和

题意：

Sol：

这其实就是一个推式子的题目， 其次他给的模数不是素数！！！所以求逆元要用扩展欧几里得算法来求解

来推式子:
假设$n\le m$(大于的话交换就好了)
利用容斥可以得到：
$$\sum _{i=1}^n\sum _{j=1}^m(n\mathrm{m}\mathrm{o}\mathrm{d}i)\ast (m\mathrm{m}\mathrm{o}\mathrm{d}j)-\sum _{i=1}^n(n\mathrm{m}\mathrm{o}\mathrm{d}i)\ast (m\mathrm{m}\mathrm{o}\mathrm{d}i)$$
同样的我们知道：$$k\mathrm{m}\mathrm{o}\mathrm{d}i=k-\lfloor \frac{k}{i}\rfloor \ast i$$
所以上式可以写成：
$$\sum _{i=1}^n\sum _{j=1}^m(n-i\ast \lfloor \frac{n}{i}\rfloor )\ast (m-j\ast \lfloor \frac{m}{j}\rfloor )-\sum _{i=1}^n(nm-mi\ast \lfloor \frac{n}{i}\rfloor -ni\ast \lfloor \frac{m}{i}\rfloor +i^2\ast \lfloor \frac{n}{i}\rfloor \ast \lfloor \frac{m}{i}\rfloor )$$
稍微化简一下：
$$\sum _{i=1}^n(n-i\ast \lfloor \frac{n}{i}\rfloor )\ast \sum _{j=1}^m(m-j\ast \lfloor \frac{m}{j}\rfloor )-\sum _{i=1}^n(nm-mi\ast \lfloor \frac{n}{i}\rfloor -ni\ast \lfloor \frac{m}{j}\rfloor +i^2\ast \lfloor \frac{n}{i}\rfloor \ast \lfloor \frac{m}{i}\rfloor )$$
最后化简为：
$$(n^2-\sum _{i=1}^n(i\ast \lfloor \frac{n}{i}\rfloor )\ast (m^2-\sum _{j=1}^m(j\ast \lfloor \frac{m}{j}\rfloor )-\sum _{i=1}^n(nm-mi\ast \lfloor \frac{n}{i}\rfloor -ni\ast \lfloor \frac{m}{j}\rfloor +i^2\ast \lfloor \frac{n}{i}\rfloor \ast \lfloor \frac{m}{i}\rfloor )$$

好了，式子推完了，用分块来求可以实现$O(\sqrt{N})$。
其中：$$1^2+2^2+\cdot \cdot \cdot +n^2=\frac{n\ast (n+1)\ast (2n+1)}{6}$$
由于三个数相乘可能long long都会爆掉，所以要用逆元。另外，此题给的模数，不是一个素数。只有当模数是素数的时候，才可以用快速幂求出$a$的逆元$a^{mod-2}$。所以此题要用扩展欧几里得算法求逆元。
另外，注意取模。很容易爆掉（）

Code：

#include <bits/stdc++.h>
#define int long long
#define IOS ios::sync_with_stdio(false); cin.tie(0);cout.tie(0)
using namespace std;

const int mod = 19940417;
const int inv2 = 9970209 ,inv6 = 3323403;
typedef long long LL;

// 1^2 + 2^2 + ··· + x^2
LL sum(int x)
{
    return x * (x + 1) % mod * (2 * x + 1) % mod * inv6 % mod;
}

signed main(){
    IOS;
    int n, m;
    cin >> n >> m;
    if(n > m) swap(n, m);
    int n1, m1, nm;
    n1 = n * n %  mod;
    int l = 1, r;
    for( ; l <= n; l = r + 1)
    {
        r = n / (n / l);
        n1 = (n1 - ((r - l + 1) * (l + r) % mod * inv2 % mod * (n / l) % mod) + mod) % mod;
    }
    m1 = m * m % mod;
    l = 1;
    for( ; l <= m; l = r + 1)
    {
        r = m / (m / l);
        m1 = (m1 - ((r - l + 1) * (l + r) % mod * inv2 % mod * (m / l) % mod) + mod) % mod;
    }
    l = 1;
    nm = 0;
    for(; l <= n; l = r + 1)
    {
        r = min(n / (n / l), m / (m / l));
        int NM = (r - l + 1) * n % mod * m % mod;
        int MI = (r - l + 1) * (l + r) % mod * inv2 % mod * m % mod * (n / l) % mod;
        int NI = (r - l + 1) * (l + r) % mod * inv2 % mod * n % mod * (m / l) % mod;
        int I2MN = (n / l) * (m / l) % mod * (sum(r) - sum(l - 1) + mod) % mod;
        nm = ((((nm % mod + NM %mod) - MI + mod) % mod - NI + mod) % mod + I2MN + mod) % mod;
    }
    int ans = ((n1 * m1) % mod - nm + mod) % mod;
    cout << ans << endl;
    return 0;
}