Matrix decomposition: Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
en>Jacobkhed
en>Ethaniel
m →‎Example: -> Orthogonal only, not unitary
Line 1: Line 1:
In [[computer science]], a '''Fibonacci heap''' is a [[Heap (data structure)|heap data structure]] consisting of a collection of [[Tree (data structure)|tree]]s. It has a better [[amortized analysis|amortized]] running time than a [[binomial heap]]. Fibonacci heaps were developed by [[Michael Fredman|Michael L. Fredman]] and [[Robert Tarjan|Robert E. Tarjan]] in 1984 and first published in a scientific journal in 1987. The name of Fibonacci heap comes from [[Fibonacci number]]s which are used in the running time analysis.
With your offense, you might experience Gunboats which can extremely easily shoot at enemy resistance coming from a very long range and Landing Crafts which you must fill up when you train units for example Rifleman, Heavy, Zooka, Warrior and Containers. To your village defenses, you might carry structures like Mortar, Bike Gun, Sniper Tower, Cannon, Flamethrower, Mine, Tank Mine, Boom Cannon and Catapult Launcher to assist you eradicate enemies.<br><br>Interweaving social styles form an absolute strong net in which we are all encased. When The Tygers of Pan Tang sang 'It's lonely at specific top. Everybody's undertaking to do you in', these people borrowed much from clash of [http://Search.Huffingtonpost.com/search?q=clans+hack&s_it=header_form_v1 clans hack] tool no survey. A society without clash of clans hack tool no survey is really like a society with no knowledge, in your it is quite smart.<br><br>Temperance is essential in almost things, and enjoying clash of clans cheats is no special. Playing for hours on finish seriously is not good-for-you, bodily or mentally. There are some games out here which know this and can include measures to remind in order to take rests. Consider the initiative yourself, although! Place an alarm a person don't play for a lot more an hour right.<br><br>Investigate evaluations and see currently the trailers before buying another video game. Makes it one thing you will be looking at before you get in which. These video games aren't low-cost, and also you likely to get nearly as much cash whenever you commerce inside a employed game title which you have purely utilized several times.<br><br>The nfl season is here and furthermore going strong, and like many fans we incredibly long for Sunday afternoon  when the games begin. If you have held and liked Soul Caliber, you will love this game. The succeeding best is the Microchip Cell which will at random , fill in some pieces. Defeating players like that by any usually means necessary can be that this reason that pushes folks to use Words that has Friends CheatIf you liked this information and you would certainly like to get more information concerning [http://prometeu.net clash of clans hack android no survey] kindly check out the website. Your current app requires you within order to answer 40 questions using varying degrees of hard times.<br><br>Had you been aware that some software games are educational tools? If you know a children that likes to participate in video games, educational options are a fantastic tactics to combine learning because of entertaining. The World wide web can connect you with thousands of parents which have similar values and tend to be more than willing of share their reviews combined with notions with you.<br><br>Disclaimer: I aggregate the guidance on this commodity by game a lot of CoC and accomplishing some web research. To the best involving my knowledge, is it authentic along with I accept amateur costed all abstracts and data. Nevertheless, it is consistently accessible we accept fabricated a aberration about or which the most important bold has afflicted back publication. Use as part of your very own risk, Do not accommodate virtually any provides. Please get in blow if that you acquisition annihilation amiss.
 
Find-minimum is ''O''(1) amortized time.<ref name="CLRS">[[Thomas H. Cormen]], [[Charles E. Leiserson]], [[Ronald L. Rivest]], and [[Clifford Stein]]. ''[[Introduction to Algorithms]]'', Second Edition. [[MIT Press]] and McGraw-Hill, 2001. ISBN 0-262-03293-7. Chapter 20: Fibonacci Heaps, pp.476&ndash;497. Third edition p518.</ref> Operations insert, decrease key, and merge (union) work in constant amortized time.<ref name="Fredman And Tarjan"/> Operations delete and delete minimum work in [[Big O notation|''O'']](log ''n'') amortized time.<ref name="Fredman And Tarjan"/> This means that starting from an empty data structure, any sequence of ''a'' operations from the first group and ''b'' operations from the second group would take ''O''(''a''&nbsp;+&nbsp;''b''&nbsp;log&nbsp;''n'') time. In a [[binomial heap]] such a sequence of operations would take ''O''((''a''&nbsp;+&nbsp;''b'') log (''n'')) time. A Fibonacci heap is thus better than a binomial heap when ''b'' is [[asymptotic analysis|asymptotically]] smaller than ''a''.
 
Using Fibonacci heaps for [[priority queue]]s improves the asymptotic running time of important algorithms, such as [[Dijkstra's algorithm]] for computing the [[shortest path]] between two nodes in a graph.
 
== Structure ==
[[Image:Fibonacci heap.png|thumbnail|250px|Figure 1. Example of a Fibonacci heap. It has three trees of degrees 0, 1 and 3. Three vertices are marked (shown in blue). Therefore the potential of the heap is 9 (3 trees + 2 * marked-vertices).]]
A Fibonacci heap is a collection of [[tree data structure|tree]]s satisfying the [[minimum-heap property]], that is, the key of a child is always greater than or equal to the key of the parent. This implies that the minimum key is always at the root of one of the trees. Compared with binomial heaps, the structure of a Fibonacci heap is more flexible. The trees do not have a prescribed shape and in the extreme case the heap can have every element in a separate tree. This flexibility allows some operations to be executed in a "lazy" manner, postponing the work for later operations. For example merging heaps is done simply by concatenating the two lists of trees, and operation ''decrease key'' sometimes cuts a node from its parent and forms a new tree.
 
However at some point some order needs to be introduced to the heap to achieve the desired running time. In particular, degrees of nodes (here degree means the number of children) are kept quite low: every node has degree at most ''O''(log ''n'') and the size of a subtree rooted in a node of degree ''k'' is at least  ''F<sub>k</sub>''<sub>&nbsp;+&nbsp;2</sub>, where ''F<sub>k</sub>'' is the ''k''th [[Fibonacci number]]. This is achieved by the rule that we can cut at most one child of each non-root node. When a second child is cut, the node itself needs to be cut from its parent and becomes the root of a new tree (see Proof of degree bounds, below). The number of trees is decreased in the operation ''delete minimum'', where trees are linked together.
 
As a result of a relaxed structure, some operations can take a long time while others are done very quickly. For the [[amortized analysis|amortized running time]] analysis we use the [[potential method]], in that we pretend that very fast operations take a little bit longer than they actually do. This additional time is then later combined and subtracted from the actual running time of slow operations. The amount of time saved for later use is measured at any given moment by a potential function. The potential of a Fibonacci heap is given by
 
:Potential = ''t'' + 2''m''
 
where ''t'' is the number of trees in the Fibonacci heap, and ''m'' is the number of marked nodes. A node is marked if at least one of its children was cut since this node was made a child of another node (all roots are unmarked).
The amortized time for an operation is given by the sum of the actual time and ''c'' times the difference in potential, where ''c'' is a constant (chosen to match the constant factors in the ''O'' notation for the actual time).
 
Thus, the root of each tree in a heap has one unit of time stored. This unit of time can be used later to link this tree with another tree at amortized time 0. Also, each marked node has two units of time stored. One can be used to cut the node from its parent. If this happens, the node becomes a root and the second unit of time will remain stored in it as in any other root.
 
== Implementation of operations ==
To allow fast deletion and concatenation, the roots of all trees are linked using a circular, [[doubly linked list]]. The children of each node are also linked using such a list. For each node, we maintain its number of children and whether the node is marked. Moreover we maintain a pointer to the root containing the minimum key.
 
Operation '''find minimum''' is now trivial because we keep the pointer to the node containing it. It does not change the potential of the heap, therefore both actual and amortized cost is constant. As mentioned above, '''merge''' is implemented simply by concatenating the lists of tree roots of the two heaps. This can be done in constant time and the potential does not change, leading again to constant amortized time.
 
Operation '''insert''' works by creating a new heap with one element and doing merge. This takes constant time, and the potential increases by one, because the number of trees increases. The amortized cost is thus still constant.
 
[[Image:Fibonacci heap extractmin1.png|170px|thumb|right|Fibonacci heap from Figure 1 after first phase of extract minimum. Node with key 1 (the minimum) was deleted and its children were added as separate trees.]]
Operation '''extract minimum''' (same as ''delete minimum'') operates in three phases. First we take the root containing the minimum element and remove it. Its children will become roots of new trees. If the number of children was ''d'', it takes time ''O''(''d'') to process all new roots and the potential increases by ''d''−1. Therefore the amortized running time of this phase is ''O''(''d'') = ''O''(log ''n'').
 
[[Image:Fibonacci heap extractmin2.png|130px|thumb|left|Fibonacci heap from Figure 1 after extract minimum is completed. First, nodes 3 and 6 are linked together. Then the result is linked with tree rooted at node 2. Finally, the new minimum is found.]]
However to complete the extract minimum operation, we need to update the pointer to the root with minimum key. Unfortunately there may be up to ''n'' roots we need to check. In the second phase we therefore decrease the number of roots by successively linking together roots of the same degree. When two roots ''u'' and ''v'' have the same degree, we make one of them a child of the other so that the one with the smaller key remains the root. Its degree will increase by one. This is repeated until every root has a different degree. To find trees of the same degree efficiently we use an array of length ''O''(log ''n'') in which we keep a pointer to one root of each degree. When a second root is found of the same degree, the two are linked and the array is updated. The actual running time is ''O''(log ''n'' + ''m'') where ''m'' is the number of roots at the beginning of the second phase. At the end we will have at most ''O''(log ''n'') roots (because each has a different degree). Therefore the difference in the potential function from before this phase to after it is: ''O''(log ''n'') − ''m'', and the amortized running time is then at most ''O''(log ''n'' + ''m'') + ''c''(''O''(log ''n'') − ''m''). With a sufficiently large choice of ''c'', this simplifies to ''O''(log ''n'').
 
In the third phase we check each of the remaining roots and find the minimum. This takes ''O''(log ''n'') time and the potential does not change. The overall amortized running time of extract minimum is therefore ''O''(log ''n'').
 
[[Image:Fibonacci heap-decreasekey.png|250px|thumb|right|Fibonacci heap from Figure 1 after decreasing key of node 9 to 0. This node as well as its two marked ancestors are cut from the tree rooted at 1 and placed as new roots.]]
Operation '''decrease key''' will take the node, decrease the key and if the heap property becomes violated (the new key is smaller than the key of the parent), the node is cut from its parent. If the parent is not a root, it is marked. If it has been marked already, it is cut as well and its parent is marked. We continue upwards until we reach either the root or an unmarked node. In the process we create some number, say ''k'', of new trees. Each of these new trees except possibly the first one was marked originally but as a root it will become unmarked. One node can become marked. Therefore the number of marked nodes changes by &minus;(''k''&nbsp;&minus;&nbsp;1)&nbsp;+&nbsp;1&nbsp;=&nbsp;&minus;&nbsp;''k''&nbsp;+&nbsp;2. Combining these 2 changes, the potential changes by 2(&minus;''k''&nbsp;+&nbsp;2)&nbsp;+&nbsp;''k''&nbsp;=&nbsp;&minus;''k''&nbsp;+&nbsp;4. The actual time to perform the cutting was ''O''(''k''), therefore (again with a sufficiently large choice of ''c'') the amortized running time is constant.
 
Finally, operation '''delete''' can be implemented simply by decreasing the key of the element to be deleted to minus infinity, thus turning it into the minimum of the whole heap. Then we call extract minimum to remove it. The amortized running time of this operation is ''O''(log ''n'').
 
==Proof of degree bounds==
The amortized performance of a Fibonacci heap depends on the degree (number of children) of any tree root being ''O''(log ''n''), where ''n'' is the size of the heap.  Here we show that the size of the (sub)tree rooted at any node ''x'' of degree ''d'' in the heap must have size at least ''F<sub>d</sub>''<sub>+2</sub>, where ''F<sub>k</sub>'' is the ''k''th [[Fibonacci number]]The degree bound follows from this and the fact (easily proved by induction) that <math>F_{d+2} \ge \varphi^d</math> for all integers <math>d\ge 0</math>, where <math>\varphi = (1+\sqrt 5)/2 \doteq 1.618</math>.  (We then have <math>n \ge F_{d+2} \ge \varphi^d</math>, and taking the log to base <math>\varphi</math> of both sides gives <math>d\le \log_{\varphi} n</math> as required.)
 
Consider any node ''x'' somewhere in the heap (''x'' need not be the root of one of the main trees). Define '''size'''(''x'') to be the size of the tree rooted at ''x'' (the number of descendants of ''x'', including ''x'' itself). We prove by induction on the height of ''x'' (the length of a longest simple path from ''x'' to a descendant leaf), that '''size'''(''x'')&nbsp;≥&nbsp;''F<sub>d</sub>''<sub>+2</sub>, where ''d'' is the degree of ''x''.
 
'''Base case:''' If ''x'' has height 0, then ''d''&nbsp;=&nbsp;0, and '''size'''(''x'')&nbsp;=&nbsp;1&nbsp;=&nbsp;''F''<sub>2</sub>.
 
'''Inductive case:'''  Suppose ''x'' has positive height and degree ''d''&gt;0.  Let ''y''<sub>1</sub>, ''y''<sub>2</sub>, ..., ''y<sub>d</sub>'' be the children of ''x'', indexed in order of the times they were most recently made children of ''x'' (''y''<sub>1</sub> being the earliest and ''y<sub>d</sub>'' the latest), and let ''c''<sub>1</sub>, ''c''<sub>2</sub>, ..., ''c<sub>d</sub>'' be their respective degrees.  We '''claim''' that ''c<sub>i</sub>''&nbsp;≥&nbsp;''i''-2 for each ''i'' with 2≤''i''≤''d'': Just before ''y<sub>i</sub>'' was made a child of ''x'', ''y''<sub>1</sub>,...,''y<sub>i</sub>''<sub>−1</sub> were already children of ''x'', and so ''x'' had degree at least ''i''−1 at that time.  Since trees are combined only when the degrees of their roots are equal, it must have been that ''y<sub>i</sub>'' also had degree at least ''i''-1 at the time it became a child of ''x''. From that time to the present, ''y<sub>i</sub>'' can only have lost at most one child (as guaranteed by the marking process), and so its current degree ''c<sub>i</sub>'' is at least ''i''−2.  This proves the '''claim'''.
 
Since the heights of all the ''y<sub>i</sub>'' are strictly less than that of ''x'', we can apply the inductive hypothesis to them to get '''size'''(''y<sub>i</sub>'')&nbsp;≥&nbsp;''F<sub>c<sub>i</sub></sub>''<sub>+2</sub>&nbsp;≥&nbsp;''F''<sub>(''i''−2)+2</sub>&nbsp;=&nbsp;''F<sub>i</sub>''. The nodes ''x'' and ''y''<sub>1</sub> each contribute at least 1 to '''size'''(''x''), and so we have
 
<math>\textbf{size}(x) \ge 2 + \sum_{i=2}^d \textbf{size}(y_i) \ge 2 + \sum_{i=2}^d F_i = 1 + \sum_{i=0}^d F_i.</math>
 
A routine induction proves that <math>1 + \sum_{i=0}^d F_i = F_{d+2}</math> for any <math>d\ge 0</math>, which gives the desired lower bound on '''size'''(''x'').
 
==Worst case==
Although the total running time of a sequence of operations starting with an empty structure is bounded by the bounds given above, some (very few) operations in the sequence can take very long to complete (in particular delete and delete minimum have linear running time in the worst case). For this reason Fibonacci heaps and other amortized data structures may not be appropriate for [[real-time computing|real-time systems]]. It is possible to create a data structure which has the same worst case performance as the Fibonacci heap has amortized performance.<ref name="bare_url">{{Citation |id = [[CiteSeerX]]: {{url|1=citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.43.8133|2=10.1.1.43.8133}} |title=Worst-Case Efficient Priority Queues |year=1996 |author=Gerth Stølting Brodal |journal=Proc. 7th ACM-SIAM Symposium on Discrete Algorithms |publisher = [[Society for Industrial and Applied Mathematics]]|pages=52–58 |doi=10.1145/313852.313883 |isbn=0-89871-366-8}}</ref> However the resulting structure (a [[Brodal queue]]) is, in the words of the creator, "quite complicated" and "[not] applicable in practice."
 
==Summary of running times==
{{Contradict-other|Heap (data structure)|there|Contradiction of running times|date=October 2013}}
{{main|Heap (data structure)#Comparison_of_theoretic_bounds_for_variants}}
 
{| class="wikitable"
|-
! ''Common Operations''
! ''Effect''
! Unsorted [[Linked List]]
! Sorted [[Linked List]]
! [[Self-balancing binary search tree]]
! [[Binary heap]]
! [[Binomial heap]]
! Fibonacci heap
! [[Brodal queue]]<ref name="bare_url" />
! [[Pairing heap]]
|-
| insert(data,key)
| Adds ''data'' to the queue, tagged with ''key''
| O(1)
| O(n)
| O(log n)
| O(log n)
| O(log n)
| O(1)
| O(1)
| O(1)<ref name="Fredman And Tarjan">{{cite journal|first1=Michael Lawrence|last1=Fredman|authorlink1=Michael Fredman|first2=Robert E.|last2=Tarjan|authorlink2=Robert Tarjan |title=Fibonacci heaps and their uses in improved network optimization algorithms| url = http://www.cl.cam.ac.uk/~sos22/supervise/dsaa/fib_heaps.pdf | format = PDF |journal=[[Journal of the Association for Computing Machinery]]|volume=34|year=1987|pages=596&ndash;615|ref=harv|doi=10.1145/28869.28874|issue=3}}</ref>
|-
| findMin() -> key,data
| Returns ''key,data'' corresponding to min-value ''key''
| O(n)
| O(1)
| O(log n) or O(1) (**)
| O(1)
| O(log n)<ref name="CLRS_generalized">{{Introduction to Algorithms|edition=1|chapter=}}</ref> or O(1) (**)
| O(1)<ref name="CLRS"/>
| O(1)
| O(1)<ref name="Fredman And Tarjan"/>
|-
| deleteMin()
| Deletes ''data'' corresponding to min-value ''key''
| O(n)
| O(1)
| O(log n)
| O(log n)
| O(log n)
| O(log n)*
| O(log n)
| O(log n)*<ref name="Fredman And Tarjan"/>
|-
| delete(node)
| Deletes ''data'' corresponding to given ''key'', given a pointer to the node being deleted
| O(1)
| O(1)
| O(log n)
| O(log n)
| O(log n)
| O(log n)*
| O(log n)
| O(log n)*<ref name="Fredman And Tarjan"/>
|-
| decreaseKey(node)
| Decreases the ''key'' of a node, given a pointer to the node being modified
| O(1)
| O(n)
| O(log n)
| O(log n)
| O(log n)
| O(1)*
| O(1)
| Unknown but bounded: <math>\Omega(\log\log n), {O(2^{2\sqrt{\log\log n}})}</math>*<ref name="Fredman And Tarjan"/><ref>{{cite journal|last=Pettie|first=Seth|title=Towards a Final Analysis of Pairing Heaps|journal=Max Planck Institut f&uuml;r Informatik|year=2005|url=http://web.eecs.umich.edu/~pettie/papers/focs05.pdf}}</ref>
|-
| merge(heap1,heap2) -> heap3
| Merges two heaps into a third
| O(1)
| O(m + n)
| O(m log(n+m))
| O(m log(n+m))
| O(log (m+n))
| O(1)*
| O(1)
| O(1)<ref name="Fredman And Tarjan"/>
|}
 
(*)Amortized time<br />
(**)With trivial modification to store an additional pointer to the minimum element
 
== References ==
{{reflist}}
 
== External links ==
* [http://www.cs.yorku.ca/~aaw/Jason/FibonacciHeapAnimation.html Java applet simulation of a Fibonacci heap]
* [http://www.mathworks.com/matlabcentral/fileexchange/30072-fibonacci-heap MATLAB implementation of Fibonacci heap]
* [http://www.labri.fr/perso/pelegrin/code/#fibonacci De-recursived and memory efficient C implementation of Fibonacci heap] (free/libre software, [http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html CeCILL-B license])
* [http://ideone.com/9jYnv C++ template Fibonacci heap, with demonstration]
* [http://github.com/evansenter/f_heap Ruby implementation of the Fibonacci heap (with tests)]
* [http://www.cs.princeton.edu/~wayne/cs423/fibonacci/FibonacciHeapAlgorithm.html Pseudocode of the Fibonacci heap algorithm]
* [http://stackoverflow.com/q/504823/194609 Efficient C++ Fibonacci heap]
* [http://stackoverflow.com/q/6273833/194609 Various Java Implementations for Fibonacci heap]
 
{{Data structures}}
 
{{DEFAULTSORT:Fibonacci Heap}}
[[Category:Fibonacci numbers]]
[[Category:Heaps (data structures)]]

Revision as of 18:21, 3 March 2014

With your offense, you might experience Gunboats which can extremely easily shoot at enemy resistance coming from a very long range and Landing Crafts which you must fill up when you train units for example Rifleman, Heavy, Zooka, Warrior and Containers. To your village defenses, you might carry structures like Mortar, Bike Gun, Sniper Tower, Cannon, Flamethrower, Mine, Tank Mine, Boom Cannon and Catapult Launcher to assist you eradicate enemies.

Interweaving social styles form an absolute strong net in which we are all encased. When The Tygers of Pan Tang sang 'It's lonely at specific top. Everybody's undertaking to do you in', these people borrowed much from clash of clans hack tool no survey. A society without clash of clans hack tool no survey is really like a society with no knowledge, in your it is quite smart.

Temperance is essential in almost things, and enjoying clash of clans cheats is no special. Playing for hours on finish seriously is not good-for-you, bodily or mentally. There are some games out here which know this and can include measures to remind in order to take rests. Consider the initiative yourself, although! Place an alarm a person don't play for a lot more an hour right.

Investigate evaluations and see currently the trailers before buying another video game. Makes it one thing you will be looking at before you get in which. These video games aren't low-cost, and also you likely to get nearly as much cash whenever you commerce inside a employed game title which you have purely utilized several times.

The nfl season is here and furthermore going strong, and like many fans we incredibly long for Sunday afternoon when the games begin. If you have held and liked Soul Caliber, you will love this game. The succeeding best is the Microchip Cell which will at random , fill in some pieces. Defeating players like that by any usually means necessary can be that this reason that pushes folks to use Words that has Friends Cheat. If you liked this information and you would certainly like to get more information concerning clash of clans hack android no survey kindly check out the website. Your current app requires you within order to answer 40 questions using varying degrees of hard times.

Had you been aware that some software games are educational tools? If you know a children that likes to participate in video games, educational options are a fantastic tactics to combine learning because of entertaining. The World wide web can connect you with thousands of parents which have similar values and tend to be more than willing of share their reviews combined with notions with you.

Disclaimer: I aggregate the guidance on this commodity by game a lot of CoC and accomplishing some web research. To the best involving my knowledge, is it authentic along with I accept amateur costed all abstracts and data. Nevertheless, it is consistently accessible we accept fabricated a aberration about or which the most important bold has afflicted back publication. Use as part of your very own risk, Do not accommodate virtually any provides. Please get in blow if that you acquisition annihilation amiss.