hdu 2899

mxy终于学会求函数极值了。

先写一道板子。

#include <bits/stdc++.h>
using namespace std;
typedef double db;
const db eps = 1e-6;
const db INF = 1e100;
const db delta = 0.98;
int sign(db k){
    if (k>eps) return 1; else if (k<-eps) return -1; return 0;
}
int cmp(db k1,db k2){return sign(k1-k2);}
int t;db y;
db get(db x){
    return 6*x*x*x*x*x*x*x+8*x*x*x*x*x*x+7*x*x*x+5*x*x-y*x;
}
db x[10];
db random(){
    return (rand()&1?1:-1)*rand()*1.0/32767;
}
void init(){
    for(int i=0;i<10;i++) x[i]=abs(random())*100;
}
db sa(){
    db ans = INF;
    db t=100;
    while (t>eps){
        for(int i=0;i<10;i++){
            db tmp = get(x[i]);
            for(int j=0;j<10;j++){
                db _x=x[i]+random()*t;
                if(cmp(_x,0)>=0&&cmp(_x,100)<=0){
                    db f = get(_x);
                    if(tmp>f)
                        x[i]=_x;
                }
            }
        }
        t*=delta;
    }
    for(int i=0;i<10;i++)
        ans=min(ans,get(x[i]));
    return ans;
}

int main(){
    ios::sync_with_stdio(false);
    scanf("%d",&t);
    while (t--){
        scanf("%lf",&y);
        init();
        printf("%.4f
",sa());
    }
}