标签：tin set 个数 else 输出 测试数据 mil limit blog
题目描述
江爷爷给你出了一道题：
给你一个图，保证每个点最多属于一个简单环，每个点度数最多为3，求这个图有多少“眼镜图形个数”
保证图联通哦~
其中“眼镜图形个数”，定义为三元组(x,y,S)，其中x和y表示图上的两个点，S表示一条x到y的简单路径，而且必
须满足：
1.x和y分别在两个不同的简单环上
2.x所在的简单环与路径S的所有交点仅有x，y所在的简单环与路径S的所有交点仅有y。
(x,y,S)与(y,x,S)算同一个眼镜
如果你无法理解，可以参考样例。
保证图是联通的
输入输出格式
输入格式：
第一行两个数n和m
之后m行，每行两个数x，y表示x和y之间有一条边。
输出格式：
输出一个数，表示眼镜的个数对19260817取膜的结果
输入输出样例
输入样例#1：
说明
样例#3，#4，#5，#6见下发的文件
非常抱歉，出了点小锅，sample5.out好像是空文件，应该是6734568
【子任务】
11 12
1 2
2 3
3 4
4 5
5 1
4 6
6 7
7 8
8 9
9 10
10 11
11 7
输出样例#1：
1
输入样例#2：
14 16
1 2
2 3
3 4
4 1
3 5
5 6
6 7
7 8
8 9
9 6
9 13
13 14
13 10
10 11
11 12
12 10
输出样例#2：
4
子任务会给出部分测试数据的特点。
如果你在解决题目中遇到了困难, 可以尝试只解决一部分测试数据。
测试点编号 n的范围 m的范围 特殊性质
测试点1 n <= 10 m <= 20
测试点2 n <= 20 m <= 40
测试点3 n <= 20 m <= 40
测试点4 n <= 2000 m <= 4000
测试点5 n <= 2000 m <= 4000
测试点6 n <= 1000000 m <= 2000000 简单环个数 <= 2000
测试点7 n <= 1000000 m <= 2000000 简单环个数 <= 2000
测试点8 n <= 1000000 m <= 2000000
测试点9 n <= 1000000 m <= 2000000
测试点10 n <= 1000000 m <= 2000000
分析：有环在是不好处理的，先把所有环给缩成一个点.观察样例可以发现，如果环与环中间还有x个环，那么实际上这两个点可以构成2^x个眼镜，因为中间的每个环既可以走上面也可以走下面.于是dfs缩点似乎可以拿70分.
其实一个眼镜中间的圆点可以看做根节点，两个镜框可以看做是树，那么我们可以以1号点为根节点来建一棵树来进行树形dp.记录的状态不能是以i为根的子树的眼镜数，因为这样不好转移，应该将状态定义为半眼镜数.对于环和圆点要分别讨论，每个环可以看做是两条路径，所以在统计答案或者更新f的时候都要*2,因为每个子树是独立的，最后要把它们合并到一起，这里利用乘法原理来解决.可以得到f[x] = Σf[son[x]],if x为环 then f[x] = f[x] * 2 + 1,因为子树的可以走两个方向，自己也可以当做一个半眼镜的端点，合并子树的时候ans += f[x] * f[son[x]],if x是环 then ans += f[x] * f[son[x]].因为f[x]随着子节点的处理越来越大，所以答案不会遗漏.
#include <cstdio> #include <cstring> #include <iostream> #include <algorithm> using namespace std; const int maxn = 1000010, maxm = 4000010, mod = 19260817; int n, m, head[maxn], nextt[maxm], to[maxm], tot = 1, cnt, fa[maxn], scc[maxn], limit; int head2[maxn], to2[maxm], nextt2[maxm], tot2 = 1; bool vis[maxn]; long long ans, f[maxn]; void add(int x, int y) { to[tot] = y; nextt[tot] = head[x]; head[x] = tot++; } void add2(int x, int y) { to2[tot2] = y; nextt2[tot2] = head2[x]; head2[x] = tot2++; } void find(int x, int y) { scc[x] = cnt; if (x == y) return; find(fa[x], y); } void dfs(int u) { vis[u] = 1; for (int i = head[u]; i; i = nextt[i]) { int v = to[i]; if (v == fa[u] || scc[v]) continue; if (vis[v]) { cnt++; find(u, v); } else { fa[v] = u; dfs(v); } } } void build(int u) { vis[u] = 1; for (int i = head[u]; i; i = nextt[i]) { int v = to[i]; if (v == fa[u] || vis[v]) continue; if (scc[u] != scc[v]) { add2(scc[u], scc[v]); add2(scc[v], scc[u]); } build(v); } } void solve(int u, int from) { f[u] = 0; for (int i = head2[u]; i; i = nextt2[i]) { int v = to2[i]; if (v == from) continue; solve(v, u); long long temp = f[u] * f[v] % mod; if (u <= limit) temp *= 2; ans = (ans + temp) % mod; f[u] = (f[u] + f[v]) % mod; } if (u <= limit) { ans = (ans + f[u]) % mod; f[u] = (f[u] * 2 + 1) % mod; } } int main() { scanf("%d%d", &n, &m); for (int i = 1; i <= m; i++) { int x, y; scanf("%d%d", &x, &y); add(x, y); add(y, x); } dfs(1); limit = cnt; for (int i = 1; i <= n; i++) if (!scc[i]) scc[i] = ++cnt; memset(vis, 0, sizeof(vis)); build(1); solve(1, 0); printf("%lld\n", ans); return 0; }
标签：tin set 个数 else 输出 测试数据 mil limit blog
原文地址：http://www.cnblogs.com/zbtrs/p/7664023.html
