月度归档: 2018 年 7 月

[ZJOI2008]泡泡堂 题解

[ZJOI2008]泡泡堂 题解

题目地址:洛谷:【P2587】[ZJOI2008]泡泡堂 – 洛谷、BZOJ:Problem 1034. — [ZJOI2008]泡泡堂BNB

题目描述

第XXXX届NOI期间,为了加强各省选手之间的交流,组委会决定组织一场省际电子竞技大赛,每一个省的代表队由n名选手组成,比赛的项目是老少咸宜的网络游戏泡泡堂。每一场比赛前,对阵双方的教练向组委会提交一份参赛选手的名单,决定了选手上场的顺序,一经确定,不得修改。比赛中,双方的一号选手,二号选手……,n号选手捉对厮杀,共进行n场比赛。每胜一场比赛得2分,平一场得1分,输一场不得分。最终将双方的单场得分相加得出总分,总分高的队伍晋级(总分相同抽签决定)。
作为浙江队的领队,你已经在事先将各省所有选手的泡泡堂水平了解的一清二楚,并将其用一个实力值来衡量。为简化问题,我们假定选手在游戏中完全不受任何外界因素干扰,即实力强的选手一定可以战胜实力弱的选手,而两个实力相同的选手一定会战平。由于完全不知道对手会使用何种策略来确定出场顺序,所以所有的队伍都采取了这样一种策略,就是完全随机决定出场顺序。
当然你不想这样不明不白的进行比赛。你想事先了解一下在最好与最坏的情况下,浙江队最终分别能得到多少分。

题意简述

你和对手都有一队人,每边$n$个人。现在要求让你和对手的队伍两两对决,能力值高的一方获胜,得2分,若能力值相等平局,得1分,否则不得分。求出所有对决的设计中,你的最高可能得分与最低可能得分是多少。

输入输出格式

输入格式:
输入文件的第一行为一个整数n,表示每支代表队的人数。
接下来n行,每行一个整数,描述了n位浙江队的选手的实力值。
接下来n行,每行一个整数,描述了你的对手的n位选手的实力值。

输出格式:
输入文件中包括两个用空格隔开的整数,分别表示浙江队在最好与最坏的情况下分别能得多少分。不要在行末输出多余的空白字符。

输入输出样例

输入样例#1:

2
1
3
2
4

输出样例#1:

2 0

输入样例#2:

6
10000000
10000000
10000000
10000000
10000000
10000000
0
0
0
0
0
0

输出样例#2:

12 12

说明

样例说明
1: 我们分别称4位选手为A,B,C,D。则可能出现以下4种对战方式,最好情况下可得2分,最坏情况下得0分。
一 二 三 四
浙江 ??? 结果 浙江 ??? 结果 浙江 ??? 结果 浙江 ??? 结果
一号选手 A C 负 A D 负 B C 胜 B D 负
二号选手 B D 负 B C 胜 A D 负 A C 负
总得分 0 2 2 0
2: 对手都是认真学习的好孩子,不会打游戏。无论如何排列出场顺序都无法改变被蹂躏的结果。浙江队总能取得全胜的结果。

20%的数据中,1<=n<=10;
40%的数据中,1<=n<=100;
60%的数据中,1<=n<=1000;
100%的数据中,1<=n<=100000,且所有选手的实力值在0到10000000之间。

题解

贪心,我们想到了田忌赛马的模型,即用自己强的打对手强的,如果打不过就用一个最弱的送死。但是这种策略是有漏洞的,有可能拿去送死的最弱的可以将对手的最弱的打败,从而多获得2分,因此我们把这种情况也加入判断。
修改后的策略是,如果我方最强能赢对方最强,则进行这种对决;否则我方最弱如果能赢对方最弱,则以弱对弱;前两种情况都处理完了,再把一个最弱的拿去送死。
最低分数可以直接让对方跑一遍这个策略,然后用总得分减掉对方的最优得分即可。

代码

// Code by KSkun, 2018/7
#include <cstdio>
#include <cctype>
#include <cstring>

#include <algorithm>

typedef long long LL;

inline char fgc() {
    static char buf[100000], *p1 = buf, *p2 = buf;
    return p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 100000, stdin), p1 == p2)
        ? EOF : *p1++;
}

