BZOJ5020: [THUWC 2017]在美妙的数学王国中畅游(LCT,泰勒展开,二项式定理)
Description
Input
Output
Sample Input
1 1 0
3 0.5 0.5
3 -0.5 0.7
appear 0 1
travel 0 1 0.3
appear 0 2
travel 1 2 0.5
disappear 0 1
appear 1 2
travel 1 2 0.5
Sample Output
1.67942554e+000
1.20000000e+000
解题思路:
题目描述如此毒瘤。
从操作3得到的启发,将多项式展开对应项相加。
这道题可以将sin(ax+b),eax+b泰勒展开。
精度的话16位肯定够。剩下的就是裸的LCT了。
听霉霉的歌写泰勒展开不容易错
代码:
1 #include<cstdio> 2 #include<cstring> 3 #include<algorithm> 4 #define lll tr[spc].ch[0] 5 #define rrr tr[spc].ch[1] 6 #define ls ch[0] 7 #define rs ch[1] 8 const int N=100010; 9 const int oo=16; 10 struct trnt{ 11 int ch[2]; 12 int fa; 13 int lzt; 14 int type; 15 bool anc; 16 double a,b; 17 double C[oo]; 18 double f[oo]; 19 double val(double x) 20 { 21 double ans=f[0]; 22 double t=x; 23 for(int i=1;i<oo;i++,t*=x) 24 ans+=f[i]*t; 25 return ans; 26 } 27 void Insert(void) 28 { 29 scanf("%d",&type); 30 scanf("%lf%lf",&a,&b); 31 return ; 32 } 33 void Taylor(double *fac) 34 { 35 double at[oo],bt[oo]; 36 for(int i=0;i<oo;i++) 37 C[i]=at[i]=bt[i]=0; 38 at[0]=1; 39 bt[0]=1; 40 for(int i=1;i<oo;i++) 41 at[i]=at[i-1]*a,bt[i]=bt[i-1]*b; 42 if(type==1) 43 {//sin(ax+b) 44 double tmp=1; 45 for(int i=1;i<oo;i+=2) 46 { 47 for(int j=0;j<=i;j++) 48 C[j]+=tmp*at[j]*bt[i-j]/fac[j]/fac[i-j]; 49 tmp*=-1.00; 50 } 51 return ; 52 } 53 if(type==2) 54 {//e^(ax+b) 55 for(int i=0;i<oo;i++) 56 { 57 for(int j=0;j<=i;j++) 58 C[j]+=fac[i]/fac[j]/fac[i-j]*at[j]*bt[i-j]/fac[i]; 59 } 60 return ; 61 } 62 if(type==3) 63 { 64 C[0]=b; 65 C[1]=a; 66 return ; 67 } 68 } 69 }tr[N]; 70 int n,m; 71 double fac[50]; 72 char tmp[10000]; 73 bool whc(int spc) 74 { 75 return tr[tr[spc].fa].rs==spc; 76 } 77 void pushup(int spc) 78 { 79 for(int i=0;i<oo;i++) 80 tr[spc].f[i]=tr[spc].C[i]; 81 if(lll) 82 for(int i=0;i<oo;i++) 83 tr[spc].f[i]+=tr[lll].f[i]; 84 if(rrr) 85 for(int i=0;i<oo;i++) 86 tr[spc].f[i]+=tr[rrr].f[i]; 87 return ; 88 } 89 void trr(int spc) 90 { 91 if(!spc) 92 return ; 93 std::swap(lll,rrr); 94 tr[spc].lzt^=1; 95 return ; 96 } 97 void pushdown(int spc) 98 { 99 if(tr[spc].lzt) 100 { 101 trr(lll); 102 trr(rrr); 103 tr[spc].lzt=0; 104 } 105 return ; 106 } 107 void recal(int spc) 108 { 109 if(!tr[spc].anc) 110 recal(tr[spc].fa); 111 pushdown(spc); 112 return ; 113 } 114 void rotate(int spc) 115 { 116 int f=tr[spc].fa; 117 bool k=whc(spc); 118 tr[f].ch[k]=tr[spc].ch[!k]; 119 tr[spc].ch[!k]=f; 120 if(tr[f].anc) 121 { 122 tr[f].anc=0; 123 tr[spc].anc=1; 124 }else 125 tr[tr[f].fa].ch[whc(f)]=spc; 126 tr[spc].fa=tr[f].fa; 127 tr[f].fa=spc; 128 tr[tr[f].ch[k]].fa=f; 129 pushup(f); 130 pushup(spc); 131 return ; 132 } 133 void splay(int spc) 134 { 135 recal(spc); 136 while(!tr[spc].anc) 137 { 138 int f=tr[spc].fa; 139 if(tr[f].anc) 140 { 141 rotate(spc); 142 return ; 143 } 144 if(whc(spc)^whc(f)) 145 rotate(spc); 146 else 147 rotate(f); 148 rotate(spc); 149 } 150 return ; 151 } 152 void access(int spc) 153 { 154 int lst=0; 155 while(spc) 156 { 157 splay(spc); 158 tr[rrr].anc=1; 159 tr[lst].anc=0; 160 rrr=lst; 161 lst=spc; 162 pushup(spc); 163 spc=tr[spc].fa; 164 } 165 return ; 166 } 167 void Mtr(int spc) 168 { 169 access(spc); 170 splay(spc); 171 trr(spc); 172 return ; 173 } 174 void split(int x,int y) 175 { 176 Mtr(x); 177 access(y); 178 splay(y); 179 return ; 180 } 181 void link(int x,int y) 182 { 183 Mtr(x); 184 tr[x].fa=y; 185 return ; 186 } 187 bool together(int x,int y) 188 { 189 split(x,y); 190 while(tr[y].ls) 191 y=tr[y].ls; 192 return x==y; 193 } 194 void cut(int x,int y) 195 { 196 split(x,y); 197 tr[x].fa=0; 198 tr[x].anc=true; 199 tr[y].ls=0; 200 pushup(y); 201 return ; 202 } 203 int main() 204 { 205 scanf("%d%d",&n,&m); 206 scanf("%s",tmp); 207 fac[0]=1; 208 for(int i=1;i<oo;i++) 209 { 210 double x=i; 211 fac[i]=fac[i-1]*x; 212 } 213 for(int i=1;i<=n;i++) 214 { 215 tr[i].Insert(); 216 tr[i].Taylor(fac); 217 tr[i].anc=1; 218 } 219 while(m--) 220 { 221 scanf("%s",tmp+1); 222 if(tmp[1]=='a') 223 { 224 int a,b; 225 scanf("%d%d",&a,&b); 226 a++,b++; 227 link(a,b); 228 }else if(tmp[1]=='d') 229 { 230 int a,b; 231 scanf("%d%d",&a,&b); 232 a++,b++; 233 cut(a,b); 234 }else if(tmp[1]=='m') 235 { 236 int x; 237 scanf("%d",&x); 238 x++; 239 splay(x); 240 tr[x].Insert(); 241 tr[x].Taylor(fac); 242 }else{ 243 int a,b; 244 scanf("%d%d",&a,&b); 245 double x; 246 scanf("%lf",&x); 247 a++,b++; 248 if(!together(a,b)) 249 puts("unreachable"); 250 else{ 251 double ret=tr[b].val(x); 252 printf("%.8e\n",ret); 253 } 254 } 255 } 256 return 0; 257 }
转载于:https://www.cnblogs.com/blog-Dr-J/p/10116746.html
BZOJ5020: [THUWC 2017]在美妙的数学王国中畅游(LCT,泰勒展开,二项式定理)相关推荐
- BZOJ5020: [THUWC 2017]在美妙的数学王国中畅游
BZOJ5020: [THUWC 2017]在美妙的数学王国中畅游 其实题面好像有点不全,建议去洛谷: P4546 [THUWC2017]在美妙的数学王国中畅游 这里还是$BZOJ$的题面. Desc ...
- [BZOJ5020] [THUWC 2017]在美妙的数学王国中畅游
Description 数字和数学规律主宰着这个世界. 机器的运转, 生命的消长, 宇宙的进程, 这些神秘而又美妙的过程无不可以用数学的语言展现出来. 这印证了一句古老的名言: "学好数理化 ...
- [bzoj5020][THUWC 2017]在美妙的数学王国中畅游
Orz w_yqts lct+泰勒展开 #include <bits/stdc++.h> using namespace std; #define db double #define N ...
- bzoj5020[THUWC 2017]在美妙的数学王国中畅游(LCT)
传送门:https://loj.ac/problem/2289 https://www.lydsy.com/JudgeOnline/problem.php?id=5020 solutionsoluti ...
- [BZOJ5020][THUWC 2017]在美妙的数学王国中畅游(LCT + 一点数学知识)
Address 洛谷 P4546 BZOJ 5020 LOJ #2289 Solution 如果只有一次函数 ax+bax+bax+b ,那么这是非常裸的 LCT ,维护 aaa 之和与 bbb 之和 ...
- 【BZOJ5020】[THUWC 2017]在美妙的数学王国中畅游 泰勒展开+LCT
[BZOJ5020][THUWC 2017]在美妙的数学王国中畅游 Description 数字和数学规律主宰着这个世界. 机器的运转, 生命的消长, 宇宙的进程, 这些神秘而又美妙的过程无不可以用数 ...
- [THUWC2017]在美妙的数学王国中畅游 LCT+泰勒展开+求导
[THUWC2017]在美妙的数学王国中畅游 LCT+泰勒展开+求导 Code: #include<bits/stdc++.h> using namespace std; #define ...
- 【bzoj5020】[THUWC 2017]在美妙的数学王国中畅游 泰勒展开+LCT
题目描述 学渣小R被大学的数学课程虐得生活不能自理,微积分的成绩曾是他在教室里上的课的最低分.然而他的某位陈姓室友却能轻松地在数学考试中得到满分.为了提升自己的数学课成绩,有一天晚上(在他睡觉的时候) ...
- 解题:THUWC 2017 在美妙的数学王国中畅游
题面 _"数字和数学规律主宰着这个世界."_ 在 @i207M 帮助下折腾了半天终于搞懂了导数和泰勒展开,引用某学长在考场上的感受:感觉整个人都泰勒展开了 显然是个奇奇怪怪的东西套 ...
最新文章
- python flask 设置 header 响应体、响应头、状态码
- Spring中@Autowired和@Resource区别
- Android中网络编程出错
- 成功解决ValueError: Expected 2D array, got 1D array instead: Reshape your data either using array.reshap
- 第四十三节,文件、文件夹、压缩包、处理模块shutil
- Qt之QLineEdit详解(附源码)
- caffe dataset
- VC++ 6.0如何创建与调用动态链接库
- windows11条件下将yafu路径添加到环境变量
- 大学学习历程简单总结
- AI:全球人工智能领域代表性学者简介及大佬们的关系圈—向大佬们学习(国外内分开,持续更新/建议收藏)
- 用 “类别混搭” 法实现移动游戏的创新
- Stata:因子变量全攻略-双向固定效应-超越对数生产函数
- Matplotlib-Python-绘制基础饼形图,分裂饼形图,环形饼形图
- mysql假死_win7系统假死的5种情况和处理方法
- python xgb模型 预测_如何使用XGBoost模型进行时间序列预测
- 百度百科计算机科学与技术,徐明伟(清华大学计算机科学与技术系教授)_百度百科...
- JVM垃圾回收GarbageCollect
- Qt学习笔记(三)——记事本
- inno setup打包脚本总结
热门文章
- 微信h5网页关闭分享以及关闭当前页面
- [线性代数]Note4--A的LU分解转置-置换-向量空间
- linux重启memcache_Linux中安装Memcached服务的方法
- 优麒麟桌面闪烁_稳定性持续增强,优麒麟 19.10.1 发布
- ThreadLocal系列(二)-InheritableThreadLocal的使用及原理解析
- java EE : tomacat 基础
- 配置HAProxy支持https协议
- cdoj 1246 每周一题 拆拆拆~ 分解质因数
- Java中文字符所占的字节数
- Windows SDK编程之一 窗口示例程序