代码源每日一题div1 环的数量

最新推荐文章于 2024-09-16 15:03:26 发布

学不会数据库

最新推荐文章于 2024-09-16 15:03:26 发布

阅读量206

点赞数 1

分类专栏：代码源每日一题div1 文章标签： c++ 算法

本文链接：https://blog.csdn.net/qq_53287227/article/details/123980155

版权

代码源每日一题div1 专栏收录该内容

16 篇文章 1 订阅

订阅专栏

这篇博客探讨了一种利用动态规划和位运算解决图论问题的方法，具体是通过状压DP和lowbit函数计算经过特定点的路径数量。在有向图中，博主详细解释了如何初始化状态、进行状态转移并修正环的计数，最终计算出不重复的路径数目。文章深入浅出，适合对图论和动态规划有兴趣的读者。

摘要由CSDN通过智能技术生成

Daimayuan Online Judge

思路：状压dp
根据 $n$ 的范围，很容易判断要用状压 $d p$ ，我们可以用一个二进制数来表示当前这 $n$ 个点，哪些点被访问，例如： $11001$ 表示点 $1, 4, 5$ 被访问
那么我们定义 $d p [i] [u]$ 表示经过的点的状态为 $i$ 且最后一个点为 $u$ 的方案数，那我们考虑起点就是状态 $i$ 这个二进制数的最低位的 $1$ 代表的点
如何求这个起点呢？我们想到 $l o w b i t$ 函数求的是一个二进制数的最低位的 $1$ 代表的值，在这基础上取对数即为这个最低位的 $1$ 的位数，即：

    int lowbit(int x){
	    return (int)log2((x)&(-x));
    }

那在 $d p [i] [u]$ 的基础上，求出起点 $k = l o w b i t (i)$ ，如果当前点 $u$ 与起点 $k$ 连通，那么 $a n s + = d p [i] [u]$
接下来枚举 $v\in a[u]$ ， $v$ 不在状态 $i$ 中且 $v$ 不会为起点( $v > k$ )，进行状态转移： $d p [i ∣ (1 < < v)] [v] + = d p [i] [u]$
由于在计算过程中，我们把图当作了有向图，一个环被计算了两次，并且每条边也被当作了环，也要减去

#include <bits/stdc++.h>
using namespace std;
#define IOS ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);
#define endl '\n';
typedef long long ll;
typedef pair<ll,ll> PII;
const int N=1e6+10;
int dp[1<<20][20],n,m,a[20][20];
int lowbit(int x){
	return (int)log2((x)&(-x));
}
int main(){
	IOS;
	cin>>n>>m;
	for(int i=1;i<=m;i++){
		int u,v;cin>>u>>v;
		a[u-1][v-1]=a[v-1][u-1]=1;
	}
	for(int i=0;i<n;i++) dp[1<<i][i]=1;
	ll ans=0;
	for(int i=1;i<(1<<n);i++){
		int k=lowbit(i);
		for(int u=0;u<n;u++){
			if(!((i>>u)&1)) continue;
			if(a[u][k]) ans+=dp[i][u];
			for(int v=k+1;v<n;v++){
				if(!a[u][v]||((i>>v)&1)) continue;
				dp[i|(1<<v)][v]+=dp[i][u];
			}
		}
	}
	cout<<(ans-m)/2<<endl;
	return 0;
}