graph.h

 

 
  
  1. #pragma once 
  2. #include"UFset.h" 
  3. #define wuqiong -1 
  4. class bian; 
  5. class graph{ 
  6.     int** juzhen; 
  7.     int num; 
  8.     string *name; 
  9. public: 
  10.     graph(int); 
  11.     ~graph(); 
  12.     void show()const; 
  13.     void printEdge(int)const; 
  14.     void insert(int,int,int); 
  15.     void del(int ,int); 
  16.     void DFS(int)const; 
  17.     void BMS(int)const; 
  18.     bian* prim(); 
  19.     bian* Kruskal(); 
  20.     void Dijstra(int); 
  21.     void floyd(); 
  22.  
  23. }; 
  24. graph::graph(int n):num(n){ 
  25.     juzhen=new int*[num]; 
  26.     for(int i=0;i!=num;i++) 
  27.         juzhen[i]=new int[num]; 
  28.     for(int i=0;i<n;i++)  
  29.         for(int j=0;j!=num;j++) 
  30.             juzhen[i][j]=0; 
  31.     name=new string[num]; 
  32.     cout<<"输入各节点名字"<<endl; 
  33.     for( int i=0; i<num;i++ ){ 
  34.         cin>>name[i]; 
  35.         for(int j=0;j<i;j++) 
  36.             if(name[i]==name[j] ){ 
  37.                 cout<<"重名了，重输!"; 
  38.                 i--; 
  39.                 break; 
  40.             } 
  41.     } 
  42. } 
  43. void graph::show()const{ 
  44.     for( int i=0; i!=num;i++ ){ 
  45.         for(int j=0;j!=num;j++) 
  46.             printf("%6d",juzhen[i][j]); 
  47.         cout<<endl; 
  48.     } 
  49.     cout<<endl; 
  50. } 
  51. void graph::printEdge(int k)const{ 
  52.     for( int i=0; i!=num;i++ ) 
  53.         if(juzhen[k][i] ) 
  54.             cout<<k<<"->"<<i<<" "<<juzhen[k][i]; 
  55.     cout<<endl; 
  56. } 
  57. void graph::insert(int f,int l,int q){ 
  58.     if(f>=num||l>=num ){ 
  59.         cout<<"输入越界"<<endl; 
  60.         return ; 
  61.     } 
  62.     juzhen[f][l]=q; 
  63.     juzhen[l][f]=q; 
  64. } 
  65. graph::~graph(){ 
  66.     int *l=(int*)juzhen; 
  67.     delete []l; 
  68.     juzhen=0; 
  69.     l=0; 
  70. } 
  71. void graph::del(int f ,int l){ 
  72.     if(f>=num||l>=num ){ 
  73.         cout<<"输入越界"<<endl; 
  74.         return ; 
  75.     } 
  76.     if(juzhen[f][l]==0){ 
  77.         cout<<"没有删除意义"<<endl; 
  78.         return ; 
  79.     } 
  80.  
  81.     juzhen[f][l]=0; 
  82.     juzhen[l][f]=0; 
  83. } 
  84. void graph::DFS(int f )const{ 
  85.     if(f>=num||f<0 ){ 
  86.         cout<<"输入错误！"<<endl; 
  87.         return ; 
  88.     } 
  89.     bool * mark=new bool[num]; 
  90.     for(int i=0;i<num;i++)  
  91.         mark[i]=false; 
  92.     stack<int> zhan; 
  93.     zhan.push(f); 
  94.     int now=f; 
  95.     cout<<name[f]<<" "; 
  96.     mark[f]=true; 
  97.     for( int i=0;i<num;i++ ){ 
  98.         if(juzhen[now][i] && mark[i]==false){ 
  99.             cout<<name[i]<<" "; 
  100.             mark[i]=true; 
  101.             zhan.push(i); 
  102.             now=i; 
  103.             i=-1; 
  104.             continue; 
  105.         } 
  106.         if(i==num-1 ){ 
  107.             if(zhan.size()==0){ 
  108.                 for( i=0;i!=num ;i++ ){ 
  109.                     if(mark[i]==false ){ 
  110.                         now=i; 
  111.                         cout<<name[i]<<" "; 
  112.                         mark[i]=true; 
  113.                         i=0; 
  114.                         zhan.push(i); 
  115.                         break; 
  116.                     } 
  117.                 } 
  118.                 if(zhan.size()!=0) 
  119.                     continue; 
  120.                 break; 
  121.             } 
  122.             now=zhan.top(); 
  123.             zhan.pop(); 
  124.             i=-1; 
  125.         } 
  126.     } 
  127.     cout<<endl; 
  128. } 
  129. void graph::BMS(int f)const{ 
  130.     queue<int> duilie; 
  131.     int now=f; 
  132.     cout<<name[now]<<" "; 
  133.     bool *mark=new bool[this->num]; 
  134.     memset(mark,0,num); 
  135.     mark[now]=true; 
  136.     while(now!=-1){ 
  137.         for( int i=0;i< this->num ;i++ ){ 
  138.             if( juzhen[now][i]!=0 && mark[i]==false ){ 
  139.                 duilie.push(i ); 
  140.                 cout<<name[i]<<" "; 
  141.                 mark[i]=true; 
  142.             } 
  143.         } 
  144.         if( ! duilie.empty()){ 
  145.             now=duilie.front(); 
  146.             duilie.pop(); 
  147.         } 
  148.         else{ 
  149.             for(int j=0 ; j<num;j++ ){ 
  150.                 if( ! mark[j] ){ 
  151.                     duilie.push(j); 
  152.                     cout<<name[j]<<" "; 
  153.                     now=j; 
  154.                     mark[j]=true; 
  155.                     break; 
  156.                 } 
  157.                 now=-1; 
  158.             } 
  159.              
  160.         } 
  161.  
  162.     } 
  163.     cout<<endl; 
  164. } 
  165. bian* graph::prim( ){ 
  166.     bian* minTree=new bian[num-1]; 
  167.     int treeNum=0; 
  168.     vector<bian> next; 
  169.     int now=0; 
  170.     bool *mark=new bool[num]; 
  171.     memset(mark,true,num ); 
  172.     mark[0]=false; 
  173.     while(treeNum!=num-1){ 
  174.         for( int i=0;i<num ;i++ ) 
  175.             if(juzhen[now][i] && mark[i] ) 
  176.                 next.push_back( bian(now,i,juzhen[now][i]) ); 
  177.         vector<bian>::iterator p=next.begin(); 
  178.         for( vector<bian>::iterator bianli=next.begin()+1 ; bianli!=next.end() ;bianli++ ){ 
  179.             if( p->quanzhi > bianli->quanzhi && mark[bianli->jieshu] ){ 
  180.                 p=bianli; 
  181.                 now=bianli->jieshu; 
  182.             } 
  183.         } 
  184.         minTree[treeNum++]=*p; 
  185.         next.erase(p); 
  186.         mark[now]=0; 
  187.     } 
  188.     return minTree; 
  189. } 
  190. bian* graph::Kruskal(){ 
  191.     UFset *p; 
  192.     p=new UFset(this->num); 
  193.     bian* k=new bian[num*(num+1)/2 ]; 
  194.     int geshu=0; 
  195.     for( int i=0; i<num;i++ ){ 
  196.         for( int j=0; j+i<num;j++ ){ 
  197.             if(juzhen[ i][j+i ]!=0 ) 
  198.                 k[geshu++]=bian(i,j+i,juzhen[ i][ j+i] ); 
  199.         } 
  200.     } 
  201.     qsort(k,geshu,sizeof(bian),aa); 
  202.     int h=0; 
  203.     bian *re=new bian[num]; 
  204.     for( int i=0;i<geshu && h<num ;i++ ){ 
  205.         if( p->isWith(k[i].qishi,k[i].jieshu) ) 
  206.             continue; 
  207.         re[h++]=k[i]; 
  208.         p->Union(k[i].jieshu,k[i].qishi); 
  209.     } 
  210.     if(h==num ) 
  211.         return re; 
  212.     else 
  213.         return 0; 
  214. } 
  215. //算法关键在于什么时候更新path 
  216. void  graph::Dijstra(int qidian ){ 
  217.     int *path=new int[num]; 
  218.     int now =qidian; 
  219.     int *mark=new int[num]; 
  220.     int *data=new int[num]; 
  221.     memset(mark,0,sizeof(int)*num); 
  222.     for(int i=0;i<this->num;i++) 
  223.         data[i]=wuqiong; 
  224.     data[now]=0; 
  225.     mark[now]=true; 
  226.     for(int i=0;i<this->num;i++) 
  227.         path[i]=now; 
  228.     for( int ii=0; ii<this->num-1;ii++ ){ 
  229.         //更新data 
  230.         for( int j=0; j<num;j++ ){ 
  231.             if(this->juzhen[now][j]!=0&&mark[j]==false){ 
  232.                 //什么时候更新path是重点！ 
  233.                 if(data[j]==wuqiong ){ 
  234.                     path[j]=now; 
  235.                     data[j]=juzhen[now][j]; 
  236.                 } 
  237.                 if((data[now]+juzhen[now][j])<data[j] ){ 
  238.                     data[j]=data[now]+juzhen[now][j]; 
  239.                     path[j]=now; 
  240.                 } 
  241.             } 
  242.         } 
  243.         int y=0; 
  244.         while(mark[y])y++; 
  245.         for( int i=y+1; i<num;i++ ) 
  246.             if(mark[i]==false && data[i]<data[y] && data[i]!=wuqiong ) 
  247.                 y=i; 
  248.         mark[y]=true; 
  249.         now=y; 
  250.     } 
  251.     //print 
  252.     for(int i=0;i<num;i++){ 
  253.         int j=i; 
  254.         cout<<qidian<<"->"<<i<<" :"<<i<<" "; 
  255.         while(j!=qidian ){ 
  256.             cout<<path[j]<<" "; 
  257.             j=path[j]; 
  258.         } 
  259.         cout<<endl; 
  260.     } 
  261. } 
  262. void graph::floyd(){ 
  263.     //初始化 
  264.     int** adj; 
  265.     adj=new int*[this->num]; 
  266.     for(int i=0 ;i<num;i++) 
  267.         adj[i]=new int[num]; 
  268.     for(int i=0;i<num;i++) 
  269.         for( int j=0;j<num;j++){ 
  270.             if(juzhen[i][j]) 
  271.                 adj[i][j]=juzhen[i][j]; 
  272.             else 
  273.                 adj[i][j]=wuqiong; 
  274.         } 
  275.     int **path; 
  276.     path=new int*[this->num]; 
  277.     for(int i=0 ;i<num;i++) 
  278.         path[i]=new int[num]; 
  279.     for(int i=0;i<num;i++) 
  280.         for( int j=0;j<num;j++) 
  281.             path[i][j]=i; 
  282.     //算法执行 
  283.     for(int i=0;i<num;i++ ){ 
  284.         for(int j=0;j<num;j++){ 
  285.             if(i==j) 
  286.                 continue; 
  287.             for( int y=0; y<num;y++){ 
  288.                 if( adj[j][i]!=wuqiong && adj[i][y]!=wuqiong && (adj[j][i]+adj[i][y]<adj[j][y] ||adj[j][y]==wuqiong)){ 
  289.                     path[j][y]=i; 
  290.                     adj[j][y]=adj[j][i]+adj[i][y]; 
  291.                 } 
  292.             } 
  293.         } 
  294.     } 
  295.     for( int i=0; i<num;i++ ){ 
  296.         cout<<"从"<<i<<"开始"<<endl; 
  297.         for(int j=0;j<num;j++){ 
  298.              
  299.             int u=j; 
  300.             cout<<i<<"->"<<j<<" :"<<" "; 
  301.             while( u!=i ){ 
  302.                 cout<<u<<" "; 
  303.                 u=path[i][u]; 
  304.             } 
  305.             cout<<i<<endl; 
  306.         } 
  307.     } 
  308.  
  309. } 
 
  
