1. 程式人生 > 其它 >CF11 D. A Simple Task

CF11 D. A Simple Task

題目傳送門:https://codeforces.com/problemset/problem/11/D

題目大意:
給定一個簡單無向圖,求圖中簡單環的個數(\(n\leqslant 19\)


\(n\)​很小的話,考慮狀壓

\(F[S][i]\)表示點集\(S\)中遍歷到\(i\)的方案數

轉移的話,我們列舉下一個點\(j\),如果\(j\notin S\),則有\(F[S|2^j][j]+=F[S][i]\);如果\(j\in S\),則成環,有\(Ans+=F[S][i]\)

但這樣存在一個問題,由於起點不固定,每個大小為\(Size\)的環都會被算\(Size\)

所以我們可以固定點集\(S\)

中,編號最小的點為起點,這樣每個環僅會被算兩次(順時針逆時針)

此外,由於沒有記錄邊是否被使用過,因此兩點(有邊相連)也會被算入答案。又因為兩點不存在順逆時針差異,故僅會被算一次

綜上所述,答案為\(\frac{Ans-m}{2}\)

/*program from Wolfycz*/
#include<map>
#include<cmath>
#include<cstdio>
#include<vector>
#include<cstring>
#include<iostream>
#include<algorithm>
#define Fi first
#define Se second
#define ll_inf 1e18
#define MK make_pair
#define sqr(x) ((x)*(x))
#define pii pair<int,int>
#define int_inf 0x7f7f7f7f
using namespace std;
typedef long long ll;
typedef unsigned int ui;
typedef unsigned long long ull;
inline char gc(){
	static char buf[1000000],*p1=buf,*p2=buf;
	return p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin),p1==p2)?EOF:*p1++;
}
template<typename T>inline T frd(T x){
	int f=1; char ch=gc();
	for (;ch<'0'||ch>'9';ch=gc())	if (ch=='-')    f=-1;
	for (;ch>='0'&&ch<='9';ch=gc())	x=(x<<1)+(x<<3)+ch-'0';
	return x*f;
}
template<typename T>inline T read(T x){
	int f=1; char ch=getchar();
	for (;ch<'0'||ch>'9';ch=getchar())	if (ch=='-')	f=-1;
	for (;ch>='0'&&ch<='9';ch=getchar())	x=(x<<1)+(x<<3)+ch-'0';
	return x*f;
}
inline void print(int x){
	if (x<0)	putchar('-'),x=-x;
	if (x>9)	print(x/10);
	putchar(x%10+'0');
}
const int N=20;
bool Map[N+10][N+10];
ll F[(1<<N)+10][N+10];
int lowbit(int x){for (int i=0;;i++)	if ((1<<i)&x)	return i;}
int main(){
//	freopen(".in","r",stdin);
//	freopen(".out","w",stdout);
	int n=read(0),m=read(0);
	ll Ans=0;
	for (int i=1;i<=m;i++){
		int x=read(0)-1,y=read(0)-1;
		Map[x][y]=Map[y][x]=1;
	}
	for (int i=0;i<n;i++)	F[1<<i][i]=1;
	for (int sta=1;sta<1<<n;sta++){
		for (int i=0;i<n;i++){
			if (!F[sta][i])	continue;
			int p=lowbit(sta);
			for (int j=p;j<n;j++){
				if (!Map[i][j])	continue;
				if (sta&(1<<j)){
					if (j==p)
						Ans+=F[sta][i];
				}else	F[sta|(1<<j)][j]+=F[sta][i];
			}
		}
	}
	printf("%lld\n",(Ans-m)>>1);
	return 0;
}

CF

作者:Wolfycz 出處:http://wolfycz.cnblogs.com/ 本文版權歸作者和部落格園共有,歡迎轉載,但必須在文章開頭註明原文出處,否則保留追究法律責任的權利