Analysis - CS @ Purdue

Teams
Task: answer questions to get team points


Score kept during semester
Periodic immediate rewards
Divide room into 4 teams




rearrange now if you like,
sit same place in later classes,
double agents allowed from my point of view, and
citizen-exchange allowed from my point of view.
1
Teams
Choose Team Captain
Choose Team Names




A: Jason, “The Big O’s”
B: Andy, “Team Rocket”
C: Andrew, “Team C”
D: Victoria, “Algorithm Maniacs”
2
First Possible Team Point
You have 43 fish in a pond.
17 fish drown.
[TEAMS] How many fish are still alive?
43
(was told to me by my 7-year old son)
3
Analysis of Algorithms
Input
Algorithm
Output
An algorithm is a step-by-step procedure for
solving a problem in a finite amount of time.
Running Time


best case
average case
worst case
120
100
Running Time
Most algorithms transform
input objects into output
objects.
The running time of an
algorithm typically grows
with the input size.
Average case time is often
difficult to determine.
We focus on the worst case
running time.
80
60
40
20
Easier to analyze
Crucial to applications such as
games, finance and robotics
Analysis of Algorithms
0
1000
2000
3000
4000
Input Size
5
Experimental Studies
9000
8000
7000
Time (ms)
Write a program
implementing the
algorithm
Run the program with
inputs of varying size and
composition
Use a function, like the
built-in clock() function, to
get an accurate measure
of the actual running time
Plot the results
6000
5000
4000
3000
2000
1000
0
0
50
100
Input Size
Analysis of Algorithms
6
Limitations of Experiments
It is necessary to implement the
algorithm, which may be difficult
Results may not be indicative of the
running time on other inputs not included
in the experiment.
In order to compare two algorithms, the
same hardware and software
environments must be used
Analysis of Algorithms
7
Theoretical Analysis
Uses a high-level description of the
algorithm instead of an implementation
Characterizes running time as a
function of the input size, n.
Takes into account all possible inputs
Allows us to evaluate the speed of an
algorithm independent of the
hardware/software environment
Analysis of Algorithms
8
Pseudocode
Example: find max
High-level description
element of an array
of an algorithm
More structured than Algorithm arrayMax(A, n)
English prose
Input array A of n integers
Less detailed than a
Output maximum element of A
program
currentMax  A[0]
Preferred notation for
for i  1 to n  1 do
describing algorithms
if A[i]  currentMax then
Hides program design
currentMax  A[i]
issues
return currentMax
Analysis of Algorithms
9
Pseudocode Details
Method/Function call
Control flow





if … then … [else …]
while … do …
repeat … until …
for … do …
Indentation replaces braces
Method declaration
Algorithm method (arg [, arg…])
Input …
Output …
var.method (arg [, arg…])
Return value
return expression
Expressions
 Assignment
(like  in C++)
 Equality testing
(like  in C++)
n2 Superscripts and other
mathematical
formatting allowed
Analysis of Algorithms
10
The Random Access Machine
(RAM) Model
A CPU
An potentially unbounded
bank of memory cells,
each of which can hold an
arbitrary number or
character
0
2
1
Memory cells are numbered and accessing
any cell in memory takes unit time.
Analysis of Algorithms
11
Primitive Operations
Basic computations
performed by an algorithm
Identifiable in pseudocode
Largely independent from the
programming language
Exact definition not important
(we will see why later)
Assumed to take a constant
amount of time in the RAM
model
Analysis of Algorithms
Examples:





Evaluating an
expression
Assigning a value
to a variable
Indexing into an
array
Calling a method
Returning from a
method
12
Counting Primitive
Operations
By inspecting the pseudocode, we can determine the
maximum number of primitive operations executed by
an algorithm, as a function of the input size
Algorithm arrayMax(A, n)
currentMax  A[0]
for i  1 to n  1 do
if A[i]  currentMax then
currentMax  A[i]
{ increment counter i }
return currentMax
# operations
2
2+n
2(n  1)
2(n  1)
2(n  1)
1
Total
Analysis of Algorithms
7n  1
13
Estimating Running Time
Algorithm arrayMax executes 7n  1 primitive
operations in the worst case.
Define:
a = Time taken by the fastest primitive operation
b = Time taken by the slowest primitive operation
Let T(n) be worst-case time of arrayMax. Then
a (7n  1)  T(n)  b(7n  1)
Hence, the running time T(n) is bounded by two
linear functions
Analysis of Algorithms
14
Power Law
Log-log plot. Plot running time T (N) vs. input
size N using log-log scale.
Power Law
Log-log plot. Plot running time T (N) vs. input size N using loglog scale.
Let T (N) = a N b, where a = 2 c
Thus lg(T (N)) = b lg N + c
After regression:
b = 2.999
c = -33.2103
Regression. Fit straight line through data points: a N b.
Hypothesis. The running time is about 1.006 × 10 –10 × N 2.999
seconds.
Growth Rate of Running Time
Changing the hardware/ software
environment


Affects T(n) by a constant factor
But does not alter the growth rate of T(n)
The linear growth rate of the running
time T(n) is an intrinsic property of
algorithm arrayMax
Analysis of Algorithms
17
Growth Rates



