当前位置： 首页 > article >正文

C++二分查找算法：有序矩阵中的第 k 个最小数组和

article 2024/11/14 6:13:36

本文涉及的基础知识点

二分查找算法合集

本题的简化

C++二分查找算法：查找和最小的 K 对数字 十分接近m恒等于2

题目

给你一个 m * n 的矩阵 mat，以及一个整数 k ，矩阵中的每一行都以非递减的顺序排列。
你可以从每一行中选出 1 个元素形成一个数组。返回所有可能数组中的第 k 个 最小 数组和。
示例 1：
输入：mat = [[1,3,11],[2,4,6]], k = 5
输出：7
解释：从每一行中选出一个元素，前 k 个和最小的数组分别是：
[1,2], [1,4], [3,2], [3,4], [1,6]。其中第 5 个的和是 7 。
示例 2：
输入：mat = [[1,3,11],[2,4,6]], k = 9
输出：17
示例 3：
输入：mat = [[1,10,10],[1,4,5],[2,3,6]], k = 7
输出：9
解释：从每一行中选出一个元素，前 k 个和最小的数组分别是：
[1,1,2], [1,1,3], [1,4,2], [1,4,3], [1,1,6], [1,5,2], [1,5,3]。其中第 7 个的和是 9 。
示例 4：
输入：mat = [[1,1,10],[2,2,9]], k = 7
输出：12
参数范围：
m == mat.length
n == mat.length[i]
1 <= m, n <= 40
1 <= k <= min(200, n ^ m)
1 <= mat[i][j] <= 5000
mat[i] 是一个非递减数组

分析

时间复杂度

O(mlog(500040)n+mkn)。GetLessKSum被调用m次，GetLessEqualSumNum共被调用mlog(500040)次。每次调用GetLessEqualSumNum，for循环共执行m次。
vRet.emplace_back极端情况下，可能被执行kn次。

主要函数介绍

GetLessKSum两行升序数据的最小k个和
GetLessEqualSumNum两行升序数据和小于等于iSum的组合数量

注意：nums[i]为正数，所以如果pre的数量大于k，只需要保留前k小，其它的被淘汰了。

二分

寻找第一个符合条件的iSum，条件如下：
和小于等于iSum的组合数量大于等于k。

代码

核心代码

class Solution {
public:
	int kthSmallest(vector<vector<int>>& mat, int k) {
		m_c = mat.front().size();
		m_iK = k;
		vector<int> pre = mat[0];
		for (int r = 1; r < mat.size(); r++)
		{
			pre = GetLessKSum(pre, mat[r]);
		}
		return pre.back();
	}
	vector<int> GetLessKSum(const vector<int>& pre, const vector<int>& cur)
	{
		int left = 0, right = 5000 * 40;
		while (right - left > 1)
		{
			const auto mid = left + (right - left) / 2;
			if (GetLessEqualSumNum(pre, cur, mid)>= m_iK)
			{
				right = mid;
			}
			else
			{
				left = mid;
			}
		}
		vector<int> vRet;
		for (const auto& pr : pre)
		{
			for (const auto& cu : cur)
			{
				if (pr + cu <= right)
				{
					vRet.emplace_back(pr + cu);
				}
				else
				{
					break;
				}
			}
		}
		sort(vRet.begin(), vRet.end());
		if (vRet.size() > m_iK)
		{
			vRet.erase(vRet.begin() + m_iK, vRet.end());
		}
		return vRet;
	}
	int GetLessEqualSumNum(const vector<int>& pre, const vector<int>& cur,int iSum)
	{
		int iNum = 0;
		for (const auto& pr : pre)
		{
			iNum += std::upper_bound(cur.begin(), cur.end(), iSum - pr)- cur.begin();
		}
		return iNum;
	}
	int m_iK;
	int m_c;
};

测试用例

template
void Assert(const T& t1, const T& t2)
{
assert(t1 == t2);
}

template
void Assert(const vector& v1, const vector& v2)
{
if (v1.size() != v2.size())
{
assert(false);
return;
}
for (int i = 0; i < v1.size(); i++)
{
Assert(v1[i], v2[i]);
}
}

int main()
{
vector<vector> mat;
int k;
int res;
{
Solution slu;
mat = { {1,3,11},{2,4,6} };
k = 5;
res = slu.kthSmallest(mat, k);
Assert(7, res);
}
{
Solution slu;
mat = { {1,3,11},{2,4,6} };
k = 9;
res = slu.kthSmallest(mat, k);
Assert(17, res);
}
{
Solution slu;
mat = { {1,10,10},{1,4,5},{2,3,6} };
k = 7;
res = slu.kthSmallest(mat, k);
Assert(9, res);
}
{
Solution slu;
mat = { {1,1,10},{2,2,9} };
k = 7;
res = slu.kthSmallest(mat, k);
Assert(12, res);
}

//CConsole::Out(res);

}

优化增加结果

vector<int> vRet;
	for (const auto& pr : pre)
	{
		for (const auto& cu : cur)
		{
			if (pr + cu < right)
			{
				vRet.emplace_back(pr + cu);
			}
			else
			{
				break;
			}
		}
	}
	while (vRet.size() < m_iK)
	{
		vRet.emplace_back(right);
	}

和小于right的数量<=k，如果不足够，则补right。时间复杂度由O(nk)降低到O(k+n)。

直接使用封装类

