码迷,mamicode.com
首页 > 其他好文 > 详细

hdu 2899 hdu 3400 三分/几何

时间:2014-10-09 01:59:18      阅读:217      评论:0      收藏:0      [点我收藏+]

标签:blog   io   os   sp   2014   c   on   log   r   

hdu2899 ;

水提,直接三分,其实求导后二分也可以。

#include<iostream>
#include<cstdio>
using namespace std;
  double y;
  double inline f( long double 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;
}
int main()
{
    int T;scanf("%d",&T);
    while(T--)
    {
        cin>>y;
        double l,r,m1,m2;
        l=0;r=100;
        while(r-l>0.000001)
        {
            m1=l+(r-l)/3;
            m2=r-(r-l)/3;
            if(f(m1)<f(m2))
            {
                r=m2;
            }
            else
            {
                l=m1;
            }
        }
       printf("%.4lf\n",f(l));
    }
    return 0;
}


hdu3400  有俩个未知,第一次确定,则第二个函数式先见后增的,而第二个反馈出最小值后,第一个未知也是先见后增的,所以三分套三分(九分??呵呵~)

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
struct xy
{
    double x,y;
};
double inline dis( xy a,xy b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
xy A,B,C,D;    int p,q,r;
double inline s(xy a ,xy b)
{
    return dis(a,A)/p+dis(a,b)/r+dis(b,D)/q;
}
double f(xy t)
{
     xy l,r,m1,m2;
     l=C;r=D;
        while(dis(r,l)>0.0001)
        {
            m1.x=l.x+(r.x-l.x)/3;
            m1.y=l.y+(r.y-l.y)/3;
            m2.x=r.x-(r.x-l.x)/3;
            m2.y=r.y-(r.y-l.y)/3;
            if(s(t,m1)<s(t,m2))
            {
                r=m2;
            }
            else
            {
                l=m1;
            }
        }
        return s(t,l);
}

int main()
{
    int T;scanf("%d",&T);
    while(T--)
    {
        scanf("%lf%lf%lf%lf%lf%lf%lf%lf",&A.x,&A.y,&B.x,&B.y,&C.x,&C.y,&D.x,&D.y);
        scanf("%d%d%d",&p,&q,&r);
        xy l,r,m1,m2;
        l=A;r=B;
        while(dis(r,l)>0.0001)
        {
            m1.x=l.x+(r.x-l.x)/3;
            m1.y=l.y+(r.y-l.y)/3;
            m2.x=r.x-(r.x-l.x)/3;
            m2.y=r.y-(r.y-l.y)/3;
            if(f(m1)<f(m2))
            {
                r=m2;
            }
            else
            {
                l=m1;
            }
        }
       printf("%.2lf\n",f(l));
    }
    return 0;
}





hdu 2899 hdu 3400 三分/几何

标签:blog   io   os   sp   2014   c   on   log   r   

原文地址:http://blog.csdn.net/u011498819/article/details/39899503

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!