思维数学(互质数）

思维数学题

题目：对于给定的整数 l, r, k ，问有多少个数 x 满足 x + k 与 x − k 互质，且 x + k, x − k ∈ [l, r]

solution：
关于互质数的一个简单结论，如果a,b互质，则gcd(a, a- b) =gcd(a, a + b) = gcd(b, a - b) = gcd(b, a + b) = 1.
说人话就是如果a, b互质，则两数a,b与其和，其差也互质。
证明：以a 和 a - b互质为例
反证法：gcd(a, a -b) >= d > 1
则d | a, d | (a - b), -> d | b
即a, b含有大于1的公约数，这与a,b互质矛盾.

所以gcd(x + k, x - k) == 1 -> gcd(x + k, 2k) == 1 || gcd(x - k, 2k) == 1.
所以我们枚举满足条件的x - k的数字的个数就是答案。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

vector<ll>fac;
void calc(ll n){
    fac.clear();
    for(int i = 2; (ll)i * i <= n; ++i){
        if(n % i ==0){
            fac.push_back(i);
            while(n % i == 0) n /= i;
        }
    }
    if(n > 1) fac.push_back(n);
}

ll solve(ll x){
    if(x <= 0) return 0;
    int n = fac.size();
    int up = (1 << n);
    ll res = x;
    for(int i = 1; i < up; ++i){
        ll tmp = 1;
        for(int j = 0; j < n; ++j){
            if(i & (1 << j)) tmp *= fac[j];
        }
        if(__builtin_parity(i)) res -= x / tmp;
        else res += x / tmp;
    }
    return res;
}

int main(){
    ios::sync_with_stdio(0);
    cin.tie(0), cout.tie(0);
    ll l, r, k;
    cin >> l >> r >> k;
    ll x = 2 * k;
    r -= x;
    if(l > r){
        cout << 0 << endl;
        return 0;
    }   
    calc(x);
    ll ans = solve(r) - solve(l - 1);
    cout << ans << endl;
    return 0;
}