• <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>
            隨筆 - 87  文章 - 279  trackbacks - 0
            <2007年10月>
            30123456
            78910111213
            14151617181920
            21222324252627
            28293031123
            45678910

            潛心看書研究!

            常用鏈接

            留言簿(19)

            隨筆分類(81)

            文章分類(89)

            相冊

            ACM OJ

            My friends

            搜索

            •  

            積分與排名

            • 積分 - 218001
            • 排名 - 117

            最新評論

            閱讀排行榜

            評論排行榜

            Sailboat

            Problem H: Sailboat

            In the sailboat race, the contestant is requested to along with the prearrange path. Sailing ship's power comes from wind power and contestant's manpower. The wind power can completely used.

            In a competition, the contestants are requested to along with a 1/4 circles with radius R, the sailboat will goto east from south. During this process, the wind direction is straight from west to the east with constant speed and power.

            In order to maintain the travel direction, the athlete must adjust the sail to the vertical angle from movement direction in any time.

            If the speed of sailboat is proportional to the power at movement direction, the proportional factor is k. Supposes the wind power is f, the athlete manpower is h, please given the time of sailboat from the beginning to the end.

            Input

            The first line of each case consists of 4 double number, that is radius of path: R, wind power: f,athlete manpower: h and proportional factor:k. In order to avoid the floating point error, you needn't output the answer directly. The next line is a integer n, the following n lines gives a double value which is candidate answer.

            Output

            For each candidate of each case, Only "Yes" or "No" should be printed. Output "Yes" if the relative error to your answer is less than 3%, otherwise "No". For example, if the model answer is 100, and the candidate is 98 or 102, you should output "Yes". Output one blank line between neighboring case

            Sample Input

            1.0 2.0 1.0 1.0
            2
            0.35
            0.76
            

            Sample Output

            No
            Yes
            

            Problem Source: provided by skywind

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

            const int MAXN = 100;
            const double PI = acos(-1.0);

            double R, F, H, K, ans;
            int n, cas;

            double func(double x) {
                
            return R / K / (H + F * cos(x));
            }


            double romberg(double a, double b, double EPS = 1e-6{
                
            double t[MAXN][MAXN] = {0}, tmp;
                
            int i, j, k, k2, m, m4;
                t[
            0][0= (func(a) + func(b)) * (b - a) / 2;
                k 
            = 1; k2 = 1;
                
            while (1{
                    tmp  
            = 0;
                    
            for (i = 1; i <= k2; i++{
                        tmp 
            += func(a + (2 * i - 1* (b - a) / (2 * k2));
                    }

                    t[
            0][k] = (t[0][k - 1+ tmp * (b - a) / k2) / 2;
                    
            for (m = 1, m4=4; m <= k; m++, m4 *= 4{
                        t[m][k 
            - m] = (m4 * t[m - 1][k - m + 1- t[m - 1][k - m]) / (m4 - 1);
                    }

                    
            if (fabs(t[k][0- t[k - 1][0]) < EPS) break;
                    k
            ++; k2 *= 2;
                }

                
            return t[k][0];
            }


            void solve() {
                
            double tmp;
                scanf(
            "%lf"&tmp);
                
            if (fabs(tmp - ans) / ans < 0.03) printf("Yes\n");
                
            else printf("No\n");
            }


            int main() {
                freopen(
            "2457.in""r", stdin);
                
            while (scanf("%lf%lf%lf%lf%d"&R, &F, &H, &K, &n) != EOF) {
                    
            if (cas) printf("\n");
                    
            else cas++;
                    ans 
            = romberg(0, PI/2);
                    
            while (n--{
                        solve();
                    }

                }

                
            return 0;
            }
            posted on 2007-10-20 01:02 閱讀(692) 評論(0)  編輯 收藏 引用 所屬分類: 算法&ACM
            亚洲国产精品综合久久网络 | 成人精品一区二区久久| 99国产欧美久久久精品蜜芽| 精品久久久久久国产91| 国产精品女同一区二区久久| 久久亚洲精品无码VA大香大香| 久久亚洲欧美国产精品| 久久黄色视频| 久久国产色AV免费看| 久久久久综合国产欧美一区二区| 精品久久亚洲中文无码| 国产精品VIDEOSSEX久久发布| 国产aⅴ激情无码久久| 久久高清一级毛片| 国产精品美女久久久久久2018| 久久国产精品无| 国产日韩久久久精品影院首页| 久久国产亚洲高清观看| 国内精品九九久久精品| 久久久久18| 国产高清美女一级a毛片久久w| 97r久久精品国产99国产精| 思思久久99热只有频精品66| 久久久99精品成人片中文字幕| 久久香蕉一级毛片| 国产精品无码久久久久久| 浪潮AV色综合久久天堂| 中文国产成人精品久久不卡| 国产午夜电影久久| 国产精自产拍久久久久久蜜| 国产L精品国产亚洲区久久 | 99热都是精品久久久久久| 久久久久亚洲AV成人片| 日产精品久久久一区二区| 亚洲精品第一综合99久久 | 国产亚洲精品美女久久久| 色综合久久久久无码专区| 中文字幕人妻色偷偷久久| 狠狠综合久久综合88亚洲| 久久人人爽人人爽人人片av高请| 无码人妻久久一区二区三区|