With b = a/a.sum(), you can multiply any solution of the vector a by any scalar without affecting b. ie, b = [1,1,1]/(1+1+1) yields the same end value as b = [.5,.5,.5]/(.5+.5+.5). This means that the inference will drift endlessly between all and any homothetic solution for a.
Hope that helps?