2020-2021 ACM-ICPC, Asia Nanjing Regional Contes F-Fireworks

题目：

思路分析：

我们计算的是最优期望

而且我们理想的是一次性制作好 然后一次性点燃

那就 就是k-1都失败 第k个成功

符合几何分布

概率期望是f(x)=(n*k+m)/(1-p**k)

证明这是一个单峰函数 直接三分找最大点即可

代码实现：

/*
*@Author:   GuoJinlong
*@Language: C++
*/
//#include <bits/stdc++.h>
/*
 *                                                     __----~~~~~~~~~~~------___
 *                                    .  .   ~~//====......          __--~ ~~
 *                    -.            \_|//     |||\\  ~~~~~~::::... /~
 *                 ___-==_       _-~o~  \/    |||  \\            _/~~-
 *         __---~~~.==~||\=_    -_--~/_-~|-   |\\   \\        _/~
 *     _-~~     .=~    |  \\-_    '-~7  /-   /  ||    \      /
 *   .~       .~       |   \\ -_    /  /-   /   ||      \   /
 *  /  ____  /         |     \\ ~-_/  /|- _/   .||       \ /
 *  |~~    ~~|--~~~~--_ \     ~==-/   | \~--===~~        .\
 *           '         ~-|      /|    |-~\~~       __--~~
 *                       |-~~-_/ |    |   ~\_   _-~            /\
 *                            /  \     \__   \/~                \__
 *                        _--~ _/ | .-~~____--~-/                  ~~==.
 *                       ((->/~   '.|||' -_|    ~~-/ ,              . _||
 *                                  -_     ~\      ~~---l__i__i__i--~~_/
 *                                  _-~-__   ~)  \--______________--~~
 *                                //.-~~~-~_--~- |-------~~~~~~~~
 *                                       //.-~~~--\
 *                       ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 *
 *                               神兽保佑            永无BUG
 */



const int MAX=1010;
int n,m;
double p;
double f(double x){
    double up=n*x+m;
    double down=1.0-pow(p,x);
    return up/down;
}
int main(){
    int t;
    cin>>t;
    while (t--) {
        cin>>n>>m;
        cin>>p;
        p=1.0-p/10000.0;
        double l=0;
        double r=1e9+10;
        while (r-l>EPS) {
            double ll=l+(r-l)/3.0;
            double rr=r-(r-l)/3.0;
            if(f(ll)>f(rr)){
                l=ll;
            }
            else r=rr;
        }
        int ans=l;
        printf("%.10f\n",min(f(ans),f(ans+1)));
    }
}

/**
*　　┏┓　　　┏┓+ +
*　┏┛┻━━━┛┻┓ + +
*　┃　　　　　　　┃
*　┃　　　━　　　┃ ++ + + +
*  ████━████+
*  ◥██◤　◥██◤ +
*　┃　　　┻　　　┃
*　┃　　　　　　　┃ + +
*　┗━┓　　　┏━┛
*　　　┃　　　┃ + + + +Code is far away from
*　　　┃　　　┃ + bug with the animal protecting
*　　　┃　 　 ┗━━━┓ 神兽保佑,代码无bug
*　　　┃ 　　　　　　 ┣┓
*　　  ┃ 　　　　　 　┏┛
*　    ┗┓┓┏━┳┓┏┛ + + + +
*　　　　┃┫┫　┃┫┫
*　　　　┗┻┛　┗┻┛+ + + +
*/