Функция Rnknack - PullRequest
       42

Функция Rnknack

0 голосов
/ 26 ноября 2018

Я запускаю программу R из статьи , где используется функция mknapsack из пакета adagio, и все хорошо.Но если я хочу использовать случайные значения, я получаю сообщение об ошибке «Ошибка возникла».

У меня есть программа:

n=16
m=5
max=700
min = 10

planks_we_have = floor(runif(n=m, min = 100, max = max))
planks_we_want = floor(runif(n=n, min = min, max = 16))

library(adagio)
# mknapsack calling signature is: mknapsack(values, weights, capacities)
solution <- mknapsack(planks_we_want, planks_we_want, planks_we_have)
# Above I added +1 cm  to each length to compensate for the loss when sawing.
solution$ksack
# Now pretty printing what to cut so that we don't make mistakes...
assignment <- data.frame(cut_this = planks_we_have[solution$ksack], into_this = planks_we_want)
t(assignment[order(assignment[,1]), ])

Результат:

Warning
In mknapsack(planks_we_want, planks_we_want, planks_we_have) :
  Error condition raised: check input data ...!

Error
In data.frame(cut_this = planks_we_have[solution$ksack], into_this = planks_we_want) :
  Arguments imply different numbers of lines: 0, 5

Iне понимаю в чем причинаИсходный код функции рюкзака мне ничего не дает:

function (p, w, k, bck = -1) 
{
    stopifnot(is.numeric(p), is.numeric(w), is.numeric(k))
    if (any(w <= 0)) 
        stop("'weights' must be a vector of positive numbers.")
    if (any(p <= 0)) 
        stop("'profits' must be a vector of positive numbers.")
    if (any(floor(p) != ceiling(p)) || any(floor(w) != ceiling(w)) || 
        any(floor(k) != ceiling(k)) || any(p >= 2^31) || any(w >= 
        2^31) || any(k >= 2^31)) 
        stop("All inputs must be positive integers < 2^31 !")
    n <- length(p)
    m <- length(k)
    if (length(w) != n) 
        stop("Profit 'p' and weight 'w' must be vectors of equal length.")
    xstar <- vector("integer", n)
    vstar <- 0
    num <- 5 * m + 14 * n + 4 * m * n + 3
    wk <- numeric(n)
    iwk <- vector("integer", num)
    S <- .Fortran("mkp", as.integer(n), as.integer(m), as.integer(p), 
        as.integer(w), as.integer(k), bs = as.integer(bck), 
        xs = as.integer(xstar), vs = as.integer(vstar), as.numeric(wk), 
        as.integer(iwk), as.integer(num), PACKAGE = "adagio")
    if (S$vs < 0) 
        warning("Error condition raised: check input data ...!")
    return(list(ksack = S$xs, value = S$vs, btracks = S$bs))
}

Версии:

R - 3.4.1
Adagio - 0.7.1

1 Ответ

0 голосов
/ 27 ноября 2018

Пожалуйста, прочитайте сначала страницу справки, если у вас есть проблемы с функцией.Глядя на возвращенное решение, оно имеет код ошибки vs=-7, а справка говорит: «vs = -7, если массив k отсортирован неправильно».Сортировка вектора мощностей может привести к другой ошибке, например, если все предметы можно положить в один рюкзак.Конечно, все это зависит от сгенерированных случайных чисел (лучше исправить случайные числа перед тем, как спрашивать).

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...