Представлены численные подходы, с помощью которых можно получать точные решения задач о нескольких частицах, взаимодействующих как между собой, так и с несколькими ветвями бозонных возбуждений. Диаграммный метод Монте-Карло позволяет точно вычислить мацубаровскую функцию Грина, а с помощью метода стохастической оптимизации можно получить свободное от приближений аналитическое продолжение. Описаны результаты, получение которых невозможно традиционными методами: свойства возбужденных состояний при автолокализации, оптические спектры полярона во всех режимах связи, анализ пределов применимости к экситону приближений Френкеля и Ванье, особенности спектров фотоэмиссии дырки, взаимодействующей с решеткой в изоляторе Moттa.