数论学习笔记-整除分块1

整除分块

主要就是一些笔记啦！参考的博客和题解：
参考题解
$$用途：快速处理形似\sum _{i=1}^n\lfloor \frac{n}{i}\rfloor 的式子$$

显然O(n)可以求，但是当数据到了10^10以上，我们就要寻找更优秀的时间复杂度的方法。

观察性质有：
$$\begin{gathered} 当\lfloor \frac{n}{i}\rfloor =\lfloor \frac{n}{i^{\prime }}\rfloor 时，i^{\prime }的最大取值为\lfloor \frac{n}{\lfloor \frac{n}{i}\rfloor }\rfloor .......(0) \\ 但是如何证明(0)式呢？下面给出(0)式的证明： \end{gathered}$$
(0)式证明如下：
$$\begin{gathered} 不妨设\lfloor \frac{n}{i}\rfloor =k,则i^{\prime }的最大取值为\lfloor \frac{n}{k}\rfloor \\ (目前是假设i^{\prime }的最大取值是\lfloor \frac{n}{k}\rfloor ) \end{gathered}$$

$$\begin{gathered} 而\lfloor \frac{n}{i}\rfloor 单调递减,所以 \\ 只要证明：i^{\prime }=k时满足\lfloor \frac{n}{\lfloor \frac{n}{k}\rfloor }\rfloor =k......(1) \\ 再证明：i^{\prime }=k时满足\lfloor \frac{n}{\lfloor \frac{n}{k}\rfloor +1}\rfloor <k.....(2) \\ 即可严格证明i=k时满足\lfloor \frac{n}{k}\rfloor 是满足条件最大的i \\ 先证明(1)式: \\ \because \lfloor \frac{n}{k}\rfloor <=\frac{n}{k} \\ \therefore \frac{n}{\lfloor \frac{n}{k}\rfloor }>=\frac{n}{\frac{n}{k}}=k \\ \therefore \lfloor \frac{n}{\lfloor \frac{n}{k}\rfloor }\rfloor >=k \\ (1)式得证.再证明(2)式： \\ 显然只要证得\frac{n}{\lfloor \frac{n}{k}\rfloor +1}<k则(2)式得证， \\ 只需证明：k\ast \lfloor \frac{n}{k}\rfloor +k>n \\ 由带余除法：n=q\ast k+r,0<=r<k \\ \therefore k\ast q+k>n \\ 而n=q\ast k+r,且0<=r<k \\ \therefore (2)式得证 \\ 综上所述：(0)式得证 \end{gathered}$$

代码：
例题1

typedef long long ll;
ll divide_block(ll n)
{
    ll res = 0;
    for(ll i = 1,j;i<=n;i=j+1)
    {
        j = n/(n/i);
        res += (j-i+1)*(n/i);
    }
    return res;
}

例题1

CQOI2007 余数之和
$$求\sum _{i=1}^{n}kmodi$$

$$\begin{gathered} 分析：乍一看mod和n/i并没有什么关系，但是可以展开有: \\ \sum _{i=1}^{n}k-\lfloor \frac{k}{i}\rfloor \ast i=n\ast k-\sum _{i=1}^n\lfloor \frac{k}{i}\rfloor \ast i \\ 而我们已经O(\sqrt{n})可以处理\sum _{i=1}^n\frac{n}{i} \\ 这里只需要考虑如何处理\sum _{i=1}^n\lfloor \frac{k}{i}\rfloor \ast i \\ 而对于一个整除分块\sum _{i=l}^r\lfloor \frac{k}{i}\rfloor ,其中的每个值都是相同的 \\ 那么不妨设T=\lfloor \frac{k}{i}\rfloor \\ 则原式=n\ast k-(\sum _{i=l}^r(T\cdot i)) \end{gathered}$$

代码：

typedef long long ll;
ll divide_block(ll n,ll k)
{
    ll ans = n*k;//根据公式来
    for(ll i = 1,j;i<=n;i=j+1)
    {
        if(k/i!=0)
        {
            j = min(k/(k/i),n);//i到j之间最多是n项，最少是k/(k/i)项
        }
        else
        {
            j = n;
        }
        ans -= (k/i)*(j-i+1)*(i+j)/2;//由推导的公式而来
        //等价于ans-=(k/i)*(j-i+1)+(i+j)/2*(j-i+1)
        //代码只是合并注释的同类项
    }
    return ans;
}

后续更新其他例题~