码迷,mamicode.com
首页 > 其他好文 > 详细

2018-8-10考试 T3. 招募(akekure)

时间:2018-08-10 21:33:43      阅读:25      评论:0      收藏:0      [点我收藏+]

标签:考试   font   c++   [1]   str   招募   题解   else   space   

题目大意:有$n$个点和$m$条边的图($n - 1 \leq m \leq n + 5$),每个点要么黑要么白,两个黑点不可以相邻,问方案数

题解:可以发现当图为一棵树的时候只需要一个树形$DP$

$$令f_{i,j}表示在第i个点,它的状态为j(1为黑,0为白)$$
$$f_{i,0}=\prod\limits_{j为i的儿子}(f_{j,0}+f_{j,1})(因为它的儿子没有限制,可以黑可以白)$$
$$f_{i,0}=\prod\limits_{j为i的儿子}f_{j,0}(它的儿子有限制,必须为白)$$


再来考虑不为树的情况,可以发现这个图最多有$6$条返祖边,我们可以暴力枚举这几条边两端的情况,这样复杂度是$O(4^{m-n+1}\times n)$。

明显是不行的,可以发现两端为黑是一定不可能的,这样复杂度成为了$O(3^{m-n+1}\times n)$。

看来可以的,但是遍历图的常数太大,还是过不了,考虑优化。

发现其实对于一条边,可以变成只处理其中的一个点,这样就只有两种状态了,一种是这个点白点,另一个点就无限制;另一种是这个点为黑点,另一个点就是白点,复杂度就成了$O(2^{m-n+1}\times n)$,轻松$AC$

卡点:1.考试时复杂度为$O(3^{m-n+1}\times n)$

 

C++ Code:

#include <cstdio>
#include <cctype>
#define maxn 100100
using namespace std;
const long long mod = 1000000007;
const long long half = 500000004;
int n, m;
long long ans;
int fa[maxn], oth[50], tot;
int f[maxn][2];

int head[maxn], cnt;
struct Edge {
	int from, to, nxt;
	bool can;
} e[(maxn + 10) << 1];
void add(int a, int b) {
	e[++cnt] = (Edge) {a, b, head[a], true}; head[a] = cnt;
}

char *ch, op[1 << 25 | 100];
inline void read(int &x) {
	x = 0;
	while (isspace(*ch)) ch++;
	while (isdigit(*ch)) x = x * 10 + (*ch++ & 15);
}

bool vis[maxn];
void dfs(int rt) {
	vis[rt] = true;
	int v;
	for (register int i = head[rt]; i; i = e[i].nxt) {
		v = e[i].to;
		if (v != fa[rt]) {
			if (!vis[v]) {
				fa[v] = rt;
				dfs(v);
			} else {
				e[i].can = false;
				if (v > rt) oth[++tot] = i;
			}
		}
	}
}
bool ispoint[maxn], mus[maxn];
inline void real(int rt) {
	int v;
	int &_0 = f[rt][0], &_1 = f[rt][1];
	_0 = _1 = 1;
	if (ispoint[rt]) f[rt][mus[rt] ^ 1] = 0;
	for (int i = head[rt]; i; i = e[i].nxt) {
		if (e[i].can) {
			v = e[i].to;
			if (v != fa[rt]) {
				real(v);
				_0 = (1ll * _0 * (f[v][1] + f[v][0])) % mod;
				_1 = (1ll * _1 * f[v][0]) % mod;
			}
		}
	}
}
void run(int x) {
	if (x > tot) {
		real(1);
		ans = (ans + f[1][1] + f[1][0]) % mod;
		return ;
	}
	int u = e[oth[x]].from, v = e[oth[x]].to;
	bool &iu = ispoint[u], &iv = ispoint[v], &mu = mus[u], &mv = mus[v];
	if (iu && iv) {
		run(x + 1);
		return ;
	}
	if (iu || iv) {
		if (iv) u ^= v ^= u ^= v;
		if (mus[u]) {
			ispoint[v] = true;
			mus[v] = false;
			run(x + 1);
			ispoint[v] = false;
		} else {
			run(x + 1);
		}
		return ;
	}
	iu = true;
	mu = false;
	run(x + 1);
	iu = false;
	iu = iv = true;
	mu = true;
	mv = false;
	run(x + 1);
	iv = iu = false;
}
int main() {
	fread(ch = op, 1, 1 << 25, stdin);
	read(n), read(m);
	for (register int i = 0; i < m; i++) {
		int a, b;
		read(a), read(b);
		add(a, b);
		add(b, a);
	}
	dfs(1);
	run(1);
	printf("%lld\n", ans);
	return 0;
}

  

2018-8-10考试 T3. 招募(akekure)

标签:考试   font   c++   [1]   str   招募   题解   else   space   

原文地址:https://www.cnblogs.com/Memory-of-winter/p/9457065.html

(0)
(0)
   
举报
评论 一句话评论(0
0条  
登录后才能评论!
© 2014 mamicode.com 版权所有 京ICP备13008772号-2
迷上了代码!