我正在一个小项目中工作,其目标是给出Gcd的定义,它给出了两个数字的gcd以及结果正确的证明。但我无法给出Gcd的完整定义。 Idris 1.3.0中Gcd的定义是完全的,但是使用assert_total来强制总体性,这违背了我的项目的目的。有人有Gcd的总定义,不使用assert_total吗?
P.S。 - 我的代码已上传到https://github.com/anotherArka/Idris-Number-Theory.git
答案 0 :(得分:7)
我有一个使用Accessible
关系的版本,表明在每次递归调用时,您找到gcd的两个数字的总和会变小:https://gist.github.com/edwinb/1907723fbcfce2fde43a380b1faa3d2c#file-gcd-idr-L25
它依赖于此Prelude.Wellfounded
:
data Accessible : (rel : a -> a -> Type) -> (x : a) -> Type where
Access : (rec : (y : a) -> rel y x -> Accessible rel y) ->
Accessible rel x
一般的想法是,您可以通过明确说明变小的内容来进行递归调用,并在每个递归调用上提供它确实变得更小的证明。对于gcd
,它看起来像这样(gcdt
,因为gcd
是在前奏中的总版本):
gcdt : Nat -> Nat -> Nat
gcdt m n with (sizeAccessible (m + n))
gcdt m Z | acc = m
gcdt Z n | acc = n
gcdt (S m) (S n) | (Access rec)
= if m > n
then gcdt (minus m n) (S n) | rec _ (minusSmaller_1 _ _)
else gcdt (S m) (minus n m) | rec _ (minusSmaller_2 _ _)
sizeAccessible
在前奏中定义,允许您在此明确声明它是变小的输入之和。递归调用小于输入,因为rec
是Access rec
的参数。
如果您想更详细地了解正在发生的事情,可以尝试用孔替换minusSmaller_1
和minusSmaller_2
来电,看看您需要证明的内容:
gcdt : Nat -> Nat -> Nat
gcdt m n with (sizeAccessible (m + n))
gcdt m Z | acc = m
gcdt Z n | acc = n
gcdt (S m) (S n) | (Access rec)
= if m > n
then gcdt (minus m n) (S n) | rec _ ?smaller1
else gcdt (S m) (minus n m) | rec _ ?smaller2
例如:
*gcd> :t smaller1
m : Nat
n : Nat
rec : (y : Nat) ->
LTE (S y) (S (plus m (S n))) -> Accessible Smaller y
--------------------------------------
smaller1 : LTE (S (plus (minus m n) (S n))) (S (plus m (S n)))
我不知道任何文档Accessible
的详细信息,至少对于Idris(您可能会找到Coq的示例),但base
库中有更多示例Data.List.Views
1}},Data.Vect.Views
和Data.Nat.Views
。
答案 1 :(得分:0)
仅供参考:idris 1.3.0(可能是1.2.0)的实现是完全的,但是使用assert_total函数来实现这一点。
:printdef gcd
gcd : (a : Nat) ->
(b : Nat) -> {auto ok : NotBothZero a b} -> Nat
gcd a 0 = a
gcd 0 b = b
gcd a (S b) = assert_total (gcd (S b)
(modNatNZ a (S b) SIsNotZ))