HYSZB PROBLEM B（mobius反演+分块+容斥原理）

这道题和hdu1965思路基本类似，
http://blog.csdn.net/abc13068938939/article/details/52198163
mobius反演得出的结果是(过程在上面的博文中有介绍)
1< =i< =m和1< =j< =n中gcd（i,j）=k的个数
令m~=m/k, n~=n/k。且不妨设m < = n。

f(1)=Sigma(d=1 : m) u(d)* ((m~)/d) *((n~)/d)。

m/=k,n/=k;
if(m>n) swap(n,m);
LL res=0ll;
for(int i=1;i<=m;i++)
    res+=u(i)*(m/i)*(n/i);

然后a< =i< =b和c< =j< =d中gcd（i,j）=k的个数=
solve(b/k,d/k)-solve(b/k,(c-1)/k)-solve((a-1)/k,d/k)+solve((a-1)/k,(c-1)/k)。
其中solve(n,m)表示1< =i< =m和1< =j< =n中gcd（i,j）=1的个数。
但是这是会超时的，
事实上m/i它在一小块内是不会变的，比如
m=100 i=26:33，m/i均等于3，也就是这一部分可以合着计算。
分块不熟可以看这道题bzoj1257(余数之和)

#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#define msc(X) memset(X,-1,sizeof(X))
#define ms(X) memset(X,0,sizeof(X))
typedef long long LL;
using namespace std;
const int MAXN=50005;
int sum[MAXN+10],mu[MAXN+10];
void getMobius(void)
{
    for(int i=1;i<=MAXN;i++)
    {
        int target=(i==1?1:0);
        int delta=target-mu[i];
        mu[i]=delta;
        for(int j=i+i;j<=MAXN;j+=i)
            mu[j]+=delta;
    }
}
LL solve(int x,int y)
{
    LL res=0ll;
    for(int i=1,r;i<=x&&i<=y;i=r+1)
    {
        int dx=x/i,dy=y/i;
        r=min(x/dx,y/dy);
        res+=(LL)(sum[r]-sum[i-1])*dx*dy;
    }
    return res;
}
int main(int argc, char const *argv[])
{
    int n;
    cin>>n;
    ms(mu);
    getMobius();
    sum[0]=0;
    for(int i=1;i<=MAXN;i++)
        sum[i]=sum[i-1]+mu[i];
    while(n--){
        int a,b,c,d,k;
        scanf("%d %d %d %d %d",&a,&b,&c,&d,&k);
        a=(a+k-1)/k,c=(c+k-1)/k;//向上取整
        b/=k,d/=k;//向下取整
        LL ans=solve(b,d)-solve(a-1,d)-solve(b,c-1)+solve(a-1,c-1);
        printf("%lld\n",ans );
    }
    return 0;
}