HDU 1695 GCD

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1695

题意：在[a,b] 和 [c,d]区间里各找一个数x和y，使得gcd(x,y) = k

思路：题目提示a和c都是1，如果不规定是1，也可以通过区间减法来计算答案。那么这题的思路是在[1,b/k] 和 [1,d/k]范围内找出互质的数，再将两数同时乘k，那么就变成范围内且gcd为k的一组数了。同时注意不要计算重复，(x,y)和(y,x)只算一次。所以计算的时候满足x<y即可。
枚举x，在[1，b/k]范围内，找出互质的数的个数，然后减去x和区间[1,x-1]互质的个数。就是x和区间[x+1,b/k]的互质个数。注意特殊情况即可。

#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <stack>
#include <map>
#include <set>
#include <vector>
#include <sstream>
#include <queue>
#include <utility>
using namespace std;

#define rep(i,j,k) for (int i=j;i<=k;i++)
#define Rrep(i,j,k) for (int i=j;i>=k;i--)

#define Clean(x,y) memset(x,y,sizeof(x))
#define LL long long
#define ULL unsigned long long
#define inf 0x7fffffff
#define mod %100000007

int T,a,b,c,d,k;
LL ans;

int prime[100009][15];


void init()
{
    Clean(prime,0);
    rep(i,2,100000)
    {
        int temp = i;
        for(int j = 2; j*j <= temp; j++)
            if ( temp % j == 0 )
            {
                prime[i][++prime[i][0]] = j;
                while( temp % j == 0 ) temp /= j;
            }
        if ( temp > 1 ) prime[i][++prime[i][0]] = temp;
    }
}

int query(int n , int m) //容斥求互质个数
{
    int ans = 0;
    int uplim = 1<<prime[m][0];
    rep(i,1,uplim-1)
    {
        int cnt = 1;
        int num = 0;
        rep(j,0,prime[m][0]-1)
        if ( i & (1<<j) )
        {
            cnt*=prime[m][j+1];
            num++;
        }
        if ( num & 1 )
            ans += n/cnt;
        else
            ans -= n/cnt;
    }
    return n-ans;
}

int main()
{
    init();
    cin>>T;
    rep(kase,1,T)
    {

        scanf("%d %d %d %d %d",&a,&b,&c,&d,&k);
        if( k == 0 || b < k || d < k )
        {
            printf("Case %d: 0\n",kase);
            continue;
        }
        b/=k;
        d/=k;
        if ( b > d ) swap(b,d);
        if ( b > 0 )ans = d;
        rep(i,2,b)
            ans += query( d , i ) - query( i-1 , i );
        printf("Case %d: %I64d\n", kase , ans);
    }
    return 0;
}