Как ускорить работу скрипта по перебору, комбинации сумм

resultdict = {
         '13098194': 5413033.09,
         '14024911': 4799851.78,
         '15478940': 4274028,
         '12794945': 3606633.61,
         '16124010': 3083362.91,
         '17045813': 2914115.9,
         '17355657': 2873580.38,
         '17194240': 2817811.2
....
}

# -*- coding: utf8 -*-
import itertools, time
from openpyxl import load_workbook
 
start_time = time.time() #Время запуска скрипта
 
wb = load_workbook('Тест.xlsx') #
sheet_ranges = wb['Sheet1']
 
total = float(sheet_ranges['E1'].value)
 
temp1 = {} #Временный словарь для загрузки данных
resultdict = {} #Словарь без дублей и сортировкой
 
for row in sheet_ranges.iter_rows():
    if row[0].value is not None and row[1].value is not None and row[1].value <= total:
        column_A = str(row[0].value)
        column_B = row[1].value
        temp1[column_A] = column_B  #Заполняем словарь данными из Excel.
temp2 = {}
[temp2.update({k:v}) for k,v in temp1.items() if v not in temp2.values()]
 
#Сортировка от большего к меньшему, если total > 1 000 000
if int(total) > 1000000:
    for i in sorted(temp2.items(),reverse=True, key=lambda para: para[1]):
        resultdict[i[0]]=i[1]
#Сортировка от меньшего к большему, если total < 1 000 000
elif int(total) < 1000000:
    for i in sorted(temp2.items(),reverse=False, key=lambda para: para[1]):
        resultdict[i[0]]=i[1]
 
# for key, value in resultdict.items():
#     print(key, value)
 
for r in range(1, len(resultdict) + 1):
    for keys_combination in itertools.combinations(resultdict, r):
        if sum([resultdict[key] for key in keys_combination]) == total:
            with open("result.txt", "w") as f:
                f.write({key: resultdict[key] for key in keys_combination},file=f)
 
 
print()
print('Желаемая сумма:',total)
print('Количество строк из файла Excel:',len(temp1))
print('Количество строк без дублей и сортировкой:',len(resultdict))
 
 
 
 
 
end_time = time.time() #Время окончания работы скрипта
print()
print(f"На вычисление потребовалось {end_time-start_time:.2f} секунд")

import numpy as np
n = 500 # количество элементов
find_sum = 1495. # требуемая сумма
a = np.array(list(range(n)), dtype=float) # входной массив 
print(a)
np.random.shuffle(a)
print(a)
c = np.add.outer(a,a)
i, j = np.where(c==find_sum)
for k in zip(i,j):
    print(k, c[k])
 
##### если i, j "пустые", то
c[np.tril_indices(n)] = np.nan
c[c>find_sum] = np.nan
c = np.add.outer(c,a)
i, j, m = np.where(c==find_sum)
for k in zip(i,j,m):
    print(k, c[k])

from bisect import bisect_left
from operator import itemgetter
from typing import Dict, List
from time import time
 
from openpyxl import load_workbook
 
 
def get_data(file_name: str) -> Dict[str, float]:
    work_book = load_workbook(file_name)
    sheet = work_book['Sheet1']
    pairs = {row[0].value: float(row[1].value) for row in sheet.iter_rows() if row[0].value and row[1].value}
    return dict(sorted(pairs.items(), key=itemgetter(1)))
 
 
def calculate(pairs: Dict[str, float], value: float) -> List[str]:
    sorted_pairs = list(pairs.items())
    min_ = sorted_pairs[0][1]
    sum_ = round(sum(pairs.values()), 2)
    if value < min_:
        raise ValueError(f'Значение {value} меньше минимального значения в таблице({min_})!')
    if value > sum_:
        raise ValueError(f'Значение {value} больше суммы всех значений в таблице({sum_})!')
    if value == sum_:
        print(f'Нужное число({value}) это сумма всех ячеек!')
        return list(pairs.keys())
    sorted_values = list(pairs.values())
    index_ = index(sorted_values, value)
    if index_ is not None:
        return [sorted_pairs[index_][0]]
    pair = _find_sum_of_two(value, sorted_values)
    if pair:
        return [sorted_pairs[pair[0]][0], sorted_pairs[pair[1]][0]]
    triple = _find_sum_of_three(value, sorted_values)
    if triple:
        return [sorted_pairs[triple[0]][0], sorted_pairs[triple[1]][0], sorted_pairs[triple[2]][0]]
    return []
 
 
def _find_sum_of_three(value, sorted_values):
    for initial_value in sorted((e for e in sorted_values if e < value), reverse=True):
        to_find = round(value - initial_value, 2)
        pairs = _find_sum_of_two(to_find, sorted_values)
        if pairs:
            return index(sorted_values, initial_value), *pairs
    return tuple()
 
 
def _find_sum_of_two(value, sorted_values: List[float]) -> tuple:
    for initial_value in sorted((e for e in sorted_values if e < value), reverse=True):
        to_find = round(value - initial_value, 2)
        index_ = index(sorted_values, to_find)
        if index_ is not None:
            return index(sorted_values, initial_value), index_
    return tuple()
 
 
def index(a, x):
    i = bisect_left(a, x)
    if i != len(a) and a[i] == x:
        return i
    return None
 
 
if __name__ == '__main__':
    pairs = get_data('Тестирование.xlsx')
    start = time()
    print(calculate(pairs, 81810.91))  # ['16579079', '15473880', '15634093']
    print(calculate(pairs, 51502.48))  # ['11550029', '11550031']
    print(calculate(pairs, 1163258.98))  # ['16278978']
    print(f'Time elapsed: {time() - start}')