数论基础之分块

数论分块

---

（一）提出问题

${\sum }_{i=1}^n\lfloor \frac{N}{i}\rfloor$
暴力？

---

（二）做法

1.显然可以使用$O（N）$的做法暴力过去，但是数据大了肯定就凉了
2.如果取N=10的话
可以发现$N/i$的值分别是10，5，3，2，2，1，1，1，1，1
这些相同的数字岂不是可以使用区间长度×相同的数字来直接计算吗，然后遍历过程中就能把变量i直接移到下一个位置（具体是哪里后续会解释）

---

（三）右端点位置的证明

现在给出一个左端点i，那么它的最大右端点$i_1$是多少，可以使得对于任意的$(i\le x\le i_1)$,都有$\lfloor \frac{N}{x}\rfloor =\lfloor \frac{N}{i}\rfloor$

事实上：$i_1=\lfloor \frac{N}{\lfloor \frac{N}{i}\rfloor }\rfloor$

证明（显然不是自己会证明的 ）：
设$\lfloor \frac{N}{i}\rfloor =k$,则根据定义可以写成$N=ki+r,(0\le r<i)$
设$\lfloor \frac{N}{i+d}\rfloor =k$,同理可得$N=ki+kd+r_1,(0\le r_1<i)$

显然我们只需求出$d_{max}=\lfloor \frac{r-r_1}{k}\rfloor =\lfloor \frac{r}{k}\rfloor$

$i_1=i+d_{max}$

$=i+\lfloor \frac{r}{k}\rfloor$

$=i+\lfloor \frac{N\%i}{\lfloor \frac{N}{i}\rfloor }\rfloor$

$=i+\lfloor \frac{N-\lfloor \frac{N}{i}\rfloor i}{\lfloor \frac{N}{i}\rfloor }\rfloor$

$=\lfloor i+\frac{N-\lfloor \frac{N}{i}\rfloor i}{\lfloor \frac{N}{i}\rfloor }\rfloor$

合并同类项后得?

$=\lfloor \frac{N}{\lfloor \frac{N}{i}\rfloor }\rfloor$

---

（四）复杂度分析

其实循环的次数就是$\lfloor \frac{N}{i}\rfloor$的不同的值的个数，即分块的个数

那么有多少个不同的数，对于一个N，可以分成两种情况考虑：
1.当$i\le \sqrt{N}$的时候，最多也只有i的范围长度即$\sqrt{n}$种可能性
2.当$i>\sqrt{N}$时，$\lfloor \frac{N}{i}\rfloor <\sqrt{n}$，最多只有范围长度$\sqrt{n}$种情况

综上，分块个数不会超过$2\sqrt{n}$个，即时间复杂度$O(\sqrt{n})$

---

（五）核心代码

ll ans=0;
for (int l=1,r=0;l<=n;l=r+1)
//每次左端点在前一个区间的右端点的下一个位置r+1
{
     r=n/(n/l);//右端点
     ans+=(r-l+1)*(n/l);
}

---

（六）例题

1.洛谷 P1403 板子题
题意：求出1-N的每个数的因数总和
做法：实际上并不用打表算，只要算${\sum }_{i=1}^n\lfloor \frac{N}{i}\rfloor$，即找出N个数中有多少个数字是含有因数i的，例如7/3=2，所以3和6是含有因数3的，则答案就加上2

int main()
{
    ll ans=0;
    ll n;
    scanf("%lld",&n);
    for (ll l=1,r=0;l<=n;l=r+1)
    {
        r=n/(n/l);
        ans+=(r-l+1)*(n/l);
    }
    WW(ans);
    return 0;
}

---

2.[CQOI2007]余数求和

题意：给出n,k，求${\sum }_{i=1}^n(kmodi)$
做法：考虑到模数并不能直接套，于是进行如下转化

${\sum }_{i=1}^n(kmodi)$

$={\sum }_{i=1}^n(k-\lfloor \frac{k}{i}\rfloor i)$

$=nk-{\sum }_{i=1}^{n}i\lfloor \frac{k}{i}\rfloor$

这就很好做了，将$\lfloor \frac{k}{i}\rfloor$的值相同的分为一块，然后$(l\le i\le r)$就相当于是一个✖上一个定值的等差数列了

	for(ll i=1,r;i<=n;i=r+1)
	{	
		if(k/i) r=min(k/(k/i),n);//不判的话会出现除0的错误
		else r=n;
		ans-=(k/i)*(r-i+1)*(i+r)/2;
	}	

---

尚不深究……