Радиоактивная порода

t = input().split() # названия элементов 
h = {t[i]: int(t[i + 1]) for i in range(0, len(t), 2)} # периоды полураспада
v = [i for i in zip(input().split(), map(float, input().split()))] # исходные излучения элементов
m = float(input()) # искомое суммарное излучение
i = 0 # номер дня
while sum([t[1] * 0.5 ** (i / h[t[0]]) for t in v]) > m: i += 1
print(i)
print(*[t[1] * 0.5 ** (i / h[t[0]]) for t in v])

(i / h[t[0]])

(i // h[t[0]])

arr = [elem for elem in input().split()]
arr = [(arr[i * 2:(i + 1) * 2][0], int(arr[i * 2:(i + 1) * 2][1])) for i in range(len(arr) // 2)]
key_array = []
 
Dict = dict()
main_decay = dict(arr)
now_decay = main_decay.copy()
 
for elem in input().split():
    if elem not in Dict:
        Dict[elem] = []
 
    key_array.append(elem)
 
elem = [float(x) for x in input().split()]
 
for i in range(len(elem)):
 
    if [] not in Dict[key_array[i]]:
        Dict[key_array[i]].append([])
 
    Dict[key_array[i]][Dict[key_array[i]].index([])].append(elem[i])
 
normal = float(input())
nod = int(list(main_decay.values())[0])
 
for elem in list(main_decay.values())[1:]:
    a, b = int(nod), int(elem)
 
    while a != 0:
        a, b = b % a, a
 
    nod = b
 
now = 0
 
while normal < sum(sum(sum(list(Dict.values()), []), [])):
    now += nod
    for key in main_decay.keys():
        if not now % main_decay[key] and key in Dict:
            for i in range(len(Dict[key])):
                for j in range(len(Dict[key][i])):
                    Dict[key][i][j] /= 2
 
answer_array = []
 
for elem in key_array:
    answer_array.append(Dict[elem][0][0])
    del Dict[elem][0]
 
print(now)
print(" ".join([str(x) for x in answer_array]))

t = list(map(int, input().split()[1::2])) # названия элементов нам нафиг не нужны
input() # И снова они нам не нужны...
v = list(map(float, input().split())) # исходные излучения элементов
m = float(input()) # искомое суммарное излучение
i = 0 # номер дня
while True:
    if sum(map(lambda x: v[x] / 2 ** (i // t[x]) , range(len(t)))) <= m:
        print(i)
        print(*map(lambda x: v[x] / 2 ** (i // t[x]), range(len(t))))
        break
    else:
        i += 1

# Считываем данные о периодах полураспада и создаем словарь main_decay
arr = [elem for elem in input().split()]  # Ввод периодов полураспада и значений элементов
arr = [(arr[i * 2:(i + 1) * 2][0], int(arr[i * 2:(i + 1) * 2][1]))
       for i in range(len(arr) // 2)]  # Группируем элементы и их периоды полураспада в кортежи
main_decay = dict(arr)  # Создаем словарь с периодами полураспада
 
# Инициализация переменных
key_array = []  # Список для хранения названий радиоактивных элементов
Dict = dict()  # Словарь для хранения активностей радиоактивных элементов
 
# Создаем словарь Dict для хранения активностей радиоактивных элементов
for elem in input().split():
    if elem not in Dict:
        Dict[elem] = []  # Инициализируем пустой список для каждого элемента
 
    key_array.append(elem)  # Добавляем элемент в список
 
elem = [float(x) for x in input().split()]  # Вводим активности каждого элемента
 
# Заполняем словарь Dict значениями активности каждого элемента
for i in range(len(elem)):
    if [] not in Dict[key_array[i]]:
        Dict[key_array[i]].append([])  # Инициализируем пустой список для каждого значения активности
 
    Dict[key_array[i]][Dict[key_array[i]].index([])].append(elem[i])  # Добавляем активность в список
 
normal = float(input())  # Допустимый радиационный фон
 
# Находим наибольший общий делитель (НОД) для периодов полураспада
nod = int(list(main_decay.values())[0])
 
for elem in list(main_decay.values())[1:]:
    a, b = int(nod), int(elem)
 
    # Алгоритм Евклида для нахождения НОД
    while a != 0:
        a, b = b % a, a
 
    nod = b
 
now = 0  # Переменная для отслеживания времени
 
# Пока текущий радиационный фон больше допустимого, уменьшаем активность элементов
while normal < sum(sum(sum(list(Dict.values()), []), [])):
    now += nod  # Увеличиваем время на НОД периодов полураспада
 
    # Уменьшаем активность каждого элемента согласно периодам полураспада
    for key in main_decay.keys():
        if not now % main_decay[key] and key in Dict:
            for i in range(len(Dict[key])):
                for j in range(len(Dict[key][i])):
                    Dict[key][i][j] /= 2  # Уменьшаем активность вдвое
 
answer_array = []
 
# Формируем ответ в виде массива активностей элементов
for elem in key_array:
    answer_array.append(Dict[elem][0][0])
    del Dict[elem][0]
 
# Выводим результат
print(now)  # Выводим время, прошедшее с начала отсчета
print(" ".join([str(x) for x in answer_array]))  # Выводим окончательное распределение активности элементов