Linear  n
Quadratic  n2
Cubic  n3
T (n )
Growth rates of
functions:
In a log-log chart,
the slope of the line
corresponds to the
growth rate of the
function
1E+30
1E+28
1E+26
1E+24
1E+22
1E+20
1E+18
1E+16
1E+14
1E+12
1E+10
1E+8
1E+6
1E+4
1E+2
1E+0
1E+0
Cubic
Quadratic
Linear
1E+2
Analysis of Algorithms
1E+4
1E+6
1E+8
1E+10
n
18
Practical implications of order-of-growth
Practical implications of order-of-growth
Practical implications of order-of-growth
Practical implications of order-of-growth
Practical implications of order-of-growth
Practical implications of order-of-growth
Practical implications of order-of-growth
Constant Factors
Quadratic
Quadratic
Linear
Linear
T (n )
1E+26
The growth rate is 1E+24
1E+22
not affected by
1E+20
 constant factors or
1E+18
1E+16
 lower-order terms
1E+14
1E+12
Examples
1E+10
 102n + 105 is a linear
1E+8
function
1E+6
1E+4
 105n2 + 108n is a
1E+2
quadratic function
1E+0
1E+0
1E+2
1E+4
1E+6
1E+8
1E+10
n
Analysis of Algorithms
26
Big-Oh Notation
10,000
Given functions f(n) and
g(n), we say that f(n) is
O(g(n)) if there are
1,000
positive constants
c and n0 such that
100
f(n)  cg(n) for n  n0
Example: 2n + 10 is O(n)
For what c and n is cg(n) 10
 f(n)?




2n + 10  cn
(c  2) n  10
n  10/(c  2)
Pick c  3 and n0  10
3n
2n+10
n
1
1
Analysis of Algorithms
10
100
1,000
n
27
Big-Oh Example
1,000,000
n^2
Example: the function
100,000
n2 is not O(n)



n2  cn
nc
The above inequality
cannot be satisfied
since c must be a
constant
100n
10n
n
10,000
1,000
100
10
1
1
Analysis of Algorithms
10
n
100
1,000
28
More Big-Oh Examples
7n-2
e.g. Is there a c and n0 such that c•n  7n-2?
What is the big-Oh?
7n-2 is O(n)
[TEAMS] What are the values for c and n0 ?
Need c > 0 and n0  1 such that 7n-2  c•n for n  n0
This is true for c = 7 and n0 = 1
Analysis of Algorithms
29
More Big-Oh Examples
 2n2+16
2n2+16 is O(n2)
need c > 0 and n0  1 such that 2n2+16  c•n2 for n  n0
this is true for c = 3 and n0 = 4
 3n3 + 2n2 + 9
3n3 + 2n2 + 9 is O(n3)
need c > 0 and n0  1 such that 3n3 + 2n2 + 9  c•n3 for n  n0
this is true for c = 4 and n0 = 3
Analysis of Algorithms
30
Big-Oh and Growth Rate
The big-Oh notation gives an upper bound on the
growth rate of a function
The statement “f(n) is O(g(n))” means that the growth
rate of f(n) is no more than the growth rate of g(n)
We can use the big-Oh notation to rank functions
according to their growth rate
f(n) is O(g(n))
g(n) is O(f(n))
g(n) grows more
Yes
No
f(n) grows more
No
Yes
Same growth
Yes
Yes
Analysis of Algorithms
31
Big-Oh Rules
If is f(n) a polynomial of degree d, then f(n) is
O(nd), i.e.,
1.
2.
Drop lower-order terms
Drop constant factors
Use the smallest possible class of functions

Say “2n is O(n)” instead of “2n is O(n2)”
Use the simplest expression of the class

Say “3n + 5 is O(n)” instead of “3n + 5 is O(3n)”
Analysis of Algorithms
32
Asymptotic Algorithm Analysis
The asymptotic analysis of an algorithm determines
the running time in big-Oh notation
To perform the asymptotic analysis


We find the worst-case number of primitive operations
executed as a function of the input size
We express this function with big-Oh notation
Example:


We determine that algorithm arrayMax executes at most
7n  1 primitive operations
We say that algorithm arrayMax “runs in O(n) time”
Since constant factors and lower-order terms are
eventually dropped anyhow, we can disregard them
when counting primitive operations
Analysis of Algorithms
33
Computing Prefix Averages
We further illustrate
asymptotic analysis with
two algorithms for prefix
averages
The i-th prefix average of
an array X is average of the
first (i + 1) elements of X:
A[i]  (X[0] + X[1] + … + X[i])/(i+1)
Computing the array A of
prefix averages of another
array X has applications to
financial analysis
Analysis of Algorithms
35
30
X
A
25
20
15
10
5
0
1 2 3 4 5 6 7
34
Prefix Averages (Quadratic)
The following algorithm computes prefix averages in
quadratic time by applying the definition
Algorithm prefixAverages1(X, n)
Input array X of n integers
Output array A of prefix averages of X #operations
A  new array of n integers
n
for i  0 to n  1 do
n
s  X[0]
n
for j  1 to i do
1 + 2 + …+ (n  1)
s  s + X[j]
1 + 2 + …+ (n  1)
A[i]  s / (i + 1)
n
return A
1
Analysis of Algorithms
35
Arithmetic Progression
The running time of
prefixAverages1 is
O(1 + 2 + …+ n)
The sum of the first n
integers is n(n + 1) / 2

