Tuesday, March 8, 2016

UVa 10341 Solution - Solve It

#include <bits/stdc++.h>
#define LL  long long
#define SZ(x) ((int)(x).size())
#define ALL(x) (x).begin(),(x).end()
#define REP(i,n) for(int i=0;i<n;i++)
#define REV(i,n) for(int i=n-1;i>=0;i--)
#define FOR(i,a,b) for(int i=a;i<=b;i++)
#define pri(a) cout<<a<<endl
#define prii(a,b) cout<<a<<" "<<b<<endl
#define priii(a,b,c) cout<<a<<" "<<b<<" "<<c<<endl
#define hlw printf("Hello World\n");
using namespace std;

const int INF = 1<<29;
const int MX  = 1e5+10;

double p,q,r,s,t,u;

double f(double x)
{
    return p * exp(-x) + q*sin(x) + r*cos(x) + s*tan(x) + t * x * x + u;
}

double bisection()
{
    double hi = 1, lo = 0;

    while(lo+1e-7 < hi)
    {
        double mid = (lo + hi) / 2;

        if(f(lo) * f(mid) <= 0)
        {
            hi = mid;
        }
        else
        {
            lo = mid;
        }
    }
    return (lo + hi) / 2;
}

int main()
{
    while(cin>>p>>q>>r>>s>>t>>u)
    {
        if(f(0) * f(1) > 0)
        {
            printf("No solution\n");
        }
        else
        {
            printf("%.4lf\n", bisection());
        }
    }

    return 0;
}

No comments:

Post a Comment