Codeforces Beta Round #99 (Div. 1) C Mushroom Gnomes - 2（单点查询）

最新推荐文章于 2022-07-13 02:01:45 发布

shiqi_614

最新推荐文章于 2022-07-13 02:01:45 发布

阅读量1.6k

点赞数

分类专栏：线段树及数据结构

本文链接：https://blog.csdn.net/shiqi_614/article/details/8738884

版权

线段树及数据结构专栏收录该内容

82 篇文章 1 订阅

订阅专栏

题意：有n棵树，m个蘑菇，每棵树有坐标a，高度h，向左边倒的概率，向右的概率（概率用0-100表示），向左倒范围[x-h,x)内的蘑菇被破坏，向右倒范围(x,x+h]范围内的蘑菇被破坏。每个蘑菇有坐标b，及它的魔力值z。

问树倒下后，所以蘑菇的魔力值的和的期望。

离散化之后，成段更新，即树倒下范围内的区间全部乘以倒向这边的概率。最后遍历每个蘑菇的存活的概率乘以它的魔力值就可以了。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
#include <map>
#include <algorithm>
using namespace std;
#define LL(x) (x<<1)
#define RR(x) (x<<1|1)
#define MID(a,b) (a+((b-a)>>1))
const int N=1e5+5;
struct OP
{
	double flag;
	int st,ed,lft;
	OP(){}
	OP(int a,int b,double c,int d){st=a;ed=b;flag=c;lft=d;}
};
struct node
{
	double P;
	int lft,rht,flag;
	int mid(){return MID(lft,rht);}
	void fun(double tmp)
	{
		P*=tmp;
		flag=1;
	}
};
int B[N],Z[N];
vector<OP> op;
vector<int> Y;
map<int,int> H;
struct Segtree
{
	node tree[N*4*4];
	void PushDown(int ind)
	{
		if(tree[ind].flag)
		{
			tree[LL(ind)].fun(tree[ind].P);
			tree[RR(ind)].fun(tree[ind].P);
			tree[ind].flag=0;	tree[ind].P=1;
		}
	}
	void build(int lft,int rht,int ind)
	{
		tree[ind].lft=lft;	tree[ind].rht=rht;
		tree[ind].P=1;		tree[ind].flag=0;
		if(lft!=rht)
		{
			int mid=tree[ind].mid();
			build(lft,mid,LL(ind));
			build(mid+1,rht,RR(ind));
		}
	}
	void updata(int st,int ed,int ind,double flag)
	{
		int lft=tree[ind].lft,rht=tree[ind].rht;
		if(st<=lft&&rht<=ed) tree[ind].fun(flag);
		else
		{
			PushDown(ind);
			int mid=tree[ind].mid();
			if(st<=mid) updata(st,ed,LL(ind),flag);
			if(ed> mid) updata(st,ed,RR(ind),flag);
		}
	}
	double query(int pos,int ind)
	{
		if(tree[ind].lft==tree[ind].rht) return tree[ind].P;
		else
		{
			PushDown(ind);
			int mid=tree[ind].mid();
			if(pos<=mid) return query(pos,LL(ind));
			else return query(pos,RR(ind));
		}
	}
}seg;
int main()
{
	int n,m;
	while(scanf("%d%d",&n,&m)!=EOF)
	{
		Y.clear(); H.clear(); op.clear();

		for(int i=0;i<n;i++)
		{
			int a,b,c,d; scanf("%d%d%d%d",&a,&b,&c,&d);
			Y.push_back(a-b); Y.push_back(a+b);
			Y.push_back(a);
			op.push_back(OP(a-b,a,(100-c)/100.0,1));
			op.push_back(OP(a,a+b,(100-d)/100.0,0));
		}
		for(int i=0;i<m;i++)
		{
			scanf("%d%d",&B[i],&Z[i]);
			Y.push_back(B[i]);
		}
		sort(Y.begin(),Y.end());
		Y.erase(unique(Y.begin(),Y.end()),Y.end());
		for(int i=0;i<(int)Y.size();i++) H[Y[i]]=i;

		seg.build(0,(int)Y.size(),1);
		for(int i=0;i<(int)op.size();i++)
		{
			if(op[i].lft) seg.updata(H[op[i].st],H[op[i].ed]-1,1,op[i].flag);
			else seg.updata(H[op[i].st]+1,H[op[i].ed],1,op[i].flag);
		}
		double res=0;
		for(int i=0;i<m;i++)
		{
			res+=Z[i]*seg.query(H[B[i]],1);
		}
		printf("%.4lf\n",res);
	}
	return 0;
}