Форум Всеукраїнської інтернет-олімпіади NetOI

#include <iostream>
#include <cmath>
#include <iomanip>
using namespace std;
int main()
{
     double  L, a1, a2, v, t;
    cin >> L >> a1 >> a2 >> v >> t;
    
    double T, t1,t2,t3,tx,vx;
    tx = (L - v*v/2/a1 - v*t - v*v/2/a2) / v;
    if (tx < 0) 
    {
        t3 = t;
        double a,b,c;
         a=a1/2*(1.0 + a1/a2);
        b=a1*t;
        c=-L;
         double x1,x2;
     x1=(-b+sqrt(pow(b,2.0)-4*a*c))/(2*a);
     x2=(-b-sqrt(pow(b,2.0)-4*a*c))/(2*a);
         if (x1<0)
        {
            t1 = x2;
            t2 = a1*x2/a2;
        }
    
        else if (x2<0)
        {
            t1 = x1;
            t2 = a1*x1/a2;
        }
        else
        {
            if (x1<x2)
            {
                t1 = x1;
                t2 = a1*x1/a2;
            }
            else
            {
                t1 = x2;
                t2 = a1*x2/a2;
            }
        }
    }
    else
    {
    t2 = t + tx;
    t1 = sqrt((L - v*(t2)   -  v*v/2/a2) * 2 / a1);
    t3 = sqrt((L - v*(t2)   -  v*v/2/a1) * 2 / a2);
    }
 
    T = t1 + t2 + t3;
    cout << fixed << setprecision(8) << T;
 return 0;
}

#include <iostream>
#include <cmath>
#include <iomanip>
using namespace std;
int main()
{
     double  L, a1, a2, v, t;
    cin >> L >> a1 >> a2 >> v >> t;
    
    double T, t1,t2,t3,tx,vx;
//T-загальний час
//t1,t2,t3 - тривалість руху на трьох ділянках
//tx-час розгону до максимально можливої швидкості за умовою
//vx-до якої максимальної швидкості можна розігнатись реально 
    
     //Щоб проїхати найшвидше потрібно розігнатись до максимальної швидкості v
    //s1 = a1*t1*t1/2   s1 = v*v/2/a1
    //s2 = v*(t+tx)  
    //s3 = a2*t3*t3/2   s3 = v*v/2/a2
    //L= v*v/2/a1 + v*(t+tx)   + v*v/2/a2
    //L= v*v/2/a1 + v*t + v*tx + v*v/2/a2
 
    //звідки: 
    tx = (L - v*v/2/a1 - v*t - v*v/2/a2) / v;
    //cout << tx << endl;
    
    //Якщо не можна розігнатись до максимальної швидкості
    if (tx < 0) 
    {
        //тоді час рівномірного руху t3 повинен бути мінімальним
        t3 = t;
 
        //s1= vx*vx/2/a1 = a1*t1*t1/2
        //s2 = vx*t
        //s3 = vx*vx/2/a2 = a2*t3*t3/2
        //vx = a1*t1 = a2*t3 => t3 = a1*t1/a2
 
        //L = s1 + s2 + s3
 
        //знаходимо t1
        //L = a1*t1*t1/2 + a1*t1*t + a2*(a1*t1/a2)*(a1*t1/a2)/2
        //L = (a1/2)*t1*t1 + a1*t*t1 + (a1*a1/a2/2)*t1*t1
        //L = (a1/2)*t1*t1 + a1*t*t1 + (a1*a1/a2/2)*t1*t1
        //a1/2(1 + a1/a2)*t1*t1 + a1*t*t1 - L = 0
        double a,b,c;
 
        a=a1/2*(1.0 + a1/a2);
        b=a1*t;
        c=-L;
 
    
    //находим t1
        double x1,x2;
 
    x1=(-b+sqrt(pow(b,2.0)-4*a*c))/(2*a);
 
    x2=(-b-sqrt(pow(b,2.0)-4*a*c))/(2*a);
 
     
        if (x1<0)
        {
            t1 = x2;
            //находим t2
            t2 = a1*x2/a2;
        }
    
        else if (x2<0)
        {
            t1 = x1;
            //находим t2
            t2 = a1*x1/a2;
        }
        else
        {
            if (x1<x2)
            {
                t1 = x1;
                //находим t2
                t2 = a1*x1/a2;
            }
            else
            {
                t1 = x2;
                //находим t2
                t2 = a1*x2/a2;
            }
        }
    }
    //Якщо можна розігнатись до максимальної швидкості
    else
    {
    //Знаходимо t2
    t2 = t + tx;
 
    //s1 + s2 + s3 – L = 0
    
    //знаходимо t1
    //a1*t1*t1/2 + v*(t2)  + v*v/2/a2 - L = 0
    t1 = sqrt((L - v*(t2)   -  v*v/2/a2) * 2 / a1);
 
    //знаходимо t3
    //v*v/2/a1+ v*(t2)  + a2*t3*t3/2  - L = 0
    t3 = sqrt((L - v*(t2)   -  v*v/2/a1) * 2 / a2);
    }
 
 
    //знаходимо загальний час:
    T = t1 + t2 + t3;
 
    cout << fixed << setprecision(8) << T;
    
 
 
     return 0;
}