NP问题:基于无向图的最小顶点覆盖的混合贪心算法(MGA)
#include<cstdio> #include<cstring> #include<vector> #include<algorithm> #include<math.h> #include<iostream> #include<set> #define MAX 100 using namespace std; int n,m; struct Edge { int u,v; }; struct Point { int id; int degree,adj_degree; } point[MAX]; //点集 vector<Edge> E; //边集 vector<int> G[MAX]; //每个点的相邻点集(相当于一个二维变长数组) int degree[MAX],adj_degree[MAX];//每个点的度数和邻接度数 bool del[MAX],vis[MAX]; //已删除的点集、已访问过的点集 int V_cnt,Vv_cnt,E_cnt; //统计已经覆盖的顶点和边 void init_edge() { E.clear(); for(int i=1; i<=n; i++) G[i].clear(); } void add_edge(int u,int v) { E.push_back((Edge) {u,v} ); E.push_back((Edge) {v,u} ); int _size=E.size(); G[u].push_back(_size-2); G[v].push_back(_size-1); } bool cmp(Point a,Point b) { return a.adj_degree>b.adj_degree; } void pretreat() //预处理:计算每个顶点的度数和邻接度数 { for(int i=1; i<=n; i++) { point[i].degree=0; if(del[i]) continue; for(int j=0,_size=G[i].size(); j<_size; j++) { Edge& e=E[G[i][j]]; if(!del[e.v]) point[i].degree++; } //cout<<"第"<<i<<"个点的度数为:"<<point[i].degree<<endl; if(point[i].degree==0) del[i]=1; } for(int i=1; i<=n; i++) { if(del[i]) { point[i].adj_degree=0; continue; } point[i].adj_degree=point[i].degree; for(int j=0,_size=G[i].size(); j<_size; j++) point[i].adj_degree+=point[E[G[i][j]].v].degree; //cout<<"第"<<i<<"个点的邻接度数为:"<<point[i].adj_degree<<endl; } sort(point+1,point+n+1,cmp); } void MinVC_MGA(bool ans[]) //最小顶点的混合贪心算法 { memset(del,0,sizeof(del)); for(int i=1; i<=n; i++) point[i].id=i; E_cnt=0; //int flag=1; //测试子图划分次数 while(E_cnt<m) { memset(vis,0,sizeof(vis)); //每次选取最大度点删除后都需要还原访问数组 pretreat(); //重新计算每个点的度数 V_cnt=0; Vv_cnt=0; for(int i=1; i<=n; i++) if(!del[i]) Vv_cnt++; for(int i=1; i<=n; i++) { if(del[point[i].id]) continue; if(vis[point[i].id]) continue; //cout<<point[i].id<<endl; ans[point[i].id]=1; //加入到最小顶点覆盖集中 del[point[i].id]=1, Vv_cnt--; for(int j=0,_size=G[point[i].id].size(); j<_size; j++) { Edge& e=E[G[point[i].id][j]]; if(del[e.v]) continue; E_cnt++; if(!vis[e.v]) { vis[e.v]=1; V_cnt++; } } if(V_cnt>=Vv_cnt) break; //该子图中的点全部被访问过了 } //flag++; } //cout<<flag<<endl; } void fileInput(){ //从文件中读取网络信息 int countN=0,countM=0; set<int> s; FILE *fp; int u,v; char fname[100]; cout<<"\n请输入顶点网络(文件名):"; gets(fname); while((fp=fopen(fname,"r"))==NULL) {cout<<"文件名输入错误,请重新输入!\n";break;} while(fscanf(fp,"%d %d\n",&u,&v)!=EOF){ if(s.count(u)==0){s.insert(u);countN++;} if(s.count(v)==0){s.insert(v);countN++;} add_edge(u,v); countM++; } m=countM; n=countN; fclose(fp); } int main() { /* cout<<"请输入图中顶点个数:"; cin>>n; cout<<"请输入图中边数:"; cin>>m; cout<<"请输入边(u,v):"<<endl; for(int i=1,u,v; i<=m; i++) { cin>>u>>v; add_edge(u,v); } */ //int maxNum=(int)ceil(log(m)/log(2)); while(1){ init_edge(); fileInput(); bool ans[n]; //记录最小顶点覆盖集 memset(ans,0,sizeof(ans)); MinVC_MGA(ans); cout<<"最小控制顶点集为:"; for(int i=1; i<=n; i++) { if(ans[i]==true) cout<<i<<" "; } } return 0; }