F # Статические ограничения типа пользователя

Я пытаюсь определить функцию, factorize, которая использует ограничения структурного типа (требует статических членов Zero, One, + и /), аналогичных Seq.sum, так что их можно использовать с int, long, bigint, и т.д. Я не могу получить синтаксис правильно и не могу найти много ресурсов по этому вопросу. Это то, что у меня есть, пожалуйста, помогите.

let inline factorize (n:^NUM) =
    ^NUM : (static member get_Zero: unit->(^NUM))
    ^NUM : (static member get_One: unit->(^NUM))
    let rec factorize (n:^NUM) (j:^NUM) (flist: ^NUM list) = 
        if n = ^NUM.One then flist
        elif n % j = ^NUM.Zero then factorize (n/j) (^NUM.One + ^NUM.One) (j::flist)
        else factorize n (j + ^NUM.One) (flist)
    factorize n (^NUM.One + ^NUM.One) []

Ответ 1

Вот как я его напишу:

module NumericLiteralG = begin
  let inline FromZero() = LanguagePrimitives.GenericZero
  let inline FromOne() = LanguagePrimitives.GenericOne
end

let inline factorize n = 
  let rec factorize n j flist =  
    if n = 1G then flist 
    elif n % j = 0G then factorize (n/j) j (j::flist) 
    else factorize n (j + 1G) (flist) 
  factorize n (1G + 1G) [] 

Тип, предполагаемый для факторизации, является слишком общим, но функция будет работать так, как вы ожидали. Вы можете принудительно ввести более строгую подпись и набор ограничений, если хотите, добавив явные типы в некоторые из общих выражений:

let inline factorize (n:^a) : ^a list = 
  let (one : ^a) = 1G
  let (zero : ^a) = 0G
  let rec factorize n (j:^a) flist =  
    if n = one then flist 
    elif n % j = zero then factorize (n/j) j (j::flist) 
    else factorize n (j + one) (flist) 
  factorize n (one + one) []

Ответ 2

Вдохновленный ответом kvb с помощью NumericLiterals, мне было предложено разработать подход, который позволил бы нам форсировать сигнатуры типа "здорового" типа, не добавляя подробные аннотации типов.

Сначала мы определяем некоторые вспомогательные функции и тип оболочки для языковых примитивов:

let inline zero_of (target:'a) : 'a = LanguagePrimitives.GenericZero<'a>
let inline one_of (target:'a) : 'a = LanguagePrimitives.GenericOne<'a>
let inline two_of (target:'a) : 'a = one_of(target) + one_of(target)
let inline three_of (target:'a) : 'a = two_of(target) + one_of(target)
let inline negone_of (target:'a) : 'a = zero_of(target) - one_of(target)

let inline any_of (target:'a) (x:int) : 'a =
    let one:'a = one_of target
    let zero:'a = zero_of target
    let xu = if x > 0 then 1 else -1
    let gu:'a = if x > 0 then one else zero-one

    let rec get i g = 
        if i = x then g
        else get (i+xu) (g+gu)
    get 0 zero 

type G<'a> = {
    negone:'a
    zero:'a
    one:'a
    two:'a
    three:'a
    any: int -> 'a
}    

let inline G_of (target:'a) : (G<'a>) = {
    zero = zero_of target
    one = one_of target
    two = two_of target
    three = three_of target
    negone = negone_of target
    any = any_of target
}

Тогда имеем:

let inline factorizeG n = 
    let g = G_of n
    let rec factorize n j flist =  
        if n = g.one then flist 
        elif n % j = g.zero then factorize (n/j) j (j::flist) 
        else factorize n (j + g.one) (flist) 
    factorize n g.two []

[ Изменить: из-за очевидной ошибки с F # 2.0/.NET 2.0, factorizen, factorizeL и factorizeI ниже выполняется значительно медленнее, чем factorizeG при компиляции в режиме Release, но в противном случае выполняется немного быстрее, чем Ожидаемый - см. Вопрос о производительности F #: что делает компилятор?]

Или мы можем сделать это еще на несколько шагов (вдохновленный экспертом F #, стр .110):

let inline factorize (g:G<'a>) n =   //'
    let rec factorize n j flist =  
        if n = g.one then flist 
        elif n % j = g.zero then factorize (n/j) j (j::flist) 
        else factorize n (j + g.one) (flist) 
    factorize n g.two []

//identical to our earlier factorizeG
let inline factorizeG n = factorize (G_of n) n

let gn = G_of 1  //int32
let gL = G_of 1L //int64
let gI = G_of 1I //bigint

//allow us to limit to only integral numeric types
//and to reap performance gain by using pre-computed instances of G
let factorizen = factorize gn
let factorizeL = factorize gL
let factorizeI = factorize gI

Кроме того, здесь представлена ​​расширенная версия kvb NumericLiteralG, которая позволяет нам использовать "2G", "-8G" и т.д. Хотя я не мог понять, как реализовать стратегию memoization (хотя это должно быть выполнимо для G.any).

module NumericLiteralG = 
    let inline FromZero() = LanguagePrimitives.GenericZero
    let inline FromOne() = LanguagePrimitives.GenericOne
    let inline FromInt32(n:int):'a =
        let one:'a = FromOne()
        let zero:'a = FromZero()
        let nu = if n > 0 then 1 else -1
        let gu:'a = if n > 0 then one else zero-one

        let rec get i g = 
            if i = n then g
            else get (i+nu) (g+gu)
        get 0 zero 

Ответ 3

Во-первых, вот тривиальный пример, который показывает, как должен выглядеть синтаксис:

let inline zero< ^NUM when ^NUM : (static member get_Zero: unit-> ^NUM)> 
    (n:^NUM) = 
  (^NUM : (static member get_Zero : unit -> ^NUM) ())

В некоторых случаях вам не нужно явно писать ограничения (компилятор F # на самом деле предупреждает вас об этом, если вы пишете выше), потому что некоторые статические члены хорошо известны компилятору и существуют стандартные функции для их использования. Таким образом, вы можете использовать эту функцию, и компилятор выведет ограничение:

let inline zero (n:^T) = 
  LanguagePrimitives.GenericZero< ^T > 

К сожалению, это действительно вам не поможет, потому что рекурсивные функции не могут быть объявлены как inline (по понятным причинам - компилятор не может встроить функцию во время компиляции, потому что он не знает, сколько раз), поэтому статические ограничения, вероятно, недостаточно эффективны для вашей проблемы.

[ EDIT: это действительно возможно для некоторых функций (см. ответ kvb)]

Я думаю, вам понадобится NumericAssociations вместо этого, которые обсуждались в в этом вопросе (они обрабатываются во время выполнения, поэтому они медленнее - но используются для реализации, например, матричного типа F # - матрица может кэшировать динамически полученную информацию, поэтому она достаточно эффективна).