新数论分块的推广及其复杂度证明 【数论】【分块】

事情是这样的。。。
最近校内集训有这样一道题：

给定一个N（N<1e12)，问有多少个M使得 (N%M)|M.

校内绝大部分的人（其实就是除了我外）都是用质因数分解过的，大概是：

因为gcd(n,m)=gcd(n%m,m),且gcd(n%m,m)==n%m (题意），所以就有gcd(n,m)=n%m.
枚举a=n%m，则要n=ax，m=ay 满足 gcd(x,y)=1.
又因为m|n-a (n-n%m),两边除以a，有y|x-1，此时显然有gcd(x,y)=1，所以上面的限制可以直接忽略。
那么再统计x-1的因数个数即可。

复杂度大约为N的所有约数的约数个数，有点绕...

这种做法太暴力且复杂度玄学，但这做法竟然跑得飞快。。。
于是我抱着一股愤怒的心情写下这篇博客，来推广一下新数论分块。

------------------------------------------------------我是萌萌哒的分割线--------------------------------------------------------

这其实是一道新数论分块（暂且叫它这个名字）的板子题。

先说结论：
类似于数论分块中，$\lfloor \frac{n}{i}\rfloor$的有有效取值只有$O(\sqrt{n})$个，$(\lfloor \frac{n}{i}\rfloor ,\lceil \frac{i}{n\%i}\rceil )$的有效取值只有$O(\sqrt{n}\ln n)$个。
($i|n$时单独划一组，显然这最多不超过$O(\sqrt{n})$组）
(上述的取整都可以分别改成上取整和下取整，这见仁见智）

推导过程如下：
$令a=\lfloor \frac{n}{i}\rfloor ,b=\lceil \frac{i}{n\%i}\rceil$，并假定存在一个常数S.
那么对于同一个a而言：
1.若$b\le S$时，则$(a,b)$最多有S组。
2.若$b>S$时：
首先有$\lfloor \frac{n}{i}\rfloor =a\iff \frac{n}{i}\ge a\iff i\le \frac{n}{a}$；
其次有$\lceil \frac{i}{n\%i}\rceil =\lceil \frac{i}{n-ai}\rceil >S\Rightarrow \lfloor \frac{i}{n-ai}\rfloor \ge S\iff \frac{i}{n-ai}\ge S\iff i\ge S\cdot (n-ai)\iff i\ge \frac{Sn}{aS+1}=\frac{n}{a+\frac{1}{S}}$
即$\frac{n}{a}\le i\le \frac{n}{a+\frac{1}{S}}$，显然$i$的整数取值最多只有$\frac{n}{a+\frac{1}{S}}-\frac{n}{a}=\frac{n}{a+a^{2}S}\le \frac{n}{a^{2}S}$

综合1和2，取$S=\frac{\sqrt{n}}{a}$，则有$S+\frac{n}{a^{2}S}=O(\frac{\sqrt{n}}{a})$

则不同的a值的总贡献大约为${\sum }_{a=1}^{n}O(\frac{\sqrt{n}}{a})=O({\sum }_{a=1}^n\frac{\sqrt{n}}{a})=O(\sqrt{n}\cdot {\sum }_{a=1}^n\frac{1}{a})=O(\sqrt{n}\ln n)$
而实际上程序跑出来的界也十分接近$\sqrt{n}\ln n$。

回到这题，其实就是问有多少个$M$使得$\frac{M}{N\%M}=\lceil \frac{M}{N\%M}\rceil$，直接套用上面的分块即可。
同时也可以证明答案不会超过$\sqrt{N}\ln N$，然而这给直接质因数分解提供了很好的复杂度证明.

-------------------------------------------------------我还是萌萌哒的分割线--------------------------------------------------------
具体写法：
先按数论整除分块，确定出$[i,j]，(\forall x\in [i,j],\lfloor \frac{n}{x}\rfloor =a(常数))$

观察到$\forall x\in [i,j],\lceil \frac{x}{n\%x}\rceil =\lceil \frac{x}{n-ax}\rceil$是一个不减的函数。
所以对于一个$l$（令$\lceil \frac{l}{n\%l}\rceil =\lceil \frac{l}{n-al}\rceil =b$），要确定一个最大$r$使得$[l,r]，\forall x\in [l,r],\lceil \frac{x}{n-ax}\rceil =b(常数)$,
即要找到最大的一个$r$使得$\lceil \frac{r}{n-ar}\rceil \le b$。

那么剩下就是简单地推一下公式了：
$\lceil \frac{r}{n-ar}\rceil \le b\iff \frac{r}{n-ar}\le b\iff r\le b\cdot (n-ar)\iff r\le \frac{bn}{ab+1}$
即$r=\lfloor \frac{bn}{ab+1}\rfloor$

注意$x|n$的情况要特判.

for (LL i=1,a,j;i<=n;i=j+1)	{
	j=n/(a=n/i);
	for (LL l=1,b,r;l<=j;l=r+1)	{
		if (n%l==0)	{
			assert(l==j);//当x|n时一定出现在整除分块的末尾处。
			//....
		} else	{
			b=l/(n%l);
			r=b*n/(a*b+1);
			//.....
		}
	}
}