• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>

            The Fourth Dimension Space

            枯葉北風(fēng)寒,忽然年以殘,念往昔,語默心酸。二十光陰無一物,韶光賤,寐難安; 不畏形影單,道途阻且慢,哪曲折,如渡飛湍。斬浪劈波酬壯志,同把酒,共言歡! -如夢令

            HDOJ 3400 三分法

            1Y。只是不明白為什么可以套兩個三分,不過從實(shí)際情況來看,在第一條直線上三分應(yīng)該也是符合凸函數(shù)性質(zhì)的。

            #include<iostream>
            #include
            <cmath>
            #include
            <algorithm>
            using namespace std;

            #define eps 1e-8


            struct point
            {
                
            double x,y;
            }
            ;
            point a,b,c,d;
            double p,q,r;

            double dist(point a,point b)
            {
                
            return sqrt( (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y) );
            }


            double cal(double t1,double t2)//以t1,t2為參數(shù)算出時間
            {
                point tm1;
                point tm2;
                tm1.x
            =a.x+(b.x-a.x)*t1;
                tm1.y
            =a.y+(b.y-a.y)*t1;

                tm2.x
            =c.x+(d.x-c.x)*t2;
                tm2.y
            =c.y+(d.y-c.y)*t2;


                
            return dist(tm1,a)/p+dist(tm2,d)/q+dist(tm1,tm2)/r;
            }



            double sanfen(double t1)//在確定t1的基礎(chǔ)上得最小值
            {
                
            double l=0;
                
            double r=1;
                
            while(l+eps<=r)
                
            {

                    
            double mid=(l+r)/2;
                    
            double mmid=(mid+r)/2;
                    
            if(cal(t1,mid)<cal(t1,mmid))
                        r
            =mmid;
                    
            else
                        l
            =mid;
                }

                
            return cal(t1,l);
            }




            int main()
            {
                
            int ca;
                scanf(
            "%d",&ca);
                
            while(ca--)
                
            {
                    scanf(
            "%lf%lf%lf%lf",&a.x,&a.y,&b.x,&b.y);
                    scanf(
            "%lf%lf%lf%lf",&c.x,&c.y,&d.x,&d.y);
                    
            //
                    scanf("%lf%lf%lf",&p,&q,&r);

                    
            double l,r;
                    l
            =0;r=1;
                    
            while(l+eps<=r)
                    
            {

                        
            double mid=(l+r)/2;
                        
            double mmid=(mid+r)/2;
                        
            if(sanfen(mid)<sanfen(mmid))
                            r
            =mmid;
                        
            else
                            l
            =mid;
                    }

                    printf(
            "%.2lf\n",sanfen(l));
                }


                
            return 0;
            }

            posted on 2010-11-07 20:57 abilitytao 閱讀(316) 評論(0)  編輯 收藏 引用


            只有注冊用戶登錄后才能發(fā)表評論。
            網(wǎng)站導(dǎo)航: 博客園   IT新聞   BlogJava   博問   Chat2DB   管理


            成人久久久观看免费毛片| 99久久www免费人成精品| 99精品国产在热久久无毒不卡| 久久精品一区二区三区AV| 亚洲午夜久久久久久久久久| 久久综合综合久久狠狠狠97色88| 国产激情久久久久影院| 亚洲精品tv久久久久久久久久| 久久午夜羞羞影院免费观看| 99久久亚洲综合精品网站| 精品久久久久中文字幕一区| 嫩草伊人久久精品少妇AV| 精品久久久久久中文字幕人妻最新| 精品无码人妻久久久久久| 波多野结衣久久精品| 伊人久久精品影院| 久久丫精品国产亚洲av不卡| 狠狠精品干练久久久无码中文字幕 | 香蕉99久久国产综合精品宅男自 | 国产成人综合久久精品尤物| 久久久久99这里有精品10| 久久国产精品波多野结衣AV| 久久综合久久美利坚合众国 | 亚洲国产精品成人AV无码久久综合影院 | 精品久久久久久国产| 国产高潮国产高潮久久久91| 日韩精品久久久久久久电影蜜臀| www亚洲欲色成人久久精品| 久久亚洲精品中文字幕| 亚洲va久久久久| 久久精品亚洲男人的天堂| 伊人热人久久中文字幕| 亚洲av成人无码久久精品| 久久精品国产亚洲av麻豆小说 | 精品久久久久久久久久中文字幕 | 狠狠久久综合伊人不卡| 一本久久a久久精品亚洲| 欧美激情精品久久久久久久| 国产精品久久毛片完整版| 精品国产乱码久久久久软件| 中文精品99久久国产 |