算法竞赛数论杂题

menji 和 gcd

题目：

一开始以为是只有l不确定，r是确定的，这样的话我们可以枚举r的所有约数，然后对其每个约数x进行判断，判断是否满足题意，具体做法是先让l % x如果 == 0则该约数可行，如果不可行且余数是y，那么就让l 加上x - y判断这个数是否比r大。时间复杂度o(sqrt(R));

但是现在两端点都是不确定的，注意到 l <= k * g < (k + 1) * g <= r 注意到我们想要得到的gcd与k的乘积是不大于r的，因此这两个未知量的较小值是不大于sqrt(R)的，因此我们可以枚举较小值i，另一个值自然可以通过R / i得到 由于在循环中k * g <= r 是一定的，因此在判断方案是否合法时只需要判断左边界l是否满足即可

代码：

#include <bits/stdc++.h>

using namespace std;

typedef long long LL;

void solve() {
    LL l, r;
    cin >> l >> r;
    LL ans = 0;
    for(LL i = 1; i <= r / i; i ++ ) {
        LL j = r / i;
        if((i - 1) * j >= l) ans = max(ans, j);
        if((j - 1) * i >= l) ans = max(ans, i); 
    }
    cout << ans << endl;
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    int t;
    cin >> t;
    while(t -- ) {
        solve();
    }
    return 0;
}

模方程

题目：

显而易见，如果a = b x将会有无穷个，如果 a <  b x将为0

现在思考a > b的情况

a = k * x + b

x = (a - b) / k;

因此k一定是(a - b)的一个因数

枚举所有因数t并判断原式是否成立

代码：

#include <bits/stdc++.h>

using namespace std;

int gcd(int a, int b) {
    return b? gcd(b, a % b): a;
}

int main() {
    int a, b;
    cin >> a >> b;
    vector<int> d;
    if(b == a) cout << "infinity";
    else {
        if(b > a) cout << "0";
        else {
            int tot = a - b;
            for(int i = 1; i <= tot / i; i ++ ) {
                if(tot % i == 0) {
                    d.push_back(i);
                    if(i != tot / i) d.push_back(tot / i);
                }
            }
            int ans = 0;
            for(auto t: d) {
                int x = tot / t;
                if(x < a) ans ++;
            }
            cout << ans;
        }
    }
    return 0;
}