HDU 1695

题意:
对于给出的 n 个询问，每次求有多少个数对 (x,y) ，
满足 a ≤ x ≤ b ， c ≤ y ≤ d ，
且 gcd(x,y) = k ， gcd(x,y) 函数为 x 和 y 的最大公约数。
(x,y)和(y,x)算一对
思路:
跟BZOJ 2301差不多
减去重复的对数,就是n和m的区间交的区间

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<queue>
#include<stack>
#include<string>
#include<vector>
#include<map>
#include<set>
using namespace std;
#define lowbit(x) (x&(-x))
typedef long long LL;
const int maxn = 100005;
const int inf=(1<<28)-1;
#define maxp 100005
bool notprimes[maxp];
int primes[maxp];
int mu[maxp];
void get_mu()
{
    memset(notprimes,false,sizeof(notprimes));
    primes[0]=0;
    mu[1]=1;
    for(int i=2;i<maxp;++i)
    {
        if(!notprimes[i])
        {
            primes[++primes[0]]=i;
            mu[i]=-1;
        }
        for(int j=1;j<=primes[0];++j)
        {
            if((LL)primes[j]*i>=maxp) break;
            notprimes[i*primes[j]]=true;
            if(i%primes[j])
                mu[i*primes[j]]=-mu[i];
            else
            {
                mu[i*primes[j]]=0;
                break;//代表i不是素数,mu[i*primes[j]]必然是0 
            }
        }
    }
}

LL Pre[maxp];
LL Solve(int n,int m,int k)
{
    n/=k,m/=k;
    int t=min(n,m);
    LL Ans=0;
    int last;
    for(int i=1;i<=t;i=last+1)
    {
        last=min(n/(n/i),m/(m/i));
        Ans+=(LL)(Pre[last]-Pre[i-1])*(n/i)*(m/i);
    }
    //printf("%lld ",Ans);
    return Ans;
}

int main()
{
    get_mu();
    Pre[0]=0;
    for(int i=1;i<maxp;++i)
        Pre[i]=Pre[i-1]+mu[i];
    int T,Case=0;
    scanf("%d",&T);
    while(T--)
    {
        int a,b,c,d,k;
        scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
        if(k==0)
        {
            printf("Case %d: 0\n",++Case);
            continue;
        }
        LL Ans=Solve(b,d,k)-Solve(b,c-1,k)-Solve(a-1,d,k)+Solve(a-1,c-1,k);
        int Left=max(a,c),Right=min(b,d);
        if(Left<=Right)
        {
            int flag=0;
            if(Left<=k&&k<=Right) flag=1;
            Ans-=(Solve(Right,Right,k)-Solve(Left-1,Right,k)-Solve(Right,Left-1,k)+Solve(Left-1,Left-1,k))/2;
        }
        printf("Case %d: %lld\n",++Case,Ans);
    }
    return 0;
}