范围最小查询 - Clojure

时间:2012-02-14 09:38:37

标签: c# clojure

我正在开发一个Clojure项目,该项目将数组a作为输入,并找到每个[i,j]i范围内j的最小值,在O(n)预处理中,每个查询O(1)。 (preproccessing需要O(n*log(n)),但通过使用并发(pmap)并将数组划分为n/log n数组,我们可以在O(n)中解决此问题

因此,我决定将数组表示为向量,将矩阵表示为向量的向量。

这是C#中的一个功能:

    void build_RMQ_log(int[] a)
    {
        ST = new int[2 * a.Length][];
        for (int i = 0; i < 2 * a.Length; i++)
            ST[i] = new int[(int)Math.Log(a.Length, 2)];

        for (int i = 0; i < a.Length; i++)
        {
            ST[i][0] = i;
            ST[i + a.Length - 1][0]=0;
        }

        for (int j = 1; j < (int)Math.Log(a.Length, 2); j++)
        {
            for (int i = 0; i < a.Length; i++)
            {
                if (a[ST[i][j - 1]] <= a[ST[i + (int)Math.Pow(2, j - 1)][j - 1]])
                    ST[i][j] = ST[i][j - 1];
                else
                    ST[i][j] = ST[i + (int)Math.Pow(2, j - 1)][j - 1];
            }
        }
    }
    i

这就是我在Clojure中写的:

    ;build_RMQ_log(int[] a)

    (defn check [row1 row2 arr j]
                  (let [x1 (nth arr (nth row1 j))
                    x2 (nth arr (nth row2 (- j 1)))
                    x3 (nth arr (nth row1 (- j 1)))]

                  (if (<= x1 x2)
                     (assoc row1 j (nth row1 (- j 1)))
                     (assoc row1 j (nth row2 (- j 1))))))

    (defn apply_fn [ST a row r]
    (let [len (count a)
     cnt (/ len (log_ len 2))]
      (loop[ii 0 r 0]
        (if (= ii (- cnt 1))
          ST
         (recur (inc ii) (check row (nth ST (+ r (Math/pow 2 (- ii 1)))) a ii))))))


   (defn build_RMQ_log [a]
     (let [len (count a)
           ST (diag_n_m len (log_ len 2))
           r 0]
     (pmap  (fn [row] (apply_fn (ST a row r))) ST )))

首先,当我尝试运行它时,它会向我显示此错误:

    #<IllegalArgumentException java.lang.IllegalArgumentException: Wrong number of  args (3) passed to: PersistentVector>

此外,我写的代码没有做我想要的,因为我不知道如果r并行工作,我怎么能改变apply_fn(代表行号)的值。即喜欢它改变c#:

for (int i = 0; i < a.Length; i++)

ri中的for - c#中的循环相似

先谢谢。

1 个答案:

答案 0 :(得分:2)

如果我没有正确理解您,您希望将r递增到apply_fn的每次调用。你可以试试这个:

(defn build_RMQ_log [a]
  (let [len (count a)
        ST (diag_n_m len (log_ len 2))
        rs (range)]
    (pmap  (fn [row r] (apply_fn (ST a row r))) ST rs)))

也就是说,您将两个集合传递给pmap,其中第二个集合是增量整数的无限集合(即[0, 1, 2, 3, ...])。