pku 2989

极大团问题：

#include <stdio.h> #include <map> using namespace std; #define MAXN 129 int n, m, cnt; struct f_int { unsigned int a,b,c,d; f_int() { a=b=c=d=0; } void operator <<(int s) { switch(--s/32) { case 0: a|=1<<s; break; case 1: b|=1<<s%32; break; case 2: c|=1<<s%32; break; case 3: d|=1<<s%32; } } bool operator==(f_int s) { if(a==s.a&&b==s.b&&c==s.c&&d==s.d) return true; return false; } bool Include(f_int s) { if((a==(a&s.a))&&(b==(b&s.b))&&(c==(c&s.c))&&(d==(d&s.d))) return true; return false; } f_int operator &(f_int s) { f_int r; r.a=a&s.a; r.b=b&s.b; r.c=c&s.c; r.d=d&s.d; return r; } bool operator <(const f_int s) const { if(d!=s.d) return d<s.d; if(c!=s.c) return c<s.c; if(b!=s.b) return b<s.b; if(a!=s.a) return a<s.a; return false; } int IDIN(int s) { switch(--s/32) { case 0: return a&1<<s; case 1: return b&1<<s%32; case 2: return c&1<<s%32; case 3: return d&1<<s%32; } } }; int main() { int a, b; while(scanf("%d%d",&n,&m)==2) { f_int st[MAXN],memory[1001]; cnt=1; memory[1].a=1; f_int tp; for(int i=0;i<m;++i) { scanf("%d%d",&a,&b); { st[a]<<b; st[b]<<a; } } for(int i=2; i<=n; ++i) { int cnt2=cnt; map<f_int,bool> hash; for(int j=1; j<=cnt2; ++j) { if(memory[j].Include(st[i])) { memory[j]<<i; continue; } tp=memory[j]&st[i]; if(!hash[tp]) { hash[tp]=true; tp<<i; int sign=0; for(int k=1; k<i; ++k) if(!tp.IDIN(k)&&tp.Include(st[k])) { sign=1; break; } if(sign) continue; memory[++cnt]=tp; if(cnt>1000) break; } } if(cnt>1000) break; } for(int i=1; i<cnt; ++i) for(int j=i+1; j<=cnt; ++j) if(memory[i]==memory[j]) printf("**%d%d**/n", i,j); if(cnt<=1000) printf("%d/n",cnt); else printf("Too many maximal sets of friends./n"); } return 0; }  

我还是写一下思路吧：

  这个就是个记忆搜索。一开始考虑问题规模为一个节点，这样的极大团是一。再考虑第二个点。这就有以下几种情况了。

  1.这个点和前面没这个点的问题中的某个极大团可以融合成一个新的极大团，这样子，极大团个数不会增加。

  2.要是这个点和前面的某个极大团不完全融合，那就找到它和那个极大团能融合的部分，这样肯定得到的是个完全子图，但是不是极大就不一定了。这要判断一下。遍历之前的点，看看有没有不在此团中但是可以融合进去的，有的话，这个团我们就不要了。后面肯定会有个更大的给你。为什么呢？假设当问题规模为i+1时，含点V(i+1)的极大团S，那么当问题为i规模时，是不是一定有个极大团为S-V(i+1),这是显然的。我把一些操作写进了struct f_int里，应该没什么别的了。

  主要是这个集合他不好表示。本来应该用hash的，但是比较麻烦，我也没写过，就用了map，集合用啥呢？用四个int就可以了，128位，代表128个数字。集合的交就是&，查找就是位移和1与一下看是不是0.

800多ms。。。

标程代码：

#include <stdio.h> #include <string.h> #include <stdlib.h> #define MAXN 128 #define SETELTS 4 #define MAXS 1000 unsigned amark = 2; typedef struct { unsigned e[SETELTS]; } set; set graph[MAXN]; set sets[MAXS + 1]; #define HSIZE 3007 struct hashel { unsigned hval; unsigned mark; set *s; } htab[HSIZE]; unsigned hval(set *s) { unsigned i, v = 0; for (i = 0; i < SETELTS; i++) v ^= s->e[i]; return v % HSIZE; } void hash (set *s) { unsigned val = hval (s), idx = val; while (htab[idx].mark == amark) { idx++; if (idx == HSIZE) idx = 0; } htab[idx].hval = val; htab[idx].mark = amark; htab[idx].s = s; } int set_eq (set *s1, set *s2) { unsigned i; for (i = 0; i < SETELTS; i++) if (s1->e[i] != s2->e[i]) return 0; return 1; } int hfind (set *s) { unsigned val = hval (s), idx = val; while (amark <= htab[idx].mark) { if (htab[idx].hval == val && set_eq (htab[idx].s, s)) return idx; idx++; if (idx == HSIZE) idx = 0; } return -1; } void unhash (set *s) { int index = hfind (s); if (index == -1) abort (); htab[index].mark = amark + 1; } int present_p (set *s) { return hfind (s) != -1; } void set_new (set *s) { memset (s, 0, sizeof (set)); } void set_add (set *s, int elt) { int eltno, bitno; eltno = elt / 32; bitno = elt % 32; s->e[eltno] |= 1 << bitno; } int in_set_p (set *s, int elt) { int eltno, bitno; eltno = elt / 32; bitno = elt % 32; return (s->e[eltno] >> bitno) & 1; } int set_intersect (set *a, set *b, set *c) { int i, ret = 1, val, new_val; for (i = 0; i < SETELTS; i++) { val = a->e[i]; new_val = val & b->e[i]; c->e[i] = new_val; ret &= (val == new_val); } return ret; } int max_clique_p (set *s, int n) { int i; set tmp; for (i = 0; i < n; i++) if (!in_set_p (s, i) && set_intersect (s, &graph[i], &tmp)) return 0; return 1; } int gen_cliques(int n) { int i, k, new_k, elt = n - 1; if (n == 1) { set_new (&sets[0]); set_add (&sets[0], 0); hash (&sets[0]); return 1; } k = gen_cliques (n - 1); if (k > MAXS) return k; new_k = k; for (i = 0; i < k; i++) { if (set_intersect (&sets[i], &graph[elt], &sets[new_k])) { unhash (&sets[i]); set_add (&sets[i], elt); hash (&sets[i]); } else { set_add (&sets[new_k], elt); if (!present_p (&sets[new_k]) && max_clique_p (&sets[new_k], n)) { hash (&sets[new_k]); new_k++; } } if (new_k > MAXS) return new_k; } return new_k; } int main(void) { int n, m, i, x, y, k; while (1) { if (scanf ("%d%d", &n, &m) != 2) return 0; for (i = 0; i < n; i++) set_new (&graph[i]); for (i = 0; i < m; i++) { scanf ("%d%d", &x, &y); x--;y--; set_add (&graph[x], y); set_add (&graph[y], x); } k = gen_cliques (n); if (k > MAXS) printf ("Too many maximal sets of friends./n"); else printf ("%d/n", k); amark+=2; } }  

70多ms。。。