题目地址:
https://leetcode.com/problems/nth-magical-number/description/
给定三个正整数 n , a , b n, a, b n,a,b,求第 n n n个能被 a a a或 b b b整除的正整数。
这个数 x x x即为满足 1 ∼ x 1\sim x 1∼x中有 n n n个正整数能被 a a a或 b b b整除的最小数。显然 x x x越大, 1 ∼ x 1\sim x 1∼x能被 a a a或 b b b整除的数的个数就越多,所以可以二分,左端点 1 1 1,右端点 n min { a , b } n\min\{a,b\} nmin{a,b}。而 1 ∼ x 1\sim x 1∼x能被 a a a或 b b b整除的数的个数可以用容斥原理,即为 ⌊ x a ⌋ + ⌊ x b ⌋ − ⌊ x ( a , b ) ⌋ \lfloor\frac{x}{a}\rfloor+\lfloor\frac{x}{b}\rfloor-\lfloor\frac{x}{(a,b)}\rfloor ⌊ax⌋+⌊bx⌋−⌊(a,b)x⌋,其中 ( a , b ) (a,b) (a,b)为 a , b a,b a,b的最小公倍数,其等于 a b gcd ( a , b ) \frac{ab}{\gcd(a, b)} gcd(a,b)ab。代码如下:
class Solution {
public:
int nthMagicalNumber(int n, int a, int b) {
const int MOD = 1e9 + 7;
int g = gcd(a, b);
auto check = [&](long x) { return x / a + x / b - x / (a * b / g) >= n; };
long l = 1, r = (long)min(a, b) * n;
while (l < r) {
long mid = l + (r - l >> 1l);
if (check(mid)) r = mid;
else l = mid + 1;
}
return l % MOD;
}
int gcd(int x, int y) { return y ? gcd(y, x % y) : x; }
};
时间复杂度 O ( log n + log min { a , b } ) O(\log n+\log \min \{a,b\}) O(logn+logmin{a,b}),空间 O ( log min { a , b } ) O(\log \min \{a,b\}) O(logmin{a,b})。