bzoj 5020(洛谷4546) [THUWC 2017]在美妙的数学王国中畅游——LCT+泰勒展开

最新推荐文章于 2021-04-28 19:07:33 发布
最新推荐文章于 2021-04-28 19:07:33 发布
阅读量112 收藏
点赞数
文章标签： php
原文链接：http://www.cnblogs.com/Narh/p/10253975.html
版权

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=5020

　　　https://www.luogu.org/problemnew/show/P4546

如果保证 x=1 ，则可以用 LCT 维护每个点的函数值。不然的话就用 LCT 拿出那条链，dfs 一下 splay 现算。可以得 60 分。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define db double
using namespace std;
int rdn()
{
  int ret=0;bool fx=1;char ch=getchar();
  while(ch>'9'||ch<'0'){if(ch=='-')fx=0;ch=getchar();}
  while(ch>='0'&&ch<='9')ret=ret*10+ch-'0',ch=getchar();
  return fx?ret:-ret;
}
const int N=1e5+5;
int n,m,fa[N],c[N][2],sta[N],top; db sm[N],vl[N]; bool fx,rev[N];
struct Node{
  int f;db a,b;
  db cal(db x)
  {
    if(f==1)return sin(a*x+b);
    if(f==2)return exp(a*x+b);
    if(f==3)return a*x+b;
  }
}a[N];
bool isroot(int cr){return c[fa[cr]][0]!=cr&&c[fa[cr]][1]!=cr;}
void pshp(int cr){sm[cr]=sm[c[cr][0]]+sm[c[cr][1]]+vl[cr];}
void Rev(int cr){if(rev[cr]){rev[cr]=0;rev[c[cr][0]]^=1;rev[c[cr][1]]^=1;swap(c[cr][0],c[cr][1]);}}
void rotate(int x)
{
  int y=fa[x],z=fa[y],d=(x==c[y][1]);
  if(!isroot(y))c[z][y==c[z][1]]=x;//!isroot(y) not !isroot(z) !!
  fa[x]=z;
  fa[y]=x; fa[c[x][!d]]=y;
  c[y][d]=c[x][!d]; c[x][!d]=y;
  if(fx)pshp(y),pshp(x);
}
void splay(int x)
{
  sta[top=1]=x;
  for(int cr=x;!isroot(cr);cr=fa[cr])sta[++top]=fa[cr];
  for(int i=top;i;i--)Rev(sta[i]);
  int y,z;
  while(!isroot(x))
    {
      y=fa[x];z=fa[y];
      if(!isroot(y))
    ((x==c[y][0])^(y==c[z][0]))?rotate(x):rotate(y);
      rotate(x);
    }
}
void access(int cr)
{
  for(int t=0;cr;splay(cr),c[cr][1]=t,pshp(cr),t=cr,cr=fa[cr]);//pshp()!!
}
void mkrt(int cr)
{
  access(cr);splay(cr);rev[cr]^=1;
}
void link(int x,int y)
{
  mkrt(x);fa[x]=y;
}
void cut(int x,int y)
{
  mkrt(x);access(y);splay(y);
  fa[x]=c[y][0]=0;if(fx)sm[y]-=vl[x];
}
void mdfy(int cr,int f,db a1,db b1)
{
  if(fx)splay(cr),sm[cr]-=vl[cr];
  a[cr].f=f;a[cr].a=a1;a[cr].b=b1;
  if(fx)vl[cr]=a[cr].cal(1),sm[cr]+=vl[cr];
}
int fnd(int cr)//root of real tree
{
  while(!isroot(cr))cr=fa[cr];return cr;
}
db dfs(int cr,db x)
{
  db ret=a[cr].cal(x);
  if(c[cr][0])ret+=dfs(c[cr][0],x);
  if(c[cr][1])ret+=dfs(c[cr][1],x);
  return ret;
}
void print(db x)
{
  bool fx=0;int t=0;
  if(x<1){fx=1;while(x<1)x*=10,t++;}
  if(x>=10){while(x>=10)x/=10,t++;}
  printf("%.8fe%c%03d\n",x,fx?'-':'+',t);
}
int main()
{
  n=rdn();m=rdn();char ch[15];scanf("%s",ch);
  if(ch[1]=='0')fx=1;
  for(int i=1;i<=n;i++)
    {
      a[i].f=rdn();scanf("%lf%lf",&a[i].a,&a[i].b);
      if(fx)vl[i]=sm[i]=a[i].cal(1);
    }
  int u,v;db a,b,x;
  while(m--)
    {
      scanf("%s",ch);
      if(ch[0]=='a'){u=rdn()+1;v=rdn()+1;link(u,v);}
      if(ch[0]=='d'){u=rdn()+1;v=rdn()+1;cut(u,v);}
      if(ch[0]=='m'){u=rdn()+1;v=rdn();scanf("%lf%lf",&a,&b);mdfy(u,v,a,b);}
      if(ch[0]=='t')
    {
      u=rdn()+1;v=rdn()+1;scanf("%lf",&x);
      mkrt(u);access(v);splay(v);
      if(fnd(u)!=v)puts("unreachable");
      else print(fx?sm[v]:dfs(v,x));
    }
    }
  return 0;
}
View Code

