2017.6.4 problem b 失败总结

        只记住了反演的一个式子，还是最不常用的、、所以就推不动了

        首先再说明一下莫比乌斯函数：

              1、积性函数，但不是完全积性函数（f(x,y)=f(x)*f(y)   gcd(x,y)==1 ），因此可以线筛

              2、

             这个东西证明要用、、别的应该就没用了

                           3（求和函数的性质）、

                          

                           4、反演公式：

   ①：②：

               ps：背错公式害死人

                 反演，就是原来根据f（）求F（）函数变成了根据F（）求f（）。

                 所以F（）必须好求，f（）必须难求，，反演题的难度就在于找F（）和F（）和f（）之间的关系

                         把上面的式子琢磨透了就可以做题了

a~i~b c~j~d 中gcd(i,j)==k的数的个数    =   a/k~i~b/k  c/k~j~d/k 中gcd(i,j)==1的数的个数 

这个题用F（i）表示1~x~b，1~y~d的最大公倍数含i这个约数的个数   显然为下取整（b/i）*下取整（d/i）

f(k)表示a~i~b c~j~d 中gcd(i,j)==k的数的个数 

所以就有

答案在统计时注意要利用容斥：    + 1~b/k的和1~d/k    - 1~a/k的和1~d/k    - 1~c/k的和1~b/k   +1~a/k的和1~c/k

这样是o（n）的（枚举1的所有倍数）    但这样是过不了的、

有个优化就是利用    下取整（b/i）*下取整（d/i）的值的范围来用乘法和莫比乌斯函数前缀和优化

下一个数就是和i一个商的数的最大数+1，

就可以  q根n  出解了

码：

#include<iostream>
#include<cstdio>
using namespace std;
#include<algorithm>
int mu[50005],qsum[50005],i,su[50005],tot,j,T,a,b,c,d,k,ans;
bool he[50005];
void eular()
{
	mu[1]=1;
	qsum[1]=1;
	for(i=2;i<50000;i++)
	{
		if(!he[i])
		{
		    su[++tot]=i;
		    mu[i]=-1;
		}	
		for(j=1;j<=tot&&i*su[j]<50000;j++)
		{
			he[i*su[j]]=1;
			if(i%su[j]==0)mu[i*su[j]]=0;  
			else   	mu[i*su[j]]=-mu[i];
		}
		qsum[i]=qsum[i-1]+mu[i];
	}
	//for(i=1;i<=500;i++)cout<<qsum[i]<<" ";
	
}
int cal(int l,int r)
{
	int i,lin,daan=0;
	if(l>r)swap(l,r);
	for(i=1;i<=l;i=lin+1)
	{
	  lin=min(l/(l/i),r/(r/i));
	  	daan+=(qsum[lin]-qsum[i-1])*(l/i)*(r/i);
	}
	return daan;
}
int main()
{
	eular();
	scanf("%d",&T);
	while(T--)
	{
   scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
   ans=0;
   ans+=cal(b/k,d/k);	
   ans-=cal((a-1)/k,d/k);
   ans-=cal((c-1)/k,b/k);
   ans+=cal((a-1)/k,(c-1)/k);	
	printf("%d\n",ans);
    } 
 }