@loj - 2289@「THUWC 2017」在美妙的数学王国中畅游
@description@
n 个点编号 0 到 n-1，每个点有一个从 [0,1] 映射到 [0,1] 的函数 f(x) 作为点权，它有以下几种形式：
正弦函数：sin(ax+b) (a∈[0,1],b∈[0,π],a+b∈[0,π])
指数函数：e^(ax+b) (a∈[−1,1],b∈[−2,0],a+b∈[−2,0])
一次函数：ax+b (a∈[−1,1],b∈[0,1],a+b∈[0,1])
有 m 个事件，种类如下：
（1）加边。
（2）删边。
（3）更改点权。
（4）询问 u 到 v 这条路径上的所有函数以 x（给定）为自变量的因变量之和。
保证一开始不存在任何边。
input
第一行两个正整数 n,m 和一个字符串 type。表示 n 个点与 m 个事件，type 是用来得部分分的。1≤n≤100000，1≤m≤200000 。
接下来 n 行，第 i 行表示编号为 i 的点的初始函数。一个整数 f 表示函数的类型，两个实数 a,b 表示函数的参数：
f=1 ,则函数为 f(x)=sin(ax+b)(a∈[0,1],b∈[0,π],a+b∈[0,π])
f=2 ,则函数为 f(x)=e^(ax+b) (a∈[−1,1],b∈[−2,0],a+b∈[−2,0])
f=3 ,则函数为 f(x)=ax+b(a∈[−1,1],b∈[0,1],a+b∈[0,1])
接下来 m 行，每行描述一个事件，事件分为四类：
appear u v：表示数学王国中出现了一条连接 u 和 v 这两座城市的魔法桥 (0≤u,v<n,u ≠ v) ，保证连接前 u 和 v 这两座城市不能互相到达。
disappear u v： 表示数学王国中连接 u 和 v 这两座城市的魔法桥消失了，保证这座魔法桥是存在的。
magic c f a b：表示城市 c 的魔法球中的魔法变成了类型为 f ，参数为 a,b 的函数
travel u v x：表示询问一个智商为 x 的人从城市 u 旅行到城市 v （即经过 u 到 v 这条路径上的所有城市，包括 u 和 v ）后，他得分的总和是多少。若无法从 u 到达 v ，则输出一行一个字符串 unreachable。
output
对于每个询问，输出一行实数，表示答案。建议使用科学计数法表示。
sample input
3 7 C1
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
9.45520207e-001
1.67942554e+000
1.20000000e+000
hint
【出题人教你学数学】
若函数 \(f(x)\) 的 \(n\) 阶导数在 \([a,b]\) 区间内连续，则对 \(f(x)\) 在 \(x_0(x_0\in[a,b])\) 处使用 n 次拉格朗日中值定理可以得到带拉格朗日余项的泰勒展开式：
其中，当 \(x>x_0\) 时，\(\xi\in[x_0,x]\)。当 \(x<x_0\) 时，\(\xi\in[x,x_0]\)。
\(f^{(n)}\) 表示函数 \(f\) 的 \(n\) 阶导数
@solution@
都加边删边了，那肯定是用 LCT 来维护了。
一次函数还好说，直接链上维护一次项系数与常数系数之和（但是竟然没有这部分的部分分？）。
考虑指数函数或是三角函数的和，这个东西是没有任何实际的意义的，无法用什么统一的数据表示。
但是题目中给的泰勒展开这玩意儿，其意义就是用多项式函数去逼近这些非多项式的函数。
那我们为什么不直接用泰勒展开来逼近就可以了？取前二十项的，误差就足够小了。
而多项式就可以直接用一次函数的方法，维护每一项的系数之和。
那么怎么求导呢？显然我一个初中生是做不来的。
直接摆公式吧：
指数函数：\(f(x) = e^x, f'(x) = e^x\)。
三角函数 sin：\(f(x) = sin(x), f'(x) = cos(x)\)。三角函数 cos：\(f(x) = cos(x), f'(x) = -sin(x)\)。
幂函数： \(f(x) = x^a, f'(x) = ax^{a-1}\)。常函数：\(f(x) = k, f'(x) = 0\)。
函数和的导数：\(h(x) = f(x) + g(x), h'(x) = f'(x) + g'(x)\)。
函数积的导数：\(h(x) = f(x)*g(x), h'(x) = f(x)*g'(x) + f'(x)*g(x)\)。
复合函数的导数：\(h(x) = f(g(x)), h'(x) = f'(g(x))*g'(x)\)。
那么对于函数 \(f(x) = e^{ax+b}\)，它的 n 阶导数为 \(f^{(n)}=a^ne^{ax+b}\)。
对于函数 \(f(x) = sin(x)\)，它的 n 阶导数的系数绝对值为 \(a^n\)，剩下的以四为循环节：\(sin(x)\)，\(cos(x)\)，\(-sin(x)\)，\(-cos(x)\)。
对于函数 \(f(x) = ax + b\)，这个就不讲了。
@accepted code@
#include<cmath>
#include<cstdio>
#include<algorithm>
using namespace std;
typedef double db;
const int MAXN = 100000;
struct node{
int rev;
node *ch[2], *fa;
db k[20], s[20];
}pl[MAXN + 5], *ad[MAXN + 5], *NIL;
bool isroot(node *x) {return x->fa->ch[0] != x && x->fa->ch[1] != x;}
void setchild(node *x, node *y, int d) {
if( x != NIL ) x->ch[d] = y;
if( y != NIL ) y->fa = x;
}
void pushup(node *x) {
for(int i=0;i<20;i++)
x->s[i] = x->ch[0]->s[i] + x->ch[1]->s[i] + x->k[i];
}
void pushdown(node *x) {
if( x->rev ) {
swap(x->ch[0], x->ch[1]);
if( x->ch[0] != NIL ) x->ch[0]->rev ^= 1;
if( x->ch[1] != NIL ) x->ch[1]->rev ^= 1;
x->rev = 0;
}
}
void rotate(node *x) {
node *y = x->fa; pushdown(y), pushdown(x);
int d = (y->ch[1] == x);
if( isroot(y) ) x->fa = y->fa;
else setchild(y->fa, x, (y->fa->ch[1] == y));
setchild(y, x->ch[!d], d);
setchild(x, y, !d);
pushup(y);
}
void splay(node *x) {
pushdown(x);
while( !isroot(x) ) {
node *y = x->fa;
if( isroot(y) )
rotate(x);
else {
if( (y->fa->ch[1] == y) == (y->ch[1] == x) )
rotate(y);
else rotate(x);
rotate(x);
}
}
pushup(x);
}
void access(node *x) {
node *y = NIL;
while( x != NIL ) {
splay(x);
x->ch[1] = y;
pushup(x);
y = x, x = x->fa;
}
}
void makeroot(node *x) {
access(x);splay(x);
x->rev ^= 1;
}
node *findroot(node *x) {
access(x), splay(x);
node *ret = x;
while( ret->ch[0] != NIL ) ret = ret->ch[0];
return ret;
}
void link(node *x, node *y) {
makeroot(x); x->fa = y;
}
void cut(node *x, node *y) {
makeroot(x); access(y), splay(y);
y->ch[0] = NIL; x->fa = NIL;
pushup(y);
}
void modify(db k[], int f, db a, db b) {
if( f == 1 ) {
db m = 1, x = sin(b), y = cos(b);
for(int i=0;i<20;i+=4) {
k[i + 0] = x*m, m = m*a/(i+1);
k[i + 1] = y*m, m = m*a/(i+2);
k[i + 2] = -x*m, m = m*a/(i+3);
k[i + 3] = -y*m, m = m*a/(i+4);
}
}
else if( f == 2 ) {
db m = 1, x = exp(b);
for(int i=0;i<20;i++)
k[i] = x*m, m = m*a/(i+1);
}
else if( f == 3 ) {
for(int i=0;i<20;i++)
k[i] = 0;
k[0] = b, k[1] = a;
}
}
db calculate(db k[], db x) {
db ret = 0, p = 1;
for(int i=0;i<20;i++)
ret += p*k[i], p *= x;
return ret;
}
void init() {NIL = &pl[0], NIL->ch[0] = NIL->ch[1] = NIL->fa = NIL;}
char op[10];
int main() {
int n, m; init();
scanf("%d%d%s", &n, &m, op);
for(int i=1;i<=n;i++) {
ad[i] = &pl[i], ad[i]->ch[0] = ad[i]->ch[1] = ad[i]->fa = NIL;
for(int j=0;j<20;j++)
ad[i]->k[j] = ad[i]->s[j] = 0;
int f; db a, b; scanf("%d%lf%lf", &f, &a, &b);
modify(ad[i]->k, f, a, b), pushup(ad[i]);
}
for(int i=1;i<=m;i++) {
scanf("%s", op);
if( op[0] == 'a' ) {
int u, v; scanf("%d%d", &u, &v); u++, v++;
link(ad[u], ad[v]);
}
if( op[0] == 'd' ) {
int u, v; scanf("%d%d", &u, &v); u++, v++;
cut(ad[u], ad[v]);
}
if( op[0] == 'm' ) {
int c, f; db a, b; scanf("%d%d%lf%lf", &c, &f, &a, &b); c++;
access(ad[c]), splay(ad[c]), modify(ad[c]->k, f, a, b), pushup(ad[c]);
}
if( op[0] == 't' ) {
int u, v; db x; scanf("%d%d%lf", &u, &v, &x); u++, v++;
makeroot(ad[u]);
if( findroot(ad[v]) != ad[u] ) puts("unreachable");
else {
db ans = calculate(ad[v]->s, x); int k = 0;
if( ans < 1 ) {
while( ans < 1 ) ans *= 10, k++;
printf("%0.8lfe-%03d\n", ans, k);
}
else {
while( ans >= 10 ) ans /= 10, k++;
printf("%0.8lfe+%03d\n", ans, k);
}
}
}
}
}
@details@
想到我考前还没写过 LCT。
想到我还要这道题没做。
于是就有了这篇博客。
另外，虽然题目是建议科学计数法，但是如果不用科学计数法就会因为浮点误差 WA 掉。
原因？你想啊，假如一个五位数从 10^(-6) 开始出现误差，你把它变成科学计数法，它在 10^(-7) 之前不就都没有误差了吗？
