hdu1245+dij,堆优化
有一片100*100的湖泊，中心坐标（0,0）,即湖泊右上角的坐标是（50,50），湖泊中间有一片以（0，0）为圆心，15为直径的圆形陆地。现有一个人在陆地，湖泊中散布着一些点可以踩，这个人要利用这些点跳到岸上，求最短路径和最短路径下的最短步数。
spfa莫名其妙的超时，dij+堆优化反而能过。。。可能spfa更适合有向图吧。
使用vector<vector<node> > g这种双重结构时，最好先g.resize(N)设置一下容量，否则直接插入会出错。
1 #include<cstdio> 2 #include<cstring> 3 #include<vector> 4 #include<cmath> 5 #include<cstdlib> 6 #include<queue> 7 using namespace std; 8 const double inf=1<<30; 9 int n; 10 double limit; 11 struct node 12 { 13 double w; 14 int v; 15 int path; 16 node(double ww,int vv,int ppath):w(ww),v(vv),path(ppath){} 17 node(){} 18 bool operator <(const node &dd) const 19 { 20 if(w!=dd.w) 21 return w>dd.w; 22 else 23 return path>dd.path; 24 } 25 }; 26 vector<vector<node> > g; 27 double dis(double x1,double y1,double x2,double y2) 28 { 29 return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)); 30 } 31 double mmax(double aa,double bb) 32 { 33 if(aa>bb) return aa; 34 else return bb; 35 } 36 int vis[110]; 37 38 int main() 39 { 40 double x[110],y[110]; 41 int i,j; 42 while(scanf("%d%lf",&n,&limit)!=EOF) 43 { 44 priority_queue<node> q; 45 g.clear(); 46 g.resize(n+10); 47 for(i=1;i<=n;i++) 48 scanf("%lf%lf",&x[i],&y[i]); 49 for(i=1;i<=n;i++) 50 { 51 double dd=dis(x[i],y[i],0.0,0.0)-7.5; 52 if(dd<=limit) 53 { 54 g[0].push_back(node(dd,i,0)); 55 g[i].push_back(node(dd,0,0)); 56 } 57 } 58 for(i=1;i<=n;i++) 59 { 60 for(j=i+1;j<=n;j++) 61 { 62 double dd=dis(x[i],y[i],x[j],y[j]); 63 if(dd<=limit) 64 { 65 g[i].push_back(node(dd,j,0)); 66 g[j].push_back(node(dd,i,0)); 67 } 68 } 69 } 70 double mm; 71 for(i=1;i<=n;i++) 72 { 73 mm=mmax(fabs(x[i]),fabs(y[i])); 74 if(50-mm<=limit) 75 { 76 g[i].push_back(node(50-mm,n+1,0)); 77 g[n+1].push_back(node(50-mm,i,0)); 78 } 79 } 80 //前面都是建图 81 q.push(node(0,0,0));//放入起点 82 memset(vis,0,sizeof(vis));//标记这个点的最短路是否已经求出 83 node z; 84 bool flag=false;//判断有无解 85 while(!q.empty()) 86 { 87 z=q.top(); 88 q.pop(); 89 if(vis[z.v]) continue;//这个点最短路已经求出，不需要继续求了 90 vis[z.v]=1;//标记这个点的最短路已求出 91 if(z.v==n+1) {flag=true;break;}//因为只需要求起点到n+1的最短路，所以可以结束了 92 for(i=0;i<g[z.v].size();i++) 93 { 94 node zz; 95 zz.v=g[z.v][i].v; 96 if(vis[zz.v]) continue; 97 zz.w=z.w+g[z.v][i].w; 98 zz.path=z.path+1; 99 q.push(zz);//反正是优先队列，不用管队列中是否已经存在zz，更新最短路的时候也不用管是否能使当前的最短路变得更短 100 } 101 } 102 if(flag) 103 printf("%.2f %d\n",z.w,z.path); 104 else printf("can't be saved\n"); 105 } 106 return 0; 107 }