inline LL readint() {
    register LL res = 0, neg = 1; register char c = fgc();
    for(; !isdigit(c); c = fgc()) if(c == '-') neg = -1;
    for(; isdigit(c); c = fgc()) res = (res << 1) + (res << 3) + c - '0';
    return res * neg;
}

const int MAXN = 100005;

int n, slf[MAXN], opp[MAXN], ls, rs, lo, ro;

int main() {
    n = readint();
    for(int i = 1; i <= n; i++) {
        slf[i] = readint();
    }
    for(int i = 1; i <= n; i++) {
        opp[i] = readint();
    }
    std::sort(slf + 1, slf + n + 1);
    std::sort(opp + 1, opp + n + 1);
    int ans1 = 0, ans2 = 0;
    ls = lo = 1; rs = ro = n;
    while(ls <= rs) { 
        if(slf[rs] > opp[ro]) {
            ans1 += 2; rs--; ro--;
        } else if(slf[ls] > opp[lo]) {
            ans1 += 2; ls++; lo++;
        } else {
            ans1 += (slf[ls] == opp[ro]);
            ls++; ro--;
        }
    }
    ls = lo = 1; rs = ro = n;
    while(ls <= rs) {
        if(opp[ro] > slf[rs]) {
            ans2 += 2; rs--; ro--;
        } else if(opp[lo] > slf[ls]) {
            ans2 += 2; ls++; lo++;
        } else {
            ans2 += (opp[lo] == slf[rs]);
            lo++; rs--;
        }
    }
    printf("%d %d", ans1, 2 * n - ans2);
    return 0;
}
Codeforces Round #495 (Div. 2) 赛后总结

Codeforces Round #495 (Div. 2) 赛后总结

比赛地址:Dashboard – Codeforces Round #495 (Div. 2) – Codeforces
官方题解:Codeforces Round #495 (Div. 2) — Editorial – Codeforces

1004A Sonya and Hotels

题意简述

一维空间坐标轴上有$n$个整点位置修建了酒店,现在想在一个未建酒店的整点位置修一个新酒店,要求新酒店离最近的酒店的距离等于$d$,求新酒店可以选择的坐标数量。
数据范围:$1 \leq 100, 1 \leq d \leq 10^9$

思路

枚举每两个酒店之间的距离,$距离=2d$的时候对答案有1的贡献,$>2d$的时候有2的贡献,首尾各有1的贡献,统计一下就好。
复杂度$O(n)$。

代码

// Code by KSkun, 2018/7
#include <cstdio>
#include <cctype>
#include <cstring>

#include <algorithm>

typedef long long LL;

inline char fgc() {
    static char buf[100000], *p1 = buf, *p2 = buf;
    return p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 100000, stdin), p1 == p2)
        ? EOF : *p1++;
}

inline LL readint() {
    register LL res = 0, neg = 1; register char c = fgc();
    for(; !isdigit(c); c = fgc()) if(c == '-') neg = -1;
    for(; isdigit(c); c = fgc()) res = (res << 1) + (res << 3) + c - '0';
    return res * neg;
}

const int MAXN = 105;

int n, d, x[MAXN];

int main() {
    n = readint(); d = readint();
    for(int i = 1; i <= n; i++) {
        x[i] = readint();
    }
    int ans = 2;
    for(int i = 2; i <= n; i++) {
        if(x[i] - x[i - 1] == 2 * d) ans++;
        if(x[i] - x[i - 1] > 2 * d) ans += 2;
    }
    printf("%d", ans);
    return 0;
}

1004B Sonya and Exhibition

题意简述

要放$n$盆花成一排,每个位置有两种选择,有$m$个人参观这些花,他们会参观一个连续区间的花,参观的美丽度定义为区间两种花的个数的乘积,求安排方案使得美丽度之和最大。
数据范围:$1 \leq n, m \leq 10^3$

思路

$$ n^2 > (n+1)(n-1) > (n+2)(n-2) > \cdots $$
因此我们可以直接输出一个类似0101010101010101这样的序列就好了。
复杂度$O(n)$。

代码

// Code by KSkun, 2018/7
#include <cstdio>
#include <cctype>
#include <cstring>

#include <algorithm>

typedef long long LL;

inline char fgc() {
    static char buf[100000], *p1 = buf, *p2 = buf;
    return p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 100000, stdin), p1 == p2)
        ? EOF : *p1++;
}

