bzoj3244 [Noi2013]树的计数
Description
我们知道一棵有根树可以进行深度优先遍历(DFS)以及广度优先遍历(BFS)来生成这棵树的DFS序以及BFS序。两棵不同的树的DFS序有可能相同,并且它们的BFS序也有可能相同,例如下面两棵树的DFS序都是1 2 4 5 3,BFS序都是1 2 3 4 5
现给定一个DFS序和BFS序,我们想要知道,符合条件的有根树中,树的高度的平均值。即,假如共有K棵不同的有根树具有这组DFS序和BFS序,且他们的高度分别是h1,h2,...,hk,那么请你输出(h1+h2..+hk)/k
Input
有3行。
第一行包含1个正整数n,表示树的节点个数。
第二行包含n个正整数,是一个1~n的排列,表示树的DFS序。
第三行包含n个正整数,是一个1~n的排列,表示树的BFS序。
输入保证至少存在一棵树符合给定的两个序列。
Output
仅包含1个实数,四舍五入保留恰好三位小数,表示树高的平均值。
Sample Input
1 2 4 5 3
1 2 3 4 5
Sample Output
HINT
【评分方式】
如果输出文件的答案与标准输出的差不超过0.001,则将获得该测试点上的分数,否则不得分。
【数据规模和约定】
20%的测试数据,满足:n≤10;
40%的测试数据,满足:n≤100;
85%的测试数据,满足:n≤2000;
100%的测试数据,满足:2≤n≤200000。
【说明】
树的高度:一棵有根树如果只包含一个根节点,那么它的高度为1。否则,它的高度为根节点的所有子树的高度的最大值加1。
对于树中任意的三个节点a , b , c ,如果a, b都是c的儿子,则a, b在BFS序中和DFS序中的相对前后位置是一致的,即要么a都在b的前方,要么a都在b的后方。
正解:结论题。
一道神奇的人类智慧题。。
考虑把点按照$dfs$序编号,然后我们研究一下$bfs$序的一些性质。
如果$bfn_{i}>bfn_{i+1}$,那么很显然,$i+1$这个点一定在$i$的下一层,它对答案的贡献为$1$;
如果$bfn_{i}+1=bfn_{i+1}$,且$i+1$到$n$能形成一段连续的$dfn$区间,那么$i+1$这个点有两种可能,一种是$i$的兄弟,一种是$i$的儿子,这两种情况概率是相等的,那么它对答案的贡献就是$0.5$;
如果$bfn_{i}+1<bfn_{i+1}$,那么很显然,$i+1$对答案没有任何贡献。
然后我们就这么做完这题了。。
1 //It is made by wfj_2048~ 2 #include <algorithm> 3 #include <iostream> 4 #include <cstring> 5 #include <cstdlib> 6 #include <cstdio> 7 #include <vector> 8 #include <cmath> 9 #include <queue> 10 #include <stack> 11 #include <map> 12 #include <set> 13 #define inf (1<<30) 14 #define N (500010) 15 #define il inline 16 #define RG register 17 #define ll long long 18 #define File(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout) 19 20 using namespace std; 21 22 int dfn[N],bfn[N],l[N],r[N],n; 23 double ans; 24 25 il int gi(){ 26 RG int x=0,q=1; RG char ch=getchar(); 27 while ((ch<'0' || ch>'9') && ch!='-') ch=getchar(); 28 if (ch=='-') q=-1,ch=getchar(); 29 while (ch>='0' && ch<='9') x=x*10+ch-48,ch=getchar(); 30 return q*x; 31 } 32 33 il void work(){ 34 n=gi(); if (n==1){ puts("1.000"); return; } 35 for (RG int i=1;i<=n;++i) dfn[gi()]=i,ans=2; 36 for (RG int i=1;i<=n;++i) bfn[i]=dfn[gi()]; l[n+1]=n+1; 37 for (RG int i=n;i;--i) l[i]=min(l[i+1],bfn[i]),r[i]=max(r[i+1],bfn[i]); 38 for (RG int i=2;i<n;++i) 39 if (bfn[i]>bfn[i+1]) ++ans; 40 else if (bfn[i]+1==bfn[i+1] && r[i+1]-l[i+1]+1==n-i) ans+=0.5; 41 printf("%0.3lf\n",ans); return; 42 } 43 44 int main(){ 45 File("count"); 46 work(); 47 return 0; 48 }