为什么插入到堆中比插入到未排序的列表中更快?
why is this insertion into a heap faster than insertion into an unsorted list?
在我的堆和未排序列表中插入 100000000 个元素后,似乎堆插入实际上更快(12 秒 vs 20 秒)。为什么是这样?我相信堆插入是 O(logn)
而未排序的列表插入是 O(1)
。我还注意到我的堆插入实现实际上并没有随着输入的数量而扩展。这也让我很困惑。
这是我 运行:
的代码
int main ()
{
clock_t unsortedStart;
clock_t heapStart;
double unsortedDuration;
double heapDuration;
int num_pushes = 100000000;
int interval = 10000;
ofstream unsorted ("unsorted.txt");
ofstream heap ("heap.txt");
UnsortedPQ<int> unsortedPQ;
HeapPQ<int> heapPQ;
unsortedStart = clock();
for (int i = 0; i < num_pushes; ++i)
{
if (i % interval == 0) {
unsortedDuration = ( clock() - unsortedStart ) / (double) CLOCKS_PER_SEC;
unsorted << unsortedDuration << " " << i << endl;
}
unsortedPQ.insertItem(rand() % 100);
}
heapStart = clock();
for (int i = 0; i < num_pushes; ++i)
{
if (i % interval == 0) {
heapDuration = ( clock() - heapStart ) / (double) CLOCKS_PER_SEC;
heap << heapDuration << " " << i << endl;
}
heapPQ.insertItem(rand() % 100);
}
return 0;
}
这是插入的堆实现(使用std::vector
):
template <class T>
void HeapPQ<T>::insertItem(T data) {
//insert into back of heap (std::vector)
dataArray.push_back(data);
int i = dataArray.size() - 1;
//sifts the inserted element up
while (i != 0 && dataArray[(i - 1) / 2] > dataArray[i]) {
swap(dataArray[i], dataArray[(i - 1) / 2]);
i = (i - 1) / 2;
}
}
这是插入的未排序列表实现(使用std::list
):
//pushes element to the back of a std::list
template <class T>
void UnsortedPQ<T>::insertItem(T data) { dataList.push_back(data); }
堆中的插入是O(logn)
,这意味着每次插入最多需要O(logn)
步。这并不意味着它必须。
在您的示例中,插入元素的平均成本是 O(1)
。为什么?
为简单起见,假设您仅以随机顺序插入 0
a 和 1
s(在当前版本中,只有数字 0..99
(rand() % 100
)插入 - 计算更复杂,但行为保持不变)。插入2*n
个元素后,堆中大约有n
0
和n
1
,堆如下:
0
0 0
00 00
...............
0 0 0 0 0 0 0
11 11 11 11 11 11 11
所以基本上,1
都在最后一级 k
,0
都在 0..k-1
。
- 如果插入了
1
,则无事可做(上面没有2
)。
- 如果插入
0
,则最多有一个交换(1
s可能在最后一层以上的水平,但在2层以上)。
这意味着我们平均只需要 0.5
交换而不是 k
。
具有相同的渐近 运行 时间,这完全取决于插入向量和列表的(摊销)成本。该列表似乎更慢(我的假设是,对于每个插入,它都需要通过 new
在堆上分配一个元素,这是一个非常慢的操作。成本取决于其他因素,例如大小插入的对象,因此可能会有所不同,哪个更快)。
让我们仔细看看您的案例,其中数字是由均匀分布生成的 [0..99]
。 n>>100
次插入后会出现如下情况(有些手忙脚乱,但大意应该很清楚):
- 堆的最后一层(
k
-th)有 n/2
个元素,由数字 50..99
组成。因此,对于 50% 的可能数字(即 50..99
),不需要转换。
- 堆的倒数第二层(
k-1
-th)有 n/4
个元素,由数字 25..49
组成。这意味着对于 25% 的可能数字,正好需要 1 个班次。
- 关卡
k-2
有 n/8
个元素,由数字 13..24
. 组成
log 100/log 2
层以上只有 0
层。所以可能的最大移位数是 m=log 100/log 2
,独立于 n
- 堆中的元素数。
所以插入的最坏情况成本是 log 100/log 2
,平均成本甚至更小:
E(insertion)=0*1/2+1*1/4+2*1/8+...<=1.0
即平均而言,我们每次插入的班次少于 1 个。
注意:这并不意味着在堆中插入的摊销成本为 O(1)
- 如果您不按随机顺序插入数字,而是首先插入所有 99
,然后 98
s,...,那么 0
s 每次插入的成本为 O(log n)
。
在我的堆和未排序列表中插入 100000000 个元素后,似乎堆插入实际上更快(12 秒 vs 20 秒)。为什么是这样?我相信堆插入是 O(logn)
而未排序的列表插入是 O(1)
。我还注意到我的堆插入实现实际上并没有随着输入的数量而扩展。这也让我很困惑。
这是我 运行:
的代码int main ()
{
clock_t unsortedStart;
clock_t heapStart;
double unsortedDuration;
double heapDuration;
int num_pushes = 100000000;
int interval = 10000;
ofstream unsorted ("unsorted.txt");
ofstream heap ("heap.txt");
UnsortedPQ<int> unsortedPQ;
HeapPQ<int> heapPQ;
unsortedStart = clock();
for (int i = 0; i < num_pushes; ++i)
{
if (i % interval == 0) {
unsortedDuration = ( clock() - unsortedStart ) / (double) CLOCKS_PER_SEC;
unsorted << unsortedDuration << " " << i << endl;
}
unsortedPQ.insertItem(rand() % 100);
}
heapStart = clock();
for (int i = 0; i < num_pushes; ++i)
{
if (i % interval == 0) {
heapDuration = ( clock() - heapStart ) / (double) CLOCKS_PER_SEC;
heap << heapDuration << " " << i << endl;
}
heapPQ.insertItem(rand() % 100);
}
return 0;
}
这是插入的堆实现(使用std::vector
):
template <class T>
void HeapPQ<T>::insertItem(T data) {
//insert into back of heap (std::vector)
dataArray.push_back(data);
int i = dataArray.size() - 1;
//sifts the inserted element up
while (i != 0 && dataArray[(i - 1) / 2] > dataArray[i]) {
swap(dataArray[i], dataArray[(i - 1) / 2]);
i = (i - 1) / 2;
}
}
这是插入的未排序列表实现(使用std::list
):
//pushes element to the back of a std::list
template <class T>
void UnsortedPQ<T>::insertItem(T data) { dataList.push_back(data); }
堆中的插入是O(logn)
,这意味着每次插入最多需要O(logn)
步。这并不意味着它必须。
在您的示例中,插入元素的平均成本是 O(1)
。为什么?
为简单起见,假设您仅以随机顺序插入 0
a 和 1
s(在当前版本中,只有数字 0..99
(rand() % 100
)插入 - 计算更复杂,但行为保持不变)。插入2*n
个元素后,堆中大约有n
0
和n
1
,堆如下:
0
0 0
00 00
...............
0 0 0 0 0 0 0
11 11 11 11 11 11 11
所以基本上,1
都在最后一级 k
,0
都在 0..k-1
。
- 如果插入了
1
,则无事可做(上面没有2
)。 - 如果插入
0
,则最多有一个交换(1
s可能在最后一层以上的水平,但在2层以上)。
这意味着我们平均只需要 0.5
交换而不是 k
。
具有相同的渐近 运行 时间,这完全取决于插入向量和列表的(摊销)成本。该列表似乎更慢(我的假设是,对于每个插入,它都需要通过 new
在堆上分配一个元素,这是一个非常慢的操作。成本取决于其他因素,例如大小插入的对象,因此可能会有所不同,哪个更快)。
让我们仔细看看您的案例,其中数字是由均匀分布生成的 [0..99]
。 n>>100
次插入后会出现如下情况(有些手忙脚乱,但大意应该很清楚):
- 堆的最后一层(
k
-th)有n/2
个元素,由数字50..99
组成。因此,对于 50% 的可能数字(即50..99
),不需要转换。 - 堆的倒数第二层(
k-1
-th)有n/4
个元素,由数字25..49
组成。这意味着对于 25% 的可能数字,正好需要 1 个班次。 - 关卡
k-2
有n/8
个元素,由数字13..24
. 组成
log 100/log 2
层以上只有0
层。所以可能的最大移位数是m=log 100/log 2
,独立于n
- 堆中的元素数。
所以插入的最坏情况成本是 log 100/log 2
,平均成本甚至更小:
E(insertion)=0*1/2+1*1/4+2*1/8+...<=1.0
即平均而言,我们每次插入的班次少于 1 个。
注意:这并不意味着在堆中插入的摊销成本为 O(1)
- 如果您不按随机顺序插入数字,而是首先插入所有 99
,然后 98
s,...,那么 0
s 每次插入的成本为 O(log n)
。