[NWRRC2017]Fygon 2.0

最新推荐文章于 2022-08-09 21:16:18 发布

OneInDark

最新推荐文章于 2022-08-09 21:16:18 发布

阅读量150

点赞数

分类专栏：数学图论 C++ 文章标签： c++

本文链接：https://blog.csdn.net/qq_42101694/article/details/120876150

版权

C++ 同时被 3 个专栏收录

547 篇文章 3 订阅

订阅专栏

数学

280 篇文章 1 订阅

订阅专栏

图论

122 篇文章 0 订阅

订阅专栏

题目

传送门 to luogu

思路

刚拿到这道题，确实没有什么很好的思路。然后 $\sf HandInDevil$ 就把它秒切了……

事实上，循环变量 $j$ 在 $[l, r]$ 范围内枚举，等价于 对于任意 $j\in[1,n]$ 的情况，仅 $l\leqslant j\leqslant r$ 时统计循环次数。因为 $1\leqslant l\wedge r\leqslant n$ 嘛。

于是考虑建出大小关系的有向图， $a\leqslant b$ 则 $a$ 向 $b$ 连边。此时强连通应当取同一个值，缩点后是 $\rm DAG$ 。看到 $n\leqslant 20$ ，是不是心中有想法了呢？

还是不会诶……感觉转移比较困难。没想到这东西可以更抽象一点：渐进复杂度确实是 $n\rightarrow+\infty$ 的取值。那么 $\frac{f(n)}{n^k}\;(n\rightarrow+\infty)$ 的含义是什么？有 $k$ 个在 $[1, n]$ 以内随机取值的变量，要求它满足 $\rm DAG$ 所限定的大小关系的概率。

对于任意一种特定的大小关系 $x_1<x_2<\cdots<x_k$ ，其概率显然是 $\frac{1}{k!}$ 。而自变量取值相等的情况，在 $n\rightarrow+\infty$ 时，概率 $\frac{1}{n}\rightarrow0^+$ ，可以直接 忽略不计。

那么怎么计算所有可行的大小关系呢？其实就是求不同的拓扑排序结果，状压 $\tt dp$ 就好了！时间复杂度 $\mathcal O(n\log n)$ 。

代码

为了不写 tarjan 我进行了一些骚操作……

#include <cstdio>
#include <iostream>
#include <vector>
#include <cstring>
#include <algorithm>
using namespace std;
# define rep(i,a,b) for(int i=(a); i<=(b); ++i)
# define drep(i,a,b) for(int i=(a); i>=(b); --i)
typedef long long int_;
inline int readint(){
	int a = 0, c = getchar(), f = 1;
	for(; '0'>c||c>'9'; c=getchar())
		if(c == '-') f = -f;
	for(; '0'<=c&&c<='9'; c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}
void writeUnsigned(const unsigned &x){
	if(x > 9) writeUnsigned(x/10);
	putchar((x-x/10*10)^48);
}
inline void writeint(const int &x){
	if(x >= 0) writeUnsigned(x);
	else putchar('-'), writeUnsigned(-x);
}
inline void getMin(int &a,const int &b){ if(b < a) a = b; }
inline void getMax(int &a,const int &b){ if(a < b) a = b; }
inline int_ getGcd(int_ a,int_ b){
	for(; b; a%=b,swap(a,b));
	return a;
}

const int MAXN = 20;
int g[1<<MAXN];
bool scc[1<<MAXN];
int bel[MAXN]; /// belong to which @c scc
int f[MAXN][MAXN]; /// if two nodes are NOT connected

int_ dp[1<<MAXN], all;
int haxi[255];
char str[MAXN];
bool used[1<<MAXN];
int main(){
	int n = readint()-1;
	rep(i,0,n-1) fill(f[i],f[i]+n,1);
	rep(i,0,n-1){
		scanf("%*s %s",str);
		haxi[int(*str)] = i;
		scanf("%*s %s",str); // range(L,
		if(str[6] != '1'){
			const int &x = haxi[int(str[6])];
			g[1<<x] |= (1<<i), f[x][i] = 0;
		}
		scanf("%s",str); // R):
		if(*str != 'n'){
			const int &x = haxi[int(*str)];
			g[1<<i] |= (1<<x), f[i][x] = 0;
		}
	}
	rep(i,0,n-1) f[i][i] = 0; // self-loop
	rep(k,0,n-1) rep(i,0,n-1) rep(j,0,n-1)
		getMin(f[i][j],f[i][k]+f[k][j]);
	for(int S=1; S<(1<<n); ++S){
		g[S] = g[S&(-S)]|g[S&(S-1)];
		scc[S] = true; // check scc
		rep(i,0,n-1) if(S>>i&1){
			rep(j,i+1,n-1) if(S>>j&1)
				if(f[i][j] || f[j][i])
					scc[S] = false;
			break; // one node is enough
		}
		if(scc[S] == false) continue;
		rep(i,0,n-1) if(S>>i&1) bel[i] = S;
	}
	{ // calculate exponent
		int ddg = 0;
		rep(i,0,n-1){
			if(!used[bel[i]]) ++ ddg;
			used[bel[i]] = true;
		}
		printf("%d ",ddg);
		rep(i,0,n-1) used[bel[i]] = false;
		rep(i,int(all=1),ddg) all *= i;
	}
	dp[0] = 1; // no node deleted
	for(int S=0; S<(1<<n)-1; ++S){
		if(!dp[S]) continue;
		rep(i,0,n-1) if(!(S>>i&1)){
			int x = bel[i]; // whole scc
			if(used[x]) continue;
			used[x] = true; // visited
			if((g[x]&(S^x)) == g[x])
				dp[S^x] += dp[S];
		}
		rep(i,0,n-1)
			used[bel[i]] = false;
	}
	int_ d = getGcd(all,dp[(1<<n)-1]);
	printf("%lld/%lld\n",dp[(1<<n)-1]/d,all/d);
	return 0;
}