CometOJ-Contest#2B她的想法、他的战斗（概率＋数学）

输入描述

第一行两个正整数 $(1\leqslant l 1⩽l<r⩽2000)。$

第二行两个正整数 $(1\leqslant L 1⩽L<R⩽2000)。$

输出描述

一个实数，表示最大的期望利润。(四舍五入后保留

如果答案为

若答案为 $v+10^{-6}v+10−6 以及 max(v-10^{-6},0)max(v−10−6,0) 四舍五入后保留 44 位小数的结果不会改变。$

样例输入 1

400 1200
600 1800

样例输出 1

200.0000

样例输入 2

1999 2000
1 2

样例输出 2

0.0000


思路：
因为主人公的售价是　Ｌ～Ｒ　所以　主人公售价的期望是　（Ｌ＋Ｒ）／２
我们给买手的钱是ｐ，那么我们主人公的利润就是　（（Ｌ＋Ｒ）／２－ｐ）
我们知道　期望＝值＊概率
利润期望＝利润值＊概率

我们来看下概率，
Ｐ必须取到　ｌ和ｒ之间，因为如果ｐ大于ｒ，那么不是最优的选择，因为可以选择ｒ可以产生更好的利益。
而ｐ小于ｌ的话，买手会无法购物商品，所以也不可以。
那么ｐ就在　ｌ和ｒ之间了，
买手能买商品的时候当且仅当　商品的市场价　恰好小于等于ｐ，也就是市场价ｙ小于ｐ时，可以购买，
那么ｙ的范围时　ｌ～ｐ，总范围是ｌ～ｒ，所以概率是　（ｐ－ｌ）／（ｒ－ｌ）
期望就是

（（Ｌ＋Ｒ）／２－ｐ）＊（ｐ－ｌ）／（ｒ－ｌ）

公式就是一个关于ｐ的二次函数，并且ａ是小于等于０的，有极大值，我们判定对称轴是否在ｌ到ｒ之间，
如果在，那么函数的极值就是答案，如果不在，那么取区间的左右极限的较大值是答案。
注意处理以下０＊负数＝－０的情况，

细节见ｃｏｄｅ

#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include <set>
#include 
#include 
#define ALL(x) (x).begin(), (x).end()
#define rt return
#define dll(x) scanf("%I64d",&x)
#define xll(x) printf("%I64d\n",x)
#define sz(a) int(a.size())
#define all(a) a.begin(), a.end()
#define rep(i,x,n) for(int i=x;i#define repd(i,x,n) for(int i=x;i<=n;i++)
#define pii pair
#define pll pair
#define gbtb ios::sync_with_stdio(false),cin.tie(0),cout.tie(0)
#define MS0(X) memset((X), 0, sizeof((X)))
#define MSC0(X) memset((X), '\0', sizeof((X)))
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define eps 1e-6
#define gg(x) getInt(&x)
#define db(x) cout<<"== [ "<using namespace std;
typedef long long ll;
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
ll lcm(ll a,ll b){return a/gcd(a,b)*b;}
ll powmod(ll a,ll b,ll MOD){ll ans=1;while(b){if(b%2)ans=ans*a%MOD;a=a*a%MOD;b/=2;}return ans;}
inline void getInt(int* p);
const int maxn=1000010;
const int inf=0x3f3f3f3f;
/*** TEMPLATE CODE * * STARTS HERE ***/
double l,r;
double L,R;
int main()
{
    //freopen("D:\\common_text\\code_stream\\in.txt","r",stdin);
    //freopen("D:\\common_text\\code_stream\\out.txt","w",stdout);
    
    gbtb;
    cin>>l>>r;
    cin>>L>>R;
    double t=(L+R)/2.00;
    double k=(r-l);
    double a=-1.00/k;
    double b=(t+l)/k;
    double c=-l*t/k;
    double x=-b/(2.00*a);
    double ans;
    if(x>=l&&x<=r)
    {
        ans=(4*a*c-b*b)/(4*a);
        // db(a*x*x+b*x+c);
    }else
    {
        ans=a*l*l+b*l+c;
        ans=max(ans,a*r*r+b*r+c);
    }
    ans+=eps;
    ans=max(0.0,ans-eps);
    cout<<fixed<4)<endl;
    
    return 0;
}

inline void getInt(int* p) {
    char ch;
    do {
        ch = getchar();
    } while (ch == ' ' || ch == '\n');
    if (ch == '-') {
        *p = -(getchar() - '0');
        while ((ch = getchar()) >= '0' && ch <= '9') {
            *p = *p * 10 - ch + '0';
        }
    }
    else {
        *p = ch - '0';
        while ((ch = getchar()) >= '0' && ch <= '9') {
            *p = *p * 10 + ch - '0';
        }
    }
}