[NOI2013] 树的计数

#include<bits/stdc++.h>
const int N = 2e5 + 10;
int ri() {
    char c = getchar(); int x = 0, f = 1; for(;c < '0' || c > '9'; c = getchar()) if(c == '-') f = -1;
    for(;c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) - '0' + c; return x * f;
}
struct Data {
    double x; int e;
}g[N], f[N];
Data operator + (Data a, Data b) {
    if(a.e < b.e)
        std::swap(a, b);
    int d = a.e - b.e; double l = b.x;
    if(d > 100) l = 0;
    else {
        for(;d--;)
            l /= 10;
    }
    a.x += l;
    if(fabs(a.x) > 10) a.x /= 10, ++a.e;
    return a;
}
void operator += (Data &a, Data b) {a = a + b;}
void operator -= (Data &a, Data b) {
    b.x = -b.x;
    a = a + b;
}
int a[N], b[N], psa[N], psb[N], n, must;
int main() {
    n = ri();
    for(int i = 1;i <= n; ++i)
        a[i] = ri(), psa[a[i]] = i;
    for(int i = 1;i <= n; ++i)
        b[i] = ri(), psb[b[i]] = i;
    int l = 0, r = 0;
    f[1].x = g[1].x = 1;
    f[2].x = g[2].x = -1;
    for(int j = 1, i;j <= n; ++j) {
        f[j] += f[j - 1]; g[j] += g[j - 1];
        if(j == n) break;
        l = std::max(l, psb[a[psa[b[j]] + 1]]);
        for(r = std::max(r, j + 2);r <= n; ++r)
            if(psa[b[r]] < psa[b[r - 1]])
                break;
        --r;
        if(l <= r) {
            f[l] += f[j]; f[r + 1] -= f[j];
            g[l] += f[j] + g[j]; g[r + 1] -= f[j] + g[j];
        }
    }
    int d = g[n].e - f[n].e; double x = g[n].x / f[n].x;
    for(;d--;)
        x *= 10;
    printf("%.3lf\n", x);
    return 0;
}

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值