#pragma once

#include <cmath>
#include <functional>
#include <iostream>
#include <functional>

using namespace std;

double ybot(double R);
double ytop(double R);

double local_ratio(double y, double z);
double next_cost(double y, double z);

double revert_exp_integral(double y);

double compute_R0();

double contiguous_ratio(double R);

double dichotomy(
    const function<double(double)>& f,  
    double c,                              
    double epsilon,                        
    double a,                              
    double b,                             
    bool is_increasing = true            
);

// Ternary search for a single parameter
double ternary_search(
    const function<double(const vector<double>&)>& func,
    vector<double>& params,
    vector<double>& valid_params, 
    int param_index,
    double left,
    double right,
    int max_iter = 100,
    double ternary_precision = 1e-8
);

// Multidimensional ternary search
vector<double> multidimensional_ternary_search(
    const function<double(const vector<double>&)>& func,
    vector<pair<double, double>> ranges,
    int max_iter = 1000,
    double iter_precision = 1e-6,
    double ternary_precision = 1e-8
);