There is a simple visual
proof of this fact
Thus, algorithm
prefixAverages1 runs in
O(n2) time
7
6
5
4
3
2
1
0
Analysis of Algorithms
1
2
3
4
5
6
36
Prefix Averages
Asymptotic analysis tells us the complexity of the
original algorithm.
[TEAMS] Given such a tool, how can we make the
algorithm more efficient?
Analysis of Algorithms
37
Prefix Averages (Linear)
The following algorithm computes prefix averages in
linear time by keeping a running sum
Algorithm prefixAverages2(X, n)
Input array X of n integers
Output array A of prefix averages of X
A  new array of n integers
s0
for i  0 to n  1 do
s  s + X[i]
A[i]  s / (i + 1)
return A
#operations
n
1
n
n
n
1
Algorithm prefixAverages2 runs in O(n) time
Analysis of Algorithms
38
Math you need to Review
Summations
Logarithms and Exponents
Proof techniques
Basic probability
properties of logarithms:
logb(xy) = logbx + logby
logb (x/y) = logbx - logby
logbxa = alogbx
logba = logxa/logxb
properties of exponentials:
a(b+c) = aba c
abc = (ab)c
ab /ac = a(b-c)
b = a logab
bc = a c*logab
Analysis of Algorithms
39
Relatives of Big-Oh
big-Omega
 f(n) is (g(n)) if there is a constant c > 0
and an integer constant n0  1 such that
f(n)  c•g(n) for n  n0
big-Theta
 f(n) is (g(n)) if there are constants c’ > 0 and c’’ > 0 and an
integer constant n0  1 such that c’•g(n)  f(n)  c’’•g(n) for n  n0
little-oh
 f(n) is o(g(n)) if, for any constant c > 0, there is an integer
constant n0  0 such that f(n)  c•g(n) for n  n0
little-omega
 f(n) is (g(n)) if, for any constant c > 0, there is an integer
constant n0  0 such that f(n)  c•g(n) for n  n0
Analysis of Algorithms
40
Intuition for Asymptotic
Notation
Big-Oh
 f(n) is O(g(n)) if f(n) is asymptotically less than or equal to g(n)
big-Omega
 f(n) is (g(n)) if f(n) is asymptotically greater than or equal to g(n)
big-Theta
 f(n) is (g(n)) if f(n) is asymptotically equal to g(n)
little-oh
 f(n) is o(g(n)) if f(n) is asymptotically strictly less than g(n)
little-omega
 f(n) is (g(n)) if is asymptotically strictly greater than g(n)
Analysis of Algorithms
41
Example Uses of the
Relatives of Big-Oh

5n2 is (n2)
f(n) is (g(n)) if there is a constant c > 0 and an integer constant n0  1
such that f(n)  c•g(n) for n  n0
[TEAMS] What are values for c and n0?
Let c = 5 and n0 = 1

5n2 is (n)

f(n) is (g(n)) if there is a constant c > 0 and an integer constant n0  1
such that f(n)  c•g(n) for n  n0
Let c = 1 and n0 = 1
5n2 is (n)
f(n) is (g(n)) if, for any constant c > 0, there is an integer constant n0 
0 such that f(n)  c•g(n) for n  n0
Need 5n02  c•n0  given c, the n0 that satisfies this is n0  c/5  0
Analysis of Algorithms
42
Euclid’s Algorithm
An algorithm for computing the greatest common divisor
(GCD) of two numbers MN:
Algorithm GCD(M, N)
while (N != 0)
rem  M mod N
MN
N  rem
endwhile
return M
[TEAMS] What is the Big-Oh?
Analysis of Algorithms
43
Euclid’s Algorithm
An algorithm for computing the greatest common divisor
(GCD) of two numbers MN:
Algorithm GCD(M, N)
while (N != 0)
rem  M mod N
MN
N  rem
endwhile
return M
The algorithm GCD runs is O(logn)
Analysis of Algorithms
44
Euclid’s Algorithm
Why?


If M  N, then (M mod N) < M/2
Thus, each iteration at least halves the
value of M
Analysis of Algorithms
45
Exponentiation
A recursive algorithm to compute X to the power N:
Algorithm pow(X, N)
if (N == 0) return 1
if (N == 1) return X
if (N is even)
return pow(X*X, N/2)
else
return pow(X*X, N/2) * X
1
[TEAMS] What is the Big-Oh?
Analysis of Algorithms
46
Exponentiation
A recursive algorithm to compute X to the power N:
Algorithm pow(X, N)
if (N == 0) return 1
if (N == 1) return X
if (N is even)
return pow(X*X, N/2)
else
return pow(X*X, N/2) * X
1
The algorithm pow is O(logn)
Analysis of Algorithms
47