20200223 专题：数论分块（咕~~~）

总览：

解决一类求和问题${\sum }_{i=1}^n\lfloor \frac{n}{i}\rfloor$
发现一段$i$中$\lfloor \frac{n}{i}\rfloor$的值不变，于是可以分块求和
假设每一块的左右区间是$[l,r]$，则知道$l$后，$r=\frac{n}{n/l}$（均为整除）
证明
那么设两个指针$l,r$，每次$l=r+1,r=\frac{n}{n/l}$后再数量乘以值累加即可

模板：

for(int l=1,r=0;l<=n;l=r+1){
		if(n/l==0)	break;
		r=n/(n/l);//如果循环为k，写成min(k/(k/l),n);
		ans+=(n/l)*(r-l+1)/2;
	}

时间复杂度：$O(\sqrt{n})$

所以时间复杂度为$O(\sqrt{n})$

T1 P2261 [CQOI2007]余数求和

[CQOI2007]余数求和
题目描述
给出正整数 $n$ 和 $k$ 计算 $G(n,k)=k\mathrm{m}\mathrm{o}\mathrm{d}1+k\mathrm{m}\mathrm{o}\mathrm{d}2+k\mathrm{m}\mathrm{o}\mathrm{d}3+\cdots +k\mathrm{m}\mathrm{o}\mathrm{d}n$ 的值 其中 $k\mathrm{m}\mathrm{o}\mathrm{d}i$ 表示 $k$ 除以 $i$ 的余数。

例如 $G(10,5)=5\mathrm{m}\mathrm{o}\mathrm{d}1+5\mathrm{m}\mathrm{o}\mathrm{d}2+5\mathrm{m}\mathrm{o}\mathrm{d}3+5\mathrm{m}\mathrm{o}\mathrm{d}4+5\mathrm{m}\mathrm{o}\mathrm{d}5\cdots +5\mathrm{m}\mathrm{o}\mathrm{d}10=0+1+2+1+0+5+5+5+5+5=29$

输入格式
两个整数 $n$ ，$k$

输出格式
答案一行

输入输出样例
输入
10 5
输出
29

说明/提示
$100\%:n,k\le 10^9$

思路：
有$a\%b=a-a/b\times b$
则$G(n,k)=k\times k-\lfloor \frac{a}{b}\rfloor \times b$
数论分块，对于一个块$[l,r]$：
$\lfloor \frac{a}{b}\rfloor$为定值，$\sum b$为块中数的和，为$\frac{(l+r)\times (r-l+1)}{2}$

代码：

#include<bits/stdc++.h>
using namespace std;
#define in Read()
#define int long long

inline int in{
	int s=0,f=1;char x;
	for(x=getchar();x<'0'||x>'9';x=getchar())	if(x=='-')	f=-1;
	for( ;x>='0'&&x<='9';x=getchar())	s=(s<<1)+(s<<3)+(x&15);
	return f==1?s:-s;
}

int n,k;
int ans;

signed main(){ 
	n=in,k=in;
	for(int l=1,r=0;l<=n;l=r+1){
		if(k/l==0)	break;
		r=min(k/(k/l),n);
		ans+=(k/l)*(l+r)*(r-l+1)/2;
	}
	printf("%lld\n",n*k-ans);
	return 0;
}

咕