GCD HDU - 1695

LINK
题意:求1~n, 和 1 ~ m 中互质的本质不同的数对的数量。
设n < m, 则可以转变成求 1 ~ n 与 1 ~ n 的数对 和 1 ~ n 与 n + 1 ~ m 的 数对, 前者直接时欧拉函数的前缀和即可(phi[i] :小于 i 且与 i 互质的数的数量。
所以主要要考虑后者,很容易想到1 ~ n 和 j 互质的数 = n - (1 ~ n 与 j 不互质的数), 那怎么求不互质的数呢,
可以将j分解为质因数, 显然质因数的倍数和j不互质,
比如 2x, 3x, 所以我们要减去这部分, n / 2, n / 3, 但有一个关键点, 这样就有可能产生重复, 比如 6 既是2的倍数, 又是3的倍数, 这样就必须要用容斥定理去重 减去
n / (2 * 3)
这个方法网上有三种, 这里借鉴的是 DFS

#include <bits/stdc++.h>

using namespace std;

const int N = 1e5 + 10;
#define ll long long

int num[N] ={0};// i的素因子个数
int p[N][20]; // i的素因子
ll phi[N] = {0};



void phi_table()
{
    phi[1] = 1;
    for(int i = 2; i < N; i ++)
    {
        if(phi[i])
        {
            phi[i] += phi[i-1];
            continue;
        }
        for(int j = i; j < N; j += i)
        {
            if(!phi[j]) phi[j] = j;
            phi[j] = phi[j] * (i - 1) / i;
            p[j][num[j] ++] = i;

        }
        phi[i] += phi[i - 1];
    }

}

int DFS(int k, int b, int cur)// 求[1, b]中与cur不互质的数的个数
{
    int rel = 0;
    for(int i = k; i < num[cur]; i ++)
    {
        rel += b / p[cur][i] - DFS(i + 1, b / p[cur][i], cur);
    }
    return rel;
}
int main()
{
    phi_table();
    int t, caseno = 1;
    cin >> t;
    while(t --)
    {
        int a, b, c, d, k;
        cin >> a >> b >> c >> d >> k;
        if(k == 0)
        {
            cout << "Case " << caseno ++ << ": " << 0 << endl;
            continue;
        }
        if(b > d) swap(b, d);
        b /= k, d /= k;
        ll ans = phi[b];
        for(int i = b + 1; i <= d; i ++)
        {
            ans += b - DFS(0, b, i);
        }
        cout << "Case " << caseno ++ << ": " << ans << endl;
    }
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值