51nod 莫比乌斯函数之和

题目在这里呀！

最近学了一下杜教筛，emm随便写个题解啦，把模板放上来而已ww

我觉得这篇博客讲解的很清晰qwq

简单讲一下，首先对于1到n的每个数i，d为i的约数，将mu[d]求和，和为1。

换句话说就是

for(int i=1;i<=n;i++) ans+=mu[d] (d|i)

ans的值为1，为什么呢，μ*1=(n==1 ? 1 : 0)也就是epsilon

然后，可以发现上式也等于

for(int i=1;i<=n;i++)
    for(int j=1;j<=n/i;j++) ans+=mu[j];

为什么呢，前者是枚举每个数，然后再枚举每个约数，后者是i枚举前者的i/d,j枚举d，两者等价。

所以后者也为1.

进一步地，

for(int i=2;i<=n;i++)
    for(int j=1;j<=n/i;j++) ans+=mu[j];

这里的ans就是1-sum(mu[i])(i=1~n)

通过移项可得，

∑i=1~n μ(i)=1−∑i=2~n ∑j=1~⌊n/i⌋ μ(j)

话说这里的格式好难看啊qwq

这样我们就可以hash+递归求解啦！

最近好像状态不太好，被数组越界搞了好长时间w。

（emm幸亏有dalao指出ww

//Suplex
#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#define N 5000000+5000
#define ll long long
using namespace std;
const int mod=1000007;
int cnt,tot,prime[N],head[mod+10],Next[N];
bool vis[N];
ll a,b,tmp,mu[N],t[500000],hashx[500000];

inline void insert(ll k,ll val){
	int key=k % mod;
	for(int i=head[key];i;i=Next[i])
		if(hashx[i]==k) return;
	++cnt;Next[cnt]=head[key];head[key]=cnt;
	hashx[cnt]=k;t[cnt]=val;
}

inline ll find(ll k){
	int key=k % mod;
	for(int i=head[key];i;i=Next[i])
		if(hashx[i]==k) return t[i];
	return -1;
}

void pre(){
	mu[1]=1;
	for(int i=2;i<N;i++){
		if(!vis[i]){
			prime[++tot]=i;
			mu[i]=-1;
		}
		for(int j=1;j<=tot && (ll)i*prime[j]<N;j++){
			vis[i*prime[j]]=1;
			if(i % prime[j]==0){
				mu[i*prime[j]]=0;
				break;
			}else mu[i*prime[j]]=-mu[i];
		}
	}
	for(int i=2;i<N;i++) mu[i]+=mu[i-1];
}

ll solve(ll n){
	if(n<=N) return mu[n];
	ll k=find(n);
	if(k != -1) return k;
	ll next,ans=1;
	for(ll i=2;i<=n;i=next+1){
		ll now=n/i;
		next=n/now;
		ans-=(ll)(next-i+1)*solve(n/i);
	}
	insert(n,ans);
	return ans;
}

int main()
{
	scanf("%lld%lld",&a,&b);
	pre();
	printf("%lld\n",solve(b)-solve(a-1));
	return 0;
}