hdu 1695 GCD （容斥定理+预处理所有数的素因子）

题目链接：hdu 1695

参考链接如下：https://blog.csdn.net/yang_7_46/article/details/9072533

 

题意：分别在 [ 1, b ] 和 [ 1, d ] 的范围内找出x ，y ( ) ，满足gcd( x, y )=k,问有多少对？注意（（5,7）和（7,5））视为一对。

 

解题思路：由gcd（x，y）=k易得 gcd（x/k , y/k)=1 ,故我们可将题目看成是在 范围 [ 1, x/k ] 和 [ 1, y/k ] 中找出互质对，对此我们可以假设 x<y ，那么我们可以枚举 y ，然后在范围 [1 , min(b/k,y) ] （min来保证x<y）中找出与y互质的有多少个，我们发现直接找互质会比较难搞，所以我们就换成找不互质的有多少个，即将y分解质因子，如果w是y的素因子，那么[1 , min(b/k,y) ] 中是w的倍数有 min(b/k,y)/w个。

用容斥定理搞搞就行了。所有不与x互素的数的个数= 1个因子倍数的个数 - 2个因子乘积的倍数的个数 + 3个……-……。

 

#include <cstdio>
#include <cstring>
#include <vector>
#include <algorithm>
using namespace std;
#define N 100005
typedef long long ll;
vector<int> x[N];
bool book[N];

void prime() ///筛出每个数的素因子
{
    memset(book, false, sizeof(book));
    for (int i=0; i<N; i++) x[i].clear();

    for (int j=2; j<N; j+=2) x[j].push_back(2);
    for (int i=3; i<N; i+=2)
        if (!book[i])
        {
            for (int j=i; j<N; j+=i)
            {
                book[j] = true;
                x[j].push_back(i);
            }
        }
}
int work(int u, int s, int w)
{
    int cnt = 0, v = 1;
    for (int i=0; i<x[w].size(); i++)
    {
        if ((1<<i) & s)
        {
            cnt++;
            v *= x[w][i];///素因子
        }
    }
    int all = u/v; ///表示 是v的倍数有u/v个
    if (cnt % 2 == 0) return -all;
    else return all;
}

int main()
{

    prime();

    int T, a, b, c, d, k;
    scanf("%d", &T);
    for (int cas=1; cas<=T; cas++)
    {
        scanf("%d%d%d%d%d", &a, &b, &c, &d, &k);
        if (k == 0)
        {
            printf("Case %d: 0\n", cas);
            continue;
        }
        b /= k, d /= k; ///预处理

        if(b>d) swap(b,d); ///假设x<y
        
        long long ans = 0;
        for (int i=1; i<=d; i++)
        {
            k = min(i, b); ///表示在范围(1,k)与d互质的个数，这里是没除开与d不互质的个数
            ans += k; ///(1,k)与d互质
            for (int j=1; j<(1<<x[i].size()); j++)///容斥定理,这里巧妙的运用二进制，i的素因子有n个,
                ans -= work(k, j, i);///那么我们就用有n位的二进制数来表示，每个值的二进制代表取或不取，
        }                            ///例如101，表示取第1个和第3个因子
        printf("Case %d: %lld\n", cas, ans);
    }
    return 0;
}

 

我的标签：做个有情怀的程序员。