bian.h
 
  
 
 
  
 
  
  1. #pragma once 
  2. #include<iostream> 
  3. #include<stack> 
  4. #include<string> 
  5. #include<queue> 
  6. #include<cstdlib> 
  7. #include<vector> 
  8. #include<cstdlib> 
  9. #include<algorithm> 
  10. #include<map> 
  11. using namespace std; 
  12.  
  13. class bian{ 
  14.     friend class graph; 
  15.     int qishi; 
  16.     int jieshu; 
  17.     int quanzhi; 
  18. public: 
  19.     void set(int a,int b,int c){ 
  20.         this->jieshu=b; 
  21.         this->qishi=a; 
  22.         this->quanzhi=c; 
  23.     } 
  24.     bian(int a=0,int b=0,int c=0){ 
  25.         this->jieshu=b; 
  26.         this->qishi=a; 
  27.         this->quanzhi=c; 
  28.     } 
  29.     void show(){ 
  30.         cout<<qishi<<"->"<<jieshu<<" : "<<quanzhi<<endl; 
  31.     } 
  32.     bool operator>(bian k){ 
  33.         return quanzhi-k.quanzhi>0; 
  34.     } 
  35.     bool operator<(bian k){ 
  36.         return quanzhi-k.quanzhi<0; 
  37.     } 
  38.     bool operator==(bian k){ 
  39.         return quanzhi==k.quanzhi; 
  40.     } 
  41.     int operator-(bian k ){ 
  42.         return quanzhi-k.quanzhi; 
  43.     } 
  44. }; 
  45. int aa(const void *a,const void *b ){ 
  46.         return *(bian*)a-*(bian*)b; 
  47. } 
 
  
