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

    
            

            




            







												


        
			
			


			
            
				
					
	
            
		
            
			pku 3072
		
            
		
            2009年7月26日

            題目鏈接：PKU 3072 Robot  

            分類：角度旋轉(zhuǎn)+Dijkastra

            題目分析與算法原型
                     此題也算是最短路的一個(gè)變形，計(jì)算代價(jià)的時(shí)候加入了每3個(gè)點(diǎn)，兩條邊的夾角的度數(shù)，因此需要自己好好的推算一下角度的公式，題目已經(jīng)提示了用atan2(double,double)這個(gè)函數(shù)，所以計(jì)算的時(shí)候把握好3個(gè)點(diǎn)中以哪個(gè)點(diǎn)為旋轉(zhuǎn)點(diǎn)，關(guān)于旋轉(zhuǎn)可以有兩個(gè)方向，取小于180度的那個(gè)方向        

            Code:

              1            #include<stdio.h>
              2#include<string.h>
              3#include<math.h>
              4#define max 1000000000
              5#define len 25
              6#define pi acos(-1.0)
              7
              8int r,n,visit[len],path[len];
              9double x[len],y[len],map[len][len],dis[len];
             10bool find;
             11
             12void init()
             13{
             14    int i,j;
             15    for(i=1;i<=n;i++)
             16        for(j=1;j<=n;j++)
             17        {
             18            if(i==j)map[i][j]=0;
             19            else
             20            {
             21                double l=sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));
             22                if(l>r)map[i][j]=max;
             23                else map[i][j]=l;
             24            }
             25        }
             26}
             27double degree(int start,int mid,int end,bool change)
             28{
             29    double ds,d1,d2;
             30
             31    if(!change)d1=atan2(y[mid]-y[start],x[mid]-x[start])*180/pi;
             32    else d1=atan2(y[start]-y[mid],x[start]-x[mid])*180/pi;
             33    d2=atan2(y[end]-y[mid],x[end]-x[mid])*180/pi;
             34    if(d1<0)d1+=360;
             35    if(d2<0)d2+=360;
             36    ds=fabs(d2-d1);
             37    if(ds>180)ds=360-ds;
             38    return ds;
             39}
             40void dij(int v0,int v)
             41{
             42    int i,j,u;
             43    for(i=1;i<n;i++)
             44    {
             45        if(i!=v0&&map[v0][i]<max)
             46        {
             47            double t=degree(n,v0,i,true);
             48            dis[i]=map[v0][i]+t;
             49        }
             50        else dis[i]=map[v0][i];
             51
             52        if(i!=v0&&dis[i]<max)path[i]=v0;
             53        else path[i]=-1;
             54    }
             55    dis[n]=map[v0][n];
             56    visit[v0]=1;
             57    for(i=1;i<n;i++)
             58    {
             59        double min=max;
             60        for(j=1;j<=n;j++)
             61            if(!visit[j]&&dis[j]<min)
             62            {
             63                u=j;
             64                min=dis[j];
             65            }
             66        if(min==max)return ;
             67
             68        if(u==v)
             69        {
             70            find=true;
             71            return ;
             72        }
             73        visit[u]=1;
             74        for(j=1;j<=n;j++)
             75            if(!visit[j]&&map[u][j]<max)
             76            {
             77                double tt=degree(path[u],u,j,false);
             78                if(dis[u]+map[u][j]+tt<dis[j])
             79                {
             80                    dis[j]=dis[u]+map[u][j]+tt;
             81                    path[j]=u;
             82                }
             83            }
             84    }
             85    return ;
             86}
             87
             88int main()
             89{
             90    int i;
             91    while(scanf("%d%d",&r,&n)!=EOF)
             92    {
             93        if(r==-1&&n==-1)break;
             94        for(i=1;i<=n;i++)scanf("%lf%lf",&x[i],&y[i]);
             95        find=false;
             96        init();
             97        memset(visit,0,sizeof(visit));
             98        dij(1,n);
             99        if(find)printf("%d\n",(int)(dis[n]+0.5));
            100        else printf("impossible\n");
            101    }
            102    return 0;
            103}
            104
            

		
            
			posted on 2009-07-26 17:49 蝸牛也Coding 閱讀(273) 評(píng)論(0)  編輯 收藏 引用  
		
            
	
            
	
	


	


            
	    
    




            



            


            

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

            


            

            		
		
            
	
            	

            



            

				
				

			
            
			
            
				
					
            
	
            
		
            <2009年7月>
            
	
            日一二三四五六
            2829301234
            567891011
            12131415161718
            19202122232425
            2627282930311
            2345678
            

            

					

            導(dǎo)航
            

					

            統(tǒng)計(jì)
            
	
					
					

            常用鏈接
            


            留言簿(8)
            


		
            隨筆檔案(78)
            
		
				
			
	



            搜索
            



            積分與排名
            



            最新評(píng)論
	
            

            
	
			
		

            


            閱讀排行榜
            



            評(píng)論排行榜
            


				
			
            
			
			
			
			
			 

	
    
    







            
	   
	



久久婷婷五月综合97色直播|
色综合久久88色综合天天|
污污内射久久一区二区欧美日韩|
国产午夜精品理论片久久|
久久久久亚洲爆乳少妇无|
久久伊人精品一区二区三区|
无码人妻少妇久久中文字幕蜜桃|
国产精品女同久久久久电影院|
久久久久综合网久久|
久久久这里有精品中文字幕|
久久天天躁夜夜躁狠狠|
精品久久人妻av中文字幕|
伊人久久大香线焦综合四虎|
久久亚洲中文字幕精品一区|
亚洲日本va中文字幕久久|
久久精品国产亚洲网站|
香蕉99久久国产综合精品宅男自
|
久久免费的精品国产V∧|
青青草原综合久久大伊人精品|
天堂无码久久综合东京热|
久久亚洲精品成人AV|
国产日韩欧美久久|
日韩乱码人妻无码中文字幕久久|
久久精品国产91久久综合麻豆自制|
久久精品国产一区二区电影|
奇米影视7777久久精品|
久久久99精品成人片中文字幕|
久久亚洲中文字幕精品有坂深雪
|
精品久久久久久99人妻|
色综合久久综合中文综合网|
国产免费久久精品99久久|
三上悠亚久久精品|
日韩中文久久|
国产精品欧美久久久天天影视|
亚洲国产精品综合久久网络|
国产亚洲精品美女久久久|
久久国产亚洲精品|
久久精品亚洲男人的天堂|
.精品久久久麻豆国产精品|
超级碰碰碰碰97久久久久|
91精品婷婷国产综合久久|