最小公倍数计数，51nod 1222，公式推导+Min_25筛

最新推荐文章于 2022-01-13 18:17:08 发布

Deep_Kevin

最新推荐文章于 2022-01-13 18:17:08 发布

阅读量184

点赞数

本文链接：https://blog.csdn.net/Deep_Kevin/article/details/100970036

版权

正题

公式还是比较好推的。

$x=p_1^{a_1}*p_2^{a_2}*...*p_k^{a_k}$

$f[x]=\frac{(\prod (2*a_i+1))+1}{2}$

为什么呢？考虑每一个质因子先将 $p_1^{a_1}$ 丢到其中的一边，然后另一边有 a_1+1 种选择，但 $(p_1^{a_1},p_1^{a_1})$ 这种情况被算了两遍，所以减掉。

乘起来，发现每一种二元组都被算了两遍，除了 (x,x) 只被算了一遍，所以+1再除以2就可以了。

然后发现 $g(x)=\prod (2a_i+1)$ 是一个积性函数，而且当 x=p^k 时函数值很好算。

g(x)=3|x=p ，所以直接上Min_25筛就可以了。

#include<bits/stdc++.h>
using namespace std;

const int N=320000;
int data=0;
long long a,b,w[N<<1],tot[N<<1];
long long n;
bool vis[N];
int p[N],id1[N],id2[N],t;

void prepare(int n){
	p[0]=0;memset(vis,false,sizeof(vis));
	for(int i=2;i<=n;i++){
		if(!vis[i]) p[++p[0]]=i;
		for(int j=1;j<=p[0] && 1ll*p[j]*i<=n;j++){
			vis[i*p[j]]=true;
			if(i%p[j]==0) break;
		}
	}
}

long long S(long long x,int t){
	int id;
	if(x>=data) id=id2[n/x];
	else id=id1[x];
	long long ans=tot[id]-3*(t-1);
	for(int i=t;i<=p[0] && 1ll*p[i]*p[i]<=x;i++)
		for(long long w=1ll*p[i]*p[i],wn=p[i],op=2;w<=x;op++)
			ans+=(2*op-1)*S(x/wn,i+1)+2*op+1,wn=w,w*=p[i];
	return ans;
}

long long solve(long long x){
	n=x;
	data=ceil(sqrt(n));
	prepare(data);
	long long l=1,r;t=0;
	while(l<=n){
		r=n/(n/l);
		if(n/l>=data) id2[r]=++t;
		else id1[n/l]=++t;
		w[t]=n/l,tot[t]=n/l-1,l=r+1;
	}
	int id;
	for(int j=1;j<=p[0];j++){
		for(int i=1;i<=t && 1ll*p[j]*p[j]<=w[i];i++){
			if(w[i]/p[j]>=data) id=id2[n/(w[i]/p[j])];
			else id=id1[w[i]/p[j]];
			tot[i]=tot[i]-(tot[id]-j+1);
		}
	}
	for(int i=1;i<=t;i++) tot[i]*=3;
	return S(n,1);
}

int main(){
	scanf("%lld %lld",&a,&b);
	printf("%lld\n",(solve(b)-solve(a-1)+b-a+1)/2);
}