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

BZOJ3772 精神污染 【主席树 + dfs序】

时间:2018-01-20 15:04:43      阅读:171      评论:0      收藏:0      [点我收藏+]

标签:发展   dft   char   string   har   ons   void   date   vector   

题目

兵库县位于日本列岛的中央位置,北临日本海,南面濑户内海直通太平洋,中央部位是森林和山地,与拥有关西机场的大阪府比邻而居,是关西地区面积最大的县,是集经济和文化于一体的一大地区,是日本西部门户,海陆空交通设施发达。濑户内海沿岸气候温暖,多晴天,有日本少见的贸易良港神户港所在的神户市和曾是豪族城邑“城下町”的姬路市等大城市,还有以疗养地而闻名的六甲山地等。
兵库县官方也大力发展旅游,为了方便,他们在县内的N个旅游景点上建立了n-1条观光道,构成了一棵图论中的树。同时他们推出了M条观光线路,每条线路由两个节点x和y指定,经过的旅游景点就是树上x到y的唯一路径上的点。保证一条路径只出现一次。
你和你的朋友打算前往兵库县旅游,但旅行社还没有告知你们最终选择的观光线路是哪一条(假设是线路A)。这时候你得到了一个消息:在兵库北有一群丧心病狂的香菜蜜,他们已经选定了一条观光线路(假设是线路B),对这条路线上的所有景点都释放了【精神污染】。这个计划还有可能影响其他的线路,比如有四个景点1-2-3-4,而【精神污染】的路径是1-4,那么1-3,2-4,1-2等路径也被视为被完全污染了。
现在你想知道的是,假设随便选择两条不同的路径A和B,存在一条路径使得如果这条路径被污染,另一条路径也被污染的概率。换句话说,一条路径被另一条路径包含的概率。

输入格式

第一行两个整数N,M
接下来N-1行,每行两个数a,b,表示A和B之间有一条观光道。
接下来M行,每行两个数x,y,表示一条旅游线路。

输出格式

所求的概率,以最简分数形式输出。

输入样例

5 3

1 2

2 3

3 4

2 5

3 5

2 5

1 4

输出样例

1/3

样例解释

可以选择的路径对有(1,2),(1,3),(2,3),只有路径1完全覆盖路径2。

提示

100%的数据满足:N,M<=100000

题解

真是精神污染= =
我们求出有多少对路径存在包含关系即可

对于路径A,如果路径B的左右端点都在A的路径上,那么A包含B

我们对每个路径左端点开一个表,储存其右端点

这样我们只需要查询每个路径上点的右端点同时也在路径上的个数

可以用dfs序 + 主席树,每个点u对应的主席树代表着到根节点的右端点的信息。
在更新u时,先继承其父亲版本,然后对于u对应的所有右端点v,在v入度处+1,出序处-1
这样一来对于u和其祖先v,两点之间的有效节点数 = 两点入度之间的值之和
因为假若有一点不在这条路径上,却夹在两点入度之间,那么一定是出入度都在其中,一加一减就没了
如果在这条路径上, 一定只有入度,还没到出度