namespace NBinarySearch
{
	template<class INDEX_TYPE,class _Pr>
	INDEX_TYPE FindFrist(INDEX_TYPE left, INDEX_TYPE right, _Pr pr)
	{
		while (right - left > 1)
		{
			const auto mid = left + (right - left) / 2;
			if (pr(mid))
			{
				right = mid;
			}
			else
			{
				left = mid;
			}
		}
		return right;
	}
}


class Solution {
public:
	int kthSmallest(vector<vector<int>>& mat, int k) {
		m_c = mat.front().size();
		m_iK = k;
		vector<int> pre = mat[0];
		for (int r = 1; r < mat.size(); r++)
		{
			pre = GetLessKSum(pre, mat[r]);
		}
		return pre.back();
	}
	vector<int> GetLessKSum(const vector<int>& pre, const vector<int>& cur)
	{
		auto GetLessEqualSumNum = [&pre, &cur, this](const int iSum)-> bool
		{
			int iNum = 0;
			for (const auto& pr : pre)
			{
				iNum += std::upper_bound(cur.begin(), cur.end(), iSum - pr) - cur.begin();
			}
			return iNum >= m_iK;
		};
		const int right = NBinarySearch::FindFrist(0, 5000 * 40, GetLessEqualSumNum);		
		vector<int> vRet;
		for (const auto& pr : pre)
		{
			for (const auto& cu : cur)
			{
				if (pr + cu < right)
				{
					vRet.emplace_back(pr + cu);
				}
				else
				{
					break;
				}
			}
		}
		while (vRet.size() < m_iK)
		{
			vRet.emplace_back(right);
		}
		sort(vRet.begin(), vRet.end());
		if (vRet.size() > m_iK)
		{
			vRet.erase(vRet.begin() + m_iK, vRet.end());
		}
		return vRet;
	}
	int GetLessEqualSumNum(const vector<int>& pre, const vector<int>& cur,int iSum)
	{
		int iNum = 0;
		for (const auto& pr : pre)
		{
			iNum += std::upper_bound(cur.begin(), cur.end(), iSum - pr)- cur.begin();
		}
		return iNum;
	}
	int m_iK;
	int m_c;
};

2023年3月暴力版

直接保留前k个。时间复杂度:O(mknlogk)
class Solution {
public:
int kthSmallest(vector<vector>& mat, int k) {
m_r = mat.size();
m_c = mat[0].size();
std::priority_queue pre;
pre.push(0);
for (int r = 0; r < mat.size(); r++)
{
std::priority_queue dp;
while (pre.size())
{
int t = pre.top();
pre.pop();
for (int c = 0; c < m_c; c++)
{
dp.push(mat[r][c] + t);
if (dp.size() > k)
{
dp.pop();
}
}
}
pre.swap(dp);
}
return pre.top();
}
int m_r, m_c;
};

扩展阅读

视频课程

有效学习：明确的目标 及时的反馈 拉伸区（难度合适），可以先学简单的课程，请移步CSDN学院，听白银讲师（也就是鄙人）的讲解。
https://edu.csdn.net/course/detail/38771

如何你想快

速形成战斗了，为老板分忧，请学习C#入职培训、C++入职培训等课程
https://edu.csdn.net/lecturer/6176

相关下载

想高屋建瓴的学习算法，请下载《喜缺全书算法册》doc版
https://download.csdn.net/download/he_zhidan/88348653

我想对大家说的话
闻缺陷则喜是一个美好的愿望，早发现问题，早修改问题，给老板节约钱。
墨子曰：事无终始，无务多业。也就是我们常说的专业的人做专业的事。
如果程序是一条龙，那算法就是他的是睛

测试环境

操作系统：win7 开发环境： VS2019 C++17
或者 操作系统：win10 开发环境：

VS2022 C++17

本文涉及的基础知识点

C++算法：前缀和、前缀乘积、前缀异或的原理、源码及测试用例 包括课程视频


http://www.kler.cn/a/134333.html

相关文章：

  • MyBatisPlus 用法详解
  • DApp开发：定制化解决方案与源码部署的一站式指南
  • vue3项目中内嵌vuepress工程两种实现方式
  • Spark 的容错机制：保障数据处理的稳定性与高效性
  • 生成 Django 中文文档 PDF 版
  • 【设计模式】关联关系与依赖关系
  • 数据结构 堆
  • 配置iTerm2打开自动执行命令
  • 打开游戏提示xapofx1_5.dll丢失如何修复？xapofx1_5.dll缺失的修复教程分享
  • 从一到无穷大 #19 TagTree，倒排索引入手是否是优化时序数据库查询的通用方案？
  • 滚雪球学Java(09-5)：Java中的赋值运算符，你真的掌握了吗？
  • mac 打不开 idea 或者 pycharm 的方法
  • JVM垃圾回收相关概念
  • V100 GPU服务器安装GPU驱动教程
  • C语言二叉树的建立和遍历
  • 【论文复现】DAE：《Annotating and Modeling Fine-grained Factuality in Summarization》
  • 云原生专栏丨基于服务网格的企业级灰度发布技术
  • Go语言常用命令详解(二)
  • 【手写数据库toadb】数据库planner的整体架构，以及逻辑查询树的设计与实现流程
  • 深度学习损失函数
  • 微信小程序手写table表格
  • DevToys：开发者的多功能瑞士军刀，让编程更高效！
  • python基础 | 类与继承
  • 小迪安全笔记（2）——web应用架构搭建漏洞HTTP数据包代理服务器
  • hive sql多表练习
  • 【原创】WeChat Server搭建