• <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   管理


            人妻无码精品久久亚瑟影视| 激情综合色综合久久综合| 亚洲欧美另类日本久久国产真实乱对白 | 久久99国产精品二区不卡| 久久久久这里只有精品| 欧美久久一区二区三区| 亚洲国产另类久久久精品| 亚洲国产成人久久精品影视| 久久人人爽人人人人片av| 天天躁日日躁狠狠久久| 国产精品伊人久久伊人电影| 久久精品中文闷骚内射| 99久久免费国产精品特黄| 久久这里只有精品首页| 99久久精品免费看国产一区二区三区| 2020久久精品国产免费| 奇米影视7777久久精品人人爽| 97久久精品人人澡人人爽| 久久精品无码一区二区无码| 亚洲伊人久久综合影院| 精品国产婷婷久久久| 成人免费网站久久久| 久久久久久久久久久| 久久这里有精品视频| 中文字幕一区二区三区久久网站| 久久国产欧美日韩精品| 亚洲国产天堂久久久久久| 久久久久亚洲精品天堂久久久久久 | 9久久9久久精品| 日韩乱码人妻无码中文字幕久久 | 国产毛片久久久久久国产毛片 | 久久久久亚洲精品天堂久久久久久| 99久久99久久精品免费看蜜桃| 久久久久久精品免费看SSS| 日韩十八禁一区二区久久| 久久人人爽人人爽人人片AV东京热| 久久九九全国免费| 国产99久久久久久免费看| 国产精品99久久久久久www| 伊人久久免费视频| 欧美一级久久久久久久大|