inline LL readint() {
    register LL res = 0, neg = 1; register char c = fgc();
    for(; !isdigit(c); c = fgc()) if(c == '-') neg = -1;
    for(; isdigit(c); c = fgc()) res = (res << 1) + (res << 3) + c - '0';
    return res * neg;
}

const int MAXN = 1005;

int n, m, l[MAXN], r[MAXN];

int main() {
    n = readint(); m = readint();
    for(int i = 1; i <= m; i++) {
        l[i] = readint(); r[i] = readint();
    }
    LL sum1 = 0, sum2 = 0;
    for(int i = 1; i <= m; i++) {
        int cnt1[2] = {0, 0}, cnt2[2] = {0, 0};
        for(int j = l[i]; j <= r[i]; j++) {
            cnt1[j & 1]++;
            cnt2[(j & 1) ^ 1]++;
        }
        sum1 += 1ll * cnt1[0] * cnt1[1];
        sum2 += 1ll * cnt2[0] * cnt2[1];
    }
    if(sum1 > sum2) {
        fprintf(stderr, "%lld\n", sum1);
        for(int i = 1; i <= n; i++) {
            putchar('0' + (i & 1));
        }
    } else {
        fprintf(stderr, "%lld\n", sum2);
        for(int i = 1; i <= n; i++) {
            putchar('0' + ((i & 1) ^ 1));
        }
    }
    return 0;
}

1004C Sonya and Robots

题意简述

有一个长为$n$的序列,两个机器人从序列的两段开始相向而行,遇到第一个与机器人中存储的数字相同的数字时就在该位置停下,你可以给机器人写入不同的数字,求让它们不相遇的数字对数。
数据范围:$1 \leq n \leq 10^5$

思路

我们从左向右扫这个序列,对于一个之前没有出现过的数字,显然它会与后面的不重复数字产生贡献,那么答案就是这些贡献的和。我们只需要知道一个位置的后面是否还有与这个位置相同的数字,动态地维护后面的不重复数字种数以及前面出现过哪些数字就可以了。
复杂度$O(n)$。

代码

// Code by KSkun, 2018/7
#include <cstdio>
#include <cctype>
#include <cstring>

#include <algorithm>

typedef long long LL;

inline char fgc() {
    static char buf[100000], *p1 = buf, *p2 = buf;
    return p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 100000, stdin), p1 == p2)
        ? EOF : *p1++;
}

inline LL readint() {
    register LL res = 0, neg = 1; register char c = fgc();
    for(; !isdigit(c); c = fgc()) if(c == '-') neg = -1;
    for(; isdigit(c); c = fgc()) res = (res << 1) + (res << 3) + c - '0';
    return res * neg;
}

const int MAXN = 100005;

int n, a[MAXN], nxt[MAXN], head[MAXN];
bool vis[MAXN];

int main() {
    n = readint();
    int cnt = 0;
    LL ans = 0;
    for(int i = 1; i <= n; i++) {
        a[i] = readint();
        if(!vis[a[i]]) {
            vis[a[i]] = true;
            head[a[i]] = i;
            cnt++;
        } else {
            nxt[head[a[i]]] = i;
            head[a[i]] = i;
        }
    }
    memset(vis, 0, sizeof(vis));
    for(int i = 1; i <= n; i++) {
        if(!nxt[i]) {
            cnt--;
        }
        if(!vis[a[i]]) ans += cnt;
        vis[a[i]] = true;
    }
    printf("%lld", ans);
    return 0;
}

1004D Sonya and Matrix

题意简述

有一个矩阵,矩阵中有1个权值为0的格子,而其他格子的权值为到0权格子的曼哈顿距离。现在给你一个长为$t$的序列,是一个包含某一个矩阵里面的所有权值的乱序可重排列。现在你要求出原来矩阵的大小$n \times m$以及0权的位置$(x, y)$。
数据范围:$1 \leq t \leq 10^6$

思路

首先,我们可以用序列统计每个数字出现的次数。在完整的图形中,我们会发现每个数字形成了菱形或者说正方形的形状,如下所示。

          5
        5 4 5
      5 4 3 4 5
    5 4 3 2 3 4 5
  5 4 3 2 1 2 3 4 5
5 4 3 2 1 0 1 2 3 4 5
  5 4 3 2 1 2 3 4 5
    5 4 3 2 3 4 5
      5 4 3 4 5
        5 4 5
          5