#include<iostream>
#include<cmath>
#include<cstdio>
#include<vector>
#include<cstring>
#include<algorithm>
#define LL long long int
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define Redge(u) for (int k = h[u],to; k; k = ed[k].nxt)
#define BUG(s,n) for (int i = 1; i <= (n); i++) cout<<s[i]<<‘ ‘; puts("");
using namespace std;
const int maxn = 100005,maxm = 4000005,INF = 1000000000;
inline int read(){
    int out = 0,flag = 1; char c = getchar();
    while (c < 48 || c > 57) {if (c == ‘-‘) flag = -1; c = getchar();}
    while (c >= 48 && c <= 57) {out = (out << 3) + (out << 1) + c - ‘0‘; c = getchar();}
    return out * flag;
}
int n,m,dfn[maxn],dft[maxn],fa[maxn][18],dep[maxn],cnt,h[maxn],ne = 2;
vector<int> end[maxn];
struct Que{int u,v;}q[maxn];
inline bool operator <(const Que& a,const Que& b){
    return a.u == b.u ? a.v < b.v : a.u < b.u;
}
struct EDGE{int to,nxt;}ed[2 * maxn];
void build(int u,int v){
    ed[ne] = (EDGE){v,h[u]}; h[u] = ne++;
    ed[ne] = (EDGE){u,h[v]}; h[v] = ne++;
}
int rt[maxn],siz,sum[maxm],ls[maxm],rs[maxm];
int update(int pre,int l,int r,int pos,int v){
    int u = ++siz; ls[u] = ls[pre]; rs[u] = rs[pre];
    if (l == r) {sum[u] = sum[pre] + v; return u;}
    int mid = l + r >> 1;
    if (mid >= pos) ls[u] = update(ls[pre],l,mid,pos,v);
    else rs[u] = update(rs[pre],mid + 1,r,pos,v);
    sum[u] = sum[ls[u]] + sum[rs[u]];
    return u;
}
int query(int a,int b,int c,int d,int l,int r,int L,int R){
    if (l >= L && r <= R) return sum[a] + sum[b] - sum[c] - sum[d];
    int mid = l + r >> 1;
    if (mid >= R) return query(ls[a],ls[b],ls[c],ls[d],l,mid,L,R);
    else if (mid < L) return query(rs[a],rs[b],rs[c],rs[d],mid + 1,r,L,R);
    else return query(ls[a],ls[b],ls[c],ls[d],l,mid,L,R) + query(rs[a],rs[b],rs[c],rs[d],mid + 1,r,L,R);
}
void dfs1(int u){
    dfn[u] = ++cnt;
    REP(i,17) fa[u][i] = fa[fa[u][i - 1]][i - 1];
    Redge(u) if ((to = ed[k].to) != fa[u][0]){
        fa[to][0] = u; dep[to] = dep[u] + 1; dfs1(to);
    }
    dft[u] = ++cnt;
}
void dfs2(int u){
    rt[u] = rt[fa[u][0]];
    for (unsigned int i = 0; i < end[u].size(); i++){
        rt[u] = update(rt[u],1,cnt,dfn[end[u][i]],1);
        rt[u] = update(rt[u],1,cnt,dft[end[u][i]],-1);
    }
    Redge(u) if ((to = ed[k].to) != fa[u][0]) dfs2(to);
}
int Lca(int u,int v){
    if (dep[u] < dep[v]) swap(u,v);
    for (int i = 0,d = dep[u] - dep[v]; (1 << i) <= d; i++)
        if ((1 << i) & d) u = fa[u][i];
    if (u == v) return u;
    for (int i = 17; i >= 0; i--)
        if (fa[u][i] != fa[v][i]) u = fa[u][i],v = fa[v][i];
    return fa[u][0];
}
int solve(int u,int v,int lca){
    int o = fa[lca][0],ans = 0;
    ans += query(rt[u],rt[v],rt[lca],rt[o],1,cnt,dfn[lca],dfn[u]);
    ans += query(rt[u],rt[v],rt[lca],rt[o],1,cnt,dfn[lca],dfn[v]);
    ans -= query(rt[u],rt[v],rt[lca],rt[o],1,cnt,dfn[lca],dfn[lca]);
    return ans - 1;
}
LL gcd(LL a,LL b){return !b ? a : gcd(b,a % b);}
int main(){
    n = read(); m = read();
    REP(i,n - 1) build(read(),read());
    REP(i,m) q[i].u = read(),q[i].v = read(),end[q[i].u].push_back(q[i].v);
    sort(q + 1,q + 1 + m);
    dfs1(1); dfs2(1);
    LL ans = 0,D = (LL)m * (m - 1) >> 1;
    REP(i,m) ans += solve(q[i].u,q[i].v,Lca(q[i].u,q[i].v));
    LL d = gcd(ans,D);
    printf("%lld/%lld\n",ans / d,D / d);
    return 0;
}

BZOJ3772 精神污染 【主席树 + dfs序】

标签:发展   dft   char   string   har   ons   void   date   vector   

原文地址:https://www.cnblogs.com/Mychael/p/8320543.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!