两排序好的数组,找中位数

描述
There are two sorted arrays A and B of size m and n respectively. Find the median of the two sorted
arrays. The overall run time complexity should be O(log(m + n)).
分析
这是一道非常经典的题。这题更通用的形式是,给定两个已经排序好的数组,找到两者所有元
素中第 k 大的元素。
O(m + n) 的解法比较直观,直接 merge 两个数组,然后求第 k 大的元素。
不过我们仅仅需要第 k 大的元素,是不需要“排序”这么复杂的操作的。可以用一个计数器,
记录当前已经找到第 m 大的元素了。同时我们使用两个指针 pA 和 pB,分别指向 A 和 B 数组的第
一个元素,使用类似于 merge sort 的原理,如果数组 A 当前元素小,那么 pA++,同时 m++;如果
数组 B 当前元素小,那么 pB++,同时 m++。最终当 m 等于 k 的时候,就得到了我们的答案,O(k)
时间,O(1) 空间。但是,当 k 很接近 m + n 的时候,这个方法还是 O(m + n) 的。
有没有更好的方案呢?我们可以考虑从 k 入手。如果我们每次都能够删除一个一定在第 k 大元
素之前的元素,那么我们需要进行 k 次。但是如果每次我们都删除一半呢?由于 A 和 B 都是有序
的,我们应该充分利用这里面的信息,类似于二分查找,也是充分利用了“有序”。
假设 A 和 B 的元素个数都大于 k/2,我们将 A 的第 k/2 个元素(即 A[k/2-1])和 B 的第 k/2
个元素(即 B[k/2-1])进行比较,有以下三种情况(为了简化这里先假设 k 为偶数,所得到的结
论对于 k 是奇数也是成立的):
• A[k/2-1] == B[k/2-1]
• A[k/2-1] > B[k/2-1]
• A[k/2-1] < B[k/2-1]
如果 A[k/2-1] < B[k/2-1],意味着 A[0] 到 A[k/2-1 的肯定在 A ∪ B 的 top k 元素的范围
内,换句话说,A[k/2-1 不可能大于 A ∪ B 的第 k 大元素。留给读者证明。
因此,我们可以放心的删除 A 数组的这 k/2 个元素。同理,当 A[k/2-1] > B[k/2-1] 时,可
以删除 B 数组的 k/2 个元素。
当 A[k/2-1] == B[k/2-1] 时,说明找到了第 k 大的元素,直接返回 A[k/2-1] 或 B[k/2-1]
即可。
因此,我们可以写一个递归函数。那么函数什么时候应该终止呢?
• 当 A 或 B 是空时,直接返回 B[k-1] 或 A[k-1];
• 当 k=1 是,返回 min(A[0], B[0]);
• 当 A[k/2-1] == B[k/2-1] 时,返回 A[k/2-1] 或 B[k/2-1]

代码:

 public class Median {

     public static void main(String[] args) {
int[] A= {2,3,5};
int[] B= {4,5,6};
int m=A.length;
int n=B.length;
int total = m + n;
if(total%2==1) //判奇偶
System.out.println( find_kth(A, m, B, n, total / 2 + 1)); //奇数
else
System.out.println((find_kth(A, m, B, n, total / 2) +
find_kth(A, m, B, n, total / 2 + 1)) / 2);
//偶数,中间位如果是偶数要除以2
}
private static double find_kth(int A[], int m, int B[], int n, int k) {
if (m > n) return find_kth(B, n, A, m, k); //m,n大小相反,则发一下
if (m == 0) return B[k - 1]; //A[0]、B[0] k=2
if (k == 1) return Math.min(A[0], B[0]);
//将K分成两部分
int pa = Math.min(k / 2, m), pb = k - pa;
if (A[pa - 1] < B[pb - 1]) {
int[] Apa=new int[m-pa];
for(int i=0;i<m-pa;i++) {
Apa[i]=A[pa+i];
}
return find_kth(Apa, m - pa, B, n, k - pa);
}
else if (A[pa - 1] > B[pb - 1]) {
int[] Bpb=new int[n-pb];
for(int i=0;i<n-pb;i++) {
Bpb[i]=B[pa+i];
} return find_kth(A, m, Bpb, n - pb, k - pb);
}else
return A[pa - 1];
}
}

最新文章

  1. web socket (记录下来方便观看)
  2. 使用代理下载Unity AssetStore上资源的方法
  3. Shell 编程基础之 [ 与 [[ 的异同
  4. 精通jQuery选择器
  5. apiCloud中图片裁剪模块FNImageClip的使用
  6. 并行编译加快VS C++项目的编译速度
  7. STRUTS2获得session和request
  8. java通过抛异常来返回提示信息
  9. 【转】OSX键盘快捷键
  10. iOS核心笔记—源代码管理工具-GIT
  11. Docker 的两类存储资源 - 每天5分钟玩转 Docker 容器技术(38)
  12. SQL查询操作及子句优先级
  13. #Java学习之路——基础阶段(第七篇)
  14. @vue/cl构建得项目下,postcss.config.js配置,将px转化成rem
  15. 使用redis限制ip访问次数
  16. Linux上svn搭建
  17. [Golang] 从零开始写Socket Server(1): Socket-Client框架
  18. iOS - UITextView放在自定义cell里面-自适应高度
  19. 用delphi开发activex打印控件
  20. rook issues

热门文章

  1. idea没有绑远程地址,如何提交到github的空项目
  2. LitePal——安卓数据库library
  3. let和const命令整理
  4. C#启动或停止 计算机中“服务”
  5. 《Unity預計算即時GI》笔记:二、光照图
  6. Hyper-V 2016 配置管理系列(应用篇)
  7. ECLIPSE 取消自动更新
  8. 卸载gitlab
  9. java Vamei快速教程05 实施接口
  10. Python开发第二篇