炮兵阵地 【POJ - 1185】【最大独立集 or 补图的最大团】【讲解最大团】

题目链接

快速链接

---

  这道题，很多人用的是状压DP来写的……但我太弱了……只能写一下图论来做了。呜呜呜……

  这道题，可以是最大独立集的证明，我们知道最大独立集是NP问题，也就是所谓的我们直接去找最大独立集是困难的，时间上比验证它要大的多得多，所以，我们得反过来想，该如何去求解它，正如很多博主说的那样，我们得记住最大独立集是等同于其补图的最大团的。但是，尴尬的是，为什么就没有人详细的说明白什么是其补图？什么是最大团呢！？那我在这里讲一下我的理解，对于补图，我们把原来相互链接的边变成断开的边，把原来不链接的两点变成链接的两点，就是对应的bool类型的所有任意两点间的边取反即可。那么，什么又是最大团？我又要说一下我粗鄙的想法了，对于最大团，我的理解就是，我们要选取相互之间，任意两点均有链接的边的最大的关系图，就是说，我们得找这样的一个团子，团子里面任意的点间都是有边相连的，而我们现在所要找的就是令这个团子最大的方案，最大的团子，就是最大团。

  这道题，也确实卡了时间，我们如果按照顺序查询的话，一定是会T的，所以，我们得逆序来查询，知道对于最大团的处理问题，我们要做的就是搜索+回溯的过程，搜索的是对应节点u，其之后的节点，对于有与它有边链接（这里当然指的是处理补图的时候，所以图也是补图的图了），而这样的节点应当是从u之后去找，就应是（u+1 ~ num）其中，num是所有的节点的情况，这样子往后找也有个好处，就是后面的节点的最大团若是已经确定的话，那么如果改点已经建立起来的最大团加上后面点的最大团，若是小于等于现在ans的值，就说明是没用的了，就不需要去取了，不管就是了，然后，去遍历剩余节点的最大团，如果和它里面的最大团的所有的点都是可以匹配上的、有建边的，那么，我们就把目前的这个访问到的节点i丢进u的最大团中继续去查找有没有更优的解。我们重复上面的操作，直到终结的时候，何时终结？就是当没有比目前的ans更优的解的时候了。

  逆序遍历这样的u，从num->1这样子去，然后不断的去搜索+回溯，我们最后就会确定最大的ans的值！

---

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
#define INF 0x3f3f3f3f
using namespace std;
typedef unsigned long long ull;
typedef long long ll;
const int maxN = 1005;
int N, M, num;
struct Node
{
    int x, y;
    Node(int a=0, int b=0):x(a), y(b) {}
}a[maxN];
bool check(int i, int j)
{
    if(a[i].x == a[j].x) return abs(a[i].y - a[j].y) <= 2;
    if(a[i].y == a[j].y) return abs(a[i].x - a[j].x) <= 2;
    return false;
}
bool mp[maxN][maxN];
char s[105][15];
int ans, cnt[maxN], /*group[maxN],*/ vis[maxN];
bool dfs(int u, int pos)
{
    for(int i=u+1; i<=num; i++)
    {
        if(cnt[i] + pos <= ans) return false;
        if(mp[i][u])
        {
            int j;
            for(j=0; j<pos; j++) if(!mp[i][vis[j]]) break;
            if(j == pos)
            {
                vis[pos] = i;
                if(dfs(i, pos+1)) return true;
            }
        }
    }
    if(ans < pos)
    {
        //for(int i=0; i<pos; i++) group[i] = vis[i];
        ans = pos;
        return true;
    }
    return false;
}
void MaxClique()
{
    ans = 0;    //初始化必须为0，考虑到没有点的情况（没法建任何的阵地）
    for(int i=num; i>0; i--)
    {
        vis[0] = i; //最大团中的第一个元素
        dfs(i, 1);
        cnt[i] = ans;   //记忆化的过程
    }
}
int main()
{
    while(scanf("%d%d", &N, &M)!=EOF)
    {
        num = 0;
        for(int i=1; i<=N; i++)
        {
            scanf("%s", s[i]+1);
            for(int j=1; j<=M; j++) if(s[i][j] == 'P') a[++num] = Node(i, j);
        }
        for(int i=1; i<=num; i++)
        {
            for(int j=1; j<i; j++)
            {
                mp[i][j] = mp[j][i] = !check(i, j);
            }
        }
        MaxClique();
        printf("%d\n", ans);
    }
    return 0;
}