2019.02.14 codechef Chef at the Food Fair(线段树+泰勒展开)
传送门
题意:现在有nnn个位置,每个位置上有一个值aia_iai.
要求支持如下两种操作:
- 区间乘vvv
- 求区间的(1−ai)(1-a_i)(1−ai)之积
思路:
考虑转换式子:
Ans=∏i=lr(1−ai)=e∑i=lrln(1−ai)Ans=\prod_{i=l}^r(1-a_i)=e^{\sum_{i=l}^rln(1-a_i)}Ans=∏i=lr(1−ai)=e∑i=lrln(1−ai)
于是只需维护∑i=lrln(1−ai)\sum_{i=l}^rln(1-a_i)∑i=lrln(1−ai)
把这个式子在x0=0x_0=0x0=0出泰勒展开可以得到:ln(1−x)=−x−12x2−13x3−⋯˙ln(1-x)=-x-\frac 12x^2-\frac 13x^3-\dot\cdotsln(1−x)=−x−21x2−31x3−⋯˙
因为允许一定的精度误差。
于是我们维护这个多项式的前若干项即可。
代码:
#include<bits/stdc++.h>
#define ri register int
using namespace std;
inline int read(){
int ans=0;
char ch=getchar();
while(!isdigit(ch))ch=getchar();
while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
return ans;
}
const int N=1e5+5;
double a[N];
int n,m;
namespace SGT{
#define lc (p<<1)
#define rc (p<<1|1)
#define mid (T[p].l+T[p].r>>1)
struct Val{
double a[100];
inline double&operator[](const int&k){return a[k];}
inline const double&operator[](const int&k)const{return a[k];}
};
struct Tag{double a;};
const Tag tag_empty=(Tag){1};
Val val_empty;
inline Val operator+(const Val&a,const Val&b){
Val ret;
for(ri i=0;i<100;++i)ret[i]=a[i]+b[i];
return ret;
}
inline void operator+=(Val&a,const Val&b){a=a+b;}
inline Tag operator+(const Tag&a,const Tag&b){return (Tag){a.a*b.a};}
inline void operator+=(Tag&a,const Tag&b){a=a+b;}
inline Val operator+(const Val&a,const Tag&b){
Val ret;
double mult=b.a;
for(ri i=0;i<100;++i)ret[i]=a[i]*mult,mult*=b.a;
return ret;
}
inline void operator+=(Val&a,const Tag&b){a=a+b;}
inline Val solve1(double a){
Val ret=val_empty;
double mult=a;
for(ri i=0;i<100;++i)ret[i]=-mult/(i+1),mult*=a;
return ret;
}
inline double solve2(Val a){
double ret=0.0;
for(ri i=0;i<100;++i)ret+=a[i];
return ret;
}
struct Node{Val val;Tag tag;int l,r;}T[N<<2];
inline void pushup(int p){T[p].val=T[lc].val+T[rc].val;}
inline void pushnow(int p,Tag v){T[p].val+=v,T[p].tag+=v;}
inline bool check(Tag a){return fabs(a.a-1)<=1e-9;}
inline void pushdown(int p){
if(check(T[p].tag))return;
pushnow(lc,T[p].tag),pushnow(rc,T[p].tag);
T[p].tag=tag_empty;
}
inline void build(int p,int l,int r){
T[p]=(Node){val_empty,tag_empty,l,r};
if(l==r){T[p].val=solve1(a[l]);return;}
build(lc,l,mid),build(rc,mid+1,r),pushup(p);
}
inline void update(int p,int ql,int qr,Tag v){
if(ql<=T[p].l&&T[p].r<=qr)return pushnow(p,v);
pushdown(p);
if(ql<=mid)update(lc,ql,qr,v);
if(qr>mid)update(rc,ql,qr,v);
pushup(p);
}
inline Val query(int p,int ql,int qr){
if(ql<=T[p].l&&T[p].r<=qr)return T[p].val;
pushdown(p);
Val ret=val_empty;
if(ql<=mid)ret+=query(lc,ql,qr);
if(qr>mid)ret+=query(rc,ql,qr);
return ret;
}
#undef lc
#undef rc
#undef mid
}
int main(){
for(ri i=0;i<100;++i)SGT::val_empty[i]=0;
n=read(),m=read();
for(ri i=1;i<=n;++i)scanf("%lf",&a[i]);
SGT::build(1,1,n);
double x;
for(ri l,r,op;m;--m){
op=read(),l=read(),r=read();
if(!op)printf("%.8lf\n",exp(SGT::solve2(SGT::query(1,l,r))));
else scanf("%lf",&x),SGT::update(1,l,r,(SGT::Tag){x});
}
return 0;
}
2019.02.14 codechef Chef at the Food Fair(线段树+泰勒展开)的更多相关文章
- 2019.02.28 bzoj4199: [Noi2015]品酒大会(sam+线段树)
传送门 题意:给一个串,每个位置有一个权值,当S[s...s+len−1]=S[t...t+len−1]&&S[s...s+len]̸=S[t..t+len]S[s...s+len-1 ...
- 2022.02.27 CF811E Vladik and Entertaining Flags(线段树+并查集)
2022.02.27 CF811E Vladik and Entertaining Flags(线段树+并查集) https://www.luogu.com.cn/problem/CF811E Ste ...
- 2019.02.15 codechef Favourite Numbers(二分+数位dp+ac自动机)
传送门 题意: 给444个整数L,R,K,nL,R,K,nL,R,K,n,和nnn个数字串,L,R,K,数字串大小≤1e18,n≤65L,R,K,数字串大小\le1e18,n\le65L,R,K,数字 ...
- 【线段树 泰勒展开】Codechef April Challenge 2018 Chef at the Food Fair
第一次写泰勒展开:本地和CC差距好大 题目大意 大厨住的城市里办了一场美食节.一条街上开设了$N$个摊位,编号为$1∼N$.这天开始时,第$i$个摊位的食物会导致食物中毒的概率是$P_i$.在这一天中 ...
- CodeChef DISTNUM2 Easy Queries 节点数组线段树
Description You are given an array A consisting of N positive integers. You have to answer Q queries ...
- UVAlive7141 BombX 14年上海区域赛D题 线段树+离散化
题意:一个无限大的棋盘, 有n个小兵, 给出了n个兵的坐标, 现在有一个长为width 高为height的炸弹放在棋盘上, 炸弹只能上下左右平移, 不能旋转. 且放炸弹的区域不能含有士兵, 炸弹可以一 ...
- 2019.01.22 bzoj3333: 排队计划(逆序对+线段树)
传送门 题意简述:给出一个序列,支持把ppp~nnn中所有小于等于apa_pap的'扯出来排序之后再放回去,要求动态维护全局逆序对. 思路:我们令fif_ifi表示第iii个位置之后比它大的数的个 ...
- 2019.01.13 loj#6515. 贪玩蓝月(线段树分治+01背包)
传送门 题意简述:有一个初始为空的双端队列,每次可以在队首和队尾插入或弹出一个二元组(wi,vi)(w_i,v_i)(wi,vi),支持询问从当前队列中选取若干个元素是的他们的和对 MODMODM ...
- 2019.01.13 bzoj4137: [FJOI2015]火星商店问题(线段树分治+可持久化01trie)
传送门 题意:序列上有nnn个商店,有两种事件会发生: sss商店上进购标价为vvv的一个物品 求编号为[l,r][l,r][l,r]之间的位置买ddd天内新进购的所有物品与一个数xxx异或值的最大值 ...
随机推荐
- MySQL 之 MHA + ProxySQL + keepalived 实现读写分离,高可用(一)
准备服务器: docker network create --subnet=192.168.0.0/16 staticnetdocker run -d --privileged -v `pwd`/my ...
- Linux下安装jdk8步骤详述
作为Java开发人员,在Linux下安装一些开发工具是必备技能,本文以安装jdk为例,详细记录了每一步的操作命令,以供参考. 0.下载jdk8 登录网址:http://www.oracle.com/t ...
- 多线程之sleep和wait的区别
它们最大本质的区别是:sleep()不释放同步锁,wait()释放同步锁. 还有用法的上的不同是:sleep(milliseconds)可以用时间指定来使他自动醒过来,如果时间不到你只能调用inter ...
- How can I manage the modules for python2 when python3 installed as well. In OSX
ref: https://stackoverflow.com/questions/53385448/how-can-i-manage-the-modules-for-python2-when-pyth ...
- java——IO流01
移动文件有一种简单方法,不需要复制文件再删除文件. package com.unir.test01; import java.io.File; import java.io.IOException; ...
- 本学期c#学习总结
本学期c#学习总结 时间转瞬即逝,大一上半学期的学习生涯已经结束.虽然以前我没什么关于学习计算机的基础,但是经过了这几个月的学习我也还是有点收获的. 我发现c#语言的关键词有很多语言特性和固定的用法, ...
- Xshell 6安装与使用教程
随着xshell5出现评估期已过的问题,发现好多人不知道怎么下载免费版的Xshell,在这里我将详细告诉大家如何下载和安装最新的Xshell6远程管理工具. Xshell安装 1.进入xshell英文 ...
- Communication Model
[Communication Model] EOSIO actions operate primarily in a message-based communication architecture. ...
- 解题(TakeBusChooseLine)
题目描述 从小明家所在公交站出发有n路公交到公司,现给出每路公交的停站数(不包括起点和终点),及每次停的时间(一路车在每个站停的时间相同)和发车的间隔,先假定每辆车同时在相对时间0分开始发车,且所有车 ...
- jquery即时获取上传文件input file文件名
截图: 代码: <input type="file" id="choosefile" style="display:none"/& ...