我们可以求出序列中的最大数字以及最大的在矩阵中完整地出现了它的菱形的数字,显然,最大数字应该在角上,由于图形的对称性,四个角是等价的,我们暂且让它在左上角$(1, 1)$。
当我们发现矩形的大小并无法让0与最大数共存的时候,显然情况不合法,这种情况的判断,我们可以使用对角线曼哈顿距离(即图形中的最长曼哈顿距离)来判断,即$n+m-2 < mx$时不合法。
另外的不合法情况就是当确定了$n, m, x, y$以后,我们可以计算出最大的整个菱形都包含进来的数字,这个数字是$\min \{ x-1, y-1, n-x, m-y \}$,再把不合法的情况扔掉就好。
剩下的合法情况已经很少了,我们直接暴力$O(t)$地验证就好。

代码

// Code by KSkun, 2018/7
#include <cstdio>
#include <cctype>
#include <cstring>
#include <cmath>

#include <algorithm>

typedef long long LL;

inline char fgc() {
    static char buf[100000], *p1 = buf, *p2 = buf;
    return p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 100000, stdin), p1 == p2)
        ? EOF : *p1++;
}

inline LL readint() {
    register LL res = 0, neg = 1; register char c = fgc();
    for(; !isdigit(c); c = fgc()) if(c == '-') neg = -1;
    for(; isdigit(c); c = fgc()) res = (res << 1) + (res << 3) + c - '0';
    return res * neg;
}

const int MAXN = 1000005;

int t, mx, cnt[MAXN], cnt2[MAXN];

inline int border(int n, int m, int x, int y) {
    return std::min(std::min(x - 1, y - 1), std::min(n - x, m - y));
}

inline bool check(int n, int m, int x, int y) {
    memset(cnt2, 0, sizeof(cnt2));
    for(int i = 1; i <= n; i++) {
        for(int j = 1; j <= m; j++) {
            cnt2[abs(x - i) + abs(y - j)]++;
        }
    }
    for(int i = 1; i <= t; i++) {
        if(cnt2[i] != cnt[i]) return false;
    }
    return true;
}

int main() {
    t = readint();
    for(int i = 1; i <= t; i++) {
        int a = readint();
        cnt[a]++;
        mx = std::max(mx, a);
    }
    int lim = 0;
    for(int i = 1; i <= t; i++) {
        if(cnt[i] != i * 4) {
            lim = i - 1; break;
        }
    }
    for(int n = 1; n * n <= t; n++) {
        if(t % n) continue;
        int m = t / n;
        if(n + m - 2 < mx) continue;
        for(int j = 1; j <= n; j++) {
            int k = mx - j + 2;
            if(k > m || k < 1) continue;
            if(border(n, m, j, k) != lim) continue;
            if(check(n, m, j, k)) {
                printf("%d %d\n%d %d", n, m, j, k);
                return 0;
            }
        }
    }
    puts("-1");
    return 0;
}

1004E Sonya and Ice Cream

题意简述

给你一棵$n$个点的树,要求从中选出一条不超过$k$个点的路径,使得树的其他点到这条路径的最短距离最大值最小。
数据范围:$1 \leq k \leq n \leq 10^5$

思路

首先,我们可以证明,这条路径在包含树的中心点的时候比较优,中心点定义为到这个点最远的点距离最小的点,显然,这个点应该在直径的中心位置,可以$O(n)$求出来。
答案具有二分性质,所以我们二分答案,接着考虑如何验证。我们可以从中心点开始DFS,每次扩展出去一个点,如果发现一个点的子树中最远的点到它的距离大于答案,则这个点一定要被加入路径。如果一个点有多个儿子一定要被加入路径(中心点特殊考虑,它允许有2个儿子必须被加入路径)则不合法。如果必须被加入路径的点数多于$k$则不合法,判断一下就好。
复杂度$O(n \log n)$。

代码

// Code by KSkun, 2018/7
#include <cstdio>
#include <cctype>
#include <cstring>

#include <algorithm>
#include <vector>

typedef long long LL;

inline char fgc() {
    static char buf[100000], *p1 = buf, *p2 = buf;
    return p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 100000, stdin), p1 == p2)
        ? EOF : *p1++;
}

inline LL readint() {
    register LL res = 0, neg = 1; register char c = fgc();
    for(; !isdigit(c); c = fgc()) if(c == '-') neg = -1;
    for(; isdigit(c); c = fgc()) res = (res << 1) + (res << 3) + c - '0';
    return res * neg;
}

