julia - Julia 中的单项式向量

Question

假设你有一个向量元组 $a$，我想在 julia 中定义一个函数 p(x)=x^a。

例如，如果 a=(1,2,3)，则结果函数将为 x^1 *y^2 * z^3。

我想为任何元组提供一个通用方法，但我不知道合适的表示法。在我的代码中，我有一个元组数组，我想为数组中的每个元组定义一个单项式。

非常感谢您的合作。

Answer 1

这是你想要的吗？

julia> function genmonomial(t::Tuple{Vararg{Integer}})
           @assert !isempty(t) && all(>=(0), t)
           return (x...) -> begin
               @assert length(x) == length(t)
               return mapreduce(x -> x[1]^x[2], *, zip(x, t))
           end
       end
genmonomial (generic function with 1 method)

julia> f = genmonomial((2,3,4))
#1 (generic function with 1 method)

julia> f(2, 1, 1)
4

julia> f(1, 2, 1)
8

julia> f(1, 1, 2)
16

julia> f(2, 2, 2)
512

julia> f(1, 1)
ERROR: AssertionError: length(x) == length(t)

julia> genmonomial(())
ERROR: AssertionError: !(isempty(t)) && all((>=)(0), t)

julia> genmonomial((1.5,))
ERROR: MethodError: no method matching genmonomial(::Tuple{Float64})
Closest candidates are:
  genmonomial(::Tuple{Vararg{Integer}}) at REPL[1]:1