• <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++博客 ::
首頁 ::
新隨筆 ::
聯(lián)系 ::
聚合
 ::
管理 ::


			
			
				

            
	
            
		
            
			
			
			
				 
				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)  編輯 收藏 引用  所屬分類: 幾何 
	
            

            







            
	    
    




            



            


            

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

            


            

            		
		
            
	
            	

            



            

				

            



    
    







            
	   
	



久久精品国产亚洲沈樵|
久久精品中文字幕一区|
国产精品久久99|
久久国产影院|
久久久久成人精品无码中文字幕|
久久亚洲欧美日本精品|
久久精品国产精品亚洲下载|
97精品依人久久久大香线蕉97|
久久国产亚洲高清观看|
无码人妻久久一区二区三区蜜桃|
久久亚洲国产精品成人AV秋霞|
狠狠88综合久久久久综合网|
久久se精品一区二区影院|
久久久久久九九99精品|
久久综合视频网站|
久久99精品综合国产首页|
久久久久久伊人高潮影院|
麻豆国内精品久久久久久|
久久久久中文字幕|
久久香蕉超碰97国产精品
|
欧美伊人久久大香线蕉综合|
97精品国产91久久久久久|
久久九九兔免费精品6|
国产成人久久精品麻豆一区|
人妻精品久久无码区|
亚洲人成无码网站久久99热国产
|
久久精品人妻一区二区三区|
狠狠色丁香久久婷婷综合五月|
国产精品99久久久久久宅男小说|
久久狠狠一本精品综合网|
2022年国产精品久久久久|
欧洲成人午夜精品无码区久久
|
久久精品国产精品亚洲下载|
久久伊人精品青青草原高清|
国产一级持黄大片99久久|
国产成人久久精品区一区二区|
久久久久久久精品妇女99|
日韩十八禁一区二区久久|
久久国内免费视频|
精品国产乱码久久久久久呢|
国产成人精品综合久久久久|