const int MAXN = 100005;

struct Edge {
    int to, w;
};
std::vector<Edge> gra[MAXN];

int n, k, du, dv, ct, fa[MAXN], dis[MAXN];

inline void dfs_dis(int u) {
    for(auto e : gra[u]) {
        if(e.to == fa[u]) continue;
        dis[e.to] = dis[u] + e.w;
        fa[e.to] = u;
        dfs_dis(e.to);
    }
}

inline void diameter() {
    dfs_dis(1);
    for(int i = 1; i <= n; i++) {
        if(dis[i] > dis[du]) du = i;
    }
    dis[du] = 0;
    memset(fa, 0, sizeof(fa));
    dfs_dis(du);
    for(int i = 1; i <= n; i++) {
        if(dis[i] > dis[dv]) dv = i;
    }
}

inline void center() {
    for(int i = dv; i; i = fa[i]) {
        if(std::max(dis[ct], dis[dv] - dis[ct]) > std::max(dis[i], dis[dv] - dis[i])) {
            ct = i;
        }
    }
}

int cnt;
bool success;

inline int dfs_check(int u, int fa, int lim) {
    int res = 0, big = 0;
    for(auto e : gra[u]) {
        if(e.to == fa) continue;
        int dis = dfs_check(e.to, u, lim);
        if(dis != -1 && dis + e.w > lim) {
            cnt++; big++;
        }
        if(dis == -1) big++;
        else res = std::max(res, dis + e.w);
    }
    if((u == ct && big > 2) || (u != ct && big > 1)) {
        success = false;
    }
    if(big || u == ct) {
        cnt++; return -1;
    }
    return res;
}

inline bool check(int mid) {
    cnt = 0; success = true;
    dfs_check(ct, 0, mid);
    return success && cnt <= k;
}

int main() {
    n = readint(); k = readint();
    for(int i = 1, u, v, w; i < n; i++) {
        u = readint(); v = readint(); w = readint();
        gra[u].push_back(Edge {v, w});
        gra[v].push_back(Edge {u, w});
    }
    diameter(); center();
    int l = -1, r = 1e9, mid;
    while(r - l > 1) {
        mid = (l + r) >> 1;
        if(check(mid)) r = mid; else l = mid;
    }
    printf("%d", r);
    return 0;
}

1004F Sonya and Bitwise OR

题意简述

维护一个长为$n$的序列,有$m$次操作,支持两种操作

  1. 单点修改
  2. 区间询问有多少连续子序列的异或和不小于$x$(对于所有的询问$x$都相同)

数据范围:$1 \leq n, m \leq 10^5$

思路

咕咕咕。

代码

咕咕咕。

[BZOJ3683]Falsita 题解

[BZOJ3683]Falsita 题解

题目地址:BZOJ:Problem 3683. — Falsita

题目描述

描述
到海边了呢……
如果没有那次选择,现在是不是会好些呢……
都过去了。
仰望着星空,迎面吹过一阵阵海风,倚靠着护栏,Fine 在海边静静地伫立着,在一个个无际的长夜后,Fine 终于放下了往事的痛楚,得到了治愈。
但是作为 Fine 的另一重人格的 Falsita 就没那么幸运了。她仍然被各种繁忙的事务困扰着。
虽然同在一副躯体中,Fine 与 Falsita 的精神世界却相差甚远,Fine 可以轻易地构造出幻梦时,Falsita 却只能停留在现实的痛楚中。
但是为了生活需要,她们还是需要经常达成共识。
让我们形式化的描述一下吧。
她们所在的精神世界是一棵以 1 号节点为根的树,每个树上的节点 u 都有一个权值Wu,她们每个人分别都在一个节点上,达成共识的方法就是两个人都到达一个共识节点(即到达它们的最近公共祖先)。
一个点 u 与另外一个点 v 之间想要达到共识需要花费的代价为Wu+Wv。
有时两人的精神有所触动时,有时点的权值会改变成某个数,有时以某个点的子树中的所有点的权值会加上某个数。
Falsita 和 Fine 经常需要达成共识,每一次询问,已知达成的共识节点,求她们花费的期望代价。

题意简述

有一棵有点权的树,以1为根。有三种操作:

  1. 单点点权加一个数
  2. 子树点权加一个数
  3. 求随机从$u$子树中选择两个点$(i, j)$且它们的LCA为$u$,$w_i+w_j$的期望

