Optimal Partitioning for Dual Pivot Quicksort
|
|
- Phillip Cobb
- 5 years ago
- Views:
Transcription
1 Optimal Partitioning for Dual Pivot Quicksort Martin Aumüller, Martin Dietzfelbinger Technische Universität Ilmenau, Germany ICALP 2013 Riga, July 12, 2013 M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 1/17
2 Dual Pivot Quicksort p M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 2/17
3 Dual Pivot Quicksort p Choose two pivots p, q with p < q. M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 2/17
4 Dual Pivot Quicksort p q Choose two pivots p, q with p < q. M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 2/17
5 Dual Pivot Quicksort p q Partition, i.e., re-arrange elements via swaps. Partition: p p p... q q q M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 2/17
6 Dual Pivot Quicksort p q Partition, i.e., re-arrange elements via swaps. Partition: p p p... q q q M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 2/17
7 Dual Pivot Quicksort p q Use recursion to sort these three subarrays. Partition: p p p... q q q M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 2/17
8 Dual Pivot Quicksort p q Use recursion to sort these three subarrays. M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 2/17
9 Dual Pivot Quicksort p Done. M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 2/17
10 Dual Pivot Quicksort: History Robert Sedgewick (1975): Analyzed a dual pivot algorithm that makes much more swaps than classical quicksort on average no further investigation. Pascal Hennequin (1991): Thorough analysis of quicksort with k 1 pivots. for k = 2, no improvements found. for k 3, slight improvements, partitioning too complicated. M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 3/17
11 Dual Pivot Quicksort: History Robert Sedgewick (1975): Analyzed a dual pivot algorithm that makes much more swaps than classical quicksort on average no further investigation. Pascal Hennequin (1991): Thorough analysis of quicksort with k 1 pivots. for k = 2, no improvements found. for k 3, slight improvements, partitioning too complicated. So, why bother? M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 3/17
12 Dual Pivot Quicksort: History Robert Sedgewick (1975): Analyzed a dual pivot algorithm that makes much more swaps than classical quicksort on average no further investigation. Pascal Hennequin (1991): Thorough analysis of quicksort with k 1 pivots. for k = 2, no improvements found. for k 3, slight improvements, partitioning too complicated. So, why bother? Java 7 (2009): Classical quicksort is replaced by a dual pivot quicksort variant due to Yaroslavskiy! Rigorous analysis of the core of this algorithm by Wild and Nebel (2012). M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 3/17
13 Classical QS vs. Yaroslavskiy s Dual Pivot QS Classical QS Yaroslavskiy Average Comparison Count 2n ln n (= 1.38n log n) 1.9n ln n Experiments: Yaroslavskiy s algorithm around 10% faster! M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 4/17
14 Classical QS vs. Yaroslavskiy s Dual Pivot QS Classical QS Yaroslavskiy Optimal DP Average Comparison Count 2n ln n (= 1.38n log n) 1.9n ln n?? Experiments: Yaroslavskiy s algorithm around 10% faster! Question: What is possible? M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 4/17
15 Focus of This Talk introduce a model that captures all dual pivot algorithms, give a unified analysis of the average comparison count, identify optimal algorithms. For analysis: input is random permutation of {1,..., n} left-most and right-most elements are the two pivots analyze average sorting cost M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 5/17
16 Reduce Sorting to Partitioning Fact (Hennequin, 1991) For dual pivot quicksort, average partitioning cost of a n + O(1) leads to average sorting cost 6 5a n ln n + O(n) a n + O(1) 6 5 a n ln n + O(n) p p p... q q q M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 6/17
17 Reduce Partitioning to Classification simplification: focus on classifying elements as small, medium, or large independent of this: clever swap strategy to obtain a partition p p p... q q q M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 7/17
18 Model: Decision Tree 3 : p σ µ λ 2 : q 4 : p 2 : p σ µ λ σ µ λ σ µ λ 4 : p 4 : q 4 : q 2 : q 2 : p 2 : p 4 : p 4 : q 4 : p σ σ µ λ λ a2 : σ a3 : σ a4 : σ a2 : σ a2 : µ a2 : λ... a3 : µ a3 : µ a3 : µ... a4 : µ a4 : µ a4 : µ a2 : λ a3 : λ a4 : λ Example: M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 8/17
19 Model: Decision Tree 3 : p σ µ λ 2 : q 4 : p 2 : p σ µ λ σ µ λ σ µ λ 4 : p 4 : q 4 : q 2 : q 2 : p 2 : p 4 : p 4 : q 4 : p σ σ µ λ λ a2 : σ a3 : σ a4 : σ a2 : σ a2 : µ a2 : λ... a3 : µ a3 : µ a3 : µ... a4 : µ a4 : µ a4 : µ a2 : λ a3 : λ a4 : λ Example: p q M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 8/17
20 Model: Decision Tree 3 : p σ µ λ 2 : q 4 : p 2 : p σ µ λ σ µ λ σ µ λ 4 : p 4 : q 4 : q 2 : q 2 : p 2 : p 4 : p 4 : q 4 : p σ σ µ λ λ a2 : σ a3 : σ a4 : σ a2 : σ a2 : µ a2 : λ... a3 : µ a3 : µ a3 : µ... a4 : µ a4 : µ a4 : µ a2 : λ a3 : λ a4 : λ Example: p q M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 8/17
21 Model: Decision Tree +1 3 : p σ µ λ 2 : q 4 : p 2 : p σ µ λ σ µ λ σ µ λ 4 : p 4 : q 4 : q 2 : q 2 : p 2 : p 4 : p 4 : q 4 : p σ σ µ λ λ a2 : σ a3 : σ a4 : σ a2 : σ a2 : µ a2 : λ... a3 : µ a3 : µ a3 : µ... a4 : µ a4 : µ a4 : µ a2 : λ a3 : λ a4 : λ Example: p σ q M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 8/17
22 Model: Decision Tree +1 3 : p σ µ λ 2 : q 4 : p 2 : p σ µ λ σ µ λ σ µ λ 4 : p 4 : q 4 : q 2 : q 2 : p 2 : p 4 : p 4 : q 4 : p σ σ µ λ λ a2 : σ a3 : σ a4 : σ a2 : σ a2 : µ a2 : λ... a3 : µ a3 : µ a3 : µ... a4 : µ a4 : µ a4 : µ a2 : λ a3 : λ a4 : λ Example: p σ q M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 8/17
23 Model: Decision Tree +1 3 : p +2 σ µ λ 2 : q 4 : p 2 : p σ µ λ σ µ λ σ µ λ 4 : p 4 : q 4 : q 2 : q 2 : p 2 : p 4 : p 4 : q 4 : p σ σ µ λ λ a2 : σ a3 : σ a4 : σ a2 : σ a2 : µ a2 : λ... a3 : µ a3 : µ a3 : µ... a4 : µ a4 : µ a4 : µ a2 : λ a3 : λ a4 : λ Example: p µ σ q M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 8/17
24 Model: Decision Tree +1 3 : p +2 σ µ λ 2 : q 4 : p 2 : p σ µ λ σ µ λ σ µ λ 4 : p 4 : q 4 : q 2 : q 2 : p 2 : p 4 : p 4 : q 4 : p σ σ µ λ λ a2 : σ a3 : σ a4 : σ a2 : σ a2 : µ a2 : λ... a3 : µ a3 : µ a3 : µ... a4 : µ a4 : µ a4 : µ a2 : λ a3 : λ a4 : λ Example: p µ σ q M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 8/17
25 Model: Decision Tree +1 3 : p +2 σ µ λ 2 : q 4 : p 2 : p σ µ λ σ µ λ σ µ λ +2 4 : p 4 : q 4 : q 2 : q 2 : p 2 : p 4 : p 4 : q 4 : p σ σ µ λ λ a2 : σ a3 : σ a4 : σ a2 : σ a2 : µ a2 : λ... a3 : µ a3 : µ a3 : µ... a4 : µ a4 : µ a4 : µ a2 : λ a3 : λ a4 : λ Example: p µ σ σ q M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 8/17
26 Model: Decision Tree +1 3 : p +2 σ µ λ 2 : q 4 : p 2 : p σ µ λ σ µ λ σ µ λ +2 4 : p 4 : q 4 : q 2 : q 2 : p 2 : p 4 : p 4 : q 4 : p σ σ µ λ λ a2 : σ a3 : σ a4 : σ a2 : σ a2 : µ a2 : λ... a3 : µ a3 : µ a3 : µ... a4 : µ a4 : µ a4 : µ a2 : λ a3 : λ a4 : λ Example: p µ σ σ q 5 comparisons. M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 8/17
27 Average Cost For the average partitioning/classification cost p n we get: p n 4/3n + average number of additional comparisons For all classification algorithms: each element has to be compared at least once against a pivot each medium element has to be compared to both ( n/3 on avg.) Additional comparisons: a small element is compared with the larger pivot first a large element is compared with the smaller pivot first M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 9/17
28 Strategy 1 (unrealistic) Optimal Strategies σ λ σ µ???? Next element Let s, l denote the number of small/large elements in area (resp.). M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 10/17
29 Strategy 1 (unrealistic) Optimal Strategies σ λ σ µ???? Next element Let s, l denote the number of small/large elements in Classification l > s: compare with larger pivot first. l s: compare with smaller pivot first. area (resp.). M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 10/17
30 Optimal Strategies Strategy 1 (unrealistic) σ λ σ µ???? Strategy 2 σ λ σ µ???? Next element Next element Let s, l denote the number of small/large elements in Classification l > s: compare with larger pivot first. l s: compare with smaller pivot first. area (resp.). M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 10/17
31 Optimal Strategies Strategy 1 (unrealistic) σ λ σ µ???? Next element Strategy 2 σ λ σ µ???? Next element Let s, l denote the number of small/large elements in area (resp.). Classification l > s: compare with larger pivot first. Theorem l s: compare with smaller pivot first. These strategies (turned into dp algorithms) have the minimum possible average sorting cost: 1.8n ln n + O(n). M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 10/17
32 Optimal Strategies Strategy 1 (unrealistic) σ λ σ µ σ λ σ σ Next element Strategy 2 σ λ σ µ σ λ σ σ Next element Let s, l denote the number of small/large elements in area (resp.). Classification l > s: compare with larger pivot first. Theorem l s: compare with smaller pivot first. These strategies (turned into dp algorithms) have the minimum possible average sorting cost: 1.8n ln n + O(n). M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 10/17
33 Optimal Strategies Strategy 1 (unrealistic) σ λ σ µ σ λ σ σ Next element Strategy 2 σ λ σ µ σ λ σ σ Next element Let s, l denote the number of small/large elements in area (resp.). Classification l > s: compare with larger pivot first. Theorem l s: compare with smaller pivot first. These strategies (turned into dp algorithms) have the minimum possible average sorting cost: 1.8n ln n + O(n). M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 10/17
34 Cost of an Arbitrary Decision Tree p n 4/3n + average number of additional comparisons M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 11/17
35 Cost of an Arbitrary Decision Tree p n 4/3n + average number of additional comparisons Central: f q s,l average number of comparisons to the larger pivot first. Lemma For the average partition cost p n of a decision tree T we have: p n = 4 3 n + 1 ( ) ( n ) f q s,l s + (n 2 f q s,l ) l + o(n). 2 (n 2) s+l n 2 dependence introduces an error term of o(n). M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 11/17
36 Cost of an Arbitrary Decision Tree p n 4/3n + average number of additional comparisons Central: f q s,l average number of comparisons to the larger pivot first. Lemma For the average partition cost p n of a decision tree T we have: p n = 4 3 n + 1 ( ) ( n ) f q s,l s + (n 2 f q s,l ) l + o(n). 2 (n 2) s+l n 2 dependence introduces an error term of o(n). M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 11/17
37 Cost of an Arbitrary Decision Tree p n 4/3n + average number of additional comparisons Central: f q s,l average number of comparisons to the larger pivot first. Lemma For the average partition cost p n of a decision tree T we have: p n = 4 3 n + 1 ( ) ( n ) f q s,l s + (n 2 f q s,l ) l + o(n). 2 (n 2) s+l n 2 dependence introduces an error term of o(n). M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 11/17
38 Cost of an Arbitrary Decision Tree p n 4/3n + average number of additional comparisons Central: f q s,l average number of comparisons to the larger pivot first. Lemma For the average partition cost p n of a decision tree T we have: p n = 4 3 n + 1 ( ) ( n ) f q s,l s + (n 2 f q s,l ) l + o(n). 2 (n 2) s+l n 2 dependence introduces an error term of o(n). M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 11/17
39 Cost of an Arbitrary Decision Tree p n 4/3n + average number of additional comparisons Central: f q s,l average number of comparisons to the larger pivot first. Lemma For the average partition cost p n of a decision tree T we have: p n = 4 3 n + 1 ( ) ( n ) f q s,l s + (n 2 f q s,l ) l + o(n). 2 (n 2) s+l n 2 dependence introduces an error term of o(n). Message: f q s,l fully describes the average cost (up to lower order terms.) M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 11/17
40 Cost of an Arbitrary Decision Tree p n 4/3n + average number of additional comparisons Central: f q s,l average number of comparisons to the larger pivot first. Lemma For the average partition cost p n of a decision tree T we have: p n = 4 3 n + 1 ( ) ( n ) f q s,l s + (n 2 f q s,l ) l + o(n). 2 (n 2) s+l n 2 dependence introduces an error term of o(n). Message: f q s,l fully describes the average cost (up to lower order terms.) Example: Yaroslavskiy s algorithm: f q s,l = l. (Leads to: 1.9n ln n.) M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 11/17
41 An Optimal Strategy Using An Oracle Assume: An oracle tells us whether or not l > s for an input: Classification l > s: Compare all elements to the larger pivot first l s: Compare all elements to the smaller pivot first M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 12/17
42 An Optimal Strategy Using An Oracle Assume: An oracle tells us whether or not l > s for an input: Classification l > s: Compare all elements to the larger pivot first l s: Compare all elements to the smaller pivot first Why is it optimal? It minimizes p n = 4 3 n + 1 ( n ) 2 (n 2) s+l n 2 ( f q s,l s + (n 2 f q s,l ) l ) + o(n). M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 12/17
43 An Optimal Strategy Using An Oracle Assume: An oracle tells us whether or not l > s for an input: Classification l > s: Compare all elements to the larger pivot first (f q s,l = n 2). l s: Compare all elements to the smaller pivot first (f q s,l = 0). Why is it optimal? It minimizes p n = 4 3 n + 1 ( n ) 2 (n 2) s+l n 2 ( f q s,l s + (n 2 f q s,l ) l ) + o(n). M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 12/17
44 An Optimal Strategy Using An Oracle Assume: An oracle tells us whether or not l > s for an input: Classification l > s: Compare all elements to the larger pivot first (f q s,l = n 2). l s: Compare all elements to the smaller pivot first (f q s,l = 0). Why is it optimal? It minimizes p n = 4 3 n + 1 ( n ) 2 (n 2) s+l n 2 ( f q s,l s + (n 2 f q s,l ) l ) + o(n). Average sorting cost: 1.8n ln n + o(n ln n) M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 12/17
45 Random Sampling Just look at, e.g., n 2/3 elements. Classification Seen more large elements than small elements? Compare all elements with the larger pivot first, otherwise with the smaller pivot first. proof that this works: standard machinery (tail estimates). difference to oracle estimation is o(n). Average sorting cost: 1.8n ln n + o(n ln n). M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 13/17
46 Experiments Our Method (N) Yaroslavskiy Clever Quicksort time [ms] n/10^6 Sorting integers: About 3% slower than Yaroslavskiy. M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 14/17
47 Experiments (2) Our Method (N) Yaroslavskiy Clever Quicksort 8000 time [ms] n/10^6 Sorting strings: About 2% faster than Yaroslavskiy. M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 15/17
48 Conclusion and Open Questions So far: gave a unified model for dual pivot quicksort algorithms, showed how to calculate the average cost for general algorithms, identified three (but only two algorithmic ) optimal dual pivot quicksort algorithms. Optimum: 1.8n ln n. (QS: 2n ln n, Yaroslavskiy: 1.9n ln n.) Open questions: Cost measure that describes running time behavior accurately? Optimal algorithms for multi-pivot quicksort with k 3 pivots? (Best (known) 3-pivot algorithm: 1.846n ln n comparisons on average) M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 16/17
49 Thanks! Any questions? M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 17/17
50 Theoretical Result Classical QS Yaroslavskiy Our algorithm Average Comparison Count 2n ln n 1.9n ln n 1.8n ln n Average Swap Count n ln n 0.6n ln n n ln n LB: 0.3n ln n M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 16/17
51 p p p... q q?? q current element M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 17/17
52 p p p... q q?? q small element M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 17/17
53 p p p... q q?? q exchange M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 17/17
54 p p p... q?? q q exchange M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 17/17
55 p p p... q?? q q inspect next element M. Aumüller Optimal Partitioning for Dual Pivot Quicksort 17/17
Empirical and Average Case Analysis
Empirical and Average Case Analysis l We have discussed theoretical analysis of algorithms in a number of ways Worst case big O complexities Recurrence relations l What we often want to know is what will
More informationYao s Minimax Principle
Complexity of algorithms The complexity of an algorithm is usually measured with respect to the size of the input, where size may for example refer to the length of a binary word describing the input,
More informationModule 4: Probability
Module 4: Probability 1 / 22 Probability concepts in statistical inference Probability is a way of quantifying uncertainty associated with random events and is the basis for statistical inference. Inference
More informationThe Limiting Distribution for the Number of Symbol Comparisons Used by QuickSort is Nondegenerate (Extended Abstract)
The Limiting Distribution for the Number of Symbol Comparisons Used by QuickSort is Nondegenerate (Extended Abstract) Patrick Bindjeme 1 James Allen Fill 1 1 Department of Applied Mathematics Statistics,
More informationIntroduction to Greedy Algorithms: Huffman Codes
Introduction to Greedy Algorithms: Huffman Codes Yufei Tao ITEE University of Queensland In computer science, one interesting method to design algorithms is to go greedy, namely, keep doing the thing that
More information15-451/651: Design & Analysis of Algorithms October 23, 2018 Lecture #16: Online Algorithms last changed: October 22, 2018
15-451/651: Design & Analysis of Algorithms October 23, 2018 Lecture #16: Online Algorithms last changed: October 22, 2018 Today we ll be looking at finding approximately-optimal solutions for problems
More informationOptimal Incentive Contract with Costly and Flexible Monitoring
Optimal Incentive Contract with Costly and Flexible Monitoring Anqi Li 1 Ming Yang 2 1 Department of Economics, Washington University in St. Louis 2 Fuqua School of Business, Duke University January 2016
More information2 all subsequent nodes. 252 all subsequent nodes. 401 all subsequent nodes. 398 all subsequent nodes. 330 all subsequent nodes
¼ À ÈÌ Ê ½¾ ÈÊÇ Ä ÅË ½µ ½¾º¾¹½ ¾µ ½¾º¾¹ µ ½¾º¾¹ µ ½¾º¾¹ µ ½¾º ¹ µ ½¾º ¹ µ ½¾º ¹¾ µ ½¾º ¹ µ ½¾¹¾ ½¼µ ½¾¹ ½ (1) CLR 12.2-1 Based on the structure of the binary tree, and the procedure of Tree-Search, any
More informationHomework Assignment #3. 1 Demonstrate how mergesort works when sorting the following list of numbers:
CISC 5835 Algorithms for Big Data Fall, 2018 Homework Assignment #3 1 Demonstrate how mergesort works when sorting the following list of numbers: 6 1 4 2 3 8 7 5 2 Given the following array (list), follows
More information15-451/651: Design & Analysis of Algorithms November 9 & 11, 2015 Lecture #19 & #20 last changed: November 10, 2015
15-451/651: Design & Analysis of Algorithms November 9 & 11, 2015 Lecture #19 & #20 last changed: November 10, 2015 Last time we looked at algorithms for finding approximately-optimal solutions for NP-hard
More informationDiscrete Mathematics for CS Spring 2008 David Wagner Final Exam
CS 70 Discrete Mathematics for CS Spring 2008 David Wagner Final Exam PRINT your name:, (last) SIGN your name: (first) PRINT your Unix account login: Your section time (e.g., Tue 3pm): Name of the person
More informationDeriving the Black-Scholes Equation and Basic Mathematical Finance
Deriving the Black-Scholes Equation and Basic Mathematical Finance Nikita Filippov June, 7 Introduction In the 97 s Fischer Black and Myron Scholes published a model which would attempt to tackle the issue
More informationSmoothed Analysis of Binary Search Trees
Smoothed Analysis of Binary Search Trees Bodo Manthey and Rüdiger Reischuk Universität zu Lübeck, Institut für Theoretische Informatik Ratzeburger Allee 160, 23538 Lübeck, Germany manthey/reischuk@tcs.uni-luebeck.de
More informationApproximate Revenue Maximization with Multiple Items
Approximate Revenue Maximization with Multiple Items Nir Shabbat - 05305311 December 5, 2012 Introduction The paper I read is called Approximate Revenue Maximization with Multiple Items by Sergiu Hart
More informationContinuing Education Course #287 Engineering Methods in Microsoft Excel Part 2: Applied Optimization
1 of 6 Continuing Education Course #287 Engineering Methods in Microsoft Excel Part 2: Applied Optimization 1. Which of the following is NOT an element of an optimization formulation? a. Objective function
More informationSelf-organized criticality on the stock market
Prague, January 5th, 2014. Some classical ecomomic theory In classical economic theory, the price of a commodity is determined by demand and supply. Let D(p) (resp. S(p)) be the total demand (resp. supply)
More information6.1 Binomial Theorem
Unit 6 Probability AFM Valentine 6.1 Binomial Theorem Objective: I will be able to read and evaluate binomial coefficients. I will be able to expand binomials using binomial theorem. Vocabulary Binomial
More informationOptimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing
Optimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing Prof. Chuan-Ju Wang Department of Computer Science University of Taipei Joint work with Prof. Ming-Yang Kao March 28, 2014
More informationNOTES ON FIBONACCI TREES AND THEIR OPTIMALITY* YASUICHI HORIBE INTRODUCTION 1. FIBONACCI TREES
0#0# NOTES ON FIBONACCI TREES AND THEIR OPTIMALITY* YASUICHI HORIBE Shizuoka University, Hamamatsu, 432, Japan (Submitted February 1982) INTRODUCTION Continuing a previous paper [3], some new observations
More informationProbability without Measure!
Probability without Measure! Mark Saroufim University of California San Diego msaroufi@cs.ucsd.edu February 18, 2014 Mark Saroufim (UCSD) It s only a Game! February 18, 2014 1 / 25 Overview 1 History of
More informationTHE TRAVELING SALESMAN PROBLEM FOR MOVING POINTS ON A LINE
THE TRAVELING SALESMAN PROBLEM FOR MOVING POINTS ON A LINE GÜNTER ROTE Abstract. A salesperson wants to visit each of n objects that move on a line at given constant speeds in the shortest possible time,
More informationApproximate Composite Minimization: Convergence Rates and Examples
ISMP 2018 - Bordeaux Approximate Composite Minimization: Convergence Rates and S. Praneeth Karimireddy, Sebastian U. Stich, Martin Jaggi MLO Lab, EPFL, Switzerland sebastian.stich@epfl.ch July 4, 2018
More information1 Overview. 2 The Gradient Descent Algorithm. AM 221: Advanced Optimization Spring 2016
AM 22: Advanced Optimization Spring 206 Prof. Yaron Singer Lecture 9 February 24th Overview In the previous lecture we reviewed results from multivariate calculus in preparation for our journey into convex
More informationOptimal Investment for Worst-Case Crash Scenarios
Optimal Investment for Worst-Case Crash Scenarios A Martingale Approach Frank Thomas Seifried Department of Mathematics, University of Kaiserslautern June 23, 2010 (Bachelier 2010) Worst-Case Portfolio
More informationElif Özge Özdamar T Reinforcement Learning - Theory and Applications February 14, 2006
On the convergence of Q-learning Elif Özge Özdamar elif.ozdamar@helsinki.fi T-61.6020 Reinforcement Learning - Theory and Applications February 14, 2006 the covergence of stochastic iterative algorithms
More informationLaws of probabilities in efficient markets
Laws of probabilities in efficient markets Vladimir Vovk Department of Computer Science Royal Holloway, University of London Fifth Workshop on Game-Theoretic Probability and Related Topics 15 November
More informationBest-Reply Sets. Jonathan Weinstein Washington University in St. Louis. This version: May 2015
Best-Reply Sets Jonathan Weinstein Washington University in St. Louis This version: May 2015 Introduction The best-reply correspondence of a game the mapping from beliefs over one s opponents actions to
More informationMarshall and Hicks Understanding the Ordinary and Compensated Demand
Marshall and Hicks Understanding the Ordinary and Compensated Demand K.J. Wainwright March 3, 213 UTILITY MAXIMIZATION AND THE DEMAND FUNCTIONS Consider a consumer with the utility function =, who faces
More informationMechanisms for Matching Markets with Budgets
Mechanisms for Matching Markets with Budgets Paul Dütting Stanford LSE Joint work with Monika Henzinger and Ingmar Weber Seminar on Discrete Mathematics and Game Theory London School of Economics July
More informationLossy compression of permutations
Lossy compression of permutations The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation As Published Publisher Wang, Da, Arya Mazumdar,
More informationHow do Variance Swaps Shape the Smile?
How do Variance Swaps Shape the Smile? A Summary of Arbitrage Restrictions and Smile Asymptotics Vimal Raval Imperial College London & UBS Investment Bank www2.imperial.ac.uk/ vr402 Joint Work with Mark
More informationIEOR E4602: Quantitative Risk Management
IEOR E4602: Quantitative Risk Management Basic Concepts and Techniques of Risk Management Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationTutorial 6. Sampling Distribution. ENGG2450A Tutors. 27 February The Chinese University of Hong Kong 1/6
Tutorial 6 Sampling Distribution ENGG2450A Tutors The Chinese University of Hong Kong 27 February 2017 1/6 Random Sample and Sampling Distribution 2/6 Random sample Consider a random variable X with distribution
More informationThe Black-Scholes Model
IEOR E4706: Foundations of Financial Engineering c 2016 by Martin Haugh The Black-Scholes Model In these notes we will use Itô s Lemma and a replicating argument to derive the famous Black-Scholes formula
More informationEfficiency and Herd Behavior in a Signalling Market. Jeffrey Gao
Efficiency and Herd Behavior in a Signalling Market Jeffrey Gao ABSTRACT This paper extends a model of herd behavior developed by Bikhchandani and Sharma (000) to establish conditions for varying levels
More informationGame Theory Fall 2003
Game Theory Fall 2003 Problem Set 5 [1] Consider an infinitely repeated game with a finite number of actions for each player and a common discount factor δ. Prove that if δ is close enough to zero then
More informationChapter ML:III. III. Decision Trees. Decision Trees Basics Impurity Functions Decision Tree Algorithms Decision Tree Pruning
Chapter ML:III III. Decision Trees Decision Trees Basics Impurity Functions Decision Tree Algorithms Decision Tree Pruning ML:III-93 Decision Trees STEIN/LETTMANN 2005-2017 Overfitting Definition 10 (Overfitting)
More informationLecture 23: April 10
CS271 Randomness & Computation Spring 2018 Instructor: Alistair Sinclair Lecture 23: April 10 Disclaimer: These notes have not been subjected to the usual scrutiny accorded to formal publications. They
More informationLecture Neyman Allocation vs Proportional Allocation and Stratified Random Sampling vs Simple Random Sampling
Math 408 - Mathematical Statistics Lecture 20-21. Neyman Allocation vs Proportional Allocation and Stratified Random Sampling vs Simple Random Sampling March 8-13, 2013 Konstantin Zuev (USC) Math 408,
More informationDRAFT. 1 exercise in state (S, t), π(s, t) = 0 do not exercise in state (S, t) Review of the Risk Neutral Stock Dynamics
Chapter 12 American Put Option Recall that the American option has strike K and maturity T and gives the holder the right to exercise at any time in [0, T ]. The American option is not straightforward
More informationOn the Optimality of a Family of Binary Trees Techical Report TR
On the Optimality of a Family of Binary Trees Techical Report TR-011101-1 Dana Vrajitoru and William Knight Indiana University South Bend Department of Computer and Information Sciences Abstract In this
More informationPricing and risk of financial products
and risk of financial products Prof. Dr. Christian Weiß Riga, 27.02.2018 Observations AAA bonds are typically regarded as risk-free investment. Only examples: Government bonds of Australia, Canada, Denmark,
More informationSharpe Ratio over investment Horizon
Sharpe Ratio over investment Horizon Ziemowit Bednarek, Pratish Patel and Cyrus Ramezani December 8, 2014 ABSTRACT Both building blocks of the Sharpe ratio the expected return and the expected volatility
More informationImprovement and Efficient Implementation of a Lattice-based Signature scheme
Improvement and Efficient Implementation of a Lattice-based Signature scheme, Johannes Buchmann Technische Universität Darmstadt TU Darmstadt August 2013 Lattice-based Signatures1 Outline Introduction
More informationChapter 16. Binary Search Trees (BSTs)
Chapter 16 Binary Search Trees (BSTs) Search trees are tree-based data structures that can be used to store and search for items that satisfy a total order. There are many types of search trees designed
More informationLecture 22. Survey Sampling: an Overview
Math 408 - Mathematical Statistics Lecture 22. Survey Sampling: an Overview March 25, 2013 Konstantin Zuev (USC) Math 408, Lecture 22 March 25, 2013 1 / 16 Survey Sampling: What and Why In surveys sampling
More informationMore On λ κ closed sets in generalized topological spaces
Journal of Algorithms and Computation journal homepage: http://jac.ut.ac.ir More On λ κ closed sets in generalized topological spaces R. Jamunarani, 1, P. Jeyanthi 2 and M. Velrajan 3 1,2 Research Center,
More informationFrom Discrete Time to Continuous Time Modeling
From Discrete Time to Continuous Time Modeling Prof. S. Jaimungal, Department of Statistics, University of Toronto 2004 Arrow-Debreu Securities 2004 Prof. S. Jaimungal 2 Consider a simple one-period economy
More informationThe Multistep Binomial Model
Lecture 10 The Multistep Binomial Model Reminder: Mid Term Test Friday 9th March - 12pm Examples Sheet 1 4 (not qu 3 or qu 5 on sheet 4) Lectures 1-9 10.1 A Discrete Model for Stock Price Reminder: The
More informationInformation aggregation for timing decision making.
MPRA Munich Personal RePEc Archive Information aggregation for timing decision making. Esteban Colla De-Robertis Universidad Panamericana - Campus México, Escuela de Ciencias Económicas y Empresariales
More informationApplications of Quantum Annealing in Computational Finance. Dr. Phil Goddard Head of Research, 1QBit D-Wave User Conference, Santa Fe, Sept.
Applications of Quantum Annealing in Computational Finance Dr. Phil Goddard Head of Research, 1QBit D-Wave User Conference, Santa Fe, Sept. 2016 Outline Where s my Babel Fish? Quantum-Ready Applications
More informationLecture 17: More on Markov Decision Processes. Reinforcement learning
Lecture 17: More on Markov Decision Processes. Reinforcement learning Learning a model: maximum likelihood Learning a value function directly Monte Carlo Temporal-difference (TD) learning COMP-424, Lecture
More informationStructural Induction
Structural Induction Jason Filippou CMSC250 @ UMCP 07-05-2016 Jason Filippou (CMSC250 @ UMCP) Structural Induction 07-05-2016 1 / 26 Outline 1 Recursively defined structures 2 Proofs Binary Trees Jason
More informationCSE 417 Dynamic Programming (pt 2) Look at the Last Element
CSE 417 Dynamic Programming (pt 2) Look at the Last Element Reminders > HW4 is due on Friday start early! if you run into problems loading data (date parsing), try running java with Duser.country=US Duser.language=en
More informationBernstein Bound is Tight
Bernstein Bound is Tight Repairing Luykx-Preneel Optimal Forgeries Mridul Nandi Indian Statistical Institute, Kolkata CRYPTO 2018 Wegman-Carter-Shoup (WCS) MAC M H κ N E K T Nonce based Authenticator Initial
More informationELEMENTS OF MONTE CARLO SIMULATION
APPENDIX B ELEMENTS OF MONTE CARLO SIMULATION B. GENERAL CONCEPT The basic idea of Monte Carlo simulation is to create a series of experimental samples using a random number sequence. According to the
More informationOptimizing Portfolios
Optimizing Portfolios An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2010 Introduction Investors may wish to adjust the allocation of financial resources including a mixture
More informationOn Existence of Equilibria. Bayesian Allocation-Mechanisms
On Existence of Equilibria in Bayesian Allocation Mechanisms Northwestern University April 23, 2014 Bayesian Allocation Mechanisms In allocation mechanisms, agents choose messages. The messages determine
More informationDynamic Pricing in Ridesharing Platforms
Dynamic Pricing in Ridesharing Platforms A Queueing Approach Sid Banerjee Ramesh Johari Carlos Riquelme Cornell Stanford Stanford rjohari@stanford.edu With thanks to Chris Pouliot, Chris Sholley, and Lyft
More informationStatistical Models of Word Frequency and Other Count Data
Statistical Models of Word Frequency and Other Count Data Martin Jansche 2004-02-12 Motivation Item counts are commonly used in NLP as independent variables in many applications: information retrieval,
More informationOnline Algorithms SS 2013
Faculty of Computer Science, Electrical Engineering and Mathematics Algorithms and Complexity research group Jun.-Prof. Dr. Alexander Skopalik Online Algorithms SS 2013 Summary of the lecture by Vanessa
More informationComputational Independence
Computational Independence Björn Fay mail@bfay.de December 20, 2014 Abstract We will introduce different notions of independence, especially computational independence (or more precise independence by
More informationMultinomial Coefficient : A Generalization of the Binomial Coefficient
Multinomial Coefficient : A Generalization of the Binomial Coefficient Example: A team plays 16 games in a season. At the end of the season, the team has 8 wins, 3 ties and 5 losses. How many different
More informationMultiproduct-Firm Oligopoly: An Aggregative Games Approach
Multiproduct-Firm Oligopoly: An Aggregative Games Approach Volker Nocke 1 Nicolas Schutz 2 1 UCLA 2 University of Mannheim ASSA ES Meetings, Philadephia, 2018 Nocke and Schutz (UCLA &Mannheim) Multiproduct-Firm
More informationChange of Measure (Cameron-Martin-Girsanov Theorem)
Change of Measure Cameron-Martin-Girsanov Theorem Radon-Nikodym derivative: Taking again our intuition from the discrete world, we know that, in the context of option pricing, we need to price the claim
More informationThe Stigler-Luckock model with market makers
Prague, January 7th, 2017. Order book Nowadays, demand and supply is often realized by electronic trading systems storing the information in databases. Traders with access to these databases quote their
More informationHedging Credit Derivatives in Intensity Based Models
Hedging Credit Derivatives in Intensity Based Models PETER CARR Head of Quantitative Financial Research, Bloomberg LP, New York Director of the Masters Program in Math Finance, Courant Institute, NYU Stanford
More informationMATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS
MATH307/37 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS School of Mathematics and Statistics Semester, 04 Tutorial problems should be used to test your mathematical skills and understanding of the lecture material.
More informationOn the smallest abundant number not divisible by the first k primes
On the smallest abundant number not divisible by the first k rimes Douglas E. Iannucci Abstract We say a ositive integer n is abundant if σ(n) > 2n, where σ(n) denotes the sum of the ositive divisors of
More informationDistribution of the Sample Mean
Distribution of the Sample Mean MATH 130, Elements of Statistics I J. Robert Buchanan Department of Mathematics Fall 2018 Experiment (1 of 3) Suppose we have the following population : 4 8 1 2 3 4 9 1
More informationDiscrete Random Variables
Discrete Random Variables MATH 130, Elements of Statistics I J. Robert Buchanan Department of Mathematics Fall 2018 Objectives During this lesson we will learn to: distinguish between discrete and continuous
More informationDiscrete Random Variables
Discrete Random Variables MATH 130, Elements of Statistics I J. Robert Buchanan Department of Mathematics Fall 2017 Objectives During this lesson we will learn to: distinguish between discrete and continuous
More informationEquivalence Nucleolus for Partition Function Games
Equivalence Nucleolus for Partition Function Games Rajeev R Tripathi and R K Amit Department of Management Studies Indian Institute of Technology Madras, Chennai 600036 Abstract In coalitional game theory,
More informationPricing Volatility Derivatives with General Risk Functions. Alejandro Balbás University Carlos III of Madrid
Pricing Volatility Derivatives with General Risk Functions Alejandro Balbás University Carlos III of Madrid alejandro.balbas@uc3m.es Content Introduction. Describing volatility derivatives. Pricing and
More informationApproximation Algorithms for Stochastic Inventory Control Models
Approximation Algorithms for Stochastic Inventory Control Models Retsef Levi Martin Pal Robin Roundy David B. Shmoys Abstract We consider stochastic control inventory models in which the goal is to coordinate
More informationNotes on the symmetric group
Notes on the symmetric group 1 Computations in the symmetric group Recall that, given a set X, the set S X of all bijections from X to itself (or, more briefly, permutations of X) is group under function
More informationCSE 417 Algorithms. Huffman Codes: An Optimal Data Compression Method
CSE 417 Algorithms Huffman Codes: An Optimal Data Compression Method 1 Compression Example 100k file, 6 letter alphabet: a 45% b 13% c 12% d 16% e 9% f 5% File Size: ASCII, 8 bits/char: 800kbits 2 3 >
More informationLIBOR models, multi-curve extensions, and the pricing of callable structured derivatives
Weierstrass Institute for Applied Analysis and Stochastics LIBOR models, multi-curve extensions, and the pricing of callable structured derivatives John Schoenmakers 9th Summer School in Mathematical Finance
More informationRandom Variables and Applications OPRE 6301
Random Variables and Applications OPRE 6301 Random Variables... As noted earlier, variability is omnipresent in the business world. To model variability probabilistically, we need the concept of a random
More informationSingle-Parameter Mechanisms
Algorithmic Game Theory, Summer 25 Single-Parameter Mechanisms Lecture 9 (6 pages) Instructor: Xiaohui Bei In the previous lecture, we learned basic concepts about mechanism design. The goal in this area
More informationarxiv: v1 [cs.dm] 4 Jan 2012
COPS AND INVISIBLE ROBBERS: THE COST OF DRUNKENNESS ATHANASIOS KEHAGIAS, DIETER MITSCHE, AND PAWE L PRA LAT arxiv:1201.0946v1 [cs.dm] 4 Jan 2012 Abstract. We examine a version of the Cops and Robber (CR)
More informationLecture l(x) 1. (1) x X
Lecture 14 Agenda for the lecture Kraft s inequality Shannon codes The relation H(X) L u (X) = L p (X) H(X) + 1 14.1 Kraft s inequality While the definition of prefix-free codes is intuitively clear, we
More informationThe investment game in incomplete markets.
The investment game in incomplete markets. M. R. Grasselli Mathematics and Statistics McMaster University RIO 27 Buzios, October 24, 27 Successes and imitations of Real Options Real options accurately
More informationLecture 5 January 30
EE 223: Stochastic Estimation and Control Spring 2007 Lecture 5 January 30 Lecturer: Venkat Anantharam Scribe: aryam Kamgarpour 5.1 Secretary Problem The problem set-up is explained in Lecture 4. We review
More informationMulti-Armed Bandit, Dynamic Environments and Meta-Bandits
Multi-Armed Bandit, Dynamic Environments and Meta-Bandits C. Hartland, S. Gelly, N. Baskiotis, O. Teytaud and M. Sebag Lab. of Computer Science CNRS INRIA Université Paris-Sud, Orsay, France Abstract This
More informationCSE 21 Winter 2016 Homework 6 Due: Wednesday, May 11, 2016 at 11:59pm. Instructions
CSE 1 Winter 016 Homework 6 Due: Wednesday, May 11, 016 at 11:59pm Instructions Homework should be done in groups of one to three people. You are free to change group members at any time throughout the
More informationPricing theory of financial derivatives
Pricing theory of financial derivatives One-period securities model S denotes the price process {S(t) : t = 0, 1}, where S(t) = (S 1 (t) S 2 (t) S M (t)). Here, M is the number of securities. At t = 1,
More informationSmoothed Analysis of the Height of Binary Search Trees
Electronic Colloquium on Computational Complexity, Report No. 63 (2005) Smoothed Analysis of the Height of Binary Search Trees Bodo Manthey Rüdiger Reischuk Universität zu Lübeck Institut für Theoretische
More informationAn orderly algorithm to enumerate finite (semi)modular lattices
An orderly algorithm to enumerate finite (semi)modular lattices BLAST 23 Chapman University October 6, 23 Outline The original algorithm: Generating all finite lattices Generating modular and semimodular
More informationLog-linear Dynamics and Local Potential
Log-linear Dynamics and Local Potential Daijiro Okada and Olivier Tercieux [This version: November 28, 2008] Abstract We show that local potential maximizer ([15]) with constant weights is stochastically
More informationWorst-Case Value-at-Risk of Non-Linear Portfolios
Worst-Case Value-at-Risk of Non-Linear Portfolios Steve Zymler Daniel Kuhn Berç Rustem Department of Computing Imperial College London Portfolio Optimization Consider a market consisting of m assets. Optimal
More informationSmoothed Analysis of Binary Search Trees and Quicksort Under Additive Noise
Smoothed Analysis of Binary Search Trees a Quicksort Uer Additive Noise Bodo Manthey 1 a Till Tantau 2 1 Saarla University, Computer Science Postfach 151150, 66041 Saarbrücken, Germany manthey@cs.uni-sb.de
More informationPhilipp Moritz Lücke
Σ 1 -partition properties Philipp Moritz Lücke Mathematisches Institut Rheinische Friedrich-Wilhelms-Universität Bonn http://www.math.uni-bonn.de/people/pluecke/ Logic & Set Theory Seminar Bristol, 14.02.2017
More informationE&G, Ch. 8: Multi-Index Models & Grouping Techniques I. Multi-Index Models.
1 E&G, Ch. 8: Multi-Index Models & Grouping Techniques I. Multi-Index Models. A. The General Multi-Index Model: R i = a i + b i1 I 1 + b i2 I 2 + + b il I L + c i Explanation: 1. Let I 1 = R m ; I 2 =
More informationDiscrete Ziggurat: A Time-Memory Trade-off for Sampling from a Gaussian Distribution over the Integers
Discrete Ziggurat: A Time-Memory Trade-off for Sampling from a Gaussian Distribution over the Integers Johannes Buchmann, Daniel Cabarcas, Florian Göpfert, Andreas Hülsing, Patrick Weiden Technische Universität
More informationCONVERGENCE OF OPTION REWARDS FOR MARKOV TYPE PRICE PROCESSES MODULATED BY STOCHASTIC INDICES
CONVERGENCE OF OPTION REWARDS FOR MARKOV TYPE PRICE PROCESSES MODULATED BY STOCHASTIC INDICES D. S. SILVESTROV, H. JÖNSSON, AND F. STENBERG Abstract. A general price process represented by a two-component
More informationThe Cascade Auction A Mechanism For Deterring Collusion In Auctions
The Cascade Auction A Mechanism For Deterring Collusion In Auctions Uriel Feige Weizmann Institute Gil Kalai Hebrew University and Microsoft Research Moshe Tennenholtz Technion and Microsoft Research Abstract
More informationFE670 Algorithmic Trading Strategies. Stevens Institute of Technology
FE670 Algorithmic Trading Strategies Lecture 4. Cross-Sectional Models and Trading Strategies Steve Yang Stevens Institute of Technology 09/26/2013 Outline 1 Cross-Sectional Methods for Evaluation of Factor
More informationOption Pricing Models. c 2013 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 205
Option Pricing Models c 2013 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 205 If the world of sense does not fit mathematics, so much the worse for the world of sense. Bertrand Russell (1872 1970)
More informationThe Binomial Probability Distribution
The Binomial Probability Distribution MATH 130, Elements of Statistics I J. Robert Buchanan Department of Mathematics Fall 2017 Objectives After this lesson we will be able to: determine whether a probability
More information