BZOJ 3994 [SDOI2015]约数个数和

Description

 设d(x)为x的约数个数，给定N、M，求  

Input

输入文件包含多组测试数据。

第一行，一个整数T，表示测试数据的组数。

接下来的T行，每行两个整数N、M。

Output

 T行，每行一个整数，表示你所求的答案。

Sample Input

2
7 4
5 6

Sample Output

110
121

HINT

 1<=N, M<=50000

1<=T<=50000

Source

Round 1 感谢yts1999上传

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

莫比乌斯反演~

推导见：http://blog.csdn.net/u012288458/article/details/51033779~

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

int t,n,m,a[N];
ll mu[N],f[N];
bool b[N];

int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0' || ch>'9') {if(ch=='-') f=-1;ch=getchar();}
	while(ch>='0' && ch<='9') {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
	return x*f;
}

ll cal(int n,int m)
{
	ll las,ans=0;
	if(n>m) swap(n,m);
	for(int i=1;i<=n;i=las+1)
	{
		las=min(n/(n/i),m/(m/i));
		ans+=f[n/i]*f[m/i]*(mu[las]-mu[i-1]);
	}
	return ans;
}

int main()
{
	mu[1]=1;
	for(int i=2;i<N;i++)
	{
		if(!b[i]) a[++a[0]]=i,mu[i]=-1;
		for(int j=1;j<=a[0] && i*a[j]<N;j++)
		{
			b[i*a[j]]=1;
			if(!(i%a[j]))
			{
				mu[i*a[j]]=0;break;
			}
			mu[i*a[j]]=-mu[i];
		}
	}
	for(int i=1;i<N;i++) mu[i]+=mu[i-1];
	for(int i=1;i<N;i++)
	{
		int las;
		for(int j=1;j<=i;j=las+1)
		{
			las=i/(i/j);
			f[i]+=(i/j)*(las-j+1);
		}
	}
	t=read();
	while(t--)
	{
		n=read();m=read();
		printf("%lld\n",cal(n,m));
	}
	return 0;
}