莫比乌斯反演学习（我就不信我搞不懂它，夏令营之前将专心数学水平突破到一定层次！！！！！！！）

bzoj2301，终于对着模版A了，也算是做过了，感觉仍然是模模糊糊的，或许还需要小心的缕一缕

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define ll long long
using namespace std;
const int N=50009;
int mo[N],pri[N],sum[N];
bool bb[N];
ll a,b,c,d,k;
void getmo()//线性筛求莫比乌斯函数
{
	memset(bb,false,sizeof(bb));
	mo[1]=sum[1]=1;//注意mo【1】是1！！！！！！！！！不是-1。。。。。
	for (int i=2;i<N;i++)
	{
		if (!bb[i])
		{
			pri[++pri[0]]=i;
			mo[i]=-1;
		}
		for (int j=1;j<=pri[0]&&pri[j]*i<N;j++)
		{
			bb[pri[j]*i]=true;
			if (i%pri[j]==0)
			{
				mo[i*pri[j]]=0;break;
			}else mo[i*pri[j]]=-mo[i];
		}
		sum[i]+=sum[i-1]+mo[i];
	}
}
ll solve(ll n,ll m,ll k)
{
	if (n>m) swap(n,m);
	n/=k;m/=k;//求n，m里gcd（i，j）==k的个数就是求（n/k，m/k）中互质的数对个数
	ll ans=0;
	int r;
	for (int i=1;i<=n;i=r+1)
	{
		r=min(n/(n/i),m/(m/i));
		ans+=(sum[r]-sum[i-1])*(n/i)*(m/i);
	}
	return ans;
}
ll solve(ll a,ll b,ll c,ll d,ll k)
{
	return solve(b,d,k)-solve(a-1,d,k)-solve(b,c-1,k)+solve(a-1,c-1,k);
}
int main()
{
	getmo();
	int T;scanf("%d",&T);
	while (T--)
	{
		scanf("%lld%lld%lld%lld%lld",&a,&b,&c,&d,&k);
		printf("%lld\n",solve(a,b,c,d,k));
	}
	return 0;
}

2016.7.18更新

终于懂了这最简单的莫比乌斯反演，也算是理解了这最基本最经典的东西，、可能还是要做一些题虐智商来找找感觉

这次是理解着，自己写的代码

其实这个地方的经典是，求n，m；互质对的个数，如果不是求互质对，都一样，两边同除以k，继续求互质对即可

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define ll long long
using namespace std;
const int N=50009;
ll a,b,c,d,k,sum[N],u[N],pri[N];
bool bb[N];
void init()
{
	sum[1]=u[1]=1;
	for (int i=2;i<N;i++)
	{
		if (!bb[i])
		{
			u[i]=-1;
			pri[++pri[0]]=i;
		}
		for (int j=1;j<=pri[0]&&i*pri[j]<N;j++)
		{
			bb[i*pri[j]]=true;
			if (i%pri[j]==0)
			{
				u[i*pri[j]]=0;break;
			}
			else u[pri[j]*i]=-u[i];
		}
		sum[i]=sum[i-1]+u[i];
	}
}
ll work(ll n,ll m,ll k)
{
	if (n>m) swap(n,m);
	n=n/k;m=m/k;
	ll r,ans=0;
	for (int l=1;l<=n;l=r+1)//
	{
		r=min(n/(n/l),m/(m/l));
		ans+=(sum[r]-sum[l-1])*(n/l)*(m/l);
	}
	return ans;
}
ll work(ll a,ll b,ll c,ll d,ll k)
{
	return work(b,d,k)-work(a-1,d,k)-work(b,c-1,k)+work(a-1,c-1,k);
}
int main()
{
	int T;
	init();
	scanf("%d",&T);
	while (T--)
	{
 		scanf("%lld%lld%lld%lld%lld",&a,&b,&c,&d,&k);
		printf("%lld\n",work(a,b,c,d,k));
	}
	return 0;
}