【线段树 泰勒展开】Codechef April Challenge 2018 Chef at the Food Fair
第一次写泰勒展开;本地和CC差距好大
题目大意
大厨住的城市里办了一场美食节。一条街上开设了$N$个摊位,编号为$1∼N$。这天开始时,第$i$个摊位的食物会导致食物中毒的概率是$P_i$。在这一天中,大厨发现某些摊位可能会根据顾客的反馈提供没那么有毒的食物。你需要处理$Q$个询问,询问有以下两类:
0 L R:求出:如果要吃遍$[L,R]$内所有摊位的食物,那么不会食物中毒的概率是多少;
1 L R T:$[L,R]$中的所有摊位的食物会导致食物中毒的概率变为了原来的$T$倍。$T$是一个小于$1$的非负实数。
对于前 $20\%$ 的数据,$n,m\le 2000$
另有 $20\%$ 的数据,$T\le 0.5$
对于 $100\%$ 的数据,满足 $n,m\le 10^5,0\le T<1,P_i\le 1,1\le L\le R\le n$,保证输入数据不超过 $6$ 位小数。
题目分析
注意到维护的操作有些不同寻常。
- 第一:维护的是$\prod (1-P_i)$
- 第二:每次操作是区间乘法
对于要求支持区间乘的问题,有一种转化套路是将它取$\ln$,那么问题就变成了维护区间和。
那么这题中还需要处理$\ln (1-P_i)$,将它泰勒展开得到$\ln(1-x)=x-{1\over 2}x^2-{1\over 3}x^3-...-{1\over n}x^n+R(x)$。我们一如既往地爆精度,只需要保留这个式子的前$MAXD=100$项和。维护时则是开$MAXD$颗线段树对每类次项分别处理区间乘法。
需要注意的是,这题需要一些常数技巧。我最先是开了$f[MAXD]$颗封装好的线段树,但由于数组的第一维是更频繁访问的一维,所以实际运行效率会比较低。如果采用形如$f[maxn<<2][MAXD]$的做法,就会快非常多。
非常迷的一点是,同一份代码在本地考试的数据下,极限数据要跑个4~5s;但是交到CC上就rank2了……
1 #include<bits/stdc++.h> 2 const int maxn = 100035; 3 const double eps = 1e-12; 4 5 int n,m,MAXD; 6 double p[maxn],w[maxn],K,S; 7 double f[maxn<<2][103],tag[maxn<<2]; 8 9 int read() 10 { 11 char ch = getchar(); 12 int num = 0, fl = 1; 13 for (; !isdigit(ch); ch=getchar()) 14 if (ch=='-') fl = -1; 15 for (; isdigit(ch); ch=getchar()) 16 num = (num<<1)+(num<<3)+ch-48; 17 return num*fl; 18 } 19 void pushup(int rt) 20 { 21 for (int i=1; i<=MAXD; i++) 22 f[rt][i] = f[rt<<1][i]+f[rt<<1|1][i]; 23 } 24 void pushdown(int rt) 25 { 26 double v = tag[rt], s = tag[rt]; 27 if (fabs(1.0-v) > eps){ 28 tag[rt<<1] *= v, tag[rt<<1|1] *= v; 29 for (int i=1; i<=MAXD; i++) 30 f[rt<<1][i] *= s, f[rt<<1|1][i] *= s, s *= v; 31 tag[rt] = 1.0; 32 } 33 } 34 void build(int rt, int l, int r) 35 { 36 tag[rt] = 1.0; 37 if (l==r){ 38 for (int i=1; i<=MAXD; i++) 39 w[l] *= p[l], f[rt][i] = w[l]/i; 40 return; 41 } 42 int mid = (l+r)>>1; 43 build(rt<<1, l, mid); 44 build(rt<<1|1, mid+1, r); 45 pushup(rt); 46 } 47 double query(int rt, int L, int R, int l, int r) 48 { 49 if (L <= l&&r <= R){ 50 double ret = 0; 51 for (int i=1; i<=MAXD; i++) 52 ret += f[rt][i]; 53 return ret; 54 } 55 int mid = (l+r)>>1; 56 double ret = 0; 57 pushdown(rt); 58 if (L <= mid) ret += query(rt<<1, L, R, l, mid); 59 if (R > mid) ret += query(rt<<1|1, L, R, mid+1, r); 60 return ret; 61 } 62 void modify(int rt, int L, int R, int l, int r, double c) 63 { 64 if (L <= l&&r <= R){ 65 double s = c; 66 tag[rt] *= s; 67 for (int i=1; i<=MAXD; i++) 68 f[rt][i] *= s, s *= c; 69 return; 70 } 71 int mid = (l+r)>>1; 72 pushdown(rt); 73 if (L <= mid) modify(rt<<1, L, R, l, mid, c); 74 if (R > mid) modify(rt<<1|1, L, R, mid+1, r, c); 75 pushup(rt); 76 } 77 int main() 78 { 79 n = read(), m = read(); 80 for (int i=1; i<=n; i++) scanf("%lf",&p[i]), w[i] = 1.0; 81 MAXD = 100; 82 build(1, 1, n); 83 for (int i=1; i<=m; i++) 84 { 85 int opt = read(), l = read(), r = read(); 86 if (opt){ 87 scanf("%lf",&K); 88 modify(1, l, r, 1, n, K); 89 }else{ 90 K = query(1, l, r, 1, n); 91 printf("%.8lf\n",exp(-K)); 92 } 93 } 94 return 0; 95 }
END
转载于:https://www.cnblogs.com/antiquality/p/10385631.html
【线段树 泰勒展开】Codechef April Challenge 2018 Chef at the Food Fair相关推荐
- [题解]CodeChef APRIL Challenge 17
Similar Dishes 题意简述 给你两个大小为 4 4的字符串集合,问你并集是不是空集. 数据范围 1≤T≤2001 \leq T \leq 200 2≤len≤10 2 \leq len \ ...
- Codechef August Challenge 2018 : Coordinate Compression
传送门 外边二分,里面拿线段树维护贪心就行了. #include<cstdio> #include<vector> #include<cstring> #inclu ...
- CFCC百套计划2 CodeChef December Challenge 2017 Chef And Easy Xor Queries
https://www.codechef.com/DEC17/problems/CHEFEXQ 题意: 位置i的数改为k 询问区间[1,i]内有多少个前缀的异或和为k 分块 sum[i][j] 表示第 ...
- Codechef January Challenge 2018 - Killjee and k-th letter
题意: 给出一个的串 s,将 s 所有子串按照字典序排列好相接起来形成一个新串q次询问,每一次询问问新串中的第 k 个字符是什么,强制在线. $|s|,q \le 2*10^{5} $ 跟所有子串有关 ...
- Codechef July Challenge 2018 : Subway Ride
传送门 首先(想了很久之后)注意到一个性质:同一条边有多种颜色的话保留3种就可以了,这是因为假如最优解要求当前位置与相邻两条边都不相同,那么只要有3条边,就肯定可以满足这一点. 完事就做一个nlogn ...
- ACM-ICPC 2018 南京赛区网络预赛 Lpl and Energy-saving Lamps 线段树
目录 ACM-ICPC 2018 南京赛区网络预赛 Lpl and Energy-saving Lamps 线段树 题面 题意 思路 ACM-ICPC 2018 南京赛区网络预赛 Lpl and En ...
- BZOJ 3221: [Codechef FEB13] Obserbing the tree树上询问( 可持久化线段树 + 树链剖分 )
树链剖分+可持久化线段树....这个一眼可以看出来, 因为可持久化所以写了标记永久化(否则就是区间修改的线段树的持久化..不会), 结果就写挂了, T得飞起...和管理员拿数据调后才发现= = 做法: ...
- 2018.3.15校内互测总结-点分治-线段树
这是曾来过咱们学校集训的一位大神出的~ T1 题目大意 给出一棵带边权的无根树,求树上前$k$大的路径的长度. $1 \leq n \leq 200000$ 题解 想了一上午点分治,却发现只会$O(n ...
- BZOJ.5249.[九省联考2018]iiidx(贪心 线段树)
BZOJ LOJ 洛谷 \(d_i\)不同就不用说了,建出树来\(DFS\)一遍. 对于\(d_i\)不同的情况: Solution 1: xxy tql! 考虑如何把这些数依次填到树里. 首先对于已 ...
最新文章
- 安装openproject
- AcWing算法提高课 Level-3 第三章 图论
- Vector的一种实现(一)
- 前端学习(975):bootstrap轮播图
- 巨头垂涎却不能染指,loT 数据库风口已至
- Apache Flink 的迁移之路,2 年处理效果提升 5 倍
- Centos 8 阿里yum源配置
- 6月29日Java实训第一天
- vue cl3、vuex、vue-router、ant design vue、axios搭建一个简易的单页面应用
- linux 基础林海峰,word数学公式编辑器-官方版合集下载-多特
- Unity3D《打地鼠》学习笔记及心得
- 研究方向之智慧仓储与物流及相关算法
- 新世界交互英语视听说下载
- Go语言Web开发GoFrame框架搭建
- 计算机毕业设计jsp旅游网站
- 经典智力题:拿球问题
- 编译官方DV300的SDKv2.0
- win7一直显示正在关机_win7电脑如何取消自动关机
- 阿拇云手机群控怎么做qq营销引流
- 宁花4000买手机 不花6元买游戏
热门文章
- 页目录项和页表项——《x86汇编语言:从实模式到保护模式》读书笔记43
- JVM_03 运行时数据区 [ 方法区]
- Hi3516A开发--/etc/passwd
- Nand Flash基础知识与坏块管理机制的研究
- 汇编语言EQU伪指令
- Vitalik Buterin:Casper 权益证明与分片技术最新进展
- VirtualAPK:滴滴 Android 插件化的实践之路
- python重复字符串n次_python装饰器听了N次也没印象,读完这篇你就懂了
- 微信为什么不进入鸿蒙生态,马化腾为何迟迟不将微信, 加入任正非的华为鸿蒙生态系统中呢?...
- 仪表自定义刻度值_史上最全的仪表选型,就差你还没看了