输入输出格式

输入格式:
输入共 m + 3 行。
第一行两个整数 n, m ,表示节点个数和操作数。
第二行 n – 1 个整数Pi,表示节点 i ( i = 2 . . . n ) 的父亲节点的编号。
第三行 n 个整数Wi。
接下来 m 行,每行表示一个操作。
1. S u delta 表示将节点 u 的权值加上 delta 。
2. M u delta 表示将以节点 u 为根的子树中的所有节点的权值加上 delta。
3. Q u 表示询问共识节点为 u 时的答案。
询问保证 u 不是叶子节点。

输出格式:
对于每组询问,输出答案,答案精确到小数点后 6 位。
你的程序输出的答案需要与标准答案之差不超过10^(-5)。

输入输出样例

输入样例#1:

4 6
1 2 2 
0 -6 3 0 
S 2 -5
M 3 8
S 1 -1
M 4 7
M 3 2
Q 1

输出样例#1:

2.000000

说明

前5个操作后,四个节点的权值分别为-1,-11,13,7。最近公共祖先为1的点对有(1,2),(1,3),(1,4),权值和分别为-12,12,6,故答案为(-12+12+6)/3=2。
对于 100% 的数据,1 ≤ n, m ≤ 300000, 0≤ |delta|, |Wi|≤ 20000。

题解

首先,期望可以表示为所有可能情况的和/情况数。对于确定的点$u$,它子树中所有可能情况的数量是可以预处理出来的。这个答案就是$u$子树中不同时在$u$的某个儿子的子树内的点对数,我们可以DFS树的使用利用子树大小来计算,假如计算$u$的方案数,$u$有儿子$v$,在遍历到$v$的时候,只需要对方案数加上$siz[v] \cdot (已经遍历过的子树大小和+1)$即可。
接下来考虑期望中分子的问题。如果不带修改,我们可以DFS的时候计算出所有节点的答案。如果用$sum[u]$表示$u$子树内所有节点的权值和,那么$u$有儿子$v$,$v$子树对$u$的答案贡献是$sum[v] \cdot (siz[u] – siz[v])$。
接下来我们考虑修改的问题,单点修改时,对该点答案的影响是$delta \cdot (siz[u] – 1)$,若该点有祖先$anc$,该点到该祖先路径上,祖先的儿子节点是$anc’$,那么对于该点的该祖先的影响是$delta \cdot (siz[anc] – siz[anc’])$。我们考虑进行树剖,在线段树上来维护这些个影响。显然,我们可以维护一个标记$delta \cdot (siz[u] – siz[v])$,其中$v$是$u$在重链上的儿子,这样,单点修改$u$就可以从$u$的父亲开始跳重链改了。但是当我们发现,一条重链最深的节点的重儿子可能不是修改点到根路径上该点的儿子,因此最深节点的情况要特别处理。我们直接把这个特别处理的修改改到该点处就好了。
对于子树加的情况,我们发现,只需要对子树里所有的节点的分子加上$delta \cdot 分母$就好了。而对于子树根的祖先们,其实跟单点修改产生的影响相同,只不过这里的影响是$delta \cdot siz[u]$。
总之,我们树剖以后在线段树上把这些个标记维护好就能实现正确的复杂度$O(n \log^2 n)$。

代码

// Code by KSkun, 2018/7
#include <cstdio>
#include <cctype>
#include <cstring>

#include <algorithm>
#include <vector>

typedef long long LL;

inline char fgc() {
    static char buf[100000], *p1 = buf, *p2 = buf;
    return p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 100000, stdin), p1 == p2)
        ? EOF : *p1++;
}

inline LL readint() {
    register LL res = 0, neg = 1; register char c = fgc();
    for(; !isdigit(c); c = fgc()) if(c == '-') neg = -1;
    for(; isdigit(c); c = fgc()) res = (res << 1) + (res << 3) + c - '0';
    return res * neg;
}

inline bool isop(char c) {
    return c == 'S' || c == 'M' || c == 'Q';
}

inline char readop() {
    char c;
    while(!isop(c = fgc())) {}
    return c;
}

const int MAXN = 300005;

int n, m, w[MAXN];
std::vector<int> gra[MAXN];

int fa[MAXN], siz[MAXN], dep[MAXN], son[MAXN], dfn[MAXN], vn[MAXN], top[MAXN], clk;
LL ans1[MAXN], ans2[MAXN], sum[MAXN];

