Пффф... я за 5 минут написал процедуру, убедился, что она даёт верные результаты, посмотрел время и шлёпнул её за ненадобностью.
У вас немеряное время оттого, что Вы запоминаете ВСЕ варианты и ВСЕ промежуточные значение в таблицах - в результате у Вас фактически всё время работает чтение и присвоение колонке объекта.
Зарезервируйте сразу массив необходимого размера и адресуйте свой двумерный массив результатов расчёта в одномерную переменную-массив, пересчёт индекса будет работать на порядок быстрее обращения к методу. А уже по окончании расчёта распихаете, буде надо, результаты куда-нить в таблицу.