51nod 1222 最小公倍数计数

定义F(n)表示最小公倍数为n的二元组的数量。

即：如果存在两个数（二元组）X，Y(X <= Y)，它们的最小公倍数为N，则F(n)的计数加1。

例如：F(6) = 5，因为[2,3] [1,6] [2,6] [3,6] [6,6]的最小公倍数等于6。

给出一个区间[a,b]，求最小公倍数在这个区间的不同二元组的数量。

例如：a = 4，b = 6。符合条件的二元组包括：

[1,4] [2,4] [4,4] [1,5] [5,5] [2,3] [1,6] [2,6] [3,6] [6,6]，共10组不同的组合。

Input

输入数据包括2个数：a, b，中间用空格分隔（1 <= a <= b <= 10^11)。

Output

输出最小公倍数在这个区间的不同二元组的数量。

Input示例

4 6

Output示例

10

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

莫比乌斯反演~

推导见：http://blog.csdn.net/u014609452/article/details/54747826~

原来用的莫比乌斯函数线性筛法记岔了：mu[i*a[j]]=-mu[i]！

#include<cstdio>
#include<iostream>
using namespace std;
#define ll long long

ll l,r,ans,a[1000001],mu[10000001];
bool b[10000001];

ll cal(ll n)
{
	ll ans=0;
	for(ll a=1;a*a*a<=n;a++)
	  for(ll b=a;a*b*b<=n;b++)
	  {
	  	ll c=n/a/b;
	  	if(c<b) break;
	  	if(a==b) ans+=(c-b)*3+1;
	  	else ans+=(c-b)*6+3;
	  }
	return ans;
}

ll sol(ll n)
{
	if(!n) return 0;
	ll ans=0;
	for(ll i=1;i*i<=n;i++)
	  if(mu[i]) ans+=cal(n/i/i)*mu[i];
	ans=(ans+n)/2;
	return ans;
}

int main()
{
	mu[1]=1;
	for(int i=2;i<=10000000;i++)
	{
		if(!b[i]) a[++a[0]]=i,mu[i]=-1;
		for(int j=1;j<=a[0] && (ll)i*a[j]<=10000000;j++)
		{
			b[i*a[j]]=1;
			if(!(i%a[j]))
			{
				mu[i*a[j]]=0;break;
			}
			mu[i*a[j]]=mu[i]*mu[a[j]];
		}
	}
	scanf("%lld%lld",&l,&r);
	printf("%lld\n",sol(r)-sol(l-1));
	return 0;
}