• <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 閱讀(292) 評論(0)  編輯 收藏 引用 所屬分類: 幾何
            亚洲香蕉网久久综合影视| 亚洲国产精品久久久久网站| 无码久久精品国产亚洲Av影片| 免费精品久久久久久中文字幕| 色综合久久天天综合| 精品久久久久一区二区三区| 久久精品国产精品亚洲艾草网美妙| 国产精品99久久免费观看| 97久久超碰国产精品2021| 久久精品无码av| 精品人妻久久久久久888| 久久久精品久久久久特色影视| 久久精品国产男包| 久久夜色tv网站| 久久99久久99小草精品免视看 | 伊人久久综在合线亚洲2019| 久久青青草原国产精品免费| 久久久久久久波多野结衣高潮| 久久精品成人免费网站| 久久精品国产亚洲AV香蕉| 久久精品国产只有精品66| Xx性欧美肥妇精品久久久久久| 亚洲精品美女久久久久99| 欧美色综合久久久久久| 99久久国产主播综合精品| 久久狠狠高潮亚洲精品| 亚洲精品乱码久久久久久中文字幕| 精品一久久香蕉国产线看播放| 精品熟女少妇a∨免费久久| 无码AV波多野结衣久久| 国产精品久久久久国产A级| 精品午夜久久福利大片| 久久99国产精品久久99| 亚洲欧美久久久久9999| 亚洲精品无码专区久久久| 久久99国产精一区二区三区| 97久久精品午夜一区二区| 理论片午午伦夜理片久久 | 狠狠色综合久久久久尤物| 久久天天躁狠狠躁夜夜不卡| 91秦先生久久久久久久|