关于本题的求导和泰勒展开可以看这个：https://www.cnblogs.com/zhoushuyu/p/8148732.html

似乎还可以麦克劳林展开：https://www.cnblogs.com/Troywar/p/8982707.html

似乎展开成 12 项的多项式就差不多了，然后可以用 LCT 维护 12 项的系数和，算的时候乘上 xk 再除以阶乘就行。

那个 x0 好像取 0 就行了。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define db double
#define ls c[cr][0]
#define rs c[cr][1]
int rdn()
{
  int ret=0;bool fx=1;char ch=getchar();
  while(ch>'9'||ch<'0'){if(ch=='-')fx=0;ch=getchar();}
  while(ch>='0'&&ch<='9')ret=ret*10+ch-'0',ch=getchar();
  return fx?ret:-ret;
}
using namespace std;
const int N=1e5+5,M=12;
int n,fa[N],c[N][2],sta[N],top; bool rev[N]; db sm[N][M+5],vl[N][M+5],jc[N];
struct Node{
  int f;db a,b;
  Node(int f=0,db a=0,db b=0):f(f),a(a),b(b) {}
  void cz(int cr)
  {
    if(f==1)
      {
    db p[5],lja=1;
    p[0]=sin(b);p[1]=cos(b);p[2]=-p[0];p[3]=-p[1];
    for(int i=0;i<M;i+=4)
      for(int j=0;j<4;j++,lja*=a)sm[cr][i+j]=vl[cr][i+j]=p[j]*lja;
      }
    if(f==2)
      {
    db ml=exp(b),lja=1;
    for(int i=0;i<M;i++,lja*=a)sm[cr][i]=vl[cr][i]=ml*lja;
      }
    if(f==3)
      {
    sm[cr][0]=vl[cr][0]=b; sm[cr][1]=vl[cr][1]=a;
    for(int i=2;i<M;i++)sm[cr][i]=vl[cr][i]=0;
      }
  }
}w[N];
bool isrt(int cr){return c[fa[cr]][0]!=cr&&c[fa[cr]][1]!=cr;}
void pshp(int cr)
{
  for(int i=0;i<M;i++)sm[cr][i]=sm[ls][i]+sm[rs][i]+vl[cr][i];
}
void Rev(int cr){if(rev[cr]){rev[cr]=0;rev[ls]^=1;rev[rs]^=1;swap(ls,rs);}}
void rotate(int x)
{
  int y=fa[x],z=fa[y],d=(x==c[y][1]);
  if(!isrt(y))c[z][y==c[z][1]]=x;
  fa[x]=z;
  fa[y]=x; fa[c[x][!d]]=y;
  c[y][d]=c[x][!d]; c[x][!d]=y;
  pshp(y);pshp(x);
}
void splay(int x)
{
  sta[top=1]=x;
  for(int cr=x;!isrt(cr);cr=fa[cr])sta[++top]=fa[cr];
  for(int i=top;i;i--)Rev(sta[i]);
  int y,z;
  while(!isrt(x))
    {
      y=fa[x];z=fa[y];
      if(!isrt(y))
    ((x==c[y][0])^(y==c[z][0]))?rotate(x):rotate(y);
      rotate(x);
    }
}
void access(int cr)
{
  for(int t=0;cr;splay(cr),rs=t,pshp(cr),t=cr,cr=fa[cr]);//pshp()!!!
}
void mkrt(int cr)
{
  access(cr);splay(cr);rev[cr]^=1;
}
void link(int x,int y)
{
  mkrt(x);fa[x]=y;
}
void cut(int x,int y)
{
  mkrt(x);access(y);splay(y);
  fa[x]=c[y][0]=0;pshp(y);
}
void mdfy(int cr,int f,db a,db b)
{
  splay(cr);
  w[cr]=Node(f,a,b);w[cr].cz(cr);
}
db cal(int cr,db x)
{
  db ret=0,ml=1;
  for(int i=0;i<M;i++,ml*=x)
    ret+=sm[cr][i]*ml/jc[i];
  return ret;
}
int fnd(int cr)
{
  while(!isrt(cr))cr=fa[cr];return cr;
}
void print(db x)
{
  bool fx=0;int t=0;
  if(x<1){fx=1;while(x<1)x*=10,t++;}
  if(x>=10){while(x>=10)x/=10,t++;}
  printf("%.8fe%c%03d\n",x,fx?'-':'+',t);
}
int main()
{
  jc[0]=1;for(int i=1;i<M;i++)jc[i]=jc[i-1]*i;
  n=rdn();int m=rdn();char ch[15];scanf("%s",ch);
  for(int i=1;i<=n;i++)
    w[i].f=rdn(),scanf("%lf%lf",&w[i].a,&w[i].b),w[i].cz(i),pshp(i);//pshp for sm[]
  int u,v;db a,b;
  while(m--)
    {
      scanf("%s",ch);
      if(ch[0]=='a'){u=rdn()+1;v=rdn()+1;link(u,v);}
      if(ch[0]=='d'){u=rdn()+1;v=rdn()+1;cut(u,v);}
      if(ch[0]=='m'){u=rdn()+1;v=rdn();scanf("%lf%lf",&a,&b);mdfy(u,v,a,b);}
      if(ch[0]=='t')
    {
      u=rdn()+1;v=rdn()+1;scanf("%lf",&a);
      mkrt(u);access(v);splay(v);
      if(fnd(u)!=v)puts("unreachable");
      else print(cal(v,a));
    }
    }
  return 0;
}

 

