Codeforces Beta Round #11 D. A Simple Task(状压dp/无向图简单环的个数)

题目

给定n(n<=19)个点，m(0<=m)的无向图，求简单环的个数

简单环是指一个没有经过重复的点，也没有经过重复的边的环

思路来源

https://www.luogu.com.cn/problem/solution/CF11D

https://blog.csdn.net/fangzhenpeng/article/details/49078233（非常详细）

题解

洛谷题解，真是不知比CF题解高到哪里去了……

这样来看，眼下进步可能是一个时间问题orz，是之前没找对方法……

首先考虑，3个点、4个点、…、19个点的环的级别大致是的，

也就是说可达19!级别，故是一个NP问题，所以要用状压压去一些重复运算

我们只关心一个固定起点，访问了一些点之后，到达了哪个终点，

中间把经过的点记录下来，在没出现环之前，

只有起点、终点、中间点集是需要关心的，

而中间点集的顺序是不需要考虑的，计一个总的路径条数就行了

 

一个圆排列只被计一次，所以把圆排列的最小的点固定为起点，就减少了重复的统计次数

dp[i][j]表示从lowbit(i)出发，当前访问集合为i，且到达了j的简单路径数，

简单环的个数，就用lowbit(i)->j，j再拓展回lowbit(i)时来统计，计最后答案为ans

但一个长度大于等于3的环lowbit(i)->j->...->p->lowbit(i)仍然会被统计两次，

一次是从p扩展到lowbit(i)，另一次是从j扩展到lowbit(i)，

此外，lowbit(i)->j时，j也可以扩展到lowbit(i)，也就是无向图的反向边被误认为是长度为2的环，共m条边

所以答案为(ans-m)/2

 

注意扩展的时候，只能从一个固定起点的状态扩展，不能中途扩展了一半变更起点，

统计答案时，分v是否在已有状态i里讨论，只关心两种情况，

一种情况是，v==lowbit(i)，此时找到简单环

另一种情况是，v没访问过，lowbit(i)->j->v，路径延长

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define pb push_back
const int N=19,M=1<<19;
ll dp[M][N],ans;//dp[i][j]表示当前从lowbit(i)出发，访问了i内的点，且到达j的简单路径数
int n,m,u,v;
vector<int>e[N];
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;++i){
        scanf("%d%d",&u,&v);
        u--;v--;
        e[u].pb(v);e[v].pb(u);
    }
    for(int i=0;i<n;++i){
        dp[1<<i][i]=1;
    }
    int up=1<<n;
    for(int i=1;i<up;++i){
        for(int j=0;j<n;++j){
            if(!dp[i][j])continue;
            for(int k=0;k<e[j].size();++k){
                int v=e[j][k];
                if((i&(-i))>(1<<v))continue;//刷后继 为保证无后效性 不能中途变更起点
                if(i&(1<<v)){
                    if((i&(-i))==(1<<v)){
                        ans+=dp[i][j];
                    }
                }
                else{
                    dp[i|(1<<v)][v]+=dp[i][j];//此时v不是i|1<<v最小的 因此保证了起始点不变
                }
            }
        }
    }
    //无向图 每条边正反建 长度为2的环会被dp[12][2]在2->1拓展时统计一次 共m条边
    //此外 dp[123][3] 3->1 和dp[123][2] 2->1 是一个环 顺时针逆时针各一次 每个简单环统计两次
    printf("%lld\n",(ans-m)/2);
    return 0;
}