• <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 閱讀(313) 評論(0)  編輯 收藏 引用


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


            久久精品国产亚洲AV蜜臀色欲| 麻豆成人久久精品二区三区免费 | 久久久无码精品亚洲日韩软件| 亚洲国产天堂久久综合网站| 久久AAAA片一区二区| 亚洲国产精品高清久久久| 久久99亚洲网美利坚合众国| 久久精品国产亚洲Aⅴ香蕉 | 久久人人爽人人爽人人片AV东京热 | 天天综合久久久网| 色综合久久夜色精品国产| 精品国产一区二区三区久久| 日本亚洲色大成网站WWW久久| 国产美女久久精品香蕉69| 亚洲v国产v天堂a无码久久| 久久国产精品久久久| 久久精品国产欧美日韩99热| 丰满少妇人妻久久久久久4| 久久精品无码午夜福利理论片| 久久久久女教师免费一区| 97久久精品人妻人人搡人人玩| 99久久香蕉国产线看观香| 久久久久香蕉视频| 久久精品国产亚洲av日韩| 久久精品国产亚洲AV蜜臀色欲| 亚洲精品乱码久久久久久蜜桃| 国产精品欧美久久久久无广告| 久久不射电影网| 国产精品对白刺激久久久| 久久精品国产亚洲精品2020| 中文字幕久久精品无码| 久久久久久久波多野结衣高潮 | 久久性精品| 久久激情亚洲精品无码?V| 国产精品一区二区久久精品无码 | 亚洲国产成人久久笫一页| 国产成人久久精品麻豆一区| 18岁日韩内射颜射午夜久久成人| 久久九九亚洲精品| 久久综合久久久| 久久久久国产一级毛片高清板|