用遞迴來解（抱歉，我用javascript來做）：

<pre class="c" name="code">


<script>
var packages = {
	a:{unit:10,cost:150},
	b:{unit:5,cost:100},
	c:{unit:2,cost:60}
};

var res = {};

function compute(p,u,a) {
	for(var i in p) {
		a.push(i);
		if(u>p[i].unit) {
			compute(p,u-p[i].unit,a);
		} else {
			var r = JSON.stringify(a);
			r = JSON.parse(r);
			var c = 0;
			for(var j in r) {
				c += p[r[j]].cost;
			}
			if(res.total) {
				if(res.total>c) {
					res = {combo:r,total:c};
				}
			} else {
				res = {combo:r,total:c};
			}
		}
		a.pop();
	}
}
compute(packages, 15, []);
alert(JSON.stringify(res));
</script>

不過這樣窮舉組合（我這樣已經會排除掉無意義的組合了）的計算次數會隨著單位快速增長，並不適合在需要時動態計算。另外，php做遞迴的效率很差（跟javascript比較的話），所以建議你還是用迴圈展開來做，我只是做好玩的。（其實可以調整一下，在遞迴的條件中加入成本的評估，超過目前最低成本這條路徑就不繼續算下去，這樣除非剛好碰到最差狀況，否則應該可以改善計算的效率，你可以自己試試看。這樣算的時候，packages陣列的順序應該也會影響效率喔...）

應該任何一種程式語言都可以。

這是一個資料結構的最短路徑的問題，你上網 Google 一下應該會有詳細的解說跟範例。

建議找一個 c 語言的例子，然後你再改寫成 php

不過也很可能待會有那位大大，直接就給你答案，也說不定。

據我所知最有效率的做法應該是 Dyanmic programming (抱歉不清楚中文譯名)

http://www.avatar.se/molbioinfo2001/dynprog/dynamic.html

我的解法一點也不高深，
只用到數學的四則運算，
較深一點的數學就 +
所以會比較繁瑣些，
幾乎像是畫「九九乘法表」的結構做出來的。

直接把程式碼貼上，
雖然是 Ruby，但還蠻好辨識在寫什麼，
再自行用其他語言來改寫；
篇幅關係，部份推想過程的說明，稍後貼在 [討論](http://ithelp.ithome.com.tw/question/10045138?tab=opini<br />
on#ooa_hash) 中：

<pre class="c" name="code">class BoxPrice
  attr_reader :solutions, :every_price, :every_total, :best
  def initialize(a,b,c,d)
    @solutions = mycount(a,b,c,d)
    @prices = every_price(@solutions)
    @total = every_total(@prices)
    @min = which_min(@total)
    @best = @solutions[@min]
  end

  def maxx(a,b)
    # 以 5 * x >= 7 為例算出x 最小的 >= b 的值
    # x = 0,1,2
    # get the x Max is 2
    # here a = 5, b = 7
    int = b/a
    mod = b%a
    if mod.zero?
      return int
    else
      return int + 1
    end
  end

  def myrange(a,b)
    # 5 * x >= 7 將可能的x 列舉成陣列
    # x = 0,1,2
    # 2 is by maxx(a,b)
    arr = []
    maxx = maxx(a,b)
    0.step(maxx,1) do |x|
      arr << x
    end
    return arr
  end

  def mycount(a,b,c,d)
    # 排出x, y, z 的各可能組合
    # ax + by + cz = d
    puts "#{a}X + #{b}Y + #{c}Z = #{d}"
    boxes = []
    arange = myrange(a,d)
    arange.each do |i|
      #puts "X = #{i}"
      no_a_total = d - i * a
      brange = myrange(b, no_a_total)
        brange.each do |j|
          no_b_total = no_a_total - j*b
          crange = myrange(c, no_b_total)
            cmax = maxx(c, no_b_total)
            if cmax > 0
              puts "#{a} x #{i} + #{b} x #{j} + #{c} x #{cmax} = #{d} ( #{a*i + b*j + c*cmax} )"
              boxes << [i, j, cmax]
            else
              puts "#{a} x #{i} + #{b} x #{j} + #{c} x 0 = #{d} ( #{a*i + b*j} )"
              boxes << [i, j, 0]
            end
        end
    end
    return boxes
  end

  def every_price(arr)
    # 將各組合換算成價錢的小計
    return arr.collect{|gr| [ gr[0] * 150, gr[1] *100, gr[2] * 60]}
  end

  def every_total(arr)
    # 各組價錢的總計
    return arr.collect {|gr| sum = 0; gr.each{|x| sum+= x}; sum}
  end

  def which_min(arr)
    # 找出各組價錢的最小值
    return arr.index arr.min
  end

end

以下是執行的結果：

<pre class="c" name="code">irb -r3Box
>> b=BoxPrice.new(10,5,2,12)
10X + 5Y + 2Z = 12
10 x 0 + 5 x 0 + 2 x 6 = 12 ( 12 )
10 x 0 + 5 x 1 + 2 x 4 = 12 ( 13 )
10 x 0 + 5 x 2 + 2 x 1 = 12 ( 12 )
10 x 0 + 5 x 3 + 2 x 0 = 12 ( 15 )
10 x 1 + 5 x 0 + 2 x 1 = 12 ( 12 )
10 x 1 + 5 x 1 + 2 x 0 = 12 ( 15 )
=> #<BoxPrice:0x8cef570 @solutions=[[0, 0, 6], [0, 1, 4], [0, 2, 1], [0, 3, 0], [1, 0, 1], [1, 1, 0]], @prices=[[0, 0, 360], [0, 100, 240], [0, 200, 60], [0, 300, 0], [150, 0, 60], [150, 100, 0]], @total=[360, 340, 260, 300, 210, 250], @min=4, @best=[1, 0, 1]>
>> b.best
=> [1, 0, 1]

上述可用 b.solutions, b.prices, b.total 來分別列出
各組可能解，各組解的小計價格，各組解的總價。