MaxHeap.h
 
  
 
 
  
 
  
  1. #pragma once 
  2. #include<iostream> 
  3. #include<algorithm> 
  4. using namespace std; 
  5.  
  6. template<class T> 
  7. class MaxHeap{ 
  8.     T* heapArray; 
  9.     int CurrentSize; 
  10.     int MaxSize; 
  11. public: 
  12.     MaxHeap(T*  ,int,int  ); 
  13.     void insert(const T&); 
  14.     void del(); 
  15.     void shift(int ); 
  16.     void show(); 
  17. }; 
  18.  
  19. template<class T> 
  20. MaxHeap<T>::MaxHeap(T* a,int c,int m ){ 
  21.     this->CurrentSize=c; 
  22.     this->MaxSize=m; 
  23.     this->heapArray=new T[c]; 
  24.     for( int i=0;i<c ;i++ ) 
  25.         this->heapArray[i]=a[i]; 
  26.     int p=(this->CurrentSize-1)/2; 
  27.     while(p!=-1 ){ 
  28.         this->shift(p--); 
  29.     } 
  30. } 
  31. template<class T> 
  32. void MaxHeap<T>::shift(int num ){ 
  33.     if( num*2+1 >= this->CurrentSize ) 
  34.         return ; 
  35.     int p; 
  36.     if( num*2+2 < this->CurrentSize && this->heapArray[num*2+1] < this->heapArray[num*2+2] ) 
  37.         p=num*2+2; 
  38.     else  
  39.         p=num*2+1; 
  40.     if(this->heapArray[p]>this->heapArray[num] ){ 
  41.         swap(this->heapArray[p],this->heapArray[num]); 
  42.         this->shift(p); 
  43.     } 
  44. } 
  45. template<class T> 
  46. void MaxHeap<T>::show(){ 
  47.     for(int i=0;i<this->CurrentSize;i++) 
  48.         cout<<this->heapArray[i]<<" "; 
  49.     cout<<endl; 
  50. } 
  51. template<class T> 
  52. void MaxHeap<T>::insert(const T& p){ 
  53.     if(this->CurrentSize==this->MaxSize) 
  54.         return ; 
  55.     this->CurrentSize++; 
  56.     this->heapArray[this->CurrentSize-1]=p; 
  57.     int y=this->CurrentSize-1; 
  58.     while(y){ 
  59.         y=(y-1)/2; 
  60.         this->shift(y); 
  61.     } 
  62. } 
  63. template<class T> 
  64. void MaxHeap<T>::del(){ 
  65.     swap(this->heapArray[0],this->heapArray[this->CurrentSize-1]); 
  66.     this->CurrentSize--; 
  67.     this->shift(0);  
  68. } 
 
  
