不知你要用哪个工具实现这个过程呢？我随手用R模拟了下你描述的这个计算，你看看能用不？（未考虑代码优化）
a <- c(1:500)
b <- rnorm(500,150,18)
c <- rnorm(500,88,12)
t <- data.frame(a,b,c)
names(t) <- c("id","capacity","cost")
mydata <- t[order(t$cost),]


result <- ""
for (i in 1:500){
  temp <- sum(mydata$capacity[1:i])
  if (temp<=5000)  {
     result <- mydata[1:i,]
   }
   else break
  }

附上结果：
length(result$capacity)#输出结果的变量数目（计算到第33个符合条件）
[1] 33
> sum(mydata$capacity[1:33])#前33个之和4954
[1] 4954.459
> sum(mydata$capacity[1:34])#前34个之和5116，超5000
[1] 5115.677

可以用Stata吗？

没资料无法试试，但你可则 d = 1 的位置就是你的答案了！（当然，我不太认为会刚好等于 5000，所以你可能找不到 d = 1）

好的！感谢大神！