786. 第 K 个最小的素数分数

786. 第 K 个最小的素数分数

给你一个按递增顺序排序的数组 arr 和一个整数 k 。数组 arr 由 1 和若干 素数  组成,且其中所有整数互不相同。

对于每对满足 0 < i < j < arr.length 的 i 和 j ,可以得到分数 arr[i] / arr[j] 。

那么第 k 个最小的分数是多少呢?  以长度为 2 的整数数组返回你的答案, 这里 answer[0] == arr[i] 且 answer[1] == arr[j] 。

示例 1:

输入:arr = [1,2,3,5], k = 3
输出:[2,5]
解释:已构造好的分数,排序后如下所示: 
1/5, 1/3, 2/5, 1/2, 3/5, 2/3
很明显第三个最小的分数是 2/5

示例 2:

输入:arr = [1,7], k = 1
输出:[1,7]

提示:

  • 2 <= arr.length <= 1000
  • 1 <= arr[i] <= 3 * 1 0 4 10^4 104
  • arr[0] == 1
  • arr[i] 是一个 素数 ,i > 0
  • arr 中的所有数字 互不相同 ,且按 严格递增 排序
  • 1 <= k <= arr.length * (arr.length - 1) / 2

解题思路

题目需要求得第 k 个最小的分数,而分数定义为每对满足 0 < i < j < arr.length 的 i 和 j ,得到的arr[i] / arr[j] ,因为数组arr是按递增顺序排序的,因此分数都是小于0的数字。为了避免浮点数的精度丢失,当我们比较当我们比较两个分数  a b \dfrac{a}{b} ba和  c d \dfrac{c}{d} dc时,我们可以使用: a × d < b × c a \times d < b \times c a×d<b×c来替代 $\dfrac{a}{b} < \dfrac{c}{d} $的判断,二者是等价的,这样就可以避免计算机浮点数丢失造成的问题。因此我们在排序中需要根据上述规则进行自定义的排序,然后选出第k个最小的分数。

代码

class Solution {
public:
    vector<int> kthSmallestPrimeFraction(vector<int> &arr, int k) {

        vector<pair<int, int>> cnt;
        for (int i = 0; i < arr.size(); ++i) {

            for (int j = i - 1; j >= 0; --j) {
                cnt.emplace_back(arr[j], arr[i]);
            }
        }
        sort(cnt.begin(), cnt.end(), [&](const auto &x, const auto &y) {
            return x.first * y.second < x.second * y.first;
        });
        return {cnt[k-1].first,cnt[k-1].second};
    }
};
以下是10种JS判断一个数是否为质数的方法: 方法1:暴力枚举 ``` function isPrime(n) { if (n <= 1) { return false; } for (let i = 2; i < n; i++) { if (n % i === 0) { return false; } } return true; } ``` 方法2:优化暴力枚举 ``` function isPrime(n) { if (n <= 1) { return false; } const sqrt = Math.sqrt(n); for (let i = 2; i <= sqrt; i++) { if (n % i === 0) { return false; } } return true; } ``` 方法3:6k±1优化 ``` function isPrime(n) { if (n <= 3) { return n > 1; } if (n % 2 === 0 || n % 3 === 0) { return false; } const sqrt = Math.sqrt(n); for (let i = 5; i <= sqrt; i += 6) { if (n % i === 0 || n % (i + 2) === 0) { return false; } } return true; } ``` 方法4:埃氏筛法 ``` function getPrimes(n) { const primes = new Array(n + 1).fill(true); primes[0] = false; primes[1] = false; for (let i = 2; i * i <= n; i++) { if (primes[i]) { for (let j = i * i; j <= n; j += i) { primes[j] = false; } } } return primes; } function isPrime(n) { if (n <= 1) { return false; } const primes = getPrimes(n); return primes[n]; } ``` 方法5:线性筛法 ``` function getPrimes(n) { const primes = []; const isPrime = new Array(n + 1).fill(true); isPrime[0] = false; isPrime[1] = false; for (let i = 2; i <= n; i++) { if (isPrime[i]) { primes.push(i); } for (let j = 0; j < primes.length && i * primes[j] <= n; j++) { isPrime[i * primes[j]] = false; if (i % primes[j] === 0) { break; } } } return primes; } function isPrime(n) { if (n <= 1) { return false; } const primes = getPrimes(n); return primes.includes(n); } ``` 方法6:费马小定理 ``` function powMod(a, b, p) { let res = 1; while (b > 0) { if (b & 1) { res = (res * a) % p; } a = (a * a) % p; b >>= 1; } return res; } function isPrime(n) { if (n <= 1) { return false; } if (n <= 3) { return true; } const k = 5; for (let i = 0; i < k; i++) { const a = Math.floor(Math.random() * (n - 3) + 2); if (powMod(a, n - 1, n) !== 1) { return false; } } return true; } ``` 方法7:米勒-拉宾素性检验 ``` function powMod(a, b, p) { let res = 1; while (b > 0) { if (b & 1) { res = (res * a) % p; } a = (a * a) % p; b >>= 1; } return res; } function isPrime(n) { if (n <= 1) { return false; } if (n <= 3) { return true; } const k = 10; let d = n - 1; while (d % 2 === 0) { d /= 2; } for (let i = 0; i < k; i++) { const a = Math.floor(Math.random() * (n - 3) + 2); let x = powMod(a, d, n); if (x === 1 || x === n - 1) { continue; } let flag = false; while (d !== n - 1) { x = (x * x) % n; d *= 2; if (x === 1) { return false; } if (x === n - 1) { flag = true; break; } } if (!flag) { return false; } } return true; } ``` 方法8:埃及分数求解 ``` function gcd(a, b) { return b === 0 ? a : gcd(b, a % b); } function isPrime(n) { if (n <= 1) { return false; } if (n === 2 || n === 3) { return true; } let sum = 0; for (let i = 2; i < n; i++) { if (gcd(i, n) === 1) { sum += 1 / i; } } return sum + 1 === n; } ``` 方法9:质因数分解 ``` function isPrime(n) { if (n <= 1) { return false; } let cnt = 0; for (let i = 2; i <= n; i++) { while (n % i === 0) { n /= i; cnt += 1; } } return cnt === 1; } ``` 方法10:Lucas素性检验 ``` function powMod(a, b, p) { let res = 1; while (b > 0) { if (b & 1) { res = (res * a) % p; } a = (a * a) % p; b >>= 1; } return res; } function lucasTest(n, p) { if (n === 0) { return 0; } if (n === 1) { return 1; } const q = Math.floor(n / p); const r = n % p; return powMod(lucasTest(q, p), 2, p) * powMod(2, r, p) % p; } function isPrime(n) { if (n <= 1) { return false; } if (n === 2 || n === 3) { return true; } const primes = [2, 3, 5, 7, 11, 13, 17, 19, 23, 29]; for (const p of primes) { if (n === p) { return true; } if (n % p === 0) { return false; } if (lucasTest(n, p) !== 1) { return false; } } return true; } ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值