bzoj2301: [HAOI2011]Problem b懵逼乌斯反演

属于结果的和好求但是结果不好求的题

（轻易能得到以k的倍数为最大公约数的对数，但是不好直接求k）

所以一波反演结束

其实反演的时候完全没有反演的感觉，就是不停地恒等变形

算是懵逼乌斯反演最简单的例题

 #include <bits/stdc++.h>
 using namespace std;
 int n,m,a,b,c,d,k,mu[],p[];bool o[];
 int calc(int n,int m)
 {
     int ret=;if(n>m) swap(n,m);
     for(int i=,j;i<=n;i=j+)
     {
         j=min(n/(n/i),m/(m/i));
         ret+=(n/i)*(m/i)*(mu[j]-mu[i-]);
     }
     return ret;
 }
 int main()
 {
     scanf("%d",&n);
     m=;mu[]=;
     for(int i=;i<=;i++)
     {
         if(!o[i]) p[++m]=i,mu[i]=-;
         for(int j=;j<=m;j++)
         if(i*p[j]<=)
         {
             o[i*p[j]]=;
             if(i%p[j]==) break;
             mu[i*p[j]]=mu[i]*mu[p[j]];
         }
         else break;
     }
     for(int i=;i<=;i++)
         mu[i]+=mu[i-];
     for(int i=;i<=n;i++)
     {
         scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
         a=(a-)/k;b/=k;c=(c-)/k;d/=k;
         printf("%d\n",calc(b,d)-calc(b,c)-calc(a,d)+calc(a,c));
     }
     return ;
 } 

改天（老是拖延。。。）总结一下懵逼乌斯反演相关知识点