数值计算问题，非智力测试

请教如下的最小化AX-B怎么解效率高


解一系列 求Xi minimize ||AiXi-Bi||

A是一特定的m X n非对称稀疏矩阵(1000=0， 其中系数=0的情况很少)得来，对不同的Ai这行的缩放系数可能每次不一样。Bi是一系列n X 3矩阵。
这个系列可能大于1000。

求教适合高效的算法，如果有可能请推荐C/C++算法包。谢谢