zoukankan      html  css  js  c++  java
  • 传说中的次短路

    我不知道我是不是对了,没数据真是忧伤,我是传统的A*+SPFA算法,

    原题是求次短路,我借此学了一下这种神奇的东西。

      1 #include <iostream>
      2 #include <fstream>
      3 #include <cstdlib>
      4 #include <cstring>
      5 /* run this program using the console pauser or add your own getch, system("pause") or input loop */
      6 using namespace std;
      7 ifstream fin("secondpath.in");
      8 ofstream fout("secondpath.out");
      9 
     10 int places=0,rodes=0;
     11 int qi=0,zhong=0;
     12 int head[1005],cnt=0;
     13 struct id{
     14  int nxt;
     15  int to;
     16  int quan;    
     17 };
     18 id lian[200005];
     19 int dui[1003],jv[1003];
     20 int jv2[1003];
     21 bool used[1003];
     22 int you_dui[200005][2],gs_cnt=0;
     23 //int zhan[200003];
     24 
     25 
     26 void add(int co,int dao,int bian);
     27 void spfa( );
     28 void a_xin( );
     29 int na(int sze);
     30 void fang(int sze,int zhi); 
     31 
     32 
     33 void add(int co,int dao,int bian){
     34 lian[++cnt].nxt=head[co];    
     35 lian[cnt].to=dao;
     36 lian[cnt].quan=bian;
     37 head[co]=cnt;
     38 return;
     39 }
     40 
     41 
     42 void spfa( ){
     43 memset(jv,127/2,sizeof(jv));
     44 int tou=0,wei=1;
     45 dui[1]=zhong;jv[zhong]=0;    
     46 used[zhong]=1;
     47 while(tou!=wei){
     48  tou++;if(tou>1001)tou=1;
     49  int cong=dui[tou];
     50  for(int x=head[cong];x!=-1;x=lian[x].nxt){
     51   int dao=lian[x].to;
     52   if(jv[dao]>jv[cong]+lian[x].quan){
     53    jv[dao]=jv[cong]+lian[x].quan;    
     54    if(used[dao]==0){
     55     wei++;if(wei>1001)wei=1;
     56     dui[wei]=dao;used[dao]=1;
     57    }    
     58   }        
     59  }    
     60  used[cong]=0;    
     61 }
     62 return;
     63 }
     64 
     65 
     66 void a_xin( ){
     67 memset(you_dui,127/2,sizeof(you_dui));
     68 memset(jv2,127/2,sizeof(jv2));
     69 jv2[qi]=0;you_dui[1][0]=0;    
     70 you_dui[1][1]=qi;gs_cnt=1;
     71 int pd=0,shang=-1;
     72 while(pd<2&&gs_cnt!=0){
     73  int sze=you_dui[1][1];
     74  int jvli=na(1);
     75  //for(int x=1;x<=gs_cnt;x++)cout<<you_dui[x][0]<<" ";
     76 //  cout<<endl;
     77 //  system ("pause");
     78  if(sze==zhong&&shang!=jvli){pd++;shang=jvli;continue;}
     79  for(int hao=head[sze];hao!=-1;hao=lian[hao].nxt){
     80   int dao=lian[hao].to;
     81   if(jv2[sze]+lian[hao].quan!=jv2[dao])fang(dao,jv2[sze]+lian[hao].quan+jv[dao]);
     82   if(jv2[sze]+lian[hao].quan<jv2[dao]){jv2[dao]=jv2[sze]+lian[hao].quan;}  
     83  }
     84 // for(int x=1;x<=gs_cnt;x++)cout<<you_dui[x][0]<<" ";
     85 // cout<<endl;
     86 //  system ("pause");    
     87 }
     88 if(pd<2){cout<<"-1";fout<<"-1";return;}
     89 cout<<shang;fout<<shang;
     90 return;
     91 }
     92 
     93 
     94 int na(int sze){
     95 int a=you_dui[sze][1];
     96 you_dui[sze][1]=you_dui[gs_cnt][1];
     97 you_dui[gs_cnt][1]=a;
     98 a=you_dui[sze][0];you_dui[sze][0]=you_dui[gs_cnt][0];
     99 you_dui[gs_cnt][0]=a;gs_cnt--;
    100 int now=1,zi=0;
    101 while(now*2<=gs_cnt){
    102  zi=now*2;
    103  if(now*2<gs_cnt&&you_dui[zi+1][0]<you_dui[zi][0])zi++;
    104  if(you_dui[zi][0]>you_dui[now][0])return a;
    105  int b=you_dui[zi][0];you_dui[zi][0]=you_dui[now][0];
    106  you_dui[now][0]=b;b=you_dui[zi][1];
    107  you_dui[zi][1]=you_dui[now][1];you_dui[now][1]=b;
    108  now=zi;
    109 }    
    110 return a;    
    111 }
    112 
    113 
    114 void fang(int sze,int zhi){
    115 you_dui[++gs_cnt][1]=sze;    
    116 you_dui[gs_cnt][0]=zhi;    
    117 int now=gs_cnt,fu=0;
    118 while(now/2>0){
    119  fu=now/2;
    120  if(you_dui[fu][0]<you_dui[now][0])return;    
    121  int a=you_dui[now][0];you_dui[now][0]=you_dui[fu][0];
    122  you_dui[fu][0]=a;a=you_dui[now][1];you_dui[now][1]=you_dui[fu][1];
    123  you_dui[fu][1]=a;now=fu;    
    124 }    
    125 return;    
    126 }
    127 
    128 
    129 int main(int argc, char** argv) {
    130 fin>>places>>rodes;
    131 memset(head,-1,sizeof(head));
    132 for(int x=1;x<=rodes;x++){
    133 int a,b,c;
    134 fin>>a>>b>>c;
    135 add(a,b,c);
    136 add(b,a,c);    
    137 }    
    138 fin>>qi>>zhong;
    139 spfa( );
    140 a_xin( );
    141     
    142  return 0;
    143 }
  • 相关阅读:
    ibatis常用sql
    在eclipse中部署maven项目的问题
    如何成为一个设计师和程序员混合型人才
    一个程序员的读书笔记:程序设计的反思
    C# 中的 == 和 equals()有什么区别?
    2014百度之星资格赛解题报告:能量变换
    2014百度之星资格赛解题报告:Xor Sum
    2014百度之星资格赛解题报告:Labyrinth
    那些年我们一起追过的ACM
    最新全球排名前50网站前端开发语言统计
  • 原文地址:https://www.cnblogs.com/Ateisti/p/5053446.html
Copyright © 2011-2022 走看看