Встроенная функция FrobeniusSolve
решает случай, когда c1, c2, ..., cn являются положительными целыми числами (а правая часть не n):
In[1]:= FrobeniusSolve[{2, 3, 5, 6}, 13]
Out[1]= {{0, 1, 2, 0}, {1, 0, 1, 1}, {1, 2, 1, 0}, {2, 1, 0, 1}, {2,
3, 0, 0}, {4, 0, 1, 0}, {5, 1, 0, 0}}
Это тот случай, когда вам нужен или вам нужен отрицательный c1, c2, ..., cn также?