Data Structures - Ellis Horowitz, Sartaj Sahni - Fundamentals - PDF

Fundamentals: Table of Contents www.itdevelopteam.com Fundamentals of Data Structures by Ellis Horowitz and Sartaj Sahni PREFACE CHAPTER 1: INTRODUCTION CHAPTER 2: ARRAYS CHAPTER 3: STACKS AND...

Fundamentals: Table of Contents www.itdevelopteam.com Fundamentals of Data Structures by Ellis Horowitz and Sartaj Sahni PREFACE CHAPTER 1: INTRODUCTION CHAPTER 2: ARRAYS CHAPTER 3: STACKS AND QUEUES CHAPTER 4: LINKED LISTS CHAPTER 5: TREES CHAPTER 6: GRAPHS CHAPTER 7: INTERNAL SORTING CHAPTER 8: EXTERNAL SORTING CHAPTER 9: SYMBOL TABLES CHAPTER 10: FILES APPENDIX A: SPARKS APPENDIX B: ETHICAL CODE IN INFORMATION PROCESSING APPENDIX C: ALGORITHM INDEX BY CHAPTER file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrDobbs_Books_Algorithms_Collection2ed/books/book1/toc.htm7/3/2004 3:56:06 PM  Fundamentals: PREFACE www.itdevelopteam.com PREFACE For many years a data structures course has been taught in computer science programs. Often it is regarded as a central course of the curriculum. It is fascinating and instructive to trace the history of how the subject matter for this course has changed. Back in the middle1960's the course was not entitled Data Structures but perhaps List Processing Languages. The major subjects were systems such as SLIP (by J. Weizenbaum), IPL-V (by A. Newell, C. Shaw, and H. Simon), LISP 1.5 (by J. McCarthy) and SNOBOL (by D. Farber, R. Griswold, and I. Polonsky). Then, in 1968, volume I of the Art of Computer Programming by D. Knuth appeared. His thesis was that list processing was not a magical thing that could only be accomplished within a specially designed system. Instead, he argued that the same techniques could be carried out in almost any language and he shifted the emphasis to efficient algorithm design. SLIP and IPL-V faded from the scene, while LISP and SNOBOL moved to the programming languages course. The new strategy was to explicitly construct a representation (such as linked lists) within a set of consecutive storage locations and to describe the algorithms by using English plus assembly language. Progress in the study of data structures and algorithm design has continued. Out of this recent work has come many good ideas which we believe should be presented to students of computer science. It is our purpose in writing this book to emphasize those trends which we see as especially valuable and long lasting. The most important of these new concepts is the need to distinguish between the specification of a data structure and its realization within an available programming language. This distinction has been mostly blurred in previous books where the primary emphasis has either been on a programming language or on representational techniques. Our attempt here has been to separate out the specification of the data structure from its realization and to show how both of these processes can be successfully accomplished. The specification stage requires one to concentrate on describing the functioning of the data structure without concern for its implementation. This can be done using English and mathematical notation, but here we introduce a programming notation called axioms. The resulting implementation independent specifications valuable in two ways: (i) to help prove that a program which uses this data structure is correct and (ii) to prove that a particular implementation of the data structure is correct. To describe a data structure in a representation independent way one needs a syntax. This can be seen at the end of section 1.1 where we also precisely define the notions of data object and data structure. This book also seeks to teach the art of analyzing algorithms but not at the cost of undue mathematical sophistication. The value of an implementation ultimately relies on its resource utilization: time and space. This implies that the student needs to be capable of analyzing these factors. A great many analyses have appeared in the literature, yet from our perspective most students don't attempt to rigorously analyze their programs. The data structures course comes at an opportune time in their training to advance and promote these ideas. For every algorithm that is given here we supply a simple, yet rigorous worst case analysis of its behavior. In some cases the average computing time is also file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrDob...Books_Algorithms_Collection2ed/books/book1/preface.htm (1 of 4)7/3/2004 3:56:18 PM  Fundamentals: PREFACE www.itdevelopteam.com derived. The growth of data base systems has put a new requirement on data structures courses, namely to cover the organization of large files. Also, many instructors like to treat sorting and searching because of the richness of its examples of data structures and its practical application. The choice of our later chapters reflects this growing interest. One especially important consideration is the choice of an algorithm description language. Such a choice is often complicated by the practical matters of student background and language availability. Our decision was to use a syntax which is particularly close to ALGOL, but not to restrict ourselves to a specific language. This gives us the ability to write very readable programs but at the same time we are not tied to the idiosyncracies of a fixed language. Wherever it seemed advisable we interspersed English descriptions so as not to obscure the main pointof an algorithm. For people who have not been exposed to the IF-THEN-ELSE, WHILE, REPEAT- UNTIL and a few other basic statements, section 1.2 defines their semantics via flowcharts. For those who have only FORTRAN available, the algorithms are directly translatable by the rules given in the appendix and a translator can be obtained (see appendix A). On the other hand, we have resisted the temptation to use language features which automatically provide sophisticated data structuring facilities. We have done so on several grounds. One reason is the need to commit oneself to a syntax which makes the book especially hard to read by those as yet uninitiated. Even more importantly, these automatic featules cover up the implementation detail whose mastery remains a cornerstone of the course. The basic audience for this book is either the computer science major with at least one year of courses or a beginning graduate student with prior training in a field other than computer science. This book contains more than one semester's worth of material and several of its chapters may be skipped without harm. The following are two scenarios which may help in deciding what chapters should be covered. The first author has used this book with sophomores who have had one semester of PL/I and one semester of assembly language. He would cover chapters one through five skipping sections 2.2, 2.3, 3.2, 4.7, 4.11, and 5.8. Then, in whatever time was left chapter seven on sorting was covered. The second author has taught the material to juniors who have had one quarter of FORTRAN or PASCAL and two quarters of introductory courses which themselves contain a potpourri of topics. In the first quarter's data structure course, chapters one through three are lightly covered and chapters four through six are completely covered. The second quarter starts with chapter seven which provides an excellent survey of the techniques which were covered in the previous quarter. Then the material on external sorting, symbol tables and files is sufficient for the remaining time. Note that the material in chapter 2 is largely mathematical and can be skipped without harm. The paradigm of class presentation that we have used is to begin each new topic with a problem, usually chosen from the computer science arena. Once defined, a high level design of its solution is made and each data structure is axiomatically specified. A tentative analysis is done to determine which operations are critical. Implementations of the data structures are then given followed by an attempt at verifying file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrDob...Books_Algorithms_Collection2ed/books/book1/preface.htm (2 of 4)7/3/2004 3:56:18 PM  Fundamentals: PREFACE www.itdevelopteam.com that the representation and specifications are consistent. The finishedalgorithm in the book is examined followed by an argument concerning its correctness. Then an analysis is done by determining the relevant parameters and applying some straightforward rules to obtain the correct computing time formula. In summary, as instructors we have tried to emphasize the following notions to our students: (i) the ability to define at a sufficiently high level of abstraction the data structures and algorithms that are needed; (ii) the ability to devise alternative implementations of a data structure; (iii) the ability to synthesize a correct algorithm; and (iv) the abilityto analyze the computing time of the resultant program. In addition there are two underlying currents which, though not explicitly emphasized are covered throughout. The first is the notion of writing nicely structured programs. For all of the programs contained herein we have tried our best to structure them appropriately. We hope that by reading programs with good style the students will pick up good writing habits. A nudge on the instructor's part will also prove useful. The second current is the choice of examples. We have tried to use those examples which prove a point well, have application to computer programming, and exhibit some of the brightest accomplishments in computer science. At the close of each chapter there is a list of references and selected readings. These are not meant to be exhaustive. They are a subset of those books and papers that we found to be the most useful. Otherwise, they are either historically significant or develop the material in the text somewhat further. Many people have contributed their time and energy to improve this book. For this we would like to thank them. We wish to thank Arvind [sic], T. Gonzalez, L. Landweber, J. Misra, and D. Wilczynski, who used the book in their own classes and gave us detailed reactions. Thanks are also due to A. Agrawal, M. Cohen, A. Howells, R. Istre, D. Ledbetter, D. Musser and to our students in CS 202, CSci 5121 and 5122 who provided many insights. For administrative and secretarial help we thank M. Eul, G. Lum, J. Matheson, S. Moody, K. Pendleton, and L. Templet. To the referees for their pungent yet favorable comments we thank S. Gerhart, T. Standish, and J. Ullman. Finally, we would like to thank our institutions, the University of Southern California and the University of Minnesota, for encouraging in every way our efforts to produce this book. Ellis Horowitz Sartaj Sahni Preface to the Ninth Printing We would like to acknowledge collectively all of the individuals who have sent us comments and corrections since the book first appeared. For this printing we have made many corrections and improvements. October 198l file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrDob...Books_Algorithms_Collection2ed/books/book1/preface.htm (3 of 4)7/3/2004 3:56:18 PM  Fundamentals: PREFACE www.itdevelopteam.com Ellis Horowitz Sartaj Sahni file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrDob...Books_Algorithms_Collection2ed/books/book1/preface.htm (4 of 4)7/3/2004 3:56:18 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com CHAPTER 1: INTRODUCTION 1.1 OVERVIEW The field of computer science is so new that one feels obliged to furnish a definition before proceeding with this book. One often quoted definition views computer science as the study of algorithms. This study encompasses four distinct areas: (i) machines for executing algorithms--this area includes everything from the smallest pocket calculator to the largest general purpose digital computer. The goal is to study various forms of machine fabrication and organization so that algorithms can be effectively carried out. (ii) languages for describing algorithms--these languages can be placed on a continuum. At one end are the languages which are closest to the physical machine and at the other end are languages designed for sophisticated problem solving. One often distinguishes between two phases of this area: language design and translation. The first calls for methods for specifying the syntax and semantics of a language. The second requires a means for translation into a more basic set of commands. (iii) foundations of algorithms--here people ask and try to answer such questions as: is a particular task accomplishable by a computing device; or what is the minimum number of operations necessary for any algorithm which performs a certain function? Abstract models of computers are devised so that these properties can be studied. (iv) analysis of algorithms--whenever an algorithm can be specified it makes sense to wonder about its behavior. This was realized as far back as 1830 by Charles Babbage, the father of computers. An algorithm's behavior pattern or performance profile is measured in terms of the computing time and space that are consumed while the algorithm is processing. Questions such as the worst and average time and how often they occur are typical. We see that in this definition of computer science, "algorithm" is a fundamental notion. Thus it deserves a precise definition. The dictionary's definition "any mechanical or recursive computational procedure" is not entirely satisfying since these terms are not basic enough. Definition: An algorithm is a finite set of instructions which, if followed, accomplish a particular task. In addition every algorithm must satisfy the following criteria: (i) input: there are zero or more quantities which are externally supplied; (ii) output: at least one quantity is produced; file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrDo...Books_Algorithms_Collection2ed/books/book1/chap01.htm (1 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com (iii) definiteness: each instruction must be clear and unambiguous; (iv) finiteness: if we trace out the instructions of an algorithm, then for all cases the algorithm will terminate after a finite number of steps; (v) effectiveness: every instruction must be sufficiently basic that it can in principle be carried out by a person using only pencil and paper. It is not enough that each operation be definite as in (iii), but it must also be feasible. In formal computer science, one distinguishes between an algorithm, and a program. A program does not necessarily satisfy condition (iv). One important example of such a program for a computer is its operating system which never terminates (except for system crashes) but continues in a wait loop until more jobs are entered. In this book we will deal strictly with programs that always terminate. Hence, we will use these terms interchangeably. An algorithm can be described in many ways. A natural language such as English can be used but we must be very careful that the resulting instructions are definite (condition iii). An improvement over English is to couple its use with a graphical form of notation such as flowcharts. This form places each processing step in a "box" and uses arrows to indicate the next step. Different shaped boxes stand for different kinds of operations. All this can be seen in figure 1.1 where a flowchart is given for obtaining a Coca-Cola from a vending machine. The point is that algorithms can be devised for many common activities. Have you studied the flowchart? Then you probably have realized that it isn't an algorithm at all! Which properties does it lack? Returning to our earlier definition of computer science, we find it extremely unsatisfying as it gives us no insight as to why the computer is revolutionizing our society nor why it has made us re-examine certain basic assumptions about our own role in the universe. While this may be an unrealistic demand on a definition even from a technical point of view it is unsatisfying. The definition places great emphasis on the concept of algorithm, but never mentions the word "data". If a computer is merely a means to an end, then the means may be an algorithm but the end is the transformation of data. That is why we often hear a computer referred to as a data processing machine. Raw data is input and algorithms are used to transform it into refined data. So, instead of saying that computer science is the study of algorithms, alternatively, we might say that computer science is the study of data: (i) machines that hold data; (ii) languages for describing data manipulation; (iii) foundations which describe what kinds of refined data can be produced from raw data; file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrDo...Books_Algorithms_Collection2ed/books/book1/chap01.htm (2 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com (iv) structures for representing data. Figure 1.1: Flowchart for obtaining a Coca-Cola There is an intimate connection between the structuring of data, and the synthesis of algorithms. In fact, a data structure and an algorithm should be thought of as a unit, neither one making sense without the other. For instance, suppose we have a list of n pairs of names and phone numbers (a1,b1)(a2,b2),..., (an, bn), and we want to write a program which when given any name, prints that person's phone number. This task is called searching. Just how we would write such an algorithm critically depends upon how the names and phone numbers are stored or structured. One algorithm might just forge ahead and examine names, a1,a2,a3,... etc., until the correct name was found. This might be fine in Oshkosh, but in Los Angeles, with hundreds of thousands of names, it would not be practical. If, however, we knew that the data was structured so that the names were in alphabetical order, then we could do much better. We could make up a second list which told us for each letter in the alphabet, where the first name with that letter appeared. For a name beginning with, say, S, we would avoid having to look at names beginning with other letters. So because of this new structure, a very different algorithm is possible. Other ideas for algorithms become possible when we realize that we can organize the data as we wish. We will discuss many more searching strategies in Chapters 7 and 9. Therefore, computer science can be defined as the study of data, its representation and transformation by a digital computer. The goal of this book is to explore many different kinds of data objects. For each object, we consider the class of operations to be performed and then the way to represent this object so that these operations may be efficiently carried out. This implies a mastery of two techniques: the ability to devise alternative forms of data representation, and the ability to analyze the algorithm which operates on that structure. The pedagogical style we have chosen is to consider problems which have arisen often in computer applications. For each problem we will specify the data object or objects and what is to be accomplished. After we have decided upon a representation of the objects, we will give a complete algorithm and analyze its computing time. After reading through several of these examples you should be confident enough to try one on your own. There are several terms we need to define carefully before we proceed. These include data structure, data object, data type and data representation. These four terms have no standard meaning in computer science circles, and they are often used interchangeably. A data type is a term which refers to the kinds of data that variables may "hold" in a programming language. In FORTRAN the data types are INTEGER, REAL, LOGICAL, COMPLEX, and DOUBLE PRECISION. In PL/I there is the data type CHARACTER. The fundamental data type of SNOBOL is the character string and in LISP it is the list (or S-expression). With every programming language there is a set of built-in data types. This means that the language allows variables to name data of that type and file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrDo...Books_Algorithms_Collection2ed/books/book1/chap01.htm (3 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com provides a set of operations which meaningfully manipulates these variables. Some data types are easy to provide because they are already built into the computer's machine language instruction set. Integer and real arithmetic are examples of this. Other data types require considerably more effort to implement. In some languages, there are features which allow one to construct combinations of the built-in types. In COBOL and PL/I this feature is called a STRUCTURE while in PASCAL it is called a RECORD. However, it is not necessary to have such a mechanism. All of the data structures we will see here can be reasonably built within a conventional programming language. Data object is a term referring to a set of elements, say D. For example the data object integers refers to D = {0, 1, 2,...}. The data object alphabetic character strings of length less than thirty one implies D = {",'A','B',...,'Z','AA',...}. Thus, D may be finite or infinite and if D is very large we may need to devise special ways of representing its elements in our computer. The notion of a data structure as distinguished from a data object is that we want to describe not only the set of objects, but the way they are related. Saying this another way, we want to describe the set of operations which may legally be applied to elements of the data object. This implies that we must specify the set of operations and show how they work. For integers we would have the arithmetic operations +, -, *, / and perhaps many others such as mod, ceil, floor, greater than, less than, etc. The data object integers plus a description of how +, -, *, /, etc. behave constitutes a data structure definition. To be more precise lets examine a modest example. Suppose we want to define the data structure natural number (abbreviated natno) where natno = {0,1,2,3,...} with the three operations being a test for zero addition and equality. The following notation can be used: structure NATNO 1 declare ZERO( ) natno 2 ISZERO(natno) boolean 3 SUCC(natno) natno 4 ADD(natno, natno) natno 5 EQ(natno, natno) boolean 6 for all x, y natno let 7 ISZERO(ZERO) ::= true; ISZERO(SUCC(x)) ::= false file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrDo...Books_Algorithms_Collection2ed/books/book1/chap01.htm (4 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com 8 ADD(ZERO, y) :: = y, ADD(SUCC(x), y) :: = SUCC(ADD(x, y)) 9 EQ(x, ZERO) :: = if ISZERO(x) then true else false 10 EQ(ZERO, SUCC(y)) :: = false EQ(SUCC(x), SUCC(y)) :: = EQ(x, y) 11 end end NATNO In the declare statement five functions are defined by giving their names, inputs and outputs. ZERO is a constant function which means it takes no input arguments and its result is the natural number zero, written as ZERO. ISZERO is a boolean function whose result is either true or false. SUCC stands for successor. Using ZERO and SUCC we can define all of the natural numbers as: ZERO, l = SUCC (ZERO), 2 = SUCC(SUCC(ZERO)), 3 = SUCC(SUCC(SUCC(ZERO))),... etc. The rules on line 8 tell us exactly how the addition operation works. For example if we wanted to add two and three we would get the following sequence of expressions: ADD(SUCC(SUCC(ZERO)),SUCC(SUCC(SUCC(ZERO)))) which, by line 8 equals SUCC(ADD(SUCC(ZERO),SUCC(SUCC(SUCC(ZERO))))) which, by line 8 equals SUCC(SUCC(ADD(ZERO,SUCC(SUCC(SUCC(ZERO)))))) which by line 8 equals SUCC(SUCC(SUCC(SUCC(SUCC(ZERO))))) Of course, this is not the way to implement addition. In practice we use bit strings which is a data structure that is usually provided on our computers. But however the ADD operation is implemented, it must obey these rules. Hopefully, this motivates the following definition. Definition: A data structure is a set of domains , a designated domain , a set of functions and a file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrDo...Books_Algorithms_Collection2ed/books/book1/chap01.htm (5 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com set of axioms. The triple denotes the data structure d and it will usually be abbreviated by writing d. In the previous example The set of axioms describes the semantics of the operations. The form in which we choose to write the axioms is important. Our goal here is to write the axioms in a representation independent way. Then, we discuss ways of implementing the functions using a conventional programming language. An implementation of a data structure d is a mapping from d to a set of other data structures e. This mapping specifies how every object of d is to be represented by the objects of e. Secondly, it requires that every function of d must be written using the functions of the implementing data structures e. Thus we say that integers are represented by bit strings, boolean is represented by zero and one, an array is represented by a set of consecutive words in memory. In current parlance the triple is referred to as an abstract data type. It is called abstract precisely because the axioms do not imply a form of representation. Another way of viewing the implementation of a data structure is that it is the process of refining an abstract data type until all of the operations are expressible in terms of directly executable functions. But at the first stage a data structure should be designed so that we know what it does, but not necessarily how it will do it. This division of tasks, called specification and implementation, is useful because it helps to control the complexity of the entire process. 1.2 SPARKS The choice of an algorithm description language must be carefully made because it plays such an important role throughout the book. We might begin by considering using some existing language; some names which come immediately to mind are ALGOL, ALGOL-W, APL, COBOL, FORTRAN, LISP, PASCAL, PL/I, SNOBOL. Though some of these are more preferable than others, the choice of a specific language leaves us with many difficulties. First of all, we wish to be able to write our algorithms without dwelling on the idiosyncracies of a given language. Secondly, some languages have already provided the mechanisms we wish to discuss. Thus we would have to make pretense to build up a capability which already exists. Finally, each language has its followers and its detractors. We would rather not have any individual rule us out simply because he did not know or, more particularly, disliked to use the language X. Furthermore it is not really necessary to write programs in a language for which a compiler exists. Instead we choose to use a language which is tailored to describing the algorithms we want to write. file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrDo...Books_Algorithms_Collection2ed/books/book1/chap01.htm (6 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com Using it we will not have to define many aspects of a language that we will never use here. Most importantly, the language we use will be close enough to many of the languages mentioned before so that a hand translation will be relatively easy to accomplish. Moreover, one can easily program a translator using some existing, but more primitive higher level language as the output (see Appendix A). We call our language SPARKS. Figure 1.2 shows how a SPARKS program could be executed on any machine. Figure 1.2: Translation of SPARKS Many language designers choose a name which is an acronym. But SPARKS was not devised in that way; it just appeared one day as Athena sprang from the head of Zeus. Nevertheless, computerniks still try to attach a meaning. Several cute ideas have been suggested, such as Structured Programming: A Reasonably Komplete Set or Smart Programmers Are Required To Know SPARKS. SPARKS contains facilities to manipulate numbers, boolean values and characters. The way to assign values is by the assignment statement variable expression. In addition to the assignment statement, SPARKS includes statements for conditional testing, iteration, input-output, etc. Several such statements can be combined on a single line if they are separated by a semi-colon. Expressions can be either arithmetic, boolean or of character type. In the boolean case there can be only one of two values, true or false. In order to produce these values, the logical operators and, or, not are provided, plus the relational operators file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrDo...Books_Algorithms_Collection2ed/books/book1/chap01.htm (7 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com A conditional statement has the form if cond then S1 if cond then S1 or else S2 where cond is a boolean expression and S1, S2 are arbitrary groups of SPARKS statements. If S1 or S2 contains more than one statement, these will be enclosed in square brackets. Brackets must be used to show how each else corresponds to one if. The meaning of this statement is given by the flow charts: We will assume that conditional expressions are evaluated in "short circuit" mode; given the boolean expression (cond1 or cond2), if condl is true then cond2 is not evaluated; or, given (condl and cond2), if cond1 is false then cond2 is not evaluated. To accomplish iteration, several statements are available. One of them is while cond do S end where cond is as before, S is as S1 before and the meaning is given by It is well known that all "proper" programs can be written using only the assignment, conditional and while statements. This result was obtained by Bohm and Jacopini. Though this is very interesting from a theoretical viewpoint, we should not take it to mean that this is the way to program. On the contrary, the more expressive our languages are, the more we can accomplish easily. So we will provide other statements such as a second iteration statement, the repeat-until, repeat S file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrDo...Books_Algorithms_Collection2ed/books/book1/chap01.htm (8 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com until cond which has the meaning In contrast to the while statement, the repeat-until guarantees that the statements of S will be executed at least once. Another iteration statement is loop S forever which has the meaning As it stands, this describes an infinite loop! However, we assume that this statement is used in conjunction with some test within S which will cause an exit. One way of exiting such a loop is by using a go to label statement which transfers control to "label." Label may be anywhere in the procedure. A more restricted form of the go to is the command exit which will cause a transfer of control to the first statement after the innermost loop which contains it. This looping statement may be a while, repeat, for or a loop-forever. exit can be used either conditionally or unconditionally, for instance loop S1 if cond then exit S2 file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrDo...Books_Algorithms_Collection2ed/books/book1/chap01.htm (9 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com forever which will execute as The last statement for iteration is called the for-loop, which has the form for vble start to finish by increment do S end Vble is a variable, while start, finish and increment are arithmetic expressions. A variable or a constant is a simple form of an expression. The clause "by increment" is optional and taken as +1 if it does not occur. We can write the meaning of this statement in SPARKS as vble start fin finish incr increment while (vble - fin) * incr 0 do S vble vble + incr end Another statement within SPARKS is the case, which allows one to distinguish easily between several alternatives without using multiple if-then-else statements. It has the form where the Si, 1 i n + 1 are groups of SPARKS statements. The semantics is easily described by the file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (10 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com following flowchart: The else clause is optional. A complete SPARKS procedure has the form procedure NAME (parameter list) S end A procedure can be used as a function by using the statement return (expr) where the value of expr is delivered as the value of the procedure. The expr may be omitted in which case a return is made to the calling procedure. The execution of an end at the end of procedure implies a return. A procedure may be invoked by using a call statement call NAME (parameter list) Procedures may call themselves, direct recursion, or there may be a sequence resulting in indirect recursion. Though recursion often carries with it a severe penalty at execution time, it remains all elegant way to describe many computing processes. This penalty will not deter us from using recursion. Many such programs are easily translatable so that the recursion is removed and efficiency achieved. A complete SPARKS program is a collection of one or more procedures, the first one taken as the main program. All procedures are treated as external, which means that the only means for communication between them is via parameters. This may be somewhat restrictive in practice, but for the purpose of exposition it helps to list all variables explicitly, as either local or parameter. The association of actual to formal parameters will be handled using the call by reference rule. This means that at run time the address of each parameter is passed to the called procedure. Parameters which are constants or values of expressions are stored into internally generated words whose addresses are then passed to the procedure. For input/output we assume the existence of two functions read (argument list), print (argument list) file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (11 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com Arguments may be variables or quoted strings. We avoid the problem of defining a "format" statement as we will need only the simplest form of input and output. The command stop halts execution of the currently executing procedure. Comments may appear anywhere on a line enclosed by double slashes, e.g. //this is a comment// Finally, we note that multi-dimensional arrays are available with arbitrary integer lower and upper bounds. An n-dimensional array A with lower and upper bounds li, ui, 1 i n may be declared by using the syntax declare A(l1:u1,...,ln:un). We have avoided introducing the record or structure concept. These are often useful features and when available they should be used. However, we will persist in building up a structure from the more elementary array concept. Finally, we emphasize that all of our variables are assumed to be of type INTEGER unless stated otherwise. Since most of the SPARKS programs will be read many more times than they will be executed, we have tried to make the code readable. This is a goal which should be aimed at by everyone who writes programs. The SPARKS language is rich enough so that one can create a good looking program by applying some simple rules of style. (i) Every procedure should carefully specify its input and output variables. (ii) The meaning of variables should be defined. (iii) The flow of the program should generally be forward except for normal looping or unavoidable instances. (iv) Indentation rules should be established and followed so that computational units of program text can more easily be identified. (v) Documentation should be short, but meaningful. Avoid sentences like ''i is increased by one." (vi) Use subroutines where appropriate. See the book The Elements of Programming Style by Kernighan and Plauger for more examples of good rules of programming. 1.3 HOW TO CREATE PROGRAMS Now that you have moved beyond the first course in computer science, you should be capable of file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (12 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com developing your programs using something better than the seat-of-the-pants method. This method uses the philosophy: write something down and then try to get it working. Surprisingly, this method is in wide use today, with the result that an average programmer on an average job turns out only between five to ten lines of correct code per day. We hope your productivity will be greater. But to improve requires that you apply some discipline to the process of creating programs. To understand this process better, we consider it as broken up into five phases: requirements, design, analysis, coding, and verification. (i) Requirements. Make sure you understand the information you are given (the input) and what results you are to produce (the output). Try to write down a rigorous description of the input and output which covers all cases. You are now ready to proceed to the design phase. Designing an algorithm is a task which can be done independently of the programming language you eventually plan to use. In fact, this is desirable because it means you can postpone questions concerning how to represent your data and what a particular statement looks like and concentrate on the order of processing. (ii) Design. You may have several data objects (such as a maze, a polynomial, or a list of names). For each object there will be some basic operations to perform on it (such as print the maze, add two polynomials, or find a name in the list). Assume that these operations already exist in the form of procedures and write an algorithm which solves the problem according to the requirements. Use a notation which is natural to the way you wish to describe the order of processing. (iii) Analysis. Can you think of another algorithm? If so, write it down. Next, try to compare these two methods. It may already be possible to tell if one will be more desirable than the other. If you can't distinguish between the two, choose one to work on for now and we will return to the second version later. (iv) Refinement and coding. You must now choose representations for your data objects (a maze as a two dimensional array of zeros and ones, a polynomial as a one dimensional array of degree and coefficients, a list of names possibly as an array) and write algorithms for each of the operations on these objects. The order in which you do this may be crucial, because once you choose a representation, the resulting algorithms may be inefficient. Modern pedagogy suggests that all processing which is independent of the data representation be written out first. By postponing the choice of how the data is stored we can try to isolate what operations depend upon the choice of data representation. You should consider alternatives, note them down and review them later. Finally you produce a complete version of your first program. It is often at this point that one realizes that a much better program could have been built. Perhaps you should have chosen the second design alternative or perhaps you have spoken to a friend who has done it better. This happens to industrial programmers as well. If you have been careful about keeping track of your previous work it may not be too difficult to make changes. One of the criteria of a good design is file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (13 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com that it can absorb changes relatively easily. It is usually hard to decide whether to sacrifice this first attempt and begin again or just continue to get the first version working. Different situations call for different decisions, but we suggest you eliminate the idea of working on both at the same time. If you do decide to scrap your work and begin again, you can take comfort in the fact that it will probably be easier the second time. In fact you may save as much debugging time later on by doing a new version now. This is a phenomenon which has been observed in practice. The graph in figure 1.3 shows the time it took for the same group to build 3 FORTRAN compilers (A, B and C). For each compiler there is the time they estimated it would take them and the time it actually took. For each subsequent compiler their estimates became closer to the truth, but in every case they underestimated. Unwarrented optimism is a familiar disease in computing. But prior experience is definitely helpful and the time to build the third compiler was less than one fifth that for the first one. Figure 1.3: History of three FORTRAN compilers (v) Verification. Verification consists of three distinct aspects: program proving, testing and debugging. Each of these is an art in itself. Before executing your program you should attempt to prove it is correct. Proofs about programs are really no different from any other kinds of proofs, only the subject matter is different. If a correct proof can be obtained, then one is assured that for all possible combinations of inputs, the program and its specification agree. Testing is the art of creating sample data upon which to run your program. If the program fails to respond correctly then debugging is needed to determine what went wrong and how to correct it. One proof tells us more than any finite amount of testing, but proofs can be hard to obtain. Many times during the proving process errors are discovered in the code. The proof can't be completed until these are changed. This is another use of program proving, namely as a methodology for discovering errors. Finally there may be tools available at your computing center to aid in the testing process. One such tool instruments your source code and then tells you for every data set: (i) the number of times a statement was executed, (ii) the number of times a branch was taken, (iii) the smallest and largest values of all variables. As a minimal requirement, the test data you construct should force every statement to execute and every condition to assume the value true and false at least once. One thing you have forgotten to do is to document. But why bother to document until the program is entirely finished and correct ? Because for each procedure you made some assumptions about its input and output. If you have written more than a few procedures, then you have already begun to forget what those assumptions were. If you note them down with the code, the problem of getting the procedures to work together will be easier to solve. The larger the software, the more crucial is the need for documentation. The previous discussion applies to the construction of a single procedure as well as to the writing of a large software system. Let us concentrate for a while on the question of developing a single procedure which solves a specific task. This shifts our emphasis away from the management and integration of the file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (14 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com various procedures to the disciplined formulation of a single, reasonably small and well-defined task. The design process consists essentially of taking a proposed solution and successively refining it until an executable program is achieved. The initial solution may be expressed in English or some form of mathematical notation. At this level the formulation is said to be abstract because it contains no details regarding how the objects will be represented and manipulated in a computer. If possible the designer attempts to partition the solution into logical subtasks. Each subtask is similarly decomposed until all tasks are expressed within a programming language. This method of design is called the top-down approach. Inversely, the designer might choose to solve different parts of the problem directly in his programming language and then combine these pieces into a complete program. This is referred to as the bottom-up approach. Experience suggests that the top-down approach should be followed when creating a program. However, in practice it is not necessary to unswervingly follow the method. A look ahead to problems which may arise later is often useful. Underlying all of these strategies is the assumption that a language exists for adequately describing the processing of data at several abstract levels. For this purpose we use the language SPARKS coupled with carefully chosen English narrative. Such an algorithm might be called pseudo-SPARKS. Let us examine two examples of top-down program development. Suppose we devise a program for sorting a set of n 1 distinct integers. One of the simplest solutions is given by the following "from those integers which remain unsorted, find the smallest and place it next in the sorted list" This statement is sufficient to construct a sorting program. However, several issues are not fully specified such as where and how the integers are initially stored and where the result is to be placed. One solution is to store the values in an array in such a way that the i-th integer is stored in the i-th array position, A(i) 1 i n. We are now ready to give a second refinement of the solution: for i 1 to n do examine A(i) to A(n) and suppose the smallest integer is at A(j); then interchange A(i) and A(j). end Note how we have begun to use SPARKS pseudo-code. There now remain two clearly defined subtasks: (i) to find the minimum integer and (ii) to interchange it with A(i). This latter problem can be solved by the code file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (15 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com t A(i); A(i) A(j); A(j) t The first subtask can be solved by assuming the minimum is A (i), checking A(i) with A(i + 1), A(i + 2),... and whenever a smaller element is found, regarding it as the new minimum. Eventually A(n) is compared to the current minimum and we are done. Putting all these observations together we get procedure SORT(A,n) 1 for i 1 to n do 2 j i 3 for k j + 1 to n do 4 if A(k) < A(j) then j k 5 end 6 t A(i); A(i) A(j); A(j) t 7 end end SORT The obvious question to ask at this point is: "does this program work correctly?" Theorem: Procedure SORT (A,n) correctly sorts a set of n 1 distinct integers, the result remains in A (1:n) such that A (1) < A (2) <... < A(n). Proof: We first note that for any i, say i = q, following the execution of lines 2 thru 6, it is the case that A (q) A(r), q < r n. Also, observe that when i becomes greater than q, A(1.. q) is unchanged. Hence, following the last execution of these lines, (i.e., i = n), we have A(1) A(2)... A(n). We observe at this point that the upper limit of the for-loop in line 1 can be changed to n - 1 without damaging the correctness of the algorithm. From the standpoint of readability we can ask if this program is good. Is there a more concise way of describing this algorithm which will still be as easy to comprehend? Substituting while statements for the for loops doesn't significantly change anything. Also, extra initialization and increment statements would be required. We might consider a FORTRAN version using the ANSI language standard file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (16 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com IF (N. LE. 1) GO TO 100 NM1 = N - 1 DO 101 I = 1, NM1 J = I JP1 = J + 1 DO 102 K = JP1, N IF (A(K).LT.A(J)) J = K 102 CONTINUE T = A(I) A(I) = A(J) A(J) = T 101 CONTINUE 100 CONTINUE FORTRAN forces us to clutter up our algorithms with extra statements. The test for N = 1 is necessary because FORTRAN DO-LOOPS always insist on executing once. Variables NM1 and JP1 are needed because of the restrictions on lower and upper limits of DO-LOOPS. Let us develop another program. We assume that we have n 1 distinct integers which are already sorted and stored in the array A(1:n). Our task is to determine if the integer x is present and if so to return j such that x = A(j); otherwise return j = 0. By making use of the fact that the set is sorted we conceive of the following efficient method: "let A(mid) be the middle element. There are three possibilities. Either x < A(mid) in which case x can only occur as A(1) to A(mid - 1); or x > A(mid) in which case x can only occur as A(mid + l) to A(n): or x = A(mid) in which case set j to mid and return. Continue in this way by keeping two pointers, lower and upper, to indicate the range of elements not yet tested." file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (17 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com At this point you might try the method out on some sample numbers. This method is referred to as binary search. Note how at each stage the number of elements in the remaining set is decreased by about one half. We can now attempt a version using SPARKS pseudo code. procedure BINSRCH(A,n,x,j) initialize lower and upper while there are more elements to check do let A(mid) be the middle element case : x > A(mid): set lower to mid + 1 : x < A(mid): set upper to mid - 1 : else: found end end not found end BINSRCH The above is not the only way we might write this program. For instance we could replace the while loop by a repeat-until statement with the same English condition. In fact there are at least six different binary search programs that can be produced which are all correct. There are many more that we might produce which would be incorrect. Part of the freedom comes from the initialization step. Whichever version we choose, we must be sure we understand the relationships between the variables. Below is one complete version. procedure BINSRCH (A,n,x,j) 1 lower 1; upper n 2 while lower upper do file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (18 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com 3 mid (lower + upper) / 2 4 case 5 : x > A(mid): lower mid + 1 6 : x < A(mid): upper mid - 1 7 : else: j mid; return 8 end 9 end 10 j 0 end To prove this program correct we make assertions about the relationship between variables before and after the while loop of steps 2-9. As we enter this loop and as long as x is not found the following holds: lower upper and A (lower) x A (upper) and SORTED (A, n) Now, if control passes out of the while loop past line 9 then we know the condition of line 2 is false lower > upper. This, combined with the above assertion implies that x is not present. Unfortunately a complete proof takes us beyond our scope but for those who wish to pursue program proving they should consult our references at the end of this chapter. An analysis of the computing time for BINSRCH is carried out in section 7.1. Recursion We have tried to emphasize the need to structure a program to make it easier to achieve the goals of readability and correctness. Actually one of the most useful syntactical features for accomplishing this is the procedure. Given a set of instructions which perform a logical operation, perhaps a very complex and long operation, they can be grouped together as a procedure. The procedure name and its parameters file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (19 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com are viewed as a new instruction which can be used in other programs. Given the input-output specifications of a procedure, we don't even have to know how the task is accomplished, only that it is available. This view of the procedure implies that it is invoked, executed and returns control to the appropriate place in the calling procedure. What this fails to stress is the fact that procedures may call themselves (direct recursion) before they are done or they may call other procedures which again invoke the calling procedure (indirect recursion). These recursive mechanisms are extremely powerful, but even more importantly, many times they can express an otherwise complex process very clearly. For these reasons we introduce recursion here. Most students of computer science view recursion as a somewhat mystical technique which only is useful for some very special class of problems (such as computing factorials or Ackermann's function). This is unfortunate because any program that can be written using assignment, the if-then-else statement and the while statement can also be written using assignment, if-then-else and recursion. Of course, this does not say that the resulting program will necessarily be easier to understand. However, there are many instances when this will be the case. When is recursion an appropriate mechanism for algorithm exposition? One instance is when the problem itself is recursively defined. Factorial fits this category, also binomial coefficients where can be recursively computed by the formula Another example is reversing a character string, S = 'x1... xn' where SUBSTRING (S,i,j) is a function which returns the string xi... xj for appropriately defined i and j and S T stands for concatenation of two strings (as in PL/I). Then the operation REVERSE is easily described recursively as procedure REVERSE(S) n LENGTH(S) if n = 1 then return (S) else return (REVERSE(SUBSTRING(S,2,n)) SUBSTRING(S,1,1)) end REVERSE file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (20 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com If this looks too simple let us develop a more complex recursive procedure. Given a set of n 1 elements the problem is to print all possible permutations of this set. For example if the set is {a,b,c}, then the set of permutations is {(a, b,c), (a,c,b), (b,a,c), (b,c,a), (c,a,b), (c,b,a)}. It is easy to see that given n elements there are n ! different permutations. A simple algorithm can be achieved by looking at the case of four elements (a,b,c,d). The answer is obtained by printing (i) a followed by all permutations of (b,c,d) (ii) b followed by all permutations of (a,c,d) (iii) c followed by all permutations of (b,a,d) (iv) d followed by all permutations of (b,c,a) The expression "followed by all permutations" is the clue to recursion. It implies that we can solve the problem for a set with n elements if we had an algorithm which worked on n - 1 elements. These considerations lead to the following procedure which is invoked by call PERM(A,1,n). A is a character string e.g. A ='abcd', and INTERCHANGE (A,k,i) exchanges the k-th character of A with the i-th character of A. procedure PERM(A,k,n) if k = n then [print (A); return] B A for i k to n do call INTERCHANGE(A,k,i) call PERM(A,k + 1,n) A B end end PERM Try this algorithm out on sets of length one, two, and three to insure that you understand how it works. Then try to do one or more of the exercises at the end of this chapter which ask for recursive procedures. file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (21 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com Another time when recursion is useful is when the data structure that the algorithm is to operate on is recursively defined. We will see several important examples of such structures, especially lists in section 4.9 and binary trees in section 5.4. Another instance when recursion is invaluable is when we want to describe a backtracking procedure. But for now we will content ourselves with examining some simple, iterative programs and show how to eliminate the iteration statements and replace them by recursion. This may sound strange, but the objective is not to show that the result is simpler to understand nor more efficient to execute. The main purpose is to make one more familiar with the execution of a recursive procedure. Suppose we start with the sorting algorithm presented in this section. To rewrite it recursively the first thing we do is to remove the for loops and express the algorithm using assignment, if-then-else and the go-to statement. procedure SORT(A,n) i 1 Ll: if i n - 1 // for i 1 to n - 1 do// then [j i; k j + 1 L2: if k n //for k j + 1 to n do// then [if A(k) < A(j) then j k k k + 1; go to L2] t A(i); A(i) A(j); A(j) t i i + 1; go to L1] end SORT Now every place where we have a label we introduce a procedure whose parameters are the variables which are already assigned a value at that point. Every place where a ''go to label'' appears, we replace that statement by a call of the procedure associated with that label. This gives us the following set of three procedures. file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (22 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com procedure SORT(A,n) call SORTL1(A,n,1) end SORT procedure SORTLl(A,n,i) if i n - 1 then [j i; call MAXL2(A,n,j,i + 1) t A(i); A(i) A(j); A(j) t call SORTL1(A,n,i + 1)] end SORTL1 procedure MAXL2(A,n,j,k) if k n then [if A(k) < A(j) then j k call MAXL2(A,n,j,k + 1)] end MAXL2 We can simplify these procedures somewhat by ignoring SORT(A,n) entirely and begin the sorting operation by call SORTL1(A,n,1). Notice how SORTL1 is directly recursive while it also uses procedure MAXL2. Procedure MAXL2 is also directly reculsive. These two procedures use eleven lines while the original iterative version was expressed in nine lines; not much of a difference. Notice how in MAXL2 the fourth parameter k is being changed. The effect of increasing k by one and restarting the procedure has essentially the same effect as the for loop. Now let us trace the action of these procedures as they sort a set of five integers When a procedure is invoked an implicit branch to its beginning is made. Thus a recursive call of a file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (23 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com program can be made to simulate a go to statement. The parameter mechanism of the procedure is a form of assignment. Thus placing the argument k + 1 as the fourth parameter of MAXL2 is equivalent to the statement k k + 1. In section 4.9 we will see the first example of a recursive data structure, the list. Also in that section are several recursive procedures, followed in some cases by their iterative equivalents. Rules are also given there for eliminating recursion. 1.4 HOW TO ANALYZE PROGRAMS One goal of this book is to develop skills for making evaluative judgements about programs. There are many criteria upon which we can judge a program, for instance: (i) Does it do what we want it to do? (ii) Does it work correctly according to the original specifications of the task? (iii) Is there documentation which describes how to use it and how it works? (iv) Are subroutines created in such a way that they perform logical sub-functions? (v) Is the code readable? The above criteria are all vitally important when it comes to writing software, most especially for large systems. Though we will not be discussing how to reach these goals, we will try to achieve them throughout this book with the programs we write. Hopefully this more subtle approach will gradually infect your own program writing habits so that you will automatically strive to achieve these goals. There are other criteria for judging programs which have a more direct relationship to performance. These have to do with computing time and storage requirements of the algorithms. Performance evaluation can be loosely divided into 2 major phases: (a) a priori estimates and (b) a posteriori testing. Both of these are equally important. First consider a priori estimation. Suppose that somewhere in one of your programs is the statement x x + 1. We would like to determine two numbers for this statement. The first is the amount of time a single execution will take; the second is the number of times it is executed. The product of these numbers will be the total time taken by this statement. The second statistic is called the frequency count, and this may file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (24 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com vary from data set to data set. One of the hardest tasks in estimating frequency counts is to choose adequate samples of data. It is impossible to determine exactly how much time it takes to execute any command unless we have the following information: (i) the machine we are executing on: (ii) its machine language instruction set; (iii) the time required by each machine instruction; (iv) the translation a compiler will make from the source to the machine language. It is possible to determine these figures by choosing a real machine and an existing compiler. Another approach would be to define a hypothetical machine (with imaginary execution times), but make the times reasonably close to those of existing hardware so that resulting figures would be representative. Neither of these alternatives seems attractive. In both cases the exact times we would determine would not apply to many machines or to any machine. Also, there would be the problem of the compiler, which could vary from machine to machine. Moreover, it is often difficult to get reliable timing figures because of clock limitations and a multi-programming or time sharing environment. Finally, the difficulty of learning another machine language outweighs the advantage of finding "exact" fictitious times. All these considerations lead us to limit our goals for an a priori analysis. Instead, we will concentrate on developing only the frequency count for all statements. The anomalies of machine configuration and language will be lumped together when we do our experimental studies. Parallelism will not be considered. Consider the three examples of Figure 1.4 below.. for i 1 to n do. for i 1 to n do. for j 1 to n do x x + l x x + 1. x x + 1. end. end file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (25 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com end (a) (b) (c) Figure 1.4: Three simple programs for frequency counting. In program (a) we assume that the statement x x + 1 is not contained within any loop either explicit or implicit. Then its frequency count is one. In program (b) the same statement will be executed n times and in program (c) n2 times (assuming n 1). Now 1, n, and n2 are said to be different and increasing orders of magnitude just like 1, 10, 100 would be if we let n = 10. In our analysis of execution we will be concerned chiefly with determining the order of magnitude of an algorithm. This means determining those statements which may have the greatest frequency count. To determine the order of magnitude, formulas such as often occur. In the program segment of figure 1.4(c) the statement x x + 1 is executed Simple forms for the above three formulas are well known, namely, In general To clarify some of these ideas, let us look at a simple program for computing the n-th Fibonacci number. The Fibonacci sequence starts as 0, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55,... Each new term is obtained by taking the sum of the two previous terms. If we call the first term of the sequence F0 then F0 = 0, F1 = 1 and in general Fn = Fn-1 + Fn-2, n 2. The program on the following page takes any non-negative integer n and prints the value Fn. file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (26 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com 1 procedure FIBONACCI 2 read (n) 3-4 if n < 0 then [print ('error'); stop] 5-6 if n = 0 then [print ('0'); stop] 7-8 if n = 1 then [print ('1'); stop] 9 fnm2 0; fnm 1 1 10 for i 2 to n do 11 fn fnm1 + fnm2 12 fnm2 fnm1 13 fnm1 fn 14 end 15 print (fn) 16 end FIBONACCI The first problem in beginning an analysis is to determine some reasonable values of n. A complete set would include four cases: n < 0, n = 0, n = 1 and n > 1. Below is a table which summarizes the frequency counts for the first three cases. Step n < 0 n = 0 n = 1 -------------------------- 1 1 1 1 2 1 1 1 3 1 1 1 file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (27 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com 4 1 0 0 5 0 1 1 6 0 1 0 7 0 0 1 8 0 0 1 9-15 0 0 0 These three cases are not very interesting. None of them exercises the program very much. Notice, though, how each if statement has two parts: the if condition and the then clause. These may have different execution counts. The most interesting case for analysis comes when n > 1. At this point the for loop will actually be entered. Steps 1, 2, 3, 5, 7 and 9 will be executed once, but steps 4, 6 and 8 not at all. Both commands in step 9 are executed once. Now, for n 2 how often is step 10 executed: not n - 1 but n times. Though 2 to n is only n - 1 executions, remember that there will be a last return to step 10 where i is incremented to n + 1, the test i > n made and the branch taken to step 15. Thus, steps 11, 12, 13 and 14 will be executed n - 1 times but step 10 will be done n times. We can summarize all of this with a table. Step Frequency Step Frequency -------------------------------- 1 1 9 2 2 1 10 n 3 1 11 n-1 4 0 12 n-1 5 1 13 n-1 6 0 14 n-1 7 1 15 1 file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (28 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com 8 0 16 1 Figure 1.5: Execution Count for Computing Fn Each statement is counted once, so step 9 has 2 statements and is executed once for a total of 2. Clearly, the actual time taken by each statement will vary. The for statement is really a combination of several statements, but we will count it as one. The total count then is 5n + 5. We will often write this as O(n), ignoring the two constants 5. This notation means that the order of magnitude is proportional to n. The notation f(n) = O(g(n)) (read as f of n equals big-oh of g of n) has a precise mathematical definition. Definition: f(n) = O(g(n)) iff there exist two constants c and no such that |f(n)| c|g(n)| for all n no. f(n) will normally represent the computing time of some algorithm. When we say that the computing time of an algorithm is O(g(n)) we mean that its execution takes no more than a constant times g(n). n is a parameter which characterizes the inputs and/or outputs. For example n might be the number of inputs or the number of outputs or their sum or the magnitude of one of them. For the Fibonacci program n represents the magnitude of the input and the time for this program is written as T(FIBONACCI) = O(n). We write O(1) to mean a computing time which is a constant. O(n) is called linear, O(n2) is called quadratic, O(n3) is called cubic, and O(2n) is called exponential. If an algorithm takes time O(log n) it is faster, for sufficiently large n, than if it had taken O(n). Similarly, O(n log n) is better than O(n2) but not as good as O(n). These seven computing times, O(1), O(log n), O(n), O(n log n), O(n2), O(n3), and O (2n) are the ones we will see most often throughout the book. If we have two algorithms which perform the same task, and the first has a computing time which is O (n) and the second O(n2), then we will usually take the first as superior. The reason for this is that as n increases the time for the second algorithm will get far worse than the time for the first. For example, if the constant for algorithms one and two are 10 and 1/2 respectively, then we get the following table of computing times: n 10n n2/2 ----------------- 1 10 1/2 5 50 12-1/2 10 100 50 file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (29 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com 15 150 112-1/2 20 200 200 25 250 312-1/2 30 300 450 For n 20, algorithm two had a smaller computing time but once past that point algorithm one became better. This shows why we choose the algorithm with the smaller order of magnitude, but we emphasize that this is not the whole story. For small data sets, the respective constants must be carefully determined. In practice these constants depend on many factors, such as the language and the machine one is using. Thus, we will usually postpone the establishment of the constant until after the program has been written. Then a performance profile can be gathered using real time calculation. Figures 1.6 and 1.7 show how the computing times (counts) grow with a constant equal to one. Notice how the times O(n) and O(n log n) grow much more slowly than the others. For large data sets, algorithms with a complexity greater than O(n log n) are often impractical. An algorithm which is exponential will work only for very small inputs. For exponential algorithms, even if we improve the constant, say by 1/2 or 1/3, we will not improve the amount of data we can handle by very much. Given an algorithm, we analyze the frequency count of each statement and total the sum. This may give a polynomial P(n) = cknk + ck-1 nk-1 +... + c1n + co where the ci are constants, ck 0 and n is a parameter. Using big-oh notation, P(n) = O(nk). On the other hand, if any step is executed 2n times or more the expression c2n + P(n) = O(2n). Another valid performance measure of an algorithm is the space it requires. Often one can trade space for time, getting a faster algorithm but using more space. We will see cases of this in subsequent chapters. Figure 1.6: Rate of Growth of Common Computing Time Functions log2n n nlog2n n2 n3 2n file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (30 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com -------------------------------------------------------- 0 1 0 1 1 2 1 2 2 4 8 4 2 4 8 16 64 16 3 8 24 64 512 256 4 16 64 256 4096 65536 5 32 160 1024 32768 2, 147, 483, 648 Figure 1.7: Values for Computing Functions We end this chapter with a problem from recreational mathematics which uses many of the SPARKS features that have been discussed. A magic square is an n x n matrix of the integers 1 to n2 such that the sum of every row, column and diagonal is the same. For example, if n = 5 we have 15 8 1 24 17 16 14 7 5 23 22 20 13 6 4 3 21 19 12 10 9 2 25 18 11 where the common sum is 65. When n is odd H. Coxeter has given a simple rule for generating a magic square: "Start with 1 in the middle of the top row; then go up and left assigning numbers in increasing order to empty squares; if you fall off the square imagine the same square as tiling the plane and continue; if a square is occupied, move down instead and continue." The magic square above was formed using this rule. We now write a SPARKS program for creating an n X n magic square for n odd. file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (31 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com procedure MAGIC(square, n) //for n odd create a magic square which is declared as an array// //square (0: n - 1, 0: n - 1)// //(i,j) is a square position. 2 key n2 is integer valued.// if n is even then [print ('input error'); stop] SQUARE 0 square (0,(n - 1)/2) 1; //store 1 in middle of first row// key 2; i 0; j (n - 1)/2 //i,j are current position// while key n2 do (k,l) ((i - 1) mod n, (j - 1) mod n) //look up and left// if square (k,l) 0 then i (i + 1) mod n //square occupied, move down// else (i,j) (k,l) //square (k,l) needs to be assigned// square (i,j) key //assign it a value// key key + 1 end print (n, square) //output result// end MAGIC MAGIC is a complete SPARKS procedure. The statement (i,j) (k,l) is a shorthand way of writing i k; j l. It emphasizes that the variables are thought of as pairs and are changed as a unit. The file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (32 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com reserved word mod computes the nonnegative remainder and is a built in function. The magic square is represented using a two dimensional array having n rows and n column. For this application it is convenient to number the rows (and columns) from zero to n - 1 rather than from one to n. Thus, when the program ''falls off the square'' the mod operator sets i and/or j back to zero or n - 1. The while loop is governed by the variable key which is an integer variable initialized to 2 and increased by one each time through the loop. Thus each statement within the while loop will be executed no more than n2 - 1 times and hence the computing time for MAGIC is O(n2). Since there are n2 positions in which the algorithm must place a number, we see that O(n2) is the best bound an algorithm could have. REFERENCES For a discussion of algorithms and how to analyze them see The Art of Computer Programming: Fundamental Algorithms, by D. E. Knuth, vol. 1, chapter 1, 2-nd edition, Addison-Wesley, 1973. For a discussion of good programming techniques see Structured Programming by O. J. Dahl, E. W. Dijkstra, and C. A. R. Hoare, Academic Press, 1972. The Elements of Programming Style by B. W. Kernighan and P. J. Plauger, McGraw-Hill, 1974. ACM Computing Surveys, Special Issue: Programming, vol. 6, no. 4, December, 1974. For a discussion of tools and procedures for developing very large software systems see Practical Strategies for Developing Large Software Systems, by E. Horowitz, Addison-Wesley, May, 1975. For a discussion of the more abstract formulation of data structures see "Toward an understanding of data structures" by J. Earley, CACM, vol. 14, no. 10, October, 1971, pp. 617-627. "Another look at data," by G. Mealy, Proc. AFIPS Fall Joint Computer Conference, vol. 31, 1967, pp. 525-534. For a further discussion of program proving see file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (33 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com "Assigning meaning to programs," by R. W. Floyd, Proc. of a Symposium in Applied Mathematics, vol. 19, J. T. Schwartz, ed., American Mathematical Society, Providence, 1967, pp. 19-32. "An interactive program verification system," by D. I. Good, R. L. London, W. W. Bledsoe, IEEE Transactions on Software Engineering, SE-1, vol. 1, March, 1975, pp. 59-67. EXERCISES 1. Look up the word algorithm or its older form algorism in the dictionary. 2. Consider the two statements: (i) Is n = 2 the largest value of n for which there exists positive integers x, y and z such that xn + yn = zn has a solution; (ii) Store 5 divided by zero into X and go to statement 10. Both do not satisfy one of the five criteria of an algorithm. Which criteria do they violate? 3. Describe the flowchart in figure 1.1 by using a combination of SPARKS and English. Can you do this without using the go to? Now make it into an algorithm. 4. Discuss how you would actually represent the list of name and telephone number pairs in a real machine. How would you handle people with the same last name. 5. Write FORTRAN equivalents of the while, repeat-until, loop-forever and for statements of SPARKS. 6. Can you think of a clever meaning for S.P.A.R.K.S.? Concentrate on the letter K first. 7. Determine the frequency counts for all statements in the following two SPARKS program segments: 1 for i 1 to n 1 i 1 2 for j l to i 2 while i n do 3 for k 1 to j 3 x x + 1 4 x x + 1 4 i i +1 5 end 5 end 6 end file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (34 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com 7 end (a) (b) 8. Horner's Rule is a means for evaluating a polynomial A(x) = anxn + an-1xn-1 +... +a1 x + a0 at a point x0 using a minimum number of multiplications. The rule is: A(x) = (... ((an.x0 + an-1)x0 +... + a1)x0 + a0 Write a SPARKS program to evaluate a polynomial using Horner's Rule. Determine how many times each statement is executed. 9. Given n boolean variables x1,..., xn we wish to print all possible combinations of truth values they can assume. For instance, if n = 2, there are four possibilities: true, true; true, false; false, true; false, false. Write a SPARKS program to accomplish this and do a frequency count. 10. Compare the two functions n2 and 2n/4 for various values of n. Determine when the second becomes larger than the first. 11. Write a SPARKS program which prints out the integer values of x, y, z in nondecreasing order. What is the computing time of your method? 12. Write a SPARKS procedure which searches an array A (1: n) for the element x. If x occurs, then set j to its position in the array else set j to zero. Try writing this without using the go to statement. 13. One useful facility we might add to SPARKS is the ability to manipulate character strings. If x, y are variables of type character, then we might like to implement the procedures: (i) z CONCAT(x,y) which concatenates a copy of string y to the end of a copy of string x and assigns the resulting string to z. Strings x and y remain unchanged. (ii) z SUBSTR(x,i,j) which copies to z the i-th to the j-th character in string x with appropriate definitions for j = 0, i > j, etc. String x is unchanged. (iii) z INDEX(x,y) which searches string x for the first occurrence of string y and sets z to its starting position in x or else zero. Implement these procedures using the array facility. 14. Write a SPARKS procedure which is given an argument STRING, whose value is a character string file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (35 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com of length n. Copy STRING into the variable FILE so that every sequence of blanks is reduced to a single blank. The last character of STRING is nonblank. 15. Design a program that counts the number of occurrences of each character in the string STRING of length n. Represent your answer in the array ANS(1:k,1:2) where ANS(i,l) is the i-th character and ANS (i,2) is the number of times it occurs in STRING. 16. Trace the action of the procedure below on the elements 2, 4, 6, 8, 10, 12, 14, 16, 18, 20 searching for l, 3, 13 and 21. i 1; j n repeat k (i + j)/2 if A(k) x then i k + 1 else j k - 1 until i > j What is the computing time for this segment in terms of n? 17. Prove by induction: 18. List as many rules of style in programming that you can think of that you would be willing to follow yourself. 19. Using the notation introduced at the end of section 1.1, define the structure Boolean with operations AND, OR, NOT, IMP and EQV (equivalent) using only the if-then-else statement. e.g. NOT (X) :: = if X then false else true. 20. Give a version of a binary search procedure which initializes lower to zero and upper to n + l. 21. Take any version of binary search, express it using assignment, if-then-else and go to and then give file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (36 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com an equivalent recursive program. 22. Analyze the computing time of procedure SORT as given in section 1.3. 23. Write a recursive procedure for computing the binomial coefficient as defined in section 1.3 where. Analyze the time and space requirements of your algorithm. 24. Ackermann's function A(m,n) is defined as follows: This function is studied because it grows very fast for small values of m and n. Write a recursive procedure for computing this function. Then write a nonrecursive algorithm for computing Ackermann's function. 25. (Tower of Hanoi) There are three towers and sixty four disks of different diameters placed on the first tower. The disks are in order of decreasing diameter as one scans up the tower. Monks were reputedly supposed to move the disks from tower 1 to tower 3 obeying the rules: (i) only one disk can be moved at any time; (ii) no disk can be placed on top of a disk with smaller diameter. Write a recursive procedure which prints the sequence of moves which accomplish this task. 26. Write an equivalent recursive version of procedure MAGIC as given in section 1.4. 27. The pigeon hole principle states that if a function f has n distinct inputs but less than n distinct outputs then there exists two inputs a, b such that a b and f(a) = f(b). Give an algorithm which finds the values a, b for which the range values are equal. 28. Given n, a positive integer determine if n is the sum of all of its divisors; i.e. if n is the sum of all t such that 1 t < n and t divides n. 29. Consider the function F(x) defined by F(x) if even(x) then x/2 else F(F(3x + 1)) Prove that F(x) terminates for all integers x. (Hint: consider integers of the form (2i + 1) 2k - 1 and use induction.) 30. If S is a set of n elements the powerset of S is the set of all possible subsets of S. For example if S = (a,b,c,) then POWERSET(S) = {( ), (a), (b), (c), (a,b), (a,c), (b,c), (a,b,c)}. Write a recursive procedure to compute powerset (S). file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (37 of 38)7/3/2004 3:56:36 PM Fundamentals: CHAPTER 1: INTRODUCTION www.itdevelopteam.com Go to Chapter 2 Back to Table of Contents file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrD...ooks_Algorithms_Collection2ed/books/book1/chap01.htm (38 of 38)7/3/2004 3:56:36 PM  Fundamentals: CHAPTER 2: ARRAYS www.itdevelopteam.com CHAPTER 2: ARRAYS 2.1 AXIOMATIZATION It is appropriate that we begin our study of data structures with the array. The array is often the only means for structuring data which is provided in a programming language. Therefore it deserves a significant amount of attention. If one asks a group of programmers to define an array, the most often quoted saying is: a consecutive set of memory locations. This is unfortunate because it clearly reveals a common point of confusion, namely the distinction between a data structure and its representation. It is true that arrays are almost always implemented by using consecutive memory, but not always. Intuitively, an array is a set of pairs, index and value. For each index which is defined, there is a value associated with that index. In mathematical terms we call this a correspondence or a mapping. However, as computer scientists we want to provide a more functional definition by giving the operations which are permitted on this data structure. For arrays this means we are concerned with only two operations which retrieve and store values. Using our notation this object can be defined as: structure ARRAY(value, index) declare CREATE( ) array RETRIEVE(array,index) value STORE(array,index,value) array; for all A array, i,j index, x value let RETRIEVE(CREATE,i) :: = error RETRIEVE(STORE(A,i,x),j) :: = if EQUAL(i,j) then x else RETRIEVE(A,j) end end ARRAY The function CREATE produces a new, empty array. RETRIEVE takes as input an array and an index, and either returns the appropriate value or an error. STORE is used to enter new index-value pairs. The file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrDo...Books_Algorithms_Collection2ed/books/book1/chap02.htm (1 of 37)7/3/2004 4:01:13 PM  Fundamentals: CHAPTER 2: ARRAYS www.itdevelopteam.com second axiom is read as "to retrieve the j-th item where x has already been stored at index i in A is equivalent to checking if i and j are equal and if so, x, or search for the j-th value in the remaining array, A." This axiom was originally given by J. McCarthy. Notice how the axioms are independent of any representation scheme. Also, i and j need not necessarily be integers, but we assume only that an EQUAL function can be devised. If we restrict the index values to be integers, then assuming a conventional random access memory we can implement STORE and RETRIEVE so that they operate in a constant amount of time. If we interpret the indices to be n-dimensional, (i1,i2,...,in), then the previous axioms apply immediately and define n-dimensional arrays. In section 2.4 we will examine how to implement RETRIEVE and STORE for multi-dimensional arrays using consecutive memory locations. 2.2 ORDERED LISTS One of the simplest and most commonly found data object is the ordered or linear list. Examples are the days of the week (MONDAY, TUESDAY, WEDNESDAY, THURSDAY, FRIDAY, SATURDAY, SUNDAY) or the values in a card deck (2, 3, 4, 5, 6, 7, 8, 9, 10, Jack, Queen, King, Ace) or the floors of a building (basement, lobby, mezzanine, first, second, third) or the years the United States fought in World War II (1941, 1942, 1943, 1944, 1945). If we consider an ordered list more abstractly, we say that it is either empty or it can be written as (a1,a2,a3,...,an) where the ai are atoms from some set S. There are a variety of operations that are performed on these lists. These operations include: file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrDo...Books_Algorithms_Collection2ed/books/book1/chap02.htm (2 of 37)7/3/2004 4:01:13 PM  Fundamentals: CHAPTER 2: ARRAYS www.itdevelopteam.com (i) find the length of the list, n; (ii) read the list from left-to-right (or right-to-left); (iii) retrieve the i-th element, ; (iv) store a new value into the i-th position, ; (v) insert a new element at position causing elements numbered i,i + 1,...,n to become numbered i + 1,i + 2,...,n + 1; (vi) delete the element at position causing elements numbered i + 1,...,n to become numbered i,i + 1,...,n - 1. See exercise 24 for a set of axioms which uses these operations to abstractly define an ordered list. It is not always necessary to be able to perform all of these operations; many times a subset will suffice. In the study of data structures we are interested in ways of representing ordered lists so that these operations can be carried out efficiently. Perhaps the most common way to represent an ordered list is by an array where we associate the list element ai with the array index i. This we will refer to as a sequential mapping, because using the conventional array representation we are storing ai and ai + 1 into consecutive locations i and i + 1 of the array. This gives us the ability to retrieve or modify the values of random elements in the list in a constant amount of time, essentially because a computer memory has random access to any word. We can access the list element values in either direction by changing the subscript values in a controlled way. It is only operations (v) and (vi) which require real effort. Insertion and deletion using sequential allocation forces us to move some of the remaining elements so the sequential mapping is preserved in its proper form. It is precisely this overhead which leads us to consider nonsequential mappings of ordered lists into arrays in Chapter 4. Let us jump right into a problem requiring ordered lists which we will solve by using one dimensional arrays. This problem has become the classical example for motivating the use of list processing techniques which we will see in later chapters. Therefore, it makes sense to look at the problem and see why arrays offer only a partially adequate solution. The problem calls for building a set of subroutines which allow for the manipulation of symbolic polynomials. By "symbolic," we mean the list of coefficients and exponents which accompany a polynomial, e.g. two such polynomials are A(x) = 3x2 + 2x + 4 and B(x) = x4 + 10x3 + 3x2 + 1 file:///C|/E%20Drive%20Data/My%20Books/Algorithm/DrDo...Books_Algorithms_Collection2ed/books/book1/chap02.htm (3 of 37)7/3/2004 4:01:13 PM  Fundamentals: CHAPTER 2: ARRAYS www.itdevelopteam.com For a start, the capabilities we would include are the four basic arithmetic operations: addition, subtraction, multiplication, and division. We will also need input and output routines and some suitable format for preparing polynomials as input. The first step is to consider how to define polynomials as a computer structure. For a mathematician a polynomial is a sum of terms where each term has the form axe; x is the variable, a is the coefficient and e is the exponent. However this is not an appropriate definition for our purposes. When defining a data object