ML语言
编程范型 | 多范型:函数式,指令式 |
---|---|
設計者 | 罗宾·米尔纳及爱丁堡大学其他人 |
发行时间 | 1973年 |
型態系統 | 类型推论,静态类型,强类型 |
衍生副語言 | |
Standard ML, OCaml | |
啟發語言 | |
ISWIM[1],PAL[2],POP-2[1],GEDANKEN[1] | |
影響語言 | |
Clojure、Coq、Cyclone、C++、Elm[3]、Futhark[4]、F#、F*、Haskell、Idris、Lazy ML[5]、Miranda、Nemerle[6]、OCaml、Opa、Rust、Scala、Standard ML、Ur[7] | |
ML(Meta Language:元语言),是一个函数式、指令式的通用的编程语言,它著称于使用了多态的Hindley–Milner类型推论[8]。ML能自动的指定多数表达式的类型,不要求显式的类型标注,而且能够确保类型安全,已经正式证明了有良好类型的ML程序不会导致运行时间类型错误[8]。
ML提供了对函数实际参数的模式匹配、垃圾回收、指令式编程、传值调用和柯里化。它被大量的用于编程语言研究之中,并且是全面规定了的和使用形式语义验证了的少数语言之一。它的类型和模式匹配使得它非常适合并且经常用于在其他形式语言上进行操作,比如在编译器构造、自动化定理证明和形式验证中。
历史
[编辑]在1970年代早期,ML由爱丁堡大学的罗宾·米尔纳及他人研制出来[1],用于在LCF定理证明器中开发证明策略[9]。LCF的语言是“PPλ”,它是一阶逻辑演算与有类型的多态λ演算的结合,以ML作为元语言。ML的语法从ISWIM及其扩展实现PAL得到了启发[2],LCF ML运行在DEC-10/TOPS-10主机的Stanford LISP 1.6下[10]。
在1980年,Luca Cardelli于爱丁堡大学的VAX/VMS系统上开发了ML编译器,它被称为“VAX ML”[11],以区别于LCF版本的“DEC-10 ML”[12]。VAX ML的编译器和运行时间系统二者,都是用Pascal书写而建立在“函数抽象机器”(FAM)之上[13]。在1982年,爱丁堡大学的Kevin Mitchell,用VAX ML重写了VAX ML编译器,随即John Scott和Alan Mycroft加入了开发,在又进行一系列重写改进之后,新的编译器被称为“Edinburgh ML”[14]。
在1981年,INRIA的Gérard Huet,将最初的LCF ML适配到Multics系统的Maclisp下,并且增加了编译器[15]。这个实现被描述于INRIA内部文档“ML手册”之中[16],它被开发者自称为“Le_ML”[17]。剑桥大学的Lawrence Paulson用它开发了基于Franz Lisp的Cambridge LCF[18],进而剑桥大学的Michael J. C. Gordon又用它开发了基于Common Lisp的第一版的HOL[19]。
在1983年,Robin Milner由两个动机驱使开始重新设计ML[20]。其一是爱丁堡大学的Rod Burstall及其小组在规定上的工作,具体化为规定语言CLEAR[21],和表达可执行规定的函数式语言HOPE[22]。这项工作与ML有关的是两方面成果:首先,HOPE拥有优雅的编程特征,特别是模式匹配[23],和子句式函数定义[24];其次,是使用在接口中的签名,进行规定的模块化构造的想法。其二是Luca Cardelli在VAX ML上的工作,通过增加命名记录和可变类型,扩展了ML中数据类型的品目[25]。
在1984年,贝尔实验室的David MacQueen提出了对Standard ML的模块系统的设计[26]。在Standard ML的持续设计期间[27],Edinburgh ML被渐进的修改,充当了Standard ML的近似原型实现[28]。在1986年,普林斯顿大学的Andrew Appel和贝尔实验室的David MacQueen,以Edinburgh Standard ML作为起步开发环境[29],开始了专注于生成高质量机器代码的Standard ML of New Jersey的活跃开发[30]。
在1990年,Robin Milner、Mads Tofte和Robert Harper制定的Standard ML的正式规定《The Definition of Standard ML》最终完成[31];在1997年,这个标准规定增补了David MacQueen为作者并进行了修订[32]。在1989年,Mads Tofte、Nick Rothwell和David N. Turner于爱丁堡大学开始开发ML Kit编译器,为了强调清晰性而非高效,将标准定义直接转译成一组Standard ML模块;在1992年和1993年期间,主要通过爱丁堡大学的Nick Rothwell和哥本哈根大学计算机科学系(DIKU)的Lars Birkedal的工作[33],ML Kit第一版完成并开源发行[34]。
在1987年,INRIA的Ascánder Suárez,基于巴黎第七大学的Guy Cousineau的“范畴抽象机器”(CAM)[35],利用Le Lisp的运行时间系统重新实现了Le_ML,并正式命名为Caml[15]。在1990年和1991年,INRIA的Xavier Leroy基于用C实现的字节码解释器[36],利用Damien Doligez提供的内存管理系统重新实现了Caml,并称其为Caml Light[37]。在1995年,Xavier Leroy又增加了机器代码编译器和高层模块系统[38],这个版本也称为“Caml Special Light”。在1996年,INRIA的Didier Rémy和Jérôme Vouillon,向Caml Special Light增加了面向对象特征[39],从而形成了OCaml[40]。
今天ML家族的两个主要的方言是Standard ML和OCaml。ML的实力大多被用于语言设计和操作,比如建构编译器、分析器、定理证明器,但是它作为通用语言也被用于生物信息和财务系统等领域。ML确立了静态类型函数式编程范型,从而在编程语言历史上占有显要地位,它的思想在影响了众多的语言,例如Haskell、Nemerle[6]、ATS和Elm[3]。
解释与编译
[编辑]ML代码片段很容易通过将其录入到“顶层”来研习,它也叫作读取﹣求值﹣输出循环或REPL。这是打印结果或定义的表达式的推论类型的交互式会话。很多SML实现提供交互式REPL,比如SML/NJ:
$ sml Standard ML of New Jersey v110.79 [built: Sun Nov 6 06:43:07 2022] -
可以在提示符-
后键入代码。例如计算1 + 2 * 3
:
- 1 + 2 * 3; val it = 7 : int
顶层推论这个表达式的类型为int
并给出结果7
。如果输入不完全则打印第二提示符=
,这时经常可以用;
终结输入。it
是给未指定变量的表达式的标准变量。输入control-C可返回解释器顶层,输入control-D可退出解释器。可以使用第三方工具rlwrap
运行SML/NJ解释器,它处理用户输入并提供readline
的行编辑、持久历史和补全功能。
下面是hello, world!的程序代码在SML/NJ解释器下执行的结果:
- val _ = print "Hello, world!\n"; Hello, world!
和使用MLton编译器进行编译执行的例子:
$ echo 'print "Hello, world!\n";' > hello-world.sml $ mlton hello-world.sml $ ./hello-world Hello, world!
核心语言
[编辑]不同于纯函数式编程语言,ML是兼具一些指令式特征的函数式编程语言。ML的特征包括:传值调用的求值策略,头等函数,带有垃圾收集的自动内存管理,参数多态,静态类型,类型推论,代数数据类型,模式匹配和例外处理。不同于Haskell,ML与大多数编程语言一样使用及早求值。
用ML书写的程序构成自要被求值的表达式,而非语句或命令,尽管一些表达式返回一个平凡的unit
值并且只为其副作用而求值。以下章节介绍采用Standard ML的语法和语义。
函数
[编辑]就像所有的函数式语言一样,ML的关键特征是函数,它被用于进行抽象。例如阶乘函数用纯ML可表达为:
fun fac 0 = 1 | fac n = n * fac (n - 1)
这里将阶乘描述为递归函数,具有一个单一的终止基础情况。它类似于在数学教科书见到的阶乘描述。多数ML代码在设施和语法上类似于数学。
凭借类型推论编译器能推导出,fac
接受整数0
作为实际参数,则形式参数n
也是整数类型int
,而fac 0
的结果是整数1
,则函数fac
的结果也是整数类型。函数fac
接受一个整数的形式参数并返回一个整数结果,它作为一个整体从而有着“从整数到整数的函数”类型int -> int
。函数及其形式参数的"类型"还可以用类型标注(annotation)来描述,使用E : t
表示法,它可以被读作表达式E
有类型t
,它是可选也可忽略的。使用类型标注,这个例子可重写为如下:
fun fac 0 = 1 | fac (n : int) : int = n * fac (n - 1)
这个函数还依赖于模式匹配,这是ML语言的重要部份。注意函数形式参数不必须在圆括号内但要用空格分隔。当一个函数的实际参数是0
,它将返回整数1
。对于所有其他情况,尝试第二行。这是一个递归,并且再次执行这个函数直到达到基础情况。它可以使用case
表达式重写为:
fun fac n = case n of 0 => 1 | n => n * fac (n - 1)
这里case
介入了模式和对应表达式的序列。它还可以重写为将标识符绑定到lambda函数:
val rec fac = fn 0 => 1 | n => n * fac (n - 1)
这里的关键字val
介入了标识符到值的绑定,fn
介入了匿名函数的定义,它可以用在fun
的位置上,但使用=>
算符而非=
。绑定到递归的匿名函数需要使用rec
关键字来指示。
通过将主要工作写入尾递归风格的内部迭代函数,借助于语言编译或解释系统进行的尾调用优化,这个函数可以得以增进性能,它的调用栈不需要随函数调用数目而成比例的增长。对这个函数可以采用向内部函数增加额外的“累加器”形式参数acc
来实现:
fun fac n = let fun loop (0, acc) = acc | loop (n, acc) = loop (n - 1, n * acc) in loop (n, 1) end
let
表达式的值是在in
和end
之间表达式的值。这个递归函数的实现不保证能够终止,因为负数实际参数会导致递归调用的无穷降链条件。更健壮的实现会在递归前检查实际参数为非负数,并在有问题的情况,即n
是负数的时候,启用例外处理:
fun fac n = let fun loop (0, acc) = acc | loop (n, acc) = loop (n - 1, n * acc) in if (n < 0) then raise Fail "negative argument" else loop (n, 1) end
类型
[编辑]有一些基本类型可以当作是“内建”的,因为它们是在Standard ML标准基本库中预先定义的,并且语言为它们提供文字表示法,比如34
是整数,而"34"
是字符串。一些最常用的基本类型是:
int
整数,比如3
或~12
。注意波浪号~表示负号。real
浮点数,比如4.2
或~6.4
。Standard ML不隐含的提升整数为浮点数,因此表达式2 + 5.67
是无效的。string
字符串,比如"this is a string"
或""
(空串)。char
字符,比如#"y"
或#"\n"
(换行符)。bool
布尔值,它是要么true
要么false
。产生bool
值的有比较算符=
、<>
、>
、>=
、<
、<=
,逻辑函数not
和短路求值的中缀算符andalso
、orelse
。
包括上述基本类型的各种类型可以用多种方式组合。一种方式是元组,它是值的有序集合,比如表达式(1, 2)
的类型是int * int
,而("foo", false)
的类型是string * bool
。可以使用#1 ("foo", false)
这样的语法来提取元组的指定次序的成员。
有0元组()
,它的类型指示为unit
。但是没有1元组,或者说在例子(1)
和1
之间没有区别,都有类型int
。元组可以嵌套,(1, 2, 3)
不同于((1, 2), 3)
和(1, (2, 3))
二者。前者的类型是int * int * int
,其他两个的类型分别是(int * int) * int
和int * (int * int)
。
组合值的另一种方式是记录。记录很像元组,除了它的成员是有名字的而非有次序的,例如{a = 5.0, b = "five"}
有类型{a : real, b : string}
,这同于类型{b : string, a : real}
。可以使用#a {a = 5.0, b = "five"}
这样的语法来选取记录的指定名字的字段。
Standard ML中的函数只接受一个值作为参数,而不是参数的一个列表,可以使用上述元组模式匹配来实际上传递多个参数。函数还可以返回一个元组。例如:
fun sum (a, b) = a + b fun average pair = sum pair div 2 infix averaged_with fun a averaged_with b = average (a, b) val c = 3 averaged_with 7 fun int_pair (n : int) = (n, n)
这里第一段创建了两个类型是int * int -> int
的函数sum
和average
。第二段创建中缀算子averaged_with
,接着定义它为类型int * int -> int
的函数。最后的int_pair
函数的类型是int -> int * int
。
下列函数是多态类型的一个例子:
fun pair x = (x, x)
编译器无法推论出的pair
的特殊化了的类型,它可以是int -> int * int
、real -> real * real
甚至是(int * real -> string) -> (int * real -> string) * (int * real -> string)
。在Standard ML中,它可以简单指定为多态类型'a -> 'a * 'a
,这里的'a
(读作“alpha”)是一个类型变量,指示任何可能的类型。在上述定义下,pair 3
和pair "x"
都是有良好定义的,分别产生(3, 3)
和("x", "x")
。
SML标准基本库包括了重载标识符:+
、-
、*
、div
、mod
、/
、~
、abs
、<
、>
、<=
、>=
。标准基本库提供了多态函数:!
、:=
、o
、before
、ignore
。中缀运算符可以有从缺省最低0
到最高9
的任何运算符优先级。标准基本库提供了如下内建中缀规定:
infix 7 * / div mod infix 6 + - ^ infixr 5 :: @ infix 4 = <> > >= < <= infix 3 := o infix 0 before
等式类型
[编辑]算符=
和<>
分别被定义为多态的等式和不等式。=
它确定两个值是否相等,有着类型''a * ''a -> bool
。这意味着它的两个运算数必须有相同的类型,这个类型必须是等式类型(eqtype
)。上述基本类型中除了real
之外,int
、real
、string
、char
和bool
都是等式类型。
例如:3 = 3
、"3" = "3"
、#"3" = #"3"
和true = true
,都是有效的表达式并求值为true
;而3 = 4
是有效表达式并求值为false
,3.0 = 3.0
是无效表达式而被编译器拒绝。这是因为IEEE浮点数等式打破了ML中对等式的一些要求。特别是nan
不等于自身,所以关系不是自反的。
元组和记录类型是等式类型,当时且仅当它的每个成员类型是等式类型;例如int * string
、{b : bool, c : char}
和unit
是等式类型,而int * real
和{x : real}
不是。函数类型永远不是等式类型,因为一般情况下不可能确定两个函数是否等价。
类型声明
[编辑]类型声明或同义词(synonym)使用type
关键字来定义。下面是给在平面中的点的类型声明,计算两点间距离,和通过海伦公式计算给定角点的三角形的面积的函数。
type loc = real * real fun heron (a: loc, b: loc, c: loc) = let fun dist ((x0, y0), (x1, y1)) = let val dx = x1 - x0 val dy = y1 - y0 in Math.sqrt (dx * dx + dy * dy) end val ab = dist (a, b) val bc = dist (b, c) val ac = dist (a, c) val s = (ab + bc + ac) / 2.0 in Math.sqrt (s * (s - ab) * (s - bc) * (s - ac)) end
数据类型
[编辑]Standard ML提供了对代数数据类型(ADT)的强力支持。一个ML数据类型可以被当作是元组的不交并(积之和)。数据类型使用datatype
关键字来定义,比如:
datatype int_or_string = INT of int | STRING of string | NEITHER
这个数据类型声明建立一个全新的数据类型int_or_string
,还有一起的新构造子(一种特殊函数或值)INT
、STRING
和NEITHER
;每个这种类型的值都是要么INT
具有一个整数,要么STRING
具有一个字符串,要么NEITHER
。写为如下这样:
val i = INT 3 val s = STRING "qq" val n = NEITHER val INT j = i
这里最后的一个声明通过模式匹配,将变量j
绑定到变量i
所绑定的整数3
。val 模式 = 表达式
是绑定的一般形式,它是良好定义的,当且仅当模式和表达式有相同的类型。
数据类型可以是多态的:
datatype 'a pair = PAIR of 'a * 'a
这个数据类型声明建立一个新的类型'a pair
家族,比如int pair
,string pair
等等。
数据类型可以是递归的:
datatype int_list = EMPTY | INT_LIST of int * int_list
这个数据类型声明建立一个新类型int_list
,这个类型的每个值是要么EMPTY
(空列表),要么是一个整数和另一个int_list
的接合。
通过datatype
创建的类型是等式类型,如果它的所有变体,要么是没有参数的空构造子,要么是有等式类型参数的构造子,并且在多态类型情况下所有类型参数也都是等式类型。递归类型在有可能情况下是等式类型,否则就不是。
列表
[编辑]基础库提供的复杂数据类型之一是列表list
。它是一个递归的、多态的数据类型,可以等价的定义为:
datatype 'a list = nil | :: of 'a * 'a list
这里的::
是中缀算符,例如3 :: 4 :: 5 :: nil
是三个整数的列表。列表是ML程序中最常用的数据类型之一,语言还为生成列表提供了特殊表示法[3, 4, 5]
。
real list
当然不是等式类型。但是没有int list
不能是等式类型的理由,所以它就是等式类型。注意类型变量不同就是不同的列表类型,比如(nil : int list) = (nil : char list)
是无效的,因为两个表达式有不同的类型,即使它们有相同的值。但是nil = nil
和(nil : int list) = nil
都是有效的。
基本库rev
函数“反转”一个列表中的元素。它的类型是'a list -> 'a list
,就是说它可以接受其元素有任何类型的列表,并返回相同类型的列表。更准确的说,它返回其元素相较于给定列表是反转次序的一个新列表,比如将[ "a", "b", "c" ]
映射成[ "c", "b", "a" ]
。中缀算符@
表示两个列表的串接。
rev
和@
一般被实现为基本库函数revAppend
的简单应用:
fun revAppend ([], l) = l | revAppend (x :: r, l) = revAppend(r, x :: l) fun rev l = revAppend(l, []) fun l1 @ l2 = revAppend(rev l1, l2)
模式匹配
[编辑]Standard ML的数据类型可以轻易的定义和用于编程,在很大程度上是由它的模式匹配,还有多数Standard ML实现的模式穷尽性检查和模式冗余检查。
模式匹配可以在语法上嵌入到变量绑定之中,比如:
val ((m: int, n: int), (r: real, s: real)) = ((2, 3), (2.0, 3.0)) type hyperlink = {protocol: string, address: string, title: string} val url : hyperlink = {title="The Standard ML Basis Library", protocol="https", address="//smlfamily.github.io/Basis/overview.html"} val {protocol=port, address=addr, title=name} = url val x as (fst, snd) = (2, true)
第一个val
绑定了四个变量m
、n
、r
和s
;第二个val
绑定了一个变量url
;第三个val
绑定了三个变量port
、addr
和name
,第四个叫分层模式,绑定了三个变量x
、fst
和snd
。
模式匹配可以在语法上嵌入到函数定义中,比如:
datatype shape = Circle of loc * real (* 中心和弧度 *) | Square of loc * real (* 左上角和边长,坐标轴对齐 *) | Triangle of loc * loc * loc (* 角点 *) fun area (Circle (_, r)) = 3.14 * r * r | area (Square (_, s)) = s * s | area (Triangle (a, b, c)) = heron (a, b, c)
次序在模式匹配中是紧要的:模式按文本次序来依次进行匹配。在特定计算中,使用下划线_
,来省略不需要它的值的子成员,这也叫做通配符(wildcard)模式。所谓的“子句形式”风格的函数定义,这里的模式紧随在函数名字之后出现,只是如下形式的一种语法糖:
fun area shape = case shape of Circle (_, r) => 3.14 * r * r | Square (_, s) => s * s | Triangle (a, b, c) => heron (a, b, c)
模式穷尽性检查将确保数据类型的每个情况都已经顾及到。[a] 如果有遗漏,则产生警告。[b] 如果模式存在冗余,也会导致一个编译时间警告。[c]
高阶函数
[编辑]函数可以接受函数作为实际参数,比如:
fun applyToBoth f x y = (f x, f y)
函数可以产生函数作为返回值,比如:
fun constantFn k = (fn anything => k)
函数可以同时接受和产生函数,比如复合函数:
fun compose (f, g) = (fn x => f (g x))
基本库的函数List.map
,是在Standard ML中最常用的Curry化高阶函数,它在概念上可定义为:
fun map' _ [] = [] | map' f (x :: xs) = f x :: map f xs
在基本库中将函数复合定义为多态中缀算符(f o g)
,map
和fold
高阶函数有较高效率的实现。[d]
例外
[编辑]例外可以使用raise
关键字发起,并通过模式匹配handle
构造来处理:
exception Undefined; fun max [x] = x | max (x :: xs) = let val m = max xs in if x > m then x else m end | max [] = raise Undefined fun main xs = let val msg = (Int.toString (max xs)) handle Undefined => "empty list...there is no max!" in print (msg ^ "\n") end
这里的^
是字符串串接算符。可以利用例外系统来实现非局部退出,例如这个函数所采用技术:
exception Zero; fun listProd ns = let fun p [] = 1 | p (0 :: _) = raise Zero | p (h :: t) = h * p t in (p ns) handle Zero => 0 end
Zero
在0
情况下发起,控制从函数p
中一起出离。
引用
[编辑]初始化基础库还以引用的形式提供了可变的存储。引用ref
可以在某种意义上被认为是如下这样定义的:
datatype 'a ref = ref of 'a
还定义了内建的:=
函数来修改引用的内容,和!
函数来检索引用的内容。阶乘函数可以使用引用定义为指令式风格:
fun factImperative n = let val i = ref n and acc = ref 1 in while !i > 0 do (acc := !acc * !i; i := !i - 1); !acc end
这里使用圆括号对以;
分隔的表达式进行了顺序复合。
可变类型'a ref
是等式类型,即使它的成员类型不是。两个引用被称为是相等的,如果它们标识相同的“ref
单元”,就是说相同的对ref
构造子调用生成的同一个指针。因此(ref 1) = (ref 1)
和(ref 1.0) = (ref 1.0)
都是有效的,并且都求值为false
,因为即使两个引用碰巧指向了相同的值,引用自身是分立的,每个都可以独立于其他而改变。
模块语言
[编辑]模块是ML用于构造大型项目和库的系统。
模块
[编辑]一个模块构成自一个签名(signature)文件和一个或多个结构文件。签名文件指定要实现的API(就像C语言头文件或Java接口文件)。结构实现这个签名(就像C语言源文件或Java类文件)。解释器通过use
命令导入它们。ML的标准库被实现为这种方式的模块。
例如,下列定义一个算术签名:
signature ARITH = sig eqtype t val zero : t val one : t val fromInt: int -> t val fromIntPair : int * int -> t val fromReal : real -> t val succ : t -> t val pred : t -> t val ~ : t -> t val + : t * t -> t val - : t * t -> t val * : t * t -> t val / : t * t -> t val == : t * t -> bool val <> : t * t -> bool val > : t * t -> bool val >= : t * t -> bool val < : t * t -> bool val <= : t * t -> bool end
和这个签名使用有理数的实现:
structure Rational : ARITH = struct datatype t = Rat of int * int val zero = Rat (0, 1) val one = Rat (1, 1) fun fromInt n = Rat (n, 1) fun ineg (a, b) = (~a, b) fun fromIntPair (num, den) = let fun reduced_fraction (numerator, denominator) = let fun gcd (n, 0) = n | gcd (n, d) = gcd (d, n mod d) val d = gcd (numerator, denominator) in if d > 1 then (numerator div d, denominator div d) else (numerator, denominator) end in if num < 0 andalso den < 0 then Rat (reduced_fraction (~num, ~den)) else if num < 0 then Rat (ineg (reduced_fraction (~num, den))) else if den < 0 then Rat (ineg (reduced_fraction (num, ~den))) else Rat (reduced_fraction (num, den)) end val SOME maxInt = Int.maxInt val minPos = 1.0 / (real maxInt) fun fromReal r = let fun convergent (i, f, h_1, k_1, h_2, k_2) = if i <> 0 andalso ((h_1 > (maxInt - h_2) div i) orelse (k_1 > (maxInt - k_2) div i)) then (h_1, k_1) else if f < minPos then (i * h_1 + h_2, i * k_1 + k_2) else convergent (trunc (1.0 / f), Real.realMod (1.0 / f), i * h_1 + h_2, i * k_1 + k_2, h_1, k_1) in if r < 0.0 then Rat (ineg (convergent (trunc (~ r), Real.realMod (~ r), 1, 0, 0, 1))) else Rat (convergent (trunc r, Real.realMod r, 1, 0, 0, 1)) end fun succ (Rat (a, b)) = fromIntPair (a + b, b) fun pred (Rat (a, b)) = fromIntPair (a - b, b) fun add (Rat (a, b), Rat (c, d)) = if b = d then fromIntPair(a + c, b) else fromIntPair (a * d + c * b, b * d) fun sub (Rat (a, b), Rat (c, d)) = if b = d then fromIntPair(a - c, b) else fromIntPair (a * d - c * b, b * d) fun mul (Rat (a, b), Rat (c, d)) = fromIntPair (a * c, b * d) fun div_ (Rat (a, b), Rat (c, d)) = fromIntPair (a * d, b * c) fun gt (Rat (a, b), Rat (c, d)) = if b = d then (a > c) else (a * d) > (b * c) fun lt (Rat (a, b), Rat (c, d)) = if b = d then (a < c) else (a * d) < (b * c) fun neg (Rat (a, b)) = Rat (~a, b) fun eq (Rat (a, b), Rat (c, d)) = ((b = d) andalso (a = c)) fun ~ a = neg a fun a + b = add (a, b) fun a - b = sub (a, b) fun a * b = mul (a, b) fun a / b = div_ (a, b) fun op == (a, b) = eq (a, b) fun a <> b = not (eq (a, b)) fun a > b = gt (a, b) fun a >= b = (gt (a, b) orelse eq (a, b)) fun a < b = lt (a, b) fun a <= b = (lt (a, b) orelse eq (a, b)) end
下面是这个结构的简单用例:
infix == local open Rational in val c = let val a = fromIntPair(2, 3) val b = fromIntPair(4, 6) in a + b end end structure R = Rational R.fromReal(3.245)
Standard ML只允许通过签名函数同实现进行交互,例如不可以直接通过这个代码中的Rat
来建立数据对象。结构块对外部隐藏所有实现细节。这里的:
叫做透明归属(ascription),可以通过所绑定的变量见到此结构的数据内容,与之相对的是:>
,它叫做不透明归属,此结构的数据内容对外完全不可见。比如上面用例有结果:val c = Rat (4,3) : Rational.t
,如果改为不透明归属则有结果:val c = - : Rational.t
。
要用有理数进行实际上的数值计算,需要处理分数形式中分母快速增大导致溢出整数类型大小范围等问题。[e]
函子
[编辑]函子接受指定签名的一个结构作为参数,并返回一个结构作为结果,下面示例的函子能在ARITH
签名的结构上计算移动平均:
signature MOVINGLIST = sig type t structure Arith : sig type t end val size : t -> int val average : t -> Arith.t val fromList : Arith.t list -> t val move : t * Arith.t -> t val expand : t * Arith.t -> t val shrink : t -> t val trunc : t -> t end
functor MovingList (S: ARITH) : MOVINGLIST = struct type t = S.t list * int * S.t structure Arith = S fun size (ml : t) = #2 ml fun average (ml : t) = #3 ml fun fromList (l : S.t list) = let val n = length l val sum = foldl S.+ S.zero l local open S in val m = sum / (fromInt n) end in if (null l) then raise Empty else (l, n, m) end fun move ((l, n, m) : t, new : S.t) = let val old = List.nth (l, n - 1) local open S in val m' = m + (new - old) / (fromInt n) end in (new :: l, n, m') end fun expand ((l, n, m) : t, new : S.t) = let val n' = n + 1; local open S in val m' = m + (new - m) / (fromInt n') end in (new :: l, n', m') end fun shrink ((l, n, m) : t) = let val old = List.nth (l, n - 1) val n' = if (n > 2) then n - 1 else 1 local open S in val m' = m + (m - old) / (fromInt n') end in (l, n', m') end fun trunc ((l, n, m) : t) = (List.take (l, n), n, m) end
和这个函子的简单用例:
structure R = Rational structure MLR = MovingList (Rational) val d = MLR.fromList [R.fromIntPair (4, 5), R.fromIntPair (2, 3)] val d = MLR.expand (d, R.fromIntPair (5, 6)) val d = MLR.move (d, R.fromIntPair (7, 8)) val d = MLR.shrink d val d = MLR.trunc d
这个用例承上节示例,Rational
结构采用了透明归属,有结果如:val d = ([Rat (4,5),Rat (2,3)],2,Rat (11,15)) : MLR.t
。如果它改为不透明归属,则对应结果为:val d = ([-,-],2,-) : MLR.t
。
示例代码
[编辑]下列例子使用了Standard ML的语法和语义。
素数
[编辑]fun prime n = let fun isPrime (l, i) = let fun existsDivisor [] = false | existsDivisor (x :: xs) = if (i mod x) = 0 then true else if (x * x) > i then false else existsDivisor xs in not (existsDivisor l) end fun iterate (acc, i) = if i > n then acc else if isPrime (acc, i) then iterate (acc @ [i], i + 2) else iterate (acc, i + 2) in if n < 2 then [] else iterate ([2], 3) end
基本库find
和exists
函数在不存在符合条件元素的时候会遍历整个列表,[f] 这里采用了特殊化了的existsDivisor
,用以在后续元素都不满足条件时立即结束运算。
fun prime n = let fun dropComposite (acc, [], _, _) = rev acc | dropComposite (acc, l as x :: xs, j, i) = if j > n then List.revAppend (acc, l) else if x < j then dropComposite (x :: acc, xs, j, i) else if x > j then dropComposite (acc, l, j + i, i) else dropComposite (acc, xs, j + i, i) fun init i = let fun loop (l, i) = if i <= 2 then l else loop (i :: l, i - 2) in loop ([], i - (i + 1) mod 2) end fun iterate (acc, []) = rev acc | iterate (acc, l as x :: xs) = if x * x > n then 2 :: List.revAppend (acc, l) else iterate (x :: acc, dropComposite ([], xs, x * x, x * 2)) in if n < 2 then [] else iterate ([], init n) end
这里基于列表的筛法实现符合埃拉托斯特尼的原初算法[41]。筛法还有基于数组的直观实现。[g] 下面是其解释器下命令行运行实例:
- fun printIntList (l: int list) = = print ((String.concatWith " " (map Int.toString l)) ^ "\n"); val printIntList = fn : int list -> unit - val _ = printIntList (prime 100); 2 3 5 7 11 13 17 19 23 29 31 37 41 43 47 53 59 61 67 71 73 79 83 89 97
汉明数
[编辑]正规数是形如的整数,对于非负整数、和,它是可以整除的数。计算升序的正规数的算法经由戴克斯特拉得以流行[42],理查德·汉明最初提出了这个问题,故而这个问题被称为“汉明问题”,这个数列也因而被称为汉明数。Dijkstra计算这些数的想法如下:
- 汉明数的序列开始于数
1
。 - 要加入序列中的数有下述形式:
2h, 3h, 5h
,这里的h
是序列已有的任意的汉明数。 - 因此,可以生成最初只有一个
1
的序列H
,并接着归并序列2H, 3H, 5H
,并以此类推。
示例汉明数程序代码,一般用来展示,确使计算只在需要时进行的纯函数式编程方式[43]。
fun Hamming_number n = let fun merge (p, q) = let fun revMerge (acc, p, q) = if not (null p) andalso not (null q) then if hd p < hd q then revMerge ((hd p) :: acc, tl p, q) else if hd p > hd q then revMerge ((hd q) :: acc, p, tl q) else revMerge ((hd p) :: acc, tl p, tl q) else if null p then List.revAppend (q, acc) else List.revAppend (p, acc) in if (null p) then q else if (null q) then p else rev (revMerge ([], p, q)) end fun mul m x = if x <= (n div m) then SOME (x * m) else NONE fun mapPrefix pred l = let fun mapp ([], acc) = rev acc | mapp (x :: xs, acc) = (case (pred x) of SOME a => mapp (xs, a :: acc) | NONE => rev acc) in mapp (l, []) end fun mergeWith f (m, i) = merge (f m, i) fun generate l = let fun listMul m = mapPrefix (mul m) l in foldl (mergeWith listMul) [] [2, 3, 5] end fun iterate (acc, l) = if (hd l) > (n div 2) then merge (l, acc) else iterate (merge (l, acc), generate l) in if n > 0 then iterate ([], [1]) else [] end
产生指定范围内的汉明数需要多轮运算,后面每轮中的三个列表元素乘积运算中都可能产生超出这个范围的结果,它们不需要出现在后续的运算之中。[h] 基本库mapPartial
函数与它所映射的函数,通过基于Option
结构的SOME
和NONE
构造子的协定,可以将所映射函数认为不符合条件的元素或者结果排除掉,它会遍历整个列表。[i] 由于这个算法采用升序列表,故而这里将它改写为mapPrefix
函数,用来在特定条件不满足条件就立即结束。下面是汉明数程序在解释器下命令行运行实例:
- fun printIntList (l: int list) = = print ((String.concatWith " " (map Int.toString l)) ^ "\n"); val printIntList = fn : int list -> unit - val _ = printIntList (Hamming_number 400); 1 2 3 4 5 6 8 9 10 12 15 16 18 20 24 25 27 30 32 36 40 45 48 50 54 60 64 72 75 80 81 90 96 100 108 120 125 128 135 144 150 160 162 180 192 200 216 225 240 243 250 256 270 288 300 320 324 360 375 384 400
续体传递风格实例
[编辑]下面是续体传递风格(CPS)[44]的高阶函数foldr和map的实现,和达成一个整数列表的合计函数的简单用例:
fun foldr' f b l = let fun f2 ([], k) = k b | f2 (a :: r, k) = f2 (r, fn x => k (f (a, x))) in f2 (l, fn x => x) end fun map' f l = foldr' (fn (x, y) => (f x) :: y) [] l fun sum l = foldr' (fn (x, y) => (x + y)) 0 l
对于输入[e1, e2, ..., en]
,sum
函数等价于函数复合(fn x => x) o (fn x => e1 + x) o (fn x => e2 + x) o ... o (fn x => en + x)
,它应用于0
上得到(e1 + (e2 + (... + (en + 0)...)))
。[45]
SML/NJ支持头等对象的续体[46]。头等续体对一门语言而言是能完全控制指令执行次序的能力。它们可以用来跳转到产生对当前函数调用的那个函数,或者跳转到此前已经退出了的函数。头等续体保存了程序执行状态,它不保存程序数据,只保存执行上下文。
排序算法
[编辑]排序算法关注计算复杂度,特别是时间复杂度,基本库函数的实现细节也要考虑在内,比如串接函数@
,它被实现为fun l1 @ l2 = revAppend(rev l1, l2)
,除非必需的情况避免使用遍历整个列表的rev
和length
函数。[j] 通过比较于这些排序算法的周知过程式编程语言比如C语言的实现,可以体察到ML在控制流程和列表数据结构上的相关限制,和与之相适应的采用尾递归的特色函数式编程风格。
插入排序
[编辑]尾递归 | 普通递归 |
---|---|
fun insertSort l = let fun insert pred (ins, l) = let fun loop (acc, []) = List.revAppend (acc, [ins]) | loop (acc, l as x :: xs) = if pred (ins, x) then List.revAppend (acc, ins :: l) else loop (x :: acc, xs) in loop ([], l) end val rec ge = fn (x, y) => (x >= y) in rev (foldl (insert ge) [] l) end | fun insertSort l = let fun insert pred (ins, []) = [ins] | insert pred (ins, l as x :: xs) = if pred (ins, x) then ins :: l else x :: (insert pred (ins, xs)) val rec ge = fn (x, y) => (x >= y) in rev (foldl (insert ge) [] l) end |
x 保存在形式参数acc 对应的分配堆中 | x 保存在调用栈栈帧中 |
插入排序算法是稳定的自适应排序,它在输入列表趋于正序的时候性能最佳,在输入列表趋于反序的时候性能最差,因此在算法实现中,需要insert
函数所插入列表保持为反序,在插入都完成后经rev
函数再反转回正序。在预期输入数据趋于随机化或者预知它经过了反转的情况下,可以采用保持要插入列表为正序的变体插入排序算法实现,它在输入列表趋于反序的时候性能最佳,在输入列表趋于正序的时候性能最差,它比自适应实现少作一次全列表反转。[k] 采用foldr
函数可以相应的保持要插入列表为正序,由于fun foldr f b l = foldl f b (rev l)
,它等同于对反转过的列表应用变体插入排序。
希尔排序
[编辑]希尔排序算法是对插入排序的改进,保持了自适应性,放弃了稳定性。[l] 下面是希尔排序的实现:
fun shellSort l = let fun insert pred (ins, l) = let fun loop (acc, []) = List.revAppend (acc, [ins]) | loop (acc, l as x :: xs) = if pred (ins, x) then List.revAppend (acc, ins :: l) else loop (x :: acc, xs) in loop ([], l) end val rec lt = fn (x, y) => (x < y) fun insertSort [] = [] | insertSort [x] = [x] | insertSort [x, y] = if (y < x) then [y, x] else [x, y] | insertSort (x :: y :: z :: xs) = let val (x, y, z) = if (y < x) then if z < y then (z, y, x) else if z < x then (y, z, x) else (y, x, z) else if z < x then (z, x, y) else if z < y then (x, z, y) else (x, y, z) in foldl (insert lt) [x, y, z] xs end fun group (lol, n) = let fun dup n = let fun loop (acc, i) = if i <= 0 then acc else loop ([] :: acc, i - 1) in loop ([], n) end fun loop ([], [], accj, lol') = List.revAppend (accj, lol') | loop (acci, [], accj, []) = loop ([], rev acci, [], rev accj) | loop (acci, [], accj, lol') = loop ([], rev acci, accj, lol') | loop (acci, lol, accj, []) = loop (acci, lol, [], rev accj) | loop (acci, [] :: ls, accj, lol') = loop (acci, ls, accj, lol') | loop (acci, (x :: xs) :: ls, accj, l' :: ls') = loop (xs :: acci, ls, (x :: l') :: accj, ls') in loop ([], lol, [], dup n) end val (lol, len) = foldl (fn (x, (l, n)) => ([x] :: l, n + 1)) ([], 0) (rev l) val incs = [1, 4, 9, 20, 46, 103, 233, 525, 1182, 2660, 5985, 13467, 30301, 68178, 153401, 345152, 776591, 1747331, 3931496, 8845866, 19903198, 44782196, 100759940, 226709866, 510097200] val gap = let val v = len * 3 div 4 val thold = if (v = 0) then 1 else v fun loop (acc, h) = if (hd h) > thold then acc else loop ((hd h) :: acc, tl h) in loop ([], incs) end fun sort (h, lol) = map insertSort (group (lol, h)) in hd (foldl sort lol gap) end
这里采用的间隔序列是OEIS A108870,即,它是徳田尚之在1992年提出的[47]。这个序列用递推公式表示为:hk = ⌈h'k⌉
,这里的h'k = 2.25·h'k-1 + 1
而h'1 = 1
。假定一个列表的长度s
位于序列两个元素之间,即hk-1 < hk ≤ s < hk+1
,如果hk ≤ ·s
,这里的n ≤ m
,则选择初始间隔为hk
,否则为hk-1
。在这个阈值下,对于不同长度s
的列表和对应的初始间隔h
,每个列表的这些初始子列表的平均长度,约在
≤ < ·
范围之内。间隔序列还可以采用OEIS A102549,它是Marcin Ciura在2001年通过实验得到的[48]。[m]
快速排序
[编辑]下面是快速排序算法的自顶向下实现:
fun quickSort [] = [] | quickSort [x] = [x] | quickSort [x, y] = if x <= y then [x, y] else [y, x] | quickSort [x, y, z] = let val (x, y) = if x <= y then (x, y) else (y, x) val (y, z) = if y <= z then (y, z) else (z, y) val (x, y) = if x <= y then (x, y) else (y, x) in [x, y, z] end | quickSort (pivot :: xs) = let fun partition pred l = let fun loop ([], p, q) = (p, q) | loop (h :: t, p, q) = if (pred h) then loop (t, h :: p, q) else loop (t, p, h :: q) in loop (l, [], []) end val (le, gt) = partition (fn x => (x <= pivot)) xs in (quickSort le) @ (pivot :: (quickSort gt)) end
基本库partition
函数实现对快速排序而言有不必要的反转,[n] 这里采用了它的简化改写。在ML中快速排序应采用自底向上实现:
fun quickSort l = let fun partition pred l = let fun loop ([], p, q) = (p, q) | loop (h :: t, p, q) = if (pred h) then loop (t, h :: p, q) else loop (t, p, h :: q) in loop (l, [], []) end fun iterate (acc, []) = acc | iterate (acc, [] :: xs) = iterate (acc, xs) | iterate (acc, [x] :: xs) = iterate (x :: acc, xs) | iterate (acc, [x, y] :: xs) = let val (x, y) = if x <= y then (x, y) else (y, x) in iterate (x :: y :: acc, xs) end | iterate (acc, [x, y, z] :: xs) = let val (x, y) = if x <= y then (x, y) else (y, x) val (x, y, z) = if y <= z then (x, y, z) else if x <= z then (x, z, y) else (z, x, y) in iterate (x :: y :: z :: acc, xs) end | iterate (acc, (pivot :: d) :: xs) = let val (le, gt) = partition (fn x => (x <= pivot)) d in iterate (acc, gt :: [pivot] :: le :: xs) end in iterate ([], [l]) end
归并排序
[编辑]下面是归并排序算法的自底向上法实现:
fun mergeSort l = let fun init (acc, []) = acc | init (acc, [x]) = [x] :: acc | init (acc, [x, y]) = if x <= y then [x, y] :: acc else [y, x] :: acc | init (acc, x :: y :: z :: xs) = let val (x, y, z) = if x <= y then if y <= z then (x, y, z) else if x <= z then (x, z, y) else (z, x, y) else if x <= z then (y, x, z) else if y <= z then (y, z, x) else (z, y, x) in init ([x, y, z] :: acc, xs) end fun mergeWith _ (acc, [], []) = acc | mergeWith _ (acc, p, []) = List.revAppend (p, acc) | mergeWith _ (acc, [], q) = List.revAppend (q, acc) | mergeWith cmp (acc, p :: ps, q :: qs) = if cmp (p, q) then mergeWith cmp (p :: acc, ps, q :: qs) else mergeWith cmp (q :: acc, p :: ps, qs) val mergeRev = mergeWith (fn (x, y) => (x > y)) val revMerge = mergeWith (fn (x, y) => (x < y)) fun iterate ([], _) = [] | iterate ([x], isRev) = if isRev then rev x else x | iterate (acc, isRev) = let val merge = if isRev then mergeRev else revMerge fun loop (acci, []) = acci | loop (acci, [x]) = (rev x) :: acci | loop (acci, x :: y :: xs) = loop (merge ([], x, y) :: acci, xs) in iterate (loop ([], acc), not isRev) end in iterate (init ([], l), false) end
考虑输入列表[x1, ..., xi, a0, ..., a9, xj, ..., xn]
,这里在xi
和xj
之间的10
个a
,具有相同的值并且需要保持其下标表示的次序,这里的xi > a
并且xj < a
,并且在这个区段前后的元素总数都能被3
整除,则它被分解成子列表的列表[Xm, ..., [xj, a8, a9], [a5, a6, a7], [a2, a3, a4], [a0, a1, xi], ..., X1]
,这里有m = n div 3
;假定这4
个含有a
的子列表两两归并,在归并正序子列表的归并条件x < y
下,能得到[X1, ..., [xi, a4, ..., a0], [a9, ..., a5, xj], ..., Xk]
;继续推演下去,在归并反序子列表的归并条件x > y
下,能得到[Xh, ..., [xj, a0, ..., a9, xi], ..., X1]
。可以看出这种归并操作能够保证排序算法的稳定性,即具有相同值的元素之间的相对次序保持不变。
分解初始的子列表采用了插入排序,还可进一步增加其大小。归并排序也有自顶向下实现。[o]
堆排序
[编辑]fun heapSort l = let val h = Array.fromList l val len = Array.length h fun get i = Array.sub (h, i) fun set i v = Array.update (h, i, v) fun siftdown (i, ins, n) = let fun sift k = let val l = k * 2 + 1 val r = l + 1 in if (r < n) andalso (get r) > (get l) then r else if (l < n) then l else k end fun loop i = let val j = sift i in if j = i orelse (get j) <= ins then set i ins else (set i (get j); loop j) end in loop i end fun heapify () = let fun loop i = if i < 0 then () else (siftdown (i, get i, len); loop (i - 1)) in loop (len div 2 - 1) end fun generate () = let fun loop (acc, i) = let val t = get 0 in if i <= 0 then t :: acc else (siftdown (0, get i, i); loop (t :: acc, i - 1)) end in if len = 0 then [] else loop ([], len - 1) end in heapify (); generate () end
在数组实现中,siftdown
函数融合了插入和筛选功能,它首先在暂时位于堆顶的要插入的元素,和从堆顶节点左右子堆的两个堆顶元素中筛选出的那个元素,二者中选择出哪个适合作堆顶元素;如果要插入元素适合,则以它为新堆顶元素而结束整个过程,否则以筛选出元素为新堆顶元素,并自顶向下逐级处理提供了新堆顶元素的子堆,将要插入元素暂时作为其堆顶元素并对它继续进行siftdown
;siftdown
只要到达了某个堆底,就插入要插入的元素而结束整个过程。
在提取堆顶元素生成结果列表时,先提取走堆顶元素的内容,再摘掉最后的堆底元素将它的内容暂时放置在堆顶,这时堆的大小也相应的减一,随后的siftdown
函数,筛选出新的堆顶元素,并把原最后堆底元素插入回堆中。
在heapify
函数建造堆的时候,首先自列表中间将元素分为前后两组,后组中的元素被视为只有一个元素的堆,然后从后往前处理前组中的元素,这时它的左右子节点已经是已有的堆或者为空,在其上进行siftdown
,从而合成一个新堆。建造堆也可以采用siftup
函数来实现,它自第二个元素开始从前往后逐个处理列表元素,其前面是已有的堆,将这个新元素自堆底向上插入到这个堆中。[p]
fun heapSort l = let datatype 'a heap = Nil | Leaf of 'a | Node of 'a * int * 'a heap * 'a heap fun key Nil = let val SOME a = Int.minInt in a end | key (Leaf k) = k | key (Node (k, _, _, _)) = k fun count Nil = 0 | count (Leaf _) = 1 | count (Node (_, c, _, _)) = c fun left Nil = Nil | left (Leaf _) = Nil | left (Node (_, _, l, _)) = l fun insert (Nil, x) = Leaf x | insert (Leaf k, l) = if l >= k then Node (l, 2, Leaf k, Nil) else Node (k, 2, Leaf l, Nil) | insert (Node (k, _, Leaf l, Nil), r) = if r >= k then Node (r, 3, Leaf k, Leaf l) else if r >= l then Node (k, 3, Leaf r, Leaf l) else Node (k, 3, Leaf l, Leaf r) | insert (Node (k, c, l, r), x) = let val (k, x) = if k >= x then (k, x) else (x, k) in if (count l) <= (count r) then Node (k, c + 1, insert (l, x), r) else if x >= (key l) then Node (k, c + 1, insert (r, x), l) else Node (k, c + 1, l, insert (r, x)) end fun extract Nil = Nil | extract (Leaf _) = Nil | extract (Node (_, _, l, Nil)) = l | extract (Node (_, c, l, r)) = let val k = key l val n = left l in if n = Nil then Node (k, c - 1, r, Nil) else if (key n) >= (key r) then Node (k, c - 1, extract l, r) else Node (k, c - 1, r, extract l) end fun heapify () = let fun loop (acc, []) = acc | loop (acc, x :: xs) = loop (insert (acc, x), xs) in loop (Nil, l) end fun generate h = let fun loop (acc, Nil) = acc | loop (acc, h) = loop ((key h) :: acc, extract h) in loop ([], h) end in generate (heapify ()) end
二叉树实现不能直接访问堆底元素,从而不适宜通过摘掉它使堆的大小减一。这里的insert
函数,在原堆顶元素和要插入元素中选择适合者作为新的堆顶元素,将落选的另一个元素作为新的要插入元素,插入到利于保持这个堆平衡的那个子树之中。这里的extract
函数只筛选不插入,它将堆的大小减一。
这里的insert
和extract
函数也可以直接转写为等价的尾递归形式,与列表情况不同,只要树结构能保持良好的平衡,采用尾递归形式就没有太大的必要性。[q] 在二叉树实现下,也可以采用siftdown
函数来初始建造堆,而不需要在节点中保存关于树状态的统计信息。[r]
基数排序
[编辑]fun radixSort l = let fun distribute (l, d) = let val t0 = ([], [], [], [], [], [], [], []) fun loop (t, []) = let fun join (acc, i) = let val f = case i of 1 => (#1 t) | 2 => (#2 t) | 3 => (#3 t) | 4 => (#4 t) | 5 => (#5 t) | 6 => (#6 t) | 7 => (#7 t) | 8 => (#8 t) | _ => [] in if i <= 0 then acc else join (List.revAppend (f, acc), i - 1) end in join ([], 8) end | loop (t, x :: xs) = let val (f0, f1, f2, f3, f4, f5, f6, f7) = t val t = case ((x div d) mod 8) of 0 => (x :: f0, f1, f2, f3, f4, f5, f6, f7) | 1 => (f0, x :: f1, f2, f3, f4, f5, f6, f7) | 2 => (f0, f1, x :: f2, f3, f4, f5, f6, f7) | 3 => (f0, f1, f2, x :: f3, f4, f5, f6, f7) | 4 => (f0, f1, f2, f3, x :: f4, f5, f6, f7) | 5 => (f0, f1, f2, f3, f4, x :: f5, f6, f7) | 6 => (f0, f1, f2, f3, f4, f5, x :: f6, f7) | 7 => (f0, f1, f2, f3, f4, f5, f6, x :: f7) | _ => t0 in loop (t, xs) end in loop (t0, l) end val SOME maxInt = Int.maxInt val max = foldl (fn (x, y) => if x > y then x else y) 0 l fun iterate (l, d) = let val l' = distribute (l, d) in if d >= (maxInt div 8 + 1) orelse ((max div d) div 8) = 0 then l' else iterate (l', d * 8) end in iterate (l, 1) end
这里采用的基数是2
的3
次幂8
,代码所使用的列表元组大小与基数大小成正比,运算量与列表中元素的总数与最大数的位数的乘积成正比。
随机数生成
[编辑]编写排序算法进行测试除了使用简单的数列,[s] 测试用列表还可以使用线性同余伪随机数函数来生成[49]:
fun randList (n, seed) = let val randx = ref seed fun lcg () = (randx := (!randx * 421 + 1663) mod 7875; !randx) (* fun lcg () = (randx := (!randx * 1366 + 150889) mod 714025; !randx) *) fun iterate (acc, i) = if i <= 0 then acc else iterate (lcg () :: acc, i - 1) in iterate ([], n) end
语言解释器
[编辑]定义和处理一个小型表达式语言是相对容易的:
exception Err; datatype ty = IntTy | BoolTy datatype exp = True | False | Int of int | Not of exp | Add of exp * exp | If of exp * exp * exp fun typeOf (True) = BoolTy | typeOf (False) = BoolTy | typeOf (Int _) = IntTy | typeOf (Not e) = if typeOf e = BoolTy then BoolTy else raise Err | typeOf (Add (e1, e2)) = if (typeOf e1 = IntTy) andalso (typeOf e2 = IntTy) then IntTy else raise Err | typeOf (If (e1, e2, e3)) = if typeOf e1 <> BoolTy then raise Err else if typeOf e2 <> typeOf e3 then raise Err else typeOf e2 fun eval (True) = True | eval (False) = False | eval (Int n) = Int n | eval (Not e) = (case eval e of True => False | False => True | _ => raise Fail "type-checking is broken") | eval (Add (e1, e2)) = let val (Int n1) = eval e1 val (Int n2) = eval e2 in Int (n1 + n2) end | eval (If (e1, e2, e3)) = if eval e1 = True then eval e2 else eval e3 fun exp_repr e = let val msg = case e of True => "True" | False => "False" | Int n => Int.toString n | _ => "" in msg ^ "\n" end (* 忽略TypeOf的返回值,它在类型错误时发起Err *) fun evalPrint e = (ignore (typeOf e); print (exp_repr (eval e)));
将这段代码录入文件比如expr-lang.sml
,并在命令行下执行sml expr-lang.sml
,可以用如下在正确类型的和不正确类型上运行的例子,测试这个新语言:
- val e1 = Add (Int 1, Int 2); (* 正确的类型 *) val e1 = Add (Int 1,Int 2) : exp - val _ = evalPrint e1; 3 - val e2 = Add (Int 1, True); (* 不正确的类型 *) val e2 = Add (Int 1,True) : exp - val _ = evalPrint e2; uncaught exception Err raised at: expr-lang.sml:25.20-25.23
注释和附录代码
[编辑]- ^ 子句集合是穷尽和不冗余的函数示例:
fun hasCorners (Circle _) = false | hasCorners _ = true
如果控制通过了第一个模式
Circle
,则这个值必定是要么Square
要么Triangle
。在任何这些情况下,这个形状都有角点,所以返回true
而不用区分具体是那种情况。 - ^ 不详尽的模式示例:
fun center (Circle (c, _)) = c | center (Square ((x, y), s)) = (x + s / 2.0, y + s / 2.0)
这里在
center
函数中没有给Triangle
的模式。 编译器发起模式不详尽的一个警告,并且如果在运行时间,Triangle
被传递给这个函数,则发起例外Match
。 - ^ 存在模式冗余的(无意义)函数示例:
fun f (Circle ((x, y), r)) = x+y | f (Circle _) = 1.0 | f _ = 0.0
匹配第二个子句的模式的任何值都也匹配第一个子句的模式,所以第二个子句是不可到达的。因此这个定义整体上是冗余的。
- ^ 映射函数的实际实现:
fun map f = let fun m [] = [] | m [a] = [f a] | m [a, b] = [f a, f b] | m [a, b, c] = [f a, f b, f c] | m (a :: b :: c :: d :: r) = f a :: f b :: f c :: f d :: m r in m end
折叠函数的实际实现:
fun foldl f b l = let fun f2 ([], b) = b | f2 (a :: r, b) = f2 (r, f (a, b)) in f2 (l, b) end fun foldr f b l = foldl f b (rev l)
过滤器函数的实际实现:
fun filter pred [] = [] | filter pred (a :: rest) = if pred a then a :: (filter pred rest) else (filter pred rest)
- ^ 对分数采取修约的有理数实现:
signature ARITH = sig eqtype t val zero : t val one : t val fromInt: int -> t val fromIntPair : int * int -> t val repr : t -> unit val succ : t -> t val pred : t -> t val ~ : t -> t val + : t * t -> t val - : t * t -> t val * : t * t -> t val / : t * t -> t end
structure Rational : ARITH = struct type t = int * int val zero = (0, 1) val one = (1, 1) val maxInt = (let val SOME a = Int.maxInt in Int.toLarge a end) fun fromInt n = (n, 1) fun neg (a, b) = (~a, b) fun fromLargePair (a, b) = (Int.fromLarge a, Int.fromLarge b) fun fromIntPair (num, den) = let fun reduced_fraction (numerator, denominator) = let fun gcd (n, 0) = n | gcd (n, d) = gcd (d, n mod d) val d = gcd (numerator, denominator) in if d > 1 then (numerator div d, denominator div d) else (numerator, denominator) end in if num < 0 andalso den < 0 then reduced_fraction (~num, ~den) else if num < 0 then neg (reduced_fraction (~num, den)) else if den < 0 then neg (reduced_fraction (num, ~den)) else reduced_fraction (num, den) end fun repr (a, b) = let val m = if (b = 0) then 0 else if (a >= 0) then a div b else ~a div b val n = if (b = 0) then 1 else if (a >= 0) then a mod b else ~a mod b val ms = Int.toString m and ns = Int.toString n and bs = Int.toString b in if n <> 0 andalso m <> 0 andalso a < 0 then print ("~" ^ ms ^ " - " ^ ns ^ "/" ^ bs ^ "\n") else if n <> 0 andalso m <> 0 then print (ms ^ " + " ^ ns ^ "/" ^ bs ^ "\n") else if n <> 0 andalso m = 0 andalso a < 0 then print ("~" ^ ns ^ "/" ^ bs ^ "\n") else if n <> 0 andalso m = 0 then print (ns ^ "/" ^ bs ^ "\n") else if a < 0 then print ("~" ^ ms ^ "\n") else print (ms ^ "\n") end fun convergent (i, n, d, h_1, k_1, h_2, k_2) = if i <> 0 andalso ((h_1 > (maxInt - h_2) div i) orelse (k_1 > (maxInt - k_2) div i)) then (h_1, k_1) else if n = 0 then (i * h_1 + h_2, i * k_1 + k_2) else convergent (d div n, d mod n, n, i * h_1 + h_2, i * k_1 + k_2, h_1, k_1) fun add ((a, b), (c, d)) = let val la = Int.toLarge a and lb = Int.toLarge b val lc = Int.toLarge c and ld = Int.toLarge d val m = la * ld + lb * lc and n = lb * ld in if b = d then fromIntPair (a + c, b) else fromLargePair (convergent (m div n, m mod n, n, 1, 0, 0, 1)) end fun sub ((a, b), (c, d)) = let val la = Int.toLarge a and lb = Int.toLarge b val lc = Int.toLarge c and ld = Int.toLarge d val m = la * ld - lb * lc and n = lb * ld in if b = d then fromIntPair (a - c, b) else if m < 0 then neg (fromLargePair (convergent (~m div n, ~m mod n, n, 1, 0, 0, 1))) else if m > 0 then fromLargePair (convergent (m div n, m mod n, n, 1, 0, 0, 1)) else (0, 1) end fun mul ((a, b), (c, d)) = let val la = Int.toLarge a and lb = Int.toLarge b val lc = Int.toLarge c and ld = Int.toLarge d val m = la * lc and n = lb * ld in fromLargePair (convergent (m div n, m mod n, n, 1, 0, 0, 1)) end fun op + ((a, b), (c, d)) = if a < 0 andalso c < 0 then neg (add ((~a, b), (~c, d))) else if a < 0 then sub ((c, d), (~a, b)) else if c < 0 then sub ((a, b), (~c, d)) else add ((a, b), (c, d)) fun op - ((a, b), (c, d)) = if a < 0 andalso c < 0 then sub ((~c, d), (~a, b)) else if a < 0 then neg (add ((~a, b), (c, d))) else if c < 0 then add ((a, b), (~c, d)) else sub ((a, b), (c, d)) fun op * ((a, b), (c, d)) = if a < 0 andalso c < 0 then mul ((~a, b), (~c, d)) else if a < 0 then neg (mul ((~a, b), (c, d))) else if c < 0 then neg (mul ((a, b), (~c, d))) else mul ((a, b), (c, d)) fun op / ((a, b), (c, d)) = if a < 0 andalso c < 0 then mul ((~a, b), (d, ~c)) else if a < 0 then neg (mul ((~a, b), (d, c))) else if c < 0 then neg (mul ((a, b), (d, ~c))) else mul ((a, b), (d, c)) fun succ a = add (a, one) fun pred a = sub (a, one) fun ~ a = neg a end
- ^ 找出函数的实际实现:
fun find pred [] = NONE | find pred (a :: rest) = if pred a then SOME a else (find pred rest)
存在谓词函数的实际实现:
fun exists pred = let fun f [] = false | f (h :: t) = pred h orelse f t in f end
- ^ 筛法基于数组的实现:
fun prime n = let val sieve = Array.array (n + 1, true); fun markComposite (j, i) = if j > n then () else (Array.update (sieve, j, false); markComposite (j + i, i)) fun iterate i = if i * i > n then () else if Array.sub (sieve, i) then (markComposite (i * i, i); iterate (i + 1)) else iterate (i + 1) fun generate (acc, i) = if i < 2 then acc else if Array.sub (sieve, i) then generate (i :: acc, i - 1) else generate (acc, i - 1) in if n < 2 then [] else (iterate 2; generate ([], n)) end
- ^ 汉明数进一步性质演示代码:
fun printIntList (l: int list) = print ((String.concatWith " " (map Int.toString l)) ^ "\n"); fun diff (p, q) = let fun revDiff (acc, p, q) = if not (null p) andalso not (null q) then if hd p < hd q then revDiff ((hd p) :: acc, tl p, q) else if hd p > hd q then revDiff (acc, p, tl q) else revDiff (acc, tl p, tl q) else if null q then List.revAppend (p, acc) else acc in if (null p) then [] else if (null q) then p else rev (revDiff ([], p, q)) end; fun Hamming_number n = let fun merge (p, q) = let fun revMerge (acc, p, q) = if not (null p) andalso not (null q) then if hd p < hd q then revMerge ((hd p) :: acc, tl p, q) else if hd p > hd q then revMerge ((hd q) :: acc, p, tl q) else revMerge ((hd p) :: acc, tl p, tl q) else if null p then List.revAppend (q, acc) else List.revAppend (p, acc) in if (null p) then q else if (null q) then p else rev (revMerge ([], p, q)) end fun mergeWith f (m, i) = merge (f m, i) fun mul m x = (x * m) fun generate l = let fun listMul m = map (mul m) l in printIntList (listMul 2); printIntList (diff (listMul 3, listMul 2)); printIntList (diff (listMul 5, merge (listMul 2, listMul 3))); foldl (mergeWith listMul) [] [2, 3, 5] end val count = ref 1; fun iterate (acc, l) = if (hd l) > (n div 2) then merge (l, acc) else (print ("round " ^ (Int.toString (!count)) ^ " for " ^ (Int.toString(length l)) ^ " number(s)\n"); count := !count + 1; iterate (merge (l, acc), generate l)) in if n > 0 then iterate ([], [1]) else [] end;
val l = Hamming_number 400; val h = List.filter (fn x => (x <= 400)) l; val _ = print ((Int.toString (length h)) ^ " numbers from " ^ (Int.toString (length l)) ^ " numbers\n");
- ^ 部份映射函数的实际实现:
fun mapPartial pred l = let fun mapp ([], l) = rev l | mapp (x :: r, l) = (case (pred x) of SOME y => mapp (r, y :: l) | NONE => mapp (r, l) (* end case *)) in mapp (l, []) end
- ^ 列表长度函数的实际实现:
fun length l = let (* fast add that avoids the overflow test *) val op + = Int.fast_add fun loop (n, []) = n | loop (n, [_]) = n + 1 | loop (n, _ :: _ :: l) = loop (n + 2, l) in loop (0, l) end
- ^ 变体的插入排序算法实现:
fun insertSort l = let fun insert pred (ins, l) = let fun loop (acc, []) = List.revAppend (acc, [ins]) | loop (acc, l as x :: xs) = if pred (ins, x) then List.revAppend (acc, ins :: l) else loop (x :: acc, xs) in loop ([], l) end val rec lt = fn (x, y) => (x < y) in foldl (insert lt) [] l end
- ^ 下面是希尔排序算法的原型实现,当输入列表趋于正序的时候,
scatter
函数将其分解为一组趋于反序的子列表,经过在输入趋于反序情况下性能最佳的变体插入排序后,gather
函数再将它们合成为一个列表:fun shellSort l = let fun insert pred (ins, l) = let fun loop (acc, []) = List.revAppend (acc, [ins]) | loop (acc, l as x :: xs) = if pred (ins, x) then List.revAppend (acc, ins :: l) else loop (x :: acc, xs) in loop ([], l) end val rec lt = fn (x, y) => (x < y) fun insertSort l = foldl (insert lt) [] l fun scatter (l, n) = let fun dup n = let fun loop (acc, i) = if i <= 0 then acc else loop ([] :: acc, i - 1) in loop ([], n) end fun loop ([], acc, lol) = List.revAppend (acc, lol) | loop (l, acc, []) = loop (l, [], rev acc) | loop (x :: xs, acc, l :: ls) = loop (xs, (x :: l) :: acc, ls) in loop (l, [], dup n) end fun gather lol = let fun loop ([], [], l) = rev l | loop (acc, [], l) = loop ([], rev acc, l) | loop (acc, [] :: ls, l) = loop (acc, ls, l) | loop (acc, (x :: xs) :: ls, l) = loop (xs :: acc, ls, x :: l) in loop ([], lol, []) end val gap = let fun loop (acc, i) = let val h = (i * 5 - 1) div 11 in if i < 5 then rev (1 :: acc) else loop (h :: acc, h) end in loop ([], length l) end fun sort (h, l) = gather (map insertSort (scatter (l, h))) in foldl sort l gap end
- ^ 希尔排序还可以采用Ciura提出的间隔序列:
val incs = [1750, 701, 301, 132, 57, 23, 10, 4, 1] val gap = let fun loop (acc, i) = let val h = (i * 4 - 1) div 9 fun iter incs = if i >= (hd incs) * 4 div 3 then incs else iter (tl incs) in if i <= ((hd incs) * 3) then List.revAppend (acc, iter incs) else loop (h :: acc, h) end in if len = 0 then [1] else loop ([], len) end
- ^ 划分函数的实际实现:
fun partition pred l = let fun loop ([], trueList, falseList) = (rev trueList, rev falseList) | loop (h :: t, trueList, falseList) = if pred h then loop (t, h :: trueList, falseList) else loop (t, trueList, h :: falseList) in loop (l, [], []) end
- ^ 归并排序算法的自顶向下实现:
fun mergeSort l = let fun sort ([], _) = [] | sort ([x], _) = [x] | sort ([x, y], _) = if x <= y then [x, y] else [y, x] | sort ([x, y, z], _) = if x <= y then if y <= z then [x, y, z] else if x <= z then [x, z, y] else [z, x, y] else if x <= z then [y, x, z] else if y <= z then [y, z, x] else [z, y, x] | sort (l, n) = let val m = n div 2 fun split (l, acc, i) = if i = 0 then (rev acc, l) else split (tl l, (hd l) :: acc, i - 1) fun merge (p, q) = let fun revMerge (acc, p, q) = if not (null p) andalso not (null q) then if (hd p) <= (hd q) then revMerge ((hd p) :: acc, tl p, q) else revMerge ((hd q) :: acc, p, tl q) else if null p then List.revAppend (q, acc) else List.revAppend (p, acc) in rev (revMerge ([], p, q)) end val (ls, rs) = split (l, [], m) in merge (sort (ls, m), sort (rs, n - m)) end in sort (l, length l) end
- ^ 堆排序算法的数组实现中,堆建造函数
heapify
也可以使用siftup
函数的数组来完成,它将新节点自堆底向上插入到合适的位置,而不对途径节点左右子树顶点进行比较:fun siftup i = let val ins = get i fun loop i = let val j = (i - 1) div 2 in if i <= 0 orelse (get j) >= ins then set i ins else (set i (get j); loop j) end in loop i end fun heapify () = let fun loop i = if i > (len - 1) then () else (siftup i; loop (i + 1)) in loop 1 end
- ^ 在堆排序算法的二叉树实现中,树插入和提取函数也可以写为等价的尾递归形式代码:
exception Err; fun revLink ([], t) = t | revLink (Nil :: xs, t) = revLink (xs, t) | revLink ((Leaf k) :: xs, t) = revLink (xs, Node (k, (count t) + 1, t, Nil)) | revLink (Node (k, c, Nil, r) :: xs, t) = revLink (xs, Node (k, c, t, r)) | revLink (Node (k, c, l, Nil) :: xs, t) = revLink (xs, Node (k, c, l, t)) | revLink (_ :: xs, t) = raise Err fun insert (t, x) = let fun loop (acc, Nil, x) = revLink (acc, Leaf x) | loop (acc, Leaf k, l) = if l >= k then revLink (acc, Node (l, 2, Leaf k, Nil)) else revLink (acc, Node (k, 2, Leaf l, Nil)) | loop (acc, Node (k, _, Leaf l, Nil), r) = if r >= k then revLink (acc, Node (r, 3, Leaf k, Leaf l)) else if r >= l then revLink (acc, Node (k, 3, Leaf r, Leaf l)) else revLink (acc, Node (k, 3, Leaf l, Leaf r)) | loop (acc, Node (k, c, l, r), x) = let val (k, x) = if k >= x then (k, x) else (x, k) in if (count l) <= (count r) then loop (Node (k, c + 1, Nil, r) :: acc, l, x) else if x >= (key l) then loop (Node (k, c + 1, Nil, l) :: acc, r, x) else loop (Node (k, c + 1, l, Nil) :: acc, r, x) end in loop ([], t, x) end fun extract t = let fun loop (acc, Nil) = revLink (acc, Nil) | loop (acc, Leaf _) = revLink (acc, Nil) | loop (acc, Node (_, _, l, Nil)) = revLink (acc, l) | loop (acc, Node (_, c, l, r)) = let val k = key l val n = left l in if n = Nil then revLink (acc, Node (k, c - 1, r, Nil)) else if (key n) >= (key r) then loop (Node (k, c - 1, Nil, r) :: acc, l) else loop (Node (k, c - 1, r, Nil) :: acc, l) end in loop ([], t) end
- ^ 堆排序算法的二叉树实现中,堆建造函数也可以主要采用
siftdown
函数来完成,这里不需要在节点中保存关于树状态的统计信息:fun heapSort l = let datatype 'a heap = Nil | Node of 'a * 'a heap * 'a heap fun key Nil = let val SOME a = Int.minInt in a end | key (Node (k, _, _)) = k fun left Nil = Nil | left (Node (_, l, _)) = l fun right Nil = Nil | right (Node (_, _, r)) = r fun leaf k = Node (k, Nil, Nil) fun sift (l, r) = if l <> Nil andalso r <> Nil then if (key l) >= (key r) then l else r else if l <> Nil then l else if r <> Nil then r else Nil fun siftdown (x, l, r) = let val superior = sift (l, r) in if superior = Nil then Node (x, Nil, Nil) else if x >= (key superior) then Node (x, l, r) else if superior = l then Node (key l, siftdown (x, left l, right l), r) else Node (key r, l, siftdown (x, left r, right r)) end fun insert (Nil, x) = Node (x, Nil, Nil) | insert (Node (k, l, r), x) = let val superior = sift (l, r) in if x >= k andalso superior = l then Node (x, l, insert (r, k)) else if x >= k then Node (x, insert (l, k), r) else if superior = l then Node (k, l, insert (r, x)) else Node (k, insert (l, x), r) end fun extract Nil = Nil | extract (Node (_, l, r)) = let val superior = sift (l, r) in if superior = Nil then Nil else if superior = l then Node (key l, extract l, r) else Node (key r, l, extract r) end fun join (l, r) = extract (Node (key Nil, l, r)) fun heapify () = let fun init (hs, ls, [], _) = (hs, ls) | init (hs, ls, x :: xs, flag) = if flag then init ((leaf x) :: hs, ls, xs, false) else init (hs, x :: ls, xs, true) val (hs, ls) = init ([], [], l, true) fun loop ([], [], []) = Nil | loop ([], [h], []) = h | loop ([], [], x :: xs) = loop ([], [leaf x], xs) | loop ([], [h], x :: xs) = loop ([], [insert (h, x)], xs) | loop (acc, [], l) = loop ([], acc, l) | loop (acc, [h], l) = loop ([], h :: acc, l) | loop (acc, l :: r :: hs, []) = loop (join (l, r) :: acc, hs, []) | loop (acc, l :: r :: hs, x :: xs) = loop (siftdown (x, l, r) :: acc, hs, xs) in loop ([], hs, ls) end fun generate h = let fun loop (acc, Nil) = acc | loop (acc, h) = loop ((key h) :: acc, extract h) in loop ([], h) end in generate (heapify ()) end
稍加处理,堆建造函数也可以只用
siftdown
函数来完成:fun heapify () = let exception Err; fun split () = let val (rs, ts) = let fun loop (acci, accj, [], i, n) = if i = n then (List.revAppend (accj, acci), []) else (acci, accj) | loop (acci, accj, x :: xs, i, n) = if i = n then loop (List.revAppend (accj, acci), [x], xs, i + 1, n * 2 + 1) else loop (acci, x :: accj, xs, i + 1, n) in loop ([], [], l, 0, 1) end fun loop (hs, ls, [], _) = (hs, ls, ts) | loop (hs, ls, x :: xs, flag) = if flag then loop (x :: hs, ls, xs, false) else loop (hs, x :: ls, xs, true) in loop ([], [], rs, true) end fun init () = let val (hs, ls, ts) = split () fun loop (acc, [], []) = (acc, ls) | loop (acc, [], ts) = (acc, List.revAppend (ts, ls)) | loop (acc, k :: hs, []) = loop ((leaf k) :: acc, hs, []) | loop (acc, k :: hs, [x]) = let val (k, x) = if k >= x then (k, x) else (x, k) in loop (Node (k, leaf x, Nil) :: acc, hs, []) end | loop (acc, k :: hs, l :: r :: rs) = let val (k, l, r) = if k >= l then if l >= r then (k, l, r) else if k >= r then (k, r, l) else (r, k, l) else if k >= r then (l, k, r) else if l >= r then (l, r, k) else (r, l, k) in loop (Node (k, leaf l, leaf r) :: acc, hs, rs) end in loop ([], hs, ts) end val (hs, ls) = init () fun loop ([], [], []) = Nil | loop ([], [h], []) = h | loop ([], [], _) = raise Err | loop ([], [h], _) = raise Err | loop (acc, [], l) = loop ([], acc, l) | loop (acc, [h], l) = loop ([], h :: acc, l) | loop (acc, l :: r :: hs, []) = raise Err | loop (acc, l :: r :: hs, x :: xs) = loop (siftdown (x, l, r) :: acc, hs, xs) in loop ([], hs, ls) end
- ^ 排序算法的简单测试代码:
fun printIntList (l: int list) = print ((String.concatWith " " (map Int.toString l)) ^ "\n") fun shuffle (l, n) = let fun split (l, acc, i) = if i = 0 then (acc, l) else split (tl l, (hd l) :: acc, i - 1) fun zip (acc, p, q, flag) = if null p then List.revAppend (q, acc) else if null q then List.revAppend (p, acc) else if flag then zip ((hd p) :: acc, tl p, q, false) else zip ((hd q) :: acc, p, tl q, true) val (p, q) = split (l, [], n div 2) in if (null l) then tl [0] else zip ([], p, q, true) end fun testsort f n = let fun loop (acc, i) = if (i <= 0) then acc else loop (i :: acc, i - 1) val sl = shuffle (loop ([], n), n) val ssl = shuffle (sl, n) in print ("source list is: "); printIntList ssl; print ("result list is: "); printIntList (f ssl) end
参见
[编辑]- Standard ML和它的实现:
- SML/NJ,由普林斯顿大学和贝尔实验室联合开发的实现,它具有并发编程扩展Concurrent ML。
- MLton,严格遵循标准定义的强力的全程序优化编译器[50]。
- HaMLet[51],由马克斯·普朗克软件系统研究所(MPI-SWS)的Andreas Rossberg编写,是一个SML解释器,意图成为精确且合宜接近的标准定义参考实现。
- OCaml,由法国国家信息与自动化研究所(INRIA)维护,是一个“工业强度”的ML方言[52],演化自最初用来实现Coq定理证明器的Caml[17]。
- Alice,由萨尔兰大学设计的Alice ML,是基于Standard ML的函数式编程语言,扩展了对并发、分布式和约束式编程的丰富支持[53]。
- ATS,由波士顿大学的Hongwei Xi和卡内基·梅隆大学的Frank Pfenning提出的Dependent ML发展而来,它向ML扩展了依赖类型。
- F#,由微软研究院(MSR)开发,是一个基于OCaml的一个以.NET为目标的编程语言。
- F*,由MSR和INRIA主导开发,是一个基于ML的依赖类型函数式编程语言。
- Futhark,由哥本哈根大学计算机科学系(DIKU)开发,是属于ML家族的一个函数式、数据并行、阵列编程语言[4]。
- Ur,由麻省理工学院的Adam Chlipala开发,是语法基于Standard ML的专门用于web开发的一个函数式编程语言[7]。
- CM-Lex和CM-Yacc,由卡内基·梅隆大学的Karl Crary开发,是用于Standard ML、OCaml和Haskell的词法分析器和语法解析器[54]。
- Amulet,是一个类ML的函数式编程语言,其编译器输出Lua文件[55]。
- Alpaca,是一个受ML启发的运行在Erlang虚拟机上的函数式编程语言[56]。
延伸阅读
[编辑]- Robin Milner, Mads Tofte, Robert Harper. The Definition of Standard ML (PDF). MIT Press. 1990 [2021-03-01]. (原始内容 (PDF)存档于2021-01-14).
- Robin Milner, Mads Tofte, Robert Harper, David MacQueen. The Definition of Standard ML, Revised (PDF). MIT Press. 1997 [2021-03-01]. ISBN 0-262-63181-4. (原始内容 (PDF)存档于2022-03-09).
- Robin Milner, Mads Tofte. Commentary on Standard ML. MIT Press. 1991 [2021-08-31]. ISBN 978-0-262-63137-2. (原始内容存档于2021-08-31).
- Emden R. Gansner, John H. Reppy. The Standard ML Basis Library (PDF). Cambridge University Press. 2004 [2021-09-17]. (原始内容 (PDF)存档于2022-01-29).
- Mads Tofte. Four Lectures on Standard ML (PDF). University of Edinburgh. 1989 [2021-09-04]. (原始内容 (PDF)存档于2022-01-28). Code examples in lectures. [2021-09-11]. 原始内容存档于2016-04-02.
- Mads Tofte. Essentials of Standard ML Modules. DIKU. 1996 [2021-09-04]. (原始内容存档于2021-09-04). The SML code (uuencoded compressed tar). [2021-09-18]. 原始内容存档于2005-03-07.
- Mads Tofte. Tips for Computer Scientists on Standard ML (Revised) (PDF). ITU. 2009 [2021-09-04]. (原始内容 (PDF)存档于2021-11-27). Examples. [2022-01-03]. 原始内容存档于2017-09-15.
- Robert Harper. Programming in Standard ML (PDF). Carnegie Mellon University. 2011 [2021-02-27]. (原始内容 (PDF)存档于2020-02-15). Examples. [2021-09-12]. (原始内容存档于2021-09-12).
- Michael J. C. Gordon. Introduction to Functional Programming. Cambridge University. 1996 [2021-09-11]. (原始内容存档于2021-04-11). Lecture notes. [2021-09-11]. (原始内容存档于2006-06-23).
- Lawrence Paulson. ML for the Working Programmer, 2nd Edition. Cambridge University Press. 1996 [2021-08-31]. ISBN 0-521-56543-X. (原始内容存档于2022-02-24). Sample programs. [2021-09-12]. (原始内容存档于2022-01-19).
- David MacQueen. Should ML be Object-Oriented? (PDF). 2002 [2021-09-10]. (原始内容 (PDF)存档于2021-10-29).
- David MacQueen, Robert Harper, John Reppy. The History of Standard ML. 2020 [2021-08-31]. (原始内容存档于2021-12-01).
- Andrew W. Appel. Compiling with Continuations. Cambridge University Press. 1992 [2022-01-03]. (原始内容存档于2022-01-03).
- Andrew W. Appel. Modern Compiler Implementation in ML. Cambridge University Press. 1998 [2022-01-12]. (原始内容存档于2022-01-12). Tiger compiler modules for programming exercises. [2021-09-12]. (原始内容存档于2022-05-06).
- Jeffrey D. Ullman. Elements of ML Programming, ML97 Edition. Prentice-Hall. 1998 [2021-12-30]. ISBN 0-13-790387-1. (原始内容存档于2022-03-12). Programs from the text. [2022-01-01]. (原始内容存档于2022-02-14).
- Anthony L. Shipman. Unix System Programming with Standard ML (PDF). 2001 [2021-09-01]. (原始内容 (PDF)存档于2021-01-21).
引用
[编辑]- ^ 1.0 1.1 1.2 1.3 Michael J. Gordon, Arthur J. Milner, Christopher P. Wadsworth. Edinburgh LCF: A Mechanised Logic of Computation. Lecture Notes in Computer Science, Vol. 78. Springer-Verlag, New York, NY, USA. 1979 [2021-12-28]. (原始内容存档于2021-12-28).
ML is a general purpose programming language. It is derived in different aspects from ISWIM, POP2 and GEDANKEN, and contains perhaps two new features. First, it has an escape and escape trapping mechanism, well-adapted to programming strategies which may be (in fact usually are) inapplicable to certain goals. Second, it has a polymorphic type discipline which combines the flexibility of programming in a typeless language with the security of compile-time type checking (as in other languages, you may also define your own types, which may be abstract and/or recursive); this is what ensures that a well-typed program cannot perform faulty proofs.
Michael J. C. Gordon. From LCF to HOL: a short history. 1996 [2021-02-27]. (原始内容存档于2016-09-05).Around 1973 Milner moved to Edinburgh University and established a project to build a successor to Stanford LCF, which was subsequently dubbed Edinburgh LCF. He initially hired Lockwood Morris and Malcolm Newey (both recent PhD graduates from Stanford) as research assistants. …… Milner, ably assisted by Morris and Newey, designed the programming language ML (an abbreviation for “Meta Language”). …… In 1975, Morris and Newey took up faculty positions at Syracuse University and the Australian National University, respectively, and were replaced by Chris Wadsworth and myself. The design and implementation of ML and Edinburgh LCF was finalised and the book “Edinburgh LCF” was written and published. In 1978, the first LCF project finished, Chris Wadsworth went off trekking in the Andes (returning to a permanent position at the Rutherford Appleton Laboratory) and I remained at Edinburgh supported by a postdoctoral fellowship and with a new research interest: hardware verification.
- ^ 2.0 2.1 M. Gordon, R. Milner, L. Morris, M. Newey, C. Wadsworth. A Metalanguage for Interactive Proof in LCF (PDF). 1978 [2021-08-31]. (原始内容 (PDF)存档于2021-05-06).
ML is a higher-order functional programming language in the tradition of ISWIM, PAL, POP2 and GEDANKEN, but differs principally in its handling of failure and, more so, of types.
- ^ 3.0 3.1 Bruce A. Tate, Fred Daoud, Ian Dees, Jack Moffitt. 3. Elm. Seven More Languages in Seven Weeks (PDF) Book version: P1.0-November 2014. The Pragmatic Programmers, LLC. 2014: 97, 101 [2021-09-04]. ISBN 978-1-941222-15-7. (原始内容 (PDF)存档于2021-03-15).
On page 101, Elm creator Evan Czaplicki says: 'I tend to say "Elm is an ML-family language" to get at the shared heritage of all these languages.' ["these languages" is referring to Haskell, OCaml, SML, and F#.]
- ^ 4.0 4.1 Troels Henriksen, Niels G. W. Serup, Martin Elsman, Fritz Henglein, Cosmin Oancea. Futhark: Purely Functional GPU-Programming with Nested Parallelism and In-Place Array Updates (PDF). Proceedings of the 38th ACM SIGPLAN Conference on Programming Language Design and Implementation. PLDI 2017. ACM. 2017 [2021-09-04]. (原始内容 (PDF)存档于2020-09-20).
- ^