hdu 3118 Arbiter 2009武汉网络赛 图论 二分图 状压

题目

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3118

题目来源：2009武汉网络赛 第四题左右

简要题意：给定一个$N$点$M$边的无向图，人在图上走，每走一步会在两种状态间相互切换。
$\qquad \qquad$人从任一个点出发时的状态和跑过若干点后回到该点的状态必须永远相同。

数据范围：$0 < T \leq 10 \\ 0 \leq N \leq 15;\quad 0 \leq M \leq 300\\ for\ every\ edge(u,v)\quad 0\leq u, v < N;\quad may\ have\ multiedge$

题解

题目其实挺巧妙的，首先必须要看出满足条件的图的充要条件是无奇环。

正确性证明：

• 充分条件：沿着奇环绕一圈就有了问题，所以有奇环肯定不满足条件，其逆否命题也成立，所以满足条件一定没奇环。
• 必要条件：如果没有奇环，那图中只有偶环。跑两个偶环还是偶环，而回到起点一定构成一个环，所以一定跑了一个偶环。绕一个偶环状态肯定不变，所以没有奇环一定满足条件，得证。

无奇环的图一定可以构成一个二分图，问题转化为找到一个边数最多的二分图。

进一步分析题目可以看到$N$的范围亲切地︿(￣︶￣)︿告诉你，你该玩个指数级的算法，于是我们可以使用状态压缩枚举二分图中点分属的两个集合。

比赛中带跑了学长，结果他去研究奇环了，图论真是一大坑啊。

实现

$\quad\ \ $建图，保存两点间边的数目即可。

$\quad\ \ $状态压缩，$01$表示分属的集合，即为无奇环的二分图的点分布情况。要删除的边数即为同属一个集合的点之间的边数。枚举的总数为$2^{N}$。

$\quad\ \ $由于边只需要计算一次，所以只需要枚举低编号到高编号的就可以了。如果有自环的话直接加到答案里就行了，不过数据里是没有的。自环的问题我之前没注意到，不过题目中的确说了是没自环的。

$\quad\ \ $去掉输入，总复杂度$\Theta(N^22^N)$。

代码

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>
#include <stack>
#include <queue>
#include <string>
#include <vector>
#include <set>
#include <map>

#define pb push_back
#define mp make_pair
#define all(x) (x).begin(),(x).end()
#define sz(x) ((int)(x).size())
#define fi first
#define se second
using namespace std;
typedef long long LL;
typedef vector<int> VI;
typedef pair<int,int> PII;
LL powmod(LL a,LL b, LL MOD) {LL res=1;a%=MOD;for(;b;b>>=1){if(b&1)res=res*a%MOD;a=a*a%MOD;}return res;}
// head
int t, n, m, u, v;
int g[20][20];

int main()
{
    scanf("%d", &t);
    while (t--) {
        scanf("%d%d", &n, &m);
        for (int i = 0; i < m; i++) {
            scanf("%d%d", &u, &v);
            g[u][v]++;
            g[v][u]++;
        }
        int MX = 1<<n, ans = 305;
        for (int i = 0; i < MX; i++) {
            int res = 0;
            for (int j = 0; j < n; j++) {
                bool a = (1<<j)&i, b;
                for (int k = j+1; k < n; k++) {
                    b = (1<<k)&i;
                    if (a == b) {
                        res += g[j][k];
                    }
                }
            }
            ans = min(res, ans);
        }
        printf("%d\n", ans);
        memset(g, 0, sizeof g);
    }
    return 0;
}