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

            gzwzm06

              C++博客 :: 首頁 :: 新隨筆 :: 聯系 :: 聚合  :: 管理 ::
              1 隨筆 :: 52 文章 :: 17 評論 :: 0 Trackbacks
              1#include <cstdio>
              2#include <cmath>
              3
              4const int SIZE = 102;
              5
              6struct TPoint
              7{
              8    double x ;
              9    double y ;
             10}
             ;
             11
             12TPoint point[SIZE];
             13int N;
             14
             15double Distance( const TPoint& p1, const TPoint& p2 )
             16{
             17    return sqrt((p1.x - p2.x) * (p1.x - p2.x) + (p1.y - p2.y) * (p1.y - p2.y)) ;
             18}

             19
             20double GetALLDist( const TPoint& p, const int& n )
             21{
             22    int i;
             23    double sum = 0.0;
             24
             25    for ( i = 0; i < n; ++i )
             26        sum += Distance( p, point[i] );
             27
             28    return sum;
             29}

             30
             31double GetFermatPoint( const int& n )
             32{
             33    TPoint st = point[0];
             34    TPoint tmp, nt;
             35    double t, step = 100, fermat = GetALLDist(st, n);
             36
             37    while ( step > 0.2 )
             38    {
             39        int ok = 1;
             40
             41        while ( ok ) {
             42            ok = 0;
             43            nt = st;
             44
             45            tmp.x = st.x, tmp.y = st.y + step;
             46
             47            t = GetALLDist(tmp, n);
             48
             49            if ( t < fermat )
             50            {
             51                fermat = t;
             52                ok = 1;
             53                nt = tmp;
             54            }

             55
             56            tmp.x = st.x, tmp.y = st.y - step;
             57
             58            t = GetALLDist(tmp, n);
             59            if ( t < fermat )
             60            {
             61                fermat = t;
             62                ok = 1;
             63                nt = tmp;
             64            }

             65
             66            tmp.x = st.x + step, tmp.y = st.y;
             67
             68            t = GetALLDist(tmp, n);
             69
             70            if ( t < fermat )
             71            {
             72                fermat = t;
             73                ok = 1;
             74                nt = tmp;
             75            }

             76
             77            tmp.x = st.x - step, tmp.y = st.y;
             78
             79            t = GetALLDist(tmp, n);
             80            if ( t < fermat )
             81            {
             82                fermat = t;
             83                ok = 1;
             84                nt = tmp;
             85            }

             86            
             87            st = nt;
             88        }

             89
             90        step /= 2.0;
             91    }

             92
             93    return fermat;
             94}

             95
             96int main()
             97{
             98    int i, ans;
             99    double t;
            100
            101    scanf("%d"&N);
            102
            103    for ( i = 0; i < N; ++i )
            104    {
            105        scanf("%lf %lf"&point[i].x, &point[i].y);
            106    }

            107
            108    t = GetFermatPoint( N );
            109
            110    ans = (int)(t + 0.5* 100 / 100;
            111
            112    printf("%d\n", ans);
            113
            114    return 0;
            115}
            posted on 2009-03-31 12:08 閱讀(286) 評論(0)  編輯 收藏 引用 所屬分類: 幾何
            国产亚洲美女精品久久久2020| 青青青国产成人久久111网站| 久久综合色区| 亚洲AV无码久久精品成人| 国产精品国色综合久久| 国内精品久久久久久久久| 综合久久精品色| 久久夜色tv网站| 久久久久一级精品亚洲国产成人综合AV区 | 理论片午午伦夜理片久久| 一本色道久久88精品综合| 久久精品成人| 97热久久免费频精品99| 欧美激情一区二区久久久| 欧美精品一区二区精品久久| 欧美成人免费观看久久| 久久国产亚洲精品麻豆| 伊人久久大香线蕉av一区| 久久久亚洲精品蜜桃臀 | 亚洲成人精品久久| 色婷婷久久综合中文久久蜜桃av| 久久久精品久久久久久| 久久96国产精品久久久| 色婷婷综合久久久久中文一区二区 | 久久久久久久久久久| 99久久精品国产一区二区三区| 亚洲va中文字幕无码久久| 亚洲人成无码久久电影网站| 国产精品激情综合久久| 久久99国产精品二区不卡| 91久久精品91久久性色| 久久久免费精品re6| 伊人色综合久久天天人手人婷| 久久毛片免费看一区二区三区| 91麻精品国产91久久久久| 青青青青久久精品国产| 91精品国产91久久久久久| 91秦先生久久久久久久| 热久久这里只有精品| 999久久久免费国产精品播放| 国产高潮国产高潮久久久91 |