UFset.h
 
  
 
 
  
 
  
  1. #pragma once 
  2. #include"bian.h" 
  3.  
  4. class UFset{ 
  5.     int num; 
  6.     int *root; 
  7.     int *next; 
  8.     int *length; 
  9. public: 
  10.     UFset(int); 
  11.     bool isWith(int,int); 
  12.     void Union(int,int); 
  13. }; 
  14.  
  15. UFset::UFset(int n):num(n){ 
  16.     this->root=new int[this->num]; 
  17.     this->next=new int[this->num]; 
  18.     this->length=new int[this->num]; 
  19.     memset(this->length,0,num*4); 
  20.     for( int i=0; i<num;i++ ) 
  21.         next[i]=root[i]=i; 
  22. } 
  23.  
  24. bool UFset::isWith(int a,int b){ 
  25.     return root[a]==root[b]; 
  26. } 
  27.  
  28. void UFset::Union(int a,int b){ 
  29.     int l=this->length[a]>this->length[b]?a:b;//目标根->l 
  30.     int r=this->root[l]; 
  31.     int y=next[a]; 
  32.     next[a]=next[b]; 
  33.     next[b]=y; 
  34.     y=root[l]; 
  35.     while(root[ next[l] ]!=y ){ 
  36.         root[ next[l] ]=y; 
  37.         l=next[next[l]]; 
  38.     } 
  39. } 
 
  
