coderforces 138 CMushroom Gnomes - 2 선분 트 리

클릭 하여 링크 열기
제목: n 개의 나무 와 m 의 사람 을 정 하고 그 다음 n 줄 은 모든 나무의 위치, 높이, 왼쪽으로 넘 어 질 확률, 오른쪽으로 넘 어 질 확률 을 정 합 니 다. 그 다음 에 m 의 개인 적 인 위 치 를 제시 하여 모든 사람 이 살아 남 을 확률 을 판단 하고 Zi 를 곱 한 다음 에 더 합 니 다.
사고방식: 선분 트 리 의 노드 유지 확률 을 사용 하여 먼저 이산 화 를 한다. 원래 의 경우 성단 갱신 에 사용 해 야 할 수도 있 지만 이 문 제 는 피 할 수 있다. 다음은 주석 을 제시 할 것 이다.

#include <map>
#include <vector>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int maxn=100010;
double num[2000010];//            2 
vector<int> G;
map<int,int> pp;
void buildtree(int le,int ri,int node){
    num[node]=1;
    if(le==ri) return ;
    int t=(le+ri)>>1;
    buildtree(le,t,node<<1);
    buildtree(t+1,ri,node<<1|1);
}
void update(int l,int r,double add,int le,int ri,int node){
    if(l<=le&&ri<=r){
        num[node]*=add;
        return ;
    }
    int t=(le+ri)>>1;
    if(l<=t) update(l,r,add,le,t,node<<1);
    if(r>t) update(l,r,add,t+1,ri,node<<1|1);
}
double query(int pos,int le,int ri,int node){
    if(le==ri) return num[node];
    int t=(le+ri)>>1;
    if(pos<=t) return num[node]*query(pos,le,t,node<<1);//            ，       ，          
    else return num[node]*query(pos,t+1,ri,node<<1|1);
}
int main(){
    int L[maxn],R[maxn],X[maxn],H[maxn],P[10010],n,m,B[10010];
    while(scanf("%d%d",&n,&m)!=-1){
        G.clear();
        for(int i=0;i<n;i++){
            scanf("%d%d%d%d",&X[i],&H[i],&L[i],&R[i]);
            G.push_back(X[i]-H[i]);
            G.push_back(X[i]-1);//       ，        X[i]   ，    1；
            G.push_back(X[i]+H[i]);
            G.push_back(X[i]+1);//    ，   1；
        }
        for(int i=0;i<m;i++){
            scanf("%d%d",&B[i],&P[i]);
            G.push_back(B[i]);
        }
        sort(G.begin(),G.end());
        G.resize(unique(G.begin(),G.end())-G.begin());//     
        int len=G.size();
        buildtree(1,len,1);
        for(int i=0;i<len;i++){
            pp[G[i]]=i+1;
        }
        for(int i=0;i<n;i++){
            update(pp[X[i]-H[i]],pp[X[i]-1],1.0-double((double)L[i]/100),1,len,1);
            update(pp[X[i]+1],pp[X[i]+H[i]],1.0-double((double)R[i]/100),1,len,1);
        }
        double ans=0;
        for(int i=0;i<m;i++){
            ans+=query(pp[B[i]],1,len,1)*P[i];
        }
        printf("%.10f
",ans);
    }
    return 0;
}