转载于:https://www.cnblogs.com/Narh/p/10253975.html

weixin_34292402
关注
【BZOJ5020】[THUWC2017] 在美妙的数学王国中畅游
Yu Gi Oh!
08-07 226
题目链接 题意 一个森林,支持加边删边 每一个点上有一个关于x的函数,函数值在[0,1]内 每次给出一个x并询问一条路径上所有点的函数值的和,x也是[0,1] Sol 所有人智商为1时,直接算函数值并用LCT维护路径即可 如果只有一次函数就维护∑a∑a\sum a和∑b∑b\sum b 然后就是怎么搞这两个不是多项式的函数了 为什么不好处理,因为我们不能通过x的次数相同的把和...
BZOJ1009狗头考试——DP+KMP+矩阵乘法优化
stevensonson的博客
10-27 303
Description 阿申准备报名参加GT考试，准考证号为N位数X1X2…Xn(0&lt;=Xi&lt;=9),他不希望准考证号上出现不吉利的数字。 他的不吉利数学A1A2…Am(0&lt;=Ai&lt;=9)有M位，不出现是指X1X2…Xn中没有恰好一段等于A1A2…Am. A1和X1可以为 0 Input 第一行输入N,M,K.接下来一行输入M位的数。 N&lt;=10^9,M&lt;=20...
【BZOJ5020】【THUWC2017】在美妙的数学王国中畅游（Link-Cut Tree，组合数学）
小蒟蒻yyb的博客
12-26 472
题解Description数字和数学规律主宰着这个世界。 机器的运转， 生命的消长， 宇宙的进程， 这些神秘而又美妙的过程无不可以用数学的语言展现出来。 这印证了一句古老的名言： “学好数理化，走遍天下都不怕。” 学渣小R被大学的数学课程虐得生活不能自理，微积分的成绩曾是他在教室里上的课的最低分。然而他的某位陈姓室友却能轻松地在数学考试中得到满分。为了提升自己的数学课成绩，有一天晚上（
BZOJ 5020 [LCT][数学]
Vectorxj的博客
09-17 412
SolutionSolutionPBS还是牛逼的啊。 把这些函数泰勒展开后只要维护多项式的系数就好了。。 我是真的菜，修改的时候没有MakeRoot。 还是不会LCT啊QAQ#include <cmath> #include <cstdio> #include <cstdlib> #include <iostream> using namespace std;const int N = 301
bzoj5020[THUWC 2017]在美妙的数学王国中畅游（LCT）
a1035719430的博客
11-30 275
传送门：https://loj.ac/problem/2289 https://www.lydsy.com/JudgeOnline/problem.php?id=5020 solutionsolutionsolution: 题目给了公式就比较容易了 x0=0x_0=0x0​=0的时候比较容易求解 sin(ax+b)=sin(b)+acos(b)1!x−a2sin(b)2!x2−a3cos(b)3...
[BZOJ5020][THUWC 2017]在美妙的数学王国中畅游（LCT + 一点数学知识）
女装的你，如此好看！
12-16 337
Address 洛谷 P4546 BZOJ 5020 LOJ #2289 Solution 如果只有一次函数 ax+bax+bax+b ，那么这是非常裸的 LCT ，维护 aaa 之和与 bbb 之和即可 然后你会发现 ea1x+b1×ea2x+b2=e(a1x+b1)(a2x+b2)=ea1a2x2+(a1b2+a2b1)x+b1b2e^{a_1x+b_1}\times e^{a_2x+b...
#网络流#BZOJ 3996 洛谷 3973 线性代数
lemondinosaur的博客
10-05 149
题目 给定一个n×nn×nn×n的矩阵BBB和一个1×n1×n1×n的矩阵CCC。求出一个1×n1×n1×n的01矩阵AAA。使得D=(A×B−C)×ATD=(A×B-C)×A^{\sf T}D=(A×B−C)×AT最大,其中ATA^{\sf T}AT为AAA的转置。输出DDD。 分析 D=∑i=1n∑j=1nAiAjBi,j−∑i=1nAiCiD=\sum_{i=1}^n\sum_{j=1}^...
BZOJ 1726: [Usaco2006 Nov]Roadblocks第二短路 ——Dijkstra+玄学
Q755100802的博客
03-07 238
这个题玄学冲过，规定每个点访问次数不能超过50次，然后找优先队列中第二次到达终点t的状态返回就ok 记录一下，怕忘了 #include <cstdio> #include <cstring> #include <algorithm> #include <iostream> #include <vector> #include &lt...
LCT（动态树）（洛谷P3690）（BZOJ3282）
forezxl的博客
01-07 545
推荐这篇论文 前置技能 Splay和树链剖分 简介 动态树问题, 简称LCT，近几年在OI中兴起的一种新型问题，是一类要求维护一个有根树森林，支持对树的分割, 合并等操作的问题。 算法用途 比如维护一个数据结构，有如下操作： 1.查询路径权值和。 2.查询LCA。 3.修改单个点的权值。 4.断开两个点，使之成为两棵树。 5.合并两个点，使之成为一颗树。 当然还...
BZOJ2721: [Violet 5]樱花（洛谷P1445）
forezxl的博客
10-28 369
数学 BZOJ题目传送门 洛谷题目传送门 来因式分解。 1x+1y=1n!xy−n!(x+y)=0(x−n!)(y−n!)=(n!)2 \frac1x+\frac1y=\frac1{n!}\\ xy-n!(x+y)=0\\ (x-n!)(y-n!)=(n!)^2 x1​+y1​=n!1​xy−n!(x+y)=0(x−n!)(y−n!)=(n!)2 然后算n!n!n!的因数个数就好了。 代码： #i...
JXI5020GFN-B
02-14
一款很好用的LED驱动器（数据串行输入，并行输出），同时驱动16路LED，各路电流差异极小，保证发光亮度一致性
[bzoj5020][THUWC 2017]在美妙的数学王国中畅游
蒟蒻YZH的博客
01-09 476
Orz w_yqts lct+泰勒展开#include <bits/stdc++.h> using namespace std; #define db double #define N 666666 #define m 20 struct lct { int c[2],f,rev; db s[m],v[m]; }tr[N]; db C[m][m]; int top,stk[N];
bzoj5020 在美妙的数学王国中畅游【LCT+泰勒展开】
cdsszjj的博客
12-20 450
解题思路：题目后貌似还有一段提示： 考试的时候完全看不懂它给的泰勒展开式，只打了可以暴力走和x=1x=1的部分，就是LCT维护路径和。考试后看了题解，原来它的意思就是取x0=0x_0=0，f(x)f(x)可以化成多项式形式，每项系数就是fi(0)i!。n=11\frac{f^i(0)}{i!}。n=11精度就够了，然后Splay维护每一项系数的和，询问时每项乘以xix^i就行了。对导数不熟悉的
bzoj 5020: [THUWC 2017]在美妙的数学王国中畅游【泰勒展开+LCT】
weixin_30551963的博客
02-20 97
参考：https://www.cnblogs.com/CQzhangyu/p/7500328.html ……其实理解了泰勒展开之后就是水题呢可是我还是用了两天时间来搞懂啊 泰勒展开是到正无穷的，但是因为精度问题，所以一般展开十几项就可以（这里展开了17项）。以下是公式： \[ e^x=\sum_{i=0}^{\infty}\frac{x^i}{i!} \] \[ sin(x)=\sum_{i=0}...
WC模拟（12.20） T1 在美妙的数学王国中畅游（bzoj5020）
scarlyw的博客
12-21 490
link cut tree + 泰勒展开
【BZOJ5020】【THUWC2017】在美妙的数学王国中畅游 LCT 泰勒展开
ez_yww的博客
10-30 1011
题目大意　　给你一棵树，每个点有一个函数f(x)f(x) 正弦函数 sin(ax+b)(a∈[0,1],b∈[0,π],a+b∈[0,π])\sin(ax+b) (a\in[0,1],b\in[0,\pi],a+b\in[0,\pi]) 指数函数 eax+b(a∈[−1,1],b∈[−2,0],a+b∈[−2,0])e^{ax+b} (a\in[−1,1],b\in[−2,0],a+b\in[−2,
[LOJ#2289 && BZOJ5020][THUWC 2017]在美妙的数学王国中畅游
CE玩家
07-02 1115
终于找到一个可以交的地方了… 考场上只写了LCT的60分暴力，因为那时候并看不懂什么泰勒展开…前段时间学了微积分，学了泰勒展开，大概知道了是怎么回事 其实题目说的很清楚了…但是那时候就是看不懂，也不会求导什么的。讲sin(ax+b)\sin(ax+b),eax+be^{ax+b}和ax+bax+b在x=0处展开可以得到： sin(ax+b)=sin(b)+acos(b)x1!−a2sin(
linux文件每行都多出一行空白,Linux：去除每一行行首的空格
weixin_35785014的博客
04-28 438
【BZOJ】1225&colon; &lbrack;HNOI2001&rsqb; 求正整数http://www.lydsy.com/JudgeOnline/problem.php?id=1225 题意:给一个数n,求一个最小的有n个约数的正整数.(n<=50000) #include ...IncDec Sequence题目链接 http://www.lydsy.co...
【BZOJ2115】Xor(线性基---求无向图1->n路径上最大的边权异或值--可出现环)
Cassie_zkq的博客
07-12 763
题目地址：https://www.lydsy.com/JudgeOnline/problem.php?id=2115 题目： Description Input 第一行包含两个整数N和 M， 表示该无向图中点的数目与边的数目。 接下来M 行描述 M 条边，每行三个整数Si，Ti ，Di，表示 Si 与Ti之间存在 一条权值为 Di的无向边。 图中可能有重边或自环。 Output 仅...
洛谷 P1640 BZOJ 1854 [SCOI2010]连续攻击游戏
最新发布
05-14
这道题目需要进行博弈论的思考。我们可以使用动态规划来解决这道题目。 首先，我们可以定义 $f_{i,j}$ 表示当前剩下 $i$ 到 $j$ 这些数时，当前玩家最多能获得的分数。我们可以发现，在这个状态中，当前玩家只能从 $i$ 或者 $j$ 中选择一个数进行取走，然后将剩下的部分转化为对手的状态。因此，我们可以得到如下的转移方程： $$ f_{i,j}=\max\{a_i+\min\{f_{i+2,j},f_{i+1,j-1}\},a_j+\min\{f_{i+1,j-1},f_{i,j-2}\}\} $$ 其中 $a_i$ 表示第 $i$ 个数的分值。上述方程的含义是，当前玩家可以选择取走 $i$ 或者 $j$，如果取走 $i$，则转化为对手的状态是 $f_{i+2,j}$ 或者 $f_{i+1,j-1}$，如果取走 $j$，则转化为对手的状态是 $f_{i,j-2}$ 或者 $f_{i+1,j-1}$。当前玩家希望获得的分数是最大的，因此需要取两种情况中的较大值。 最终答案就是 $f_{1,n}$，表示从 $1$ 到 $n$ 这些数中，当前玩家最多能获得的分数。 最后，需要注意的是，我们需要使用逆序的方式来进行动态规划，即从小到大枚举区间长度 $len$，然后枚举区间起点 $i$，最终计算 $f_{1,n}$。这样可以保证在计算 $f_{i,j}$ 时，$f_{i+1,j}$ 和 $f_{i,j-1}$ 已经被计算出来了。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值