本篇是in place algorithm的学习笔记。目前学习的是in place merge与in place martrix transposition这两个算法。

1.in place merge

论文链接：

 论文讨论的是如何O(n)时间复杂度，O(1)空间复杂度合并两个相邻的有序数组。

b) 将sublist1与sublist2按sqrt(n)进行block划分, 余下的尾元素移到开头作为buffer( buffer size >= sqrt(n)，我们不关注buffer内的元素是否有序 )

c) 选择排序( O(sqrt(n)*sqrt(n)) ), 将各个block排序，使得各个block的末尾元素呈非下降序列。

 

a) 找到首个block尾元素a[i]，使得a[i] > a[i+1]，显然之前的满足a[i] <= a[i+1]则之前的已成有序序列。划分好series1和series2.

b) c) 将series1与series2通过buffer进行merge, 直到series1内的元素全部排到正确位置，如c) 所示。

　　由于按段尾非下降排序各个block，显然series1会先于series2排完。

 

 a) b) c) 将series2当作series1，series2的下一个block当作series2，重复之前的操作。以上操作的总的时间复杂度尾O(n).

 d) sqrt(n) <= buffer size < 2*sqrt(n)， 冒泡排序/选择排序，将buffer变为有序序列，时间复杂度O(sqrt(n)*sqrt(n)).

  

附: 论文中，buffer中的元素是最大的，最后一步操作sort buffer内的元素，即可。

个人想法: 我们需要保证 buffer元素 >= 非buffer元素，否则我们仍然需要insert sort。因此，本人认为在merge前需要进行处理，从sublist1和sublist2中抠出前sqrt(n)大的元素，作为buffer。但这样可能会破坏block，导致某个block元素不够。个人认为可以抠多一点，再补回去。比如先扣除前 3*sqrt(n) 大的元素。这样补回到sublist1，至少还能余下sqrt(n)个元素。

总的时间复杂度为O(n). 可能常数略大。

 

2.in place matrix transposition

wiki:

 对于非方阵，进行置换操作。对于一个Circle，记录下初始节点，每次用前驱节点替换当前节点；用初始节点替换最后一个节点。

Non-square matrices: Following the cycles

for each length>1 cycle C of the permutation    pick a starting address s in C    let D = data at s    let x = predecessor of s in the cycle    while x ≠ s        move data from x to successor of x        let x = predecessor of x    move data from D to successor of s

Improving memory locality at the cost of greater total data movement

 reference to