个人复习（四）-除法分块

  学莫比乌斯反演学着学着，发现除法分块不会，学吧学吧，反正莫比乌斯反演都学了。

除法分块

      简单来说就是将商相同的数放在一起计算，可以把时间复杂度从O（n）降到O（sqrt(n)。最经典的就是∑(i从1到n)   n/i. 因为1-n的商会有相同的部分，可以把一起计算来加快速度。

 模板（后面有链接）

//经典模板：
long long solve(int N){
    int L= 1,R;
    long long tot = 0;
    while (L <= N){
        R = N / (N / L);  //求出右边的范围
        tot += (N / L) * (R - L + 1); //求出商相同区间的值
        L = R + 1; //区间往后移动1
    }
    return tot;
}

 

举个栗子：
求∑(i从1到10) n/i

i=1         10

i=2         5

i=3         3

i=4         4

i=5         5

i=6         1

i=7         1

i=8         1

i=9         1

i=10       1

我们可以看出i=6-10的结果都是1，我们可以将他们直接相加，求出他们的范围就好了，范围怎么求呢？
我们来看 R=N/(N/L)

先求出商的大小，然后再去找商为它的最大数是多少。

这样就可以求出它的范围（范围是连续的，所以我们每次记录左边界的范围）

上面就是最简单的除法分块情况

 

LG的2007题算是我写的第一题（转了点弯），需要先将题目个的公式化简成能用除法分块的形式：n*k-∑((k/i) * i) （i<=k）的形式  

相对的 分块的公式也有点不同：

（k/left）*(right-left)*(left+right)/2;

  因为每项后面还需要乘以i，所以我们需要将一个范围的i相加，所以还需要额外的乘以 (left+right)/2;

另外：注意的就是∑的范围，可能n<K 所以需要注意范围

AC代码，不是特别简洁：
 

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<queue>
#include<algorithm>
#define Max 1005
#define max(a,b) a>b?a:b;
#define min(a,b) a>b?b:a;

using namespace std;

{
	long long left,right,n; 
} 
int main()
{
	long long n,k;
	scanf("%lld%lld",&n,&k);
	long long ans,m;
	ans=n*k;
	int l=1,r;
	if(n>=k)
	while(l<=k)   //n>k的情况 
	{
		r=k/(k/l);
		ans-=(k/l)*(r-l+1)*(r+l)/2;
		l=r+1; 
	}
	else
	while(l<=n)  //n<k的情况 
	{
		r=k/(k/l);
		if(r>n)
		r=n;
		ans-=(k/l)*(r-l+1)*(l+r)/2;
		l=r+1;
	}
	
	printf("%lld\n",ans);
	
	return 0;
}

别人简洁的代码（别人的题解）：
 

#include<bits/stdc++.h>
using namespace std;
int main()
{
    long long n,k,l;
    scanf("%lld %lld",&n,&k);
    long long ans = n * k;
    for (long long i = 1 ; i <= n ; i = l + 1)
    {
        long long a = k/i;
        l = (a != 0) ? min(k/a,n) : n;
        long long cha = a * (i+l)*(l-i+1)/2;
        ans -= cha;
    }
    cout<<ans<<endl;
    return 0;
}

 

参考博客：

除法分块

除法分块证明