inline void dfs1(int u) {
    siz[u] = 1; sum[u] = w[u];
    for(int i = 0; i < gra[u].size(); i++) {
        int v = gra[u][i];
        if(v == fa[u]) continue;
        fa[v] = u; dep[v] = dep[u] + 1;
        dfs1(v);
        ans2[u] += 1ll * siz[v] * siz[u];
        siz[u] += siz[v];
        sum[u] += sum[v];
        if(siz[v] > siz[son[u]]) son[u] = v;
    }
}

inline void dfs2(int u, int tp) {
    dfn[u] = ++clk; vn[dfn[u]] = u; top[u] = tp;
    ans1[u] += 1ll * w[u] * (siz[u] - 1);
    if(son[u]) {
        ans1[u] += 1ll * (siz[u] - siz[son[u]]) * sum[son[u]];
        dfs2(son[u], tp);
    }
    for(int i = 0; i < gra[u].size(); i++) {
        int v = gra[u][i];
        if(v == fa[u] || v == son[u]) continue;
        ans1[u] += 1ll * (siz[u] - siz[v]) * sum[v];
        dfs2(v, v);
    }
}

#define lch o << 1
#define rch o << 1 | 1
#define mid ((l + r) >> 1)

LL tag1[MAXN << 2], tag2[MAXN << 2];
// 1: ptadd 2: subtadd

inline void pushdown(int o, int l, int r) {
    if(tag1[o]) {
        tag1[lch] += tag1[o]; tag1[rch] += tag1[o];
        if(l == mid) ans1[vn[l]] += tag1[o] * ans2[vn[l]];
        if(mid + 1 == r) ans1[vn[r]] += tag1[o] * ans2[vn[r]];
        tag1[o] = 0;
    }
    if(tag2[o]) {
        tag2[lch] += tag2[o]; tag2[rch] += tag2[o];
        if(l == mid) ans1[vn[l]] += (siz[vn[l]] - siz[son[vn[l]]]) * tag2[o];
        if(mid + 1 == r) ans1[vn[r]] += (siz[vn[r]] - siz[son[vn[r]]]) * tag2[o];
        tag2[o] = 0;
    }
}

inline void modify(int o, int l, int r, int ll, int rr, LL v, int type) {
    if(l >= ll && r <= rr) {
        if(type == 1) {
            if(l == r) ans1[vn[l]] += v * ans2[vn[l]];
            else tag1[o] += v;
        } else if(type == 2) {
            if(l == r) ans1[vn[l]] += (siz[vn[l]] - siz[son[vn[l]]]) * v;
            else tag2[o] += v;
        }
        return;
    }
    pushdown(o, l, r);
    if(ll <= mid) modify(lch, l, mid, ll, rr, v, type);
    if(rr > mid) modify(rch, mid + 1, r, ll, rr, v, type);
}

inline void update(int o, int l, int r, int x) {
    if(l == r) return;
    pushdown(o, l, r);
    if(x <= mid) update(lch, l, mid, x);
    else update(rch, mid + 1, r, x);
}

inline void modify(int u, int _lu, LL v) {
    int tu = top[u], lu = _lu;
    while(u) {
        if(dfn[tu] != dfn[u]) modify(1, 1, n, dfn[tu], dfn[u] - 1, v, 2);
        ans1[u] += v * (siz[u] - siz[lu]);
        lu = u; u = fa[tu]; tu = top[u];
    }
}

int main() {
    n = readint(); m = readint();
    for(int i = 2, p; i <= n; i++) {
        p = readint();
        gra[i].push_back(p); gra[p].push_back(i);
    }
    for(int i = 1; i <= n; i++) {
        w[i] = readint();
    }
    dfs1(1);
    dfs2(1, 1);
    char op; int u, delta;
    while(m--) {
        op = readop(); u = readint();
        if(op != 'Q') delta = readint();
        if(op == 'S') {
            ans1[u] += 1ll * delta * (siz[u] - 1);
            modify(fa[u], u, delta);
        } else if(op == 'M') {
            modify(1, 1, n, dfn[u], dfn[u] + siz[u] - 1, delta * 2ll, 1);
            modify(fa[u], u, 1ll * delta * siz[u]);
        } else {
            update(1, 1, n, dfn[u]);
            printf("%.6lf\n", double(ans1[u]) / ans2[u]);
        }
    }
    return 0;
}