最优R-Forge项目

Question

我从Martello算法的背包代码开始研究垃圾箱的包装问题。它是用旧的Fortran IV或66码写的。我在results Optimist (管理员HansW.Borcher)发现了一个非常有趣的项目，通过R控制台，您可以调用它用Fortran 66编写的子例程，并运行它来检查结果。如果您希望以更现代的长度编写代码，并检查是否到达相同的结果，则这是有用的。我下载了R x64 3.3.1和优化包。我不知道如何从R中运行这个包--我是说:调用Fortran中的子例程，输入数据并查看view的结果。

有什么建议吗？

提前谢谢。爱德华多

Answer 1

原始的Martello和Toth Fortran例程可以在Optimist package项目的‘背包’包中获得(而不是在'adagio‘包中)。不幸的是，这些Fortran代码不能通过CRAN分发。原因是它们是在ACM许可下发布的，与GPL不兼容。我问Silvano Martello教授是否愿意改变许可证，但他不能也不想改变(正如他明确告诉我的那样)。

给你一个开始:我想你已经安装了R，更好的是，R和RStudio。启动R后，首先必须安装软件包一次，每次启动R时都要加载它：(您需要有一个可用的Fortran编译器，但我想您已经有了)。

> install.packages("knapsack", repos="http://R-Forge.R-project.org")
> library(knapsack)

然后您可以调用help或example on，例如背包函数。目前实现的功能是knapsack和subsetsum。帮助页将向您展示如何应用这些例程。p和w (利润和权重)必须是长度相等的整数值的向量，p[i]/w[i]是严格递减的序列：

> p = c(15, 100, 90, 60, 40, 15, 10,  1)
> w = c( 2,  20, 20, 30, 40, 30, 60, 10)
> cap = 102

现在您可以调用背包函数并显示结果：

> res = knapsack(p, w, cap)
> res
# [1]  1 2 3 4 6

还有子集和例程。Martello和Toth书中的其他代码还没有包装在“背包”中。但如果你有兴趣的话，这很容易做到。这些额外的例程是解决垃圾箱包装，分配和改变制造问题。