#include <stdio.h>
#include <stdlib.h>
#include <signal.h>
#include <unistd.h>
#include <math.h>

double f(double x);
void find_zero_x();
void ctrlc_handler(int signum);
void choose_path();
// int cnt_root();

double a, b, eps, x, *pa = &a, *pb = &b;
char choice, *ch = &choice;

int main()
{
    double tmp;
    /*
    struct sigaction sa;
    sa.sa_handler = ctrlc_handler;
    sa.sa_flags = 0;
    sigemptyset(&sa.sa_mask);
    */

    printf("Введите интервал [a, b] и точность eps: ");
    scanf("%lf %lf %lf", pa, pb, &eps);

    if (*pa > *pb) // Если пользователь перепутал границы
    {
        printf("Границы указаны неверно, запуск задачи с исправленным условием\n");
        tmp = *pa;
        *pa = *pb;
        *pb = tmp;
    }
    if (f(*pa) * f(*pb) > 0)
    {
        printf("В указанных границах нуля функции не существует.\n");
        printf("Нажмите Ctrl+C и поменяйте границы\n");
    }
    if (eps < 0) // Если пользователь указал отрицательную точность
    {
        eps = eps * (-1);
        printf("Точность отрицательна, условие приведено в порядок, запуск решения...\n");
    }
    if (*pa == *pb)
    {
        if (f(*pa) == 0)
        {
            printf("\nКорень уравнения: %lf\n", *pa);
            printf("Работа программы завершена.\n");
            exit(0);
        }
        else
        {
            printf("Нажмите Ctrl+C и поменяйте границы\n");
        }
    }

    find_zero_x();
    /*
    if (sigaction(SIGINT, &sa, NULL) == -1)
    {
        perror("sigaction");
        return 1;
    }
    */
    printf("Ожидание сигнала Ctrl+C...\n");
    signal(SIGINT, ctrlc_handler);
    while (1)
    {
        sleep(1);
    }

    return 0;
}

double f(double x) // Функция обязана быть линейной
{
    return x * x - 16;
}

void find_zero()
{
    double x;
    x = (*pa + *pb) / 2; // метод деления отрезка пополам
    if (f(*pa) * f(x) < 0)
        *pb = x;
    else
        *pa = x;
    printf("Текущее приближение: %lf\n", x);

    if (fabs(b - a) < eps)
    {
        printf("\nКорень уравнения: %lf\n", x);
        printf("Работа программы завершена.\n");
        exit(0);
    }
    else
    {
        return;
    }
}

void ctrlc_handler(int signum)
{
    printf("\nПолучен сигнал Ctrl+C\n");
    choose_path();
}

void choose_path()
{
    /*
    printf("\nПродолжить поиск корня? (C - продолжить, A - закончить работу программы, R - начать поиск на другом отрезке): \n");
    choice = 'Q';
    *ch = getchar();

    if (*ch == 'A')
    {
        printf("Работа программы завершена.\n");
        exit(0);
    }

    else if (*ch == 'C')
    {
        find_zero_x();
        printf("Ожидание сигнала Ctrl+C...\n");
    }
    else if (*ch == 'R')
    {
        printf("Введите новые значения границ a и b: ");
        scanf("%lf %lf", pa, pb);
        find_zero_x();
        printf("Ожидание сигнала Ctrl+C...\n");
    }
    else
    {
        printf("Ошибка ввода, повторная попытка...");
        choose_path();
    }
    */
    printf("\nПродолжить поиск корня?\
        (C - продолжить, A - закончить работу программы, R - начать поиск на другом отрезке): \n");
    do
    {

        *ch = getchar();
        switch (*ch)
        {
        case 'C':
            find_zero_x();
            printf("Ожидание сигнала Ctrl+C...\n");
            return;
        case 'A':
            printf("Работа программы завершена.\n");
            exit(0);
        case 'R':
            printf("Введите новые значения границ a и b: ");
            scanf("%lf %lf", pa, pb);
            if (f(*pa) * f(*pb) > 0)
            {
                printf("В указанных границах нуля функции не существует.\n");
                printf("Нажмите Ctrl+C и поменяйте границы\n");
            }
            if (*pa == *pb)
            {
                if (f(*pa) == 0)
                {
                    printf("\nКорень уравнения: %lf\n", *pa);
                    printf("Работа программы завершена.\n");
                    exit(0);
                }
                else
                {
                    printf("Нажмите Ctrl+C и поменяйте границы\n");
                }
            }
            find_zero_x();
            printf("Ожидание сигнала Ctrl+C...\n");
            return;
        }
    } while (*ch != 'A' || *ch != 'R' || *ch != 'C');
    printf("Ошибка ввода, повторная попытка...");
    choose_path();
}
/*
int cnt_root()
{
    if (fabs(*pb - *pa) < eps)
        return (f(*pa) * f(*pb) <= 0) ? 1 : 0;
    else
    {
        double m = (*pa + *pb) / 2;
        if(f(*pa)*f(m)<=0)
        {
            return cnt_root(*pa, m);
        } else {
            return cnt_root(m, *pb);
        }
    }
}
*/

void find_zero_x()
{
    float mid = (*pa + *pb) / 2;
    mid = (*pb + *pa) / 2;
    if (f(mid) > 0)
    {
        *pa = (f(*pa) > f(*pb)) ? mid : *pa;
        *pb = (f(*pa) < f(*pb)) ? mid : *pb;
    }
    else
    {
        *pa = (f(*pa) < f(*pb)) ? mid : *pa;
        *pb = (f(*pa) > f(*pb)) ? mid : *pb;
    }

    printf("Текущее приближение: %lf\n", *pa);

    if (fabs(f(mid)) < eps)
    {
        printf("\nКорень уравнения: %lf\n", *pa);
        printf("Работа программы завершена.\n");
        exit(0);
    }
    else
    {
        return;
    }
}