codeforces 712E

题目大意

有$n$个赌场，你在$i$赌场时，有$p_i$的概率走到$i+1$,有$1-p_i$的概率走到$i-1$.保证任何时候$p_i\leq p_{i+1}$

有$q$次操作，修改一个赌场的$p$值；或者询问$[l,r]$表示从第$l$个赌场走到$r$的概率,他在走的过程中不会离开区间$[l,r]$.

解题思路

设$f_i$表示从$i$能走到$n$的概率。
易得$f_i=f_{i-1}*(1-p_{i})+f_{i+1}*p_i\Rightarrow f_i-f_{i-1}=p_i*(f_{i+1}-f_{i-1})$.

令$g_i=f_i-f_{i-1}$
得$g_i=p_i*(g_i+g_{i+1})\Rightarrow g_{i+1}=g_i*\frac{1-p_i}{p_i}$

令$t_i=\frac{1-p_i}{p_i}$
则$g_{i+1}=g_i*t_i$

易得$\sum_{i=1}^n g_i=1$且$f_1=g_1$.

$$\therefore g_1+g_1*t_1+g_1*t_1*t_2+……+g_1*t_1*t_1*t_2*t_3*……*t_{n-1}=1$$
$$g_1*(1+t_1+t_1*t_2+……+t_1*t_1*t_2*t_3*……*t_{n-1})=1$$

只需维护
$(1+t_1+t_1*t_2+……+t_1*t_1*t_2*t_3*……*t_{n-1})$
即可。
这个可以使用线段树简单实现。

参考代码

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define maxn 100005
#define db double
using namespace std;

struct seg{
    db a,b;
}t[maxn * 4],now;

int n,q;

db p[maxn];

void build(int v,int l,int r){
    if (l==r) {
        t[v].a=t[v].b=(1-p[l])/p[l];
        return;
    }
    int mid=(l+r) >> 1;
    build(v << 1,l,mid);
    build(v << 1 | 1,mid+1,r);
    t[v].a=t[v << 1].a*t[v << 1 | 1].a;
    t[v].b=t[v << 1].b+t[v << 1 | 1].b*t[v << 1].a;
}

void change(int v,int l,int r,int x){
    if (l==r) {
        t[v].a=t[v].b=(1-p[l])/p[l];
        return;
    }
    int mid=(l+r) >> 1;
    if (x<=mid) change(v << 1,l,mid,x);
    else change(v << 1 | 1,mid+1,r,x);
    t[v].a=t[v << 1].a*t[v << 1 | 1].a;
    t[v].b=t[v << 1].b+t[v << 1 | 1].b*t[v << 1].a;
}

void query(int v,int l,int r,int x,int y){
    if (l==x && r==y) {
        now.b=now.b+now.a*t[v].b;
        now.a=now.a*t[v].a;
        return;
    }
    int mid=(l+r) >> 1;
    if (y<=mid) query(v << 1,l,mid,x,y);
    else if (x>mid) query(v << 1 | 1,mid+1,r,x,y);
    else query(v << 1,l,mid,x,mid),query(v << 1 | 1,mid+1,r,mid+1,y);
}

int main(){
    scanf("%d%d",&n,&q);
    fo(i,1,n) {
        int x,y;
        scanf("%d%d",&x,&y);
        p[i]=1.0* x / y;
    }
    build(1,1,n);
    while (q--) {
        int Type,x,y,z;
        scanf("%d",&Type);
        if (Type==1) {
            scanf("%d%d%d",&x,&y,&z);
            p[x]=1.0 * y / z;
            change(1,1,n,x);
        }
        else {
            scanf("%d%d",&x,&y);
            db ans;
            now.a=1;
            now.b=0;
            query(1,1,n,x,y);
            if (now.b<=1e15) ans=1/(1+now.b);
            else ans=0;
            printf("%.10lf\n",ans);
        }
    }
    return 0;
}