1 #include <fstream>
 2 #include <iostream>
 3 #include <stdlib.h>
 4 #include <conio.h>
 5 using namespace std;
 6 
 7 #define MAX_v 50 //定义一个最大顶点个数
 8 typedef struct{
 9     int a[MAX_v][MAX_v]; //无向图G的邻接矩阵
10     int v; //无向图G的顶点
11     int e; //无向图G的边
12     int x[50]; //顶点与当前团的连接，x[i]=1 表示有连接——即x[i]==1代表在当前最大团的解内
13     int bestx[50]; //当前最优解
14     int cnum; //当前团的顶点数目
15     int bestn; //最大团的顶点数目
16 }MCP;
17 
18 void Creat(MCP &G);
19 void Backtrack(MCP &G,int i);
20 
21 void Creat(MCP &G){
22     int i,j;
23     ifstream fin("data.txt");
24     if (!fin)
25     {
26         cout<<"不能打开文件:"<<"data.txt"<<endl;
27         exit(1);
28     }
29     fin>>G.v;
30     for (int i=1;i<=G.v;i++)
31         for (int j=1;j<=G.v;j++)
32             fin>>G.a[i][j];
33     for(i=1;i<=G.v;i++) //初始化
34     {
35         G.bestx[i]=0;
36         G.x[i]=0;
37         G.bestn=0;
38         G.cnum=0;
39     }
40     cout<<"———回溯法求解最大团问题———"<<endl;
41     cout<<"输入初始化无向图矩阵为:"<<endl; //初始化
42     for(i=1;i<=G.v;i++)
43     {
44         for(j=1;j<=G.v;j++)
45             cout<<G.a[i][j]<<" ";
46         cout<<endl;
47     }
48 }
49 
50 void Backtrack(MCP &G,int i){
51     if (i>G.v){        //output()阶段
52         for (int j=1; j<=G.v; j++)
53             G.bestx[j] = G.x[j];    //记录最优解
54         G.bestn =G.cnum;
55         return ;
56     }
57     //检查顶点i与当前团的连接
58     int OK = 1;
59     for (int j=1; j<=i ; j++)
60         if (G.x[j]&& G.a[i][j]==0){        //G.x[j]:顶点j在当前解的最大团内；G.a[i][j]:待考察i顶点与最大团中前i-1个顶点间边的关系
61             //i不与j相连
62             OK = 0;
63             break;
64         }
65         if (OK) {        //进入左子树
66             G.x[i] = 1;//把i加入团
67             G.cnum++;
68             Backtrack(G,i+1);
69             G.x[i]=0;
70             G.cnum-- ;
71         }
72         if (G.cnum+G.v- i>G.bestn){        //进入右子树——剪枝函数
73             G.x[i] = 0;
74             Backtrack(G,i+1);
75         }
76 }
77 
78 int main(){
79     MCP G;
80     Creat(G);
81     Backtrack(G,1);
82     cout<<"最大团包含的顶点数为:"<<G.bestn<<endl;
83     cout<<"最大团方案为:( ";
84     for (int i=1;i<=G.v;i++)
85         if(G.bestx[i]==1){
86             cout<<i<<" ";
87         }
88         cout<<")"<<endl;
89         getch();
90 }