图定义.cpp(测试代码)
 
  
 
 
  
 
  
  1. #include"graph.h" 
  2. //此图为无向图 
  3. int main(){ 
  4.     int num,f,l,p; 
  5.     cout<<"输入结点个数"<<endl; 
  6.     cin>>num; 
  7.     graph k(num); 
  8.     cout<<"下面输入节点（按ctrl +Z结束输入）"<<endl<<"f  l  p"<<endl; 
  9. /*  while( cin.good() ){ 
  10.         cin>>f; 
  11.         cin>>l; 
  12.         cin>>p; 
  13.         if(cin.good()) 
  14.             k.insert(f,l,p); 
  15.     } 
  16. */ 
  17.     k.insert(0,1,6); 
  18.     k.insert(0,2,1); 
  19.     k.insert(0,3,5); 
  20.     k.insert(1,2,5); 
  21.     k.insert(2,3,5); 
  22.     k.insert(1,4,3); 
  23.     k.insert(2,4,6); 
  24.     k.insert(3,5,2); 
  25.     k.insert(2,5,4); 
  26.     k.insert(4,5,6);  
  27.     cout<<"show():"<<endl; 
  28.     k.show(); 
  29. /*  cout<<"del(int,int)"<<endl; 
  30.     cin.clear(); 
  31.     cin>>f>>l; 
  32.     k.del(f,l); 
  33.     cout<<"show():"<<endl; 
  34.     k.show(); 
  35.     cout<<"k.DFS(0);"<<endl; 
  36.     k.DFS(0); 
  37.     cout<<"k.BMS(0);"<<endl; 
  38.     k.BMS(0); 
  39.     cout<<"k.prim()"<<endl; 
  40.     bian *op=k.prim(); 
  41.     for(int i=0;i!=num-1;i++) 
  42.         op[i].show(); 
  43.     cout<<endl; 
  44.     op=k.Kruskal(); 
  45.     if(op) 
  46.         for(int i=0;i!=num-1;i++) 
  47.             op[i].show(); 
  48.     else 
  49.         cout<<"no in!!!!"<<endl; 
  50.     */ 
  51.     k.Dijstra(0); 
  52.     k.floyd(); 
  53. }