Hdu 1695 - GCD (容斥)

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

题目大意:
给定区间 [a,b] 与区间 [c,d] ,求两区间内各取一数gcd为k的种数, a,b b,a 视为同一种

分析:
由于题目给定了一个特殊条件,a=c=1,所以其实转化为了求 [1,b/k] [1,d/k] 内互质数对的个数
首先一个区间内与x互质的数的个数我们可以用容斥快速的求得,而区间又不是很大,所以可以直接暴力整个区间,求所有的互质数对

代码:

//#include<bits/stdc++.h>
#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
typedef long long ll;

const ll mod = 1e9+7;

ll a,b,c,d,k;
const int maxn = 100000+10;
vector<int> ft[120005];
bool isprime[120000];

void init()
{
    for (int i = 2 ; i <= maxn ; i ++)
    {
        if (!isprime[i])
        {
            for (int j = i ; j <= maxn ; j += i)
            {
                isprime[j] = true;
                ft[j].push_back(i);
            }
        }
    }

}

ll solve(int x,int sta,ll n)
{
    int pos = 0;
    ll temp = 1;
    while (sta)
    {
        if (sta&1)
            temp *= ft[x][pos];
        pos ++;
        sta >>= 1;
    }

    return n/temp;
}



int main(){
    init();
    int T,t=1;
    scanf("%d",&T);
    while (T--)
    {
        scanf("%lld%lld%lld%lld%lld",&a,&b,&c,&d,&k);
        printf("Case %d: ",t++);
        if (k==0)
        {
            printf("0\n");
            continue;
        }
        b/=k;d/=k;

        if (b*d==0)
            printf("0\n");
        else
        {
            if (b>d)
                swap(b,d);
            ll ans = 0;
            for (int i = 1 ; i <= d ; i ++)
            {
                ans += min(i*1LL,b);
                for (int j = 1 ; j < (1<<ft[i].size()); j ++)
                {
                    ll temp = solve(i,j,min(i*1LL,b));int cp = j,cnt = 0;
                    //printf("cp = %d\t",cp);
                    while (cp)
                    {
                        cp -= cp&(-cp);
                        cnt++;
                    }
                    //printf("cnt = %d\n",cnt);
                    if (cnt&1)  ans -= temp;
                    else    ans += temp;
                }
            }
            printf("%lld\n",ans);

        }

    }



    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值