并不对劲的多项式求ln,exp

最新推荐文章于 2020-06-10 17:37:32 发布
最新推荐文章于 2020-06-10 17:37:32 发布
阅读量89 收藏
点赞数
原文链接：http://www.cnblogs.com/xzyf/p/10519392.html
版权
ln
解释

设\(g(x)=ln(f(x))\)，两边同时求导，则有：\(g'(x)=ln'(f(x))*f'(x)=f^{-1}(x)*f'(x)\)(1)
因为\(f(x)\)是个多项式，所以设\(f(x)=\sum_{i=0}^{n}a_i*x^i\)，则有\(f'(x)=\sum_{i=0}^{n-1}a_{i+1}*(i+1)*x^i\)
设\(h(x)=f^{-1}(x)*f'(x)\)，对(1)式两边同时求积分，得：\(g(x)=\int h(x)\space dx\)
因为\(h(x)\)是个多项式，设\(h(x)=\sum_{i=0}^{n}b_i*x^i\)就有：\(\int h(x)\space dx=\sum_{i=1}^{n}\frac{b_{i-1}}{i}x^i\)

代码
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<iomanip>
#include<iostream>
#include<map>
#include<queue>
#include<set>
#include<stack>
#include<vector>
#define rep(i,x,y) for(register int i=(x);i<=(y);++i)
#define dwn(i,x,y) for(register int i=(x);i>=(y);--i)
#define maxn 800010
#define LL long long 
#define mo(x) (x>=mod?x-mod:(x<0?x+mod:x))
using namespace std;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)&&ch!='-')ch=getchar();
    if(ch=='-')f=-1,ch=getchar();
    while(isdigit(ch))x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
    return x*f;
}
void write(int x)
{
    if(x==0){putchar('0');return;}
    int f=0;char ch[20];
    if(x<0)putchar('-'),x=-x;
    while(x)ch[++f]=x%10+'0',x/=10;
    while(f)putchar(ch[f--]);
    return;
}
const LL mod=998244353;
int f[maxn],g[maxn],h[maxn],nowlen,nown,r[maxn],tmp[maxn],n;
int mul(int x,int y){int res=1;while(y){if(y&1)res=(LL)res*(LL)x%mod;x=(LL)x*(LL)x%mod,y>>=1;}return res;}
void dnt(int * u,int fh)
{
    rep(i,0,nown-1)r[i]=(r[i>>1]>>1)|((i&1)<<(nowlen-1));
    rep(i,0,nown-1)if(i<r[i])swap(u[i],u[r[i]]);
    for(int i=1;i<nown;i<<=1)
    {
        int wn=mul(3,(mod-1)/(i<<1)); 
        if(fh==-1)wn=mul(wn,mod-2);
        for(int j=0;j<nown;j+=(i<<1))
        {
            int w=1;
            rep(k,0,i-1){int x=u[j+k],y=(LL)w*(LL)u[i+j+k]%mod;u[j+k]=mo(x+y),u[i+j+k]=mo(x-y),w=(LL)w*(LL)wn%mod;}
        }
    }
    if(fh==-1)
    {
        int inv=mul(nown,mod-2);
        rep(i,0,nown-1)u[i]=(LL)u[i]*(LL)inv%mod;
    }
}
void getny(int * u,int * v,int nlen)
{
    v[0]=mul(u[0],mod-2);
    for(int len=0,tmpn=1;tmpn<nlen;len++,tmpn<<=1)
    {
        nown=tmpn<<1,nowlen=len+1;
        rep(i,0,nown-1)tmp[i]=u[i];
        nown<<=1,nowlen++;
        rep(i,(tmpn<<1),nown-1)tmp[i]=0;
        dnt(tmp,1),dnt(v,1);
        rep(i,0,nown-1)v[i]=mo(2ll-(LL)tmp[i]*(LL)v[i]%mod)*(LL)v[i]%mod;
        dnt(v,-1);
        rep(i,(tmpn<<1),nown-1)v[i]=0;
    }
    rep(i,nlen,nown)v[i]=0;
}
void getup(int * u,int * v,int nlen)
{
    rep(i,1,nlen-1)v[i-1]=(LL)i*(LL)u[i]%mod;v[nlen-1]=0;
}
void getdx(int * u,int * v,int nlen)
{
    rep(i,1,nlen-1)v[i]=(LL)u[i-1]*(LL)mul(i,mod-2)%mod;v[0]=0;
}
void getln(int * u,int nlen)
{
    getup(u,h,nlen),getny(u,g,nlen);
    for(nowlen=0,nown=1;nown<(nlen+nlen);nowlen++,nown<<=1);
    dnt(h,1),dnt(g,1);
    rep(i,0,nown-1)h[i]=(LL)h[i]*(LL)g[i]%mod;
    rep(i,0,nown-1)g[i]=0;
    dnt(h,-1);getdx(h,g,nlen);
    rep(i,nlen,nown)g[i]=0;
}
int main()
{
    n=read();
    rep(i,0,n-1)f[i]=read();
    getln(f,n);
    rep(i,0,n-1)write(g[i]),putchar(' ');
    return 0;
}
exp
解释

这个人（点这里）讲得很清楚

代码
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<iomanip>
#include<iostream>
#include<map>
#include<queue>
#include<set>
#include<stack>
#include<vector>
#define rep(i,x,y) for(register int i=(x);i<=(y);++i)
#define dwn(i,x,y) for(register int i=(x);i>=(y);--i)
#define maxn 800010
#define LL long long 
#define mo(x) (x>=mod?x-mod:(x<0?x+mod:x))
using namespace std;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)&&ch!='-')ch=getchar();
    if(ch=='-')f=-1,ch=getchar();
    while(isdigit(ch))x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
    return x*f;
}
void write(int x)
{
    if(x==0){putchar('0');return;}
    int f=0;char ch[20];
    if(x<0)putchar('-'),x=-x;
    while(x)ch[++f]=x%10+'0',x/=10;
    while(f)putchar(ch[f--]);
    return;
}
const LL mod=998244353;
int f[maxn],g[maxn],h[maxn],ans[maxn],nowlen,nown,r[maxn],tmp[maxn],tmp2[maxn],n;
int mul(int x,int y){int res=1;while(y){if(y&1)res=(LL)res*(LL)x%mod;x=(LL)x*(LL)x%mod,y>>=1;}return res;}
void dnt(int * u,int fh)
{
    rep(i,0,nown-1)r[i]=(r[i>>1]>>1)|((i&1)<<(nowlen-1));
    rep(i,0,nown-1)if(i<r[i])swap(u[i],u[r[i]]);
    for(int i=1;i<nown;i<<=1)
    {
        int wn=mul(3,(mod-1)/(i<<1)); 
        if(fh==-1)wn=mul(wn,mod-2);
        for(int j=0;j<nown;j+=(i<<1))
        {
            int w=1;
            rep(k,0,i-1){int x=u[j+k],y=(LL)w*(LL)u[i+j+k]%mod;u[j+k]=mo(x+y),u[i+j+k]=mo(x-y),w=(LL)w*(LL)wn%mod;}
        }
    }
    if(fh==-1)
    {
        int inv=mul(nown,mod-2);
        rep(i,0,nown-1)u[i]=(LL)u[i]*(LL)inv%mod;
    }
}
void getny(int * u,int * v,int nlen)
{
    v[0]=mul(u[0],mod-2);
    for(int len=0,tmpn=1;tmpn<nlen;len++,tmpn<<=1)
    {
        nown=tmpn<<1,nowlen=len+1;
        rep(i,0,nown-1)tmp[i]=u[i];
        nown<<=1,nowlen++;
        rep(i,(tmpn<<1),nown-1)tmp[i]=0;
        dnt(tmp,1),dnt(v,1);
        rep(i,0,nown-1)v[i]=mo(2ll-(LL)tmp[i]*(LL)v[i]%mod)*(LL)v[i]%mod;
        dnt(v,-1);
        rep(i,(tmpn<<1),nown-1)v[i]=0;
    }
    rep(i,0,nown)tmp[i]=0; 
    rep(i,nlen,nown)v[i]=0;
}
void getup(int * u,int * v,int nlen)
{
    rep(i,1,nlen-1)v[i-1]=(LL)i*(LL)u[i]%mod;v[nlen-1]=0;
}
void getdx(int * u,int * v,int nlen)
{
    rep(i,1,nlen-1)v[i]=(LL)u[i-1]*(LL)mul(i,mod-2)%mod;v[0]=0;
}
void getln(int * u,int * v,int nlen)
{
    getup(u,h,nlen),getny(u,g,nlen);
    for(nowlen=0,nown=1;nown<(nlen+nlen);nowlen++,nown<<=1);
    dnt(h,1),dnt(g,1);
    rep(i,0,nown-1)h[i]=(LL)h[i]*(LL)g[i]%mod;
    rep(i,0,nown-1)g[i]=0;
    dnt(h,-1);getdx(h,v,nlen);
    rep(i,nlen,nown)v[i]=0;
}
void getexp(int * u,int * v,int nlen)
{
    rep(i,0,(nlen<<2))v[i]=0;
    v[0]=1;
    for(int len=0,tmpn=1;tmpn<nlen;len++,tmpn<<=1)
    {
        rep(i,0,(tmpn<<1))tmp2[i]=0;
        getln(v,tmp2,(tmpn<<1));
        nown=(tmpn<<2),nowlen=len+2;
        rep(i,(tmpn<<1),nown)tmp2[i]=0;
        rep(i,0,(tmpn<<1)-1)tmp[i]=u[i];
        rep(i,(tmpn<<1),nown)tmp[i]=0;
        dnt(tmp2,1),dnt(v,1),dnt(tmp,1);
        rep(i,0,nown-1)v[i]=(LL)v[i]*(LL)mo(mo(1ll-tmp2[i])+tmp[i])%mod;
        dnt(v,-1);
        rep(i,(tmpn<<1),nown)v[i]=0;
    }
    rep(i,nlen,nown)v[i]=0;
}
int main()
{
    n=read();
    rep(i,0,n-1)f[i]=read();
    getexp(f,ans,n);
    rep(i,0,n-1)write(ans[i]),putchar(' ');
    return 0;
}

转载于:https://www.cnblogs.com/xzyf/p/10519392.html

weixin_30851409
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
multinomial-ln:log（多项式系数）从对数一直向下
07-08
multinomialLn ( [ 1 , 5 , 7 , 5 , 5 , 1 , 2 , 3 ] )// 106.60573852903470 由于用于生成这些对数多项式系数的阶乘函数是近似的，因此这样的结果并不准确，尽管考虑到该值的大小 ( e^106.6 )，但实际值非常接近。...
Coefficients of Polynomial Chaos Expansion (PCE)：这个简单的 Matlab 代码计算对数正态分布不确定参数 X 的 PCE 系数，最高可达 N 阶。-matlab开发
05-29
标题中的"Polynomial Chaos Expansion (PCE)"，也称为多项式混沌展开，是一种在工程、科学和金融领域广泛应用的...通过深入理解这个代码，工程师和科学家们能够更有效地处理不确定性问题，优化设计决策，并降低风险。
多项式求ln
gigo_64 sink in my world
08-06 171
求 两边求导 所以对于g’求导，乘g的逆元，然后得到f'，再积分回去。 关于求导： 关于积分： 关于代码 #include<bits/stdc++.h> using namespace std; #define in read() #define int long long int in{ int cnt=0,f=1;char ch=0; while(!isdigit...
多项式求ln,求exp,开方,快速幂 学习总结
weixin_30737363的博客
06-21 356
按理说Po姐姐三月份来讲课的时候我就应该学了 但是当时觉得比较难加上自己比较懒，所以就QAQ了 现在不得不重新弄一遍了 首先说多项式求ln 设G(x)=lnF(x) 我们两边求导可以得到G'(x)=F‘(x)/F(x) 则G(x)就是F’(x)/F(x)的积分 我们知道多项式求导和积分是O(n)的,多项式求逆是O(nlogn)的 所以总时间复杂度O(nlogn) ...
BZOJ3456 城市规划 【多项式求ln】
weixin_30268921的博客
06-15 106
题目链接 BZOJ3456 题解 真是一道经典好题，至此已经写了分治\(NTT\)，多项式求逆，多项式求\(ln\)三种写法 我们发现我们要求的是大小为\(n\)无向联通图的数量 而\(n\)个点的无向图是由若干个无向联通图构成的 那么我们设\(F(x)\)为无向联通图数量的指数型生成函数 设\(G(x)\)为无向图数量的指数型生成函数 \(G(x)\)很好求 而 \[G(x) = \frac{F...
洛谷 4389 付公主的背包——多项式求ln、exp
weixin_33973600的博客
02-18 116
题目：https://www.luogu.org/problemnew/show/P4389 关于泰勒展开： https://blog.csdn.net/SoHardToNamed/article/details/80550935 https://www.cnblogs.com/guo-xiang/p/6662881.html 大概就是：\( f(x) = \sum\limits_{i=0...
多项式小结 （求逆）（ln）（exp）（除法）（多点求值）（快速插值）（任意模数NTT）
FSYo 的博客
11-07 644
多项式求逆 f(x)∗g(x)=1(mod xn)f(x)*g(x)=1(mod\ x^n)f(x)∗g(x)=1(mod xn) 先求出 f(x)∗g(x)′=1(mod xn2)f(x)*g(x)'=1(mod\ x^{\frac{n}{2}})f(x)∗g(x)′=1(mod x2n​) 然后有 f(x)∗g(x)′=f(x)∗g(x)(mod&nbs...
ln_exp_pow.rar_POW_ln e
09-14
"ln_exp_pow.rar_POW_ln e"这个标题暗示了我们即将探讨的主题是关于自然对数、指数函数以及它们之间的关系，特别是涉及到自然底数e的幂运算。 自然对数ln（对数以e为底）是指数函数的逆运算，而e是自然对数的基数，...
Maple常用计算命令
08-15
- **Randpoly, Randprime**：用于生成有限域内的随机多项式和首一素数多项式。 - **ithprime**：用于找到第 i 个素数。 - **nextprime, prevprime**：分别用于找到比给定数更大的下一个素数和比给定数更小的前一个...
实验四 MATLAB数值计算.pdf
07-14
总的来说，实验四通过MATLAB实践，帮助学生熟悉并掌握了数值计算的核心技术，包括数据处理、函数拟合、微积分运算以及微分方程求解等，这些都是科学研究和工程应用中不可或缺的技能。通过这样的实验训练，学生能够更...
多项式全家桶（三）：多项式的ln，exp
周道-Althen的博客
04-03 1050
      &ThinSpace;\ \ \ \ \ \ \,      对数和指数是很重要的东西了，复杂的多项式都和他们有关系，所以说掌握很重要，这里不建议光背板子，因为这两个板子都有致命的限制，而在实际操作的时候，可以通过一些方法绕过这个限制直接求解，这个就很重要了。...
[LOJ#6268] 分拆数 [多项式求ln][多项式exp]
最后一次修改
04-24 241
Link LOJ - https://loj.ac/problem/6268 http://danihao123.is-programmer.com/posts/213392.html
多项式ln,exp学习小计
Deep_Thinking的博客
06-10 743
主要作为一个模板和总结 多项式ln 求G(x)=ln(F(x))G(x)=ln(F(x))G(x)=ln(F(x)) 两边求导 G′(x)=F′(x)F(x)G'(x)=\frac{F'(x)}{F(x)}G′(x)=F(x)F′(x)​，注意这是复合函数求导。 多项式求逆，求导，乘法即可求出G′(x)G'(x)G′(x) 再积分回去就是G(x)G(x)G(x) 一般会钦定常数项为1，否则ln的时候先求导再积分（这里是不定积分）你会发现常数项丢失了，得到的结果并不是原来的多项式，而是在每一项的基础上除以了
指数型生成函数 及 多项式求ln
weixin_30515513的博客
06-15 272
指数型生成函数 我们知道普通型生成函数解决的是组合问题，而指数型生成函数解决的是排列问题 对于数列\(\{a_n\}\)，我们定义其指数型生成函数为 \[G(x) = a_0 + a_1x + a_2\frac{x^2}{2!} + a_3\frac{x^3}{3!} + a_4\frac{x^4}{4!} + \dots = \sum\limits_{i = 0}^{\infty} a_i\f...
多项式exp，ln，求逆板子
LYD729
02-16 251
题目是jzoj 5923 #include&amp;lt;cstdio&amp;gt; #include&amp;lt;cstring&amp;gt; #include&amp;lt;algorithm&amp;gt; #define fo(i,a,b) for(int i=a;i&amp;lt;=b;++i) #define fd(i,b,a) for(int i=b;i&amp;gt;=a;--i) #define max(x,y) ((x
luoguP4389 付公主的背包 多项式求逆 多项式求ln 多项式求exp 生成函数
maxtir的博客
12-29 216
luoguP4389 付公主的背包 题目传送门 分析 神仙题系列。。。 首先写出每件商品的生成函数，假设体积为VVV A(x)=∑i∞xiv=11−xvA(x)=\sum_i^{\infty}x^{iv}=\frac{1}{1-x^{v}}A(x)=∑i∞​xiv=1−xv1​ 答案就是 Ans(x)≡∏11−xvi(modxm+1)Ans(x)\equiv\prod\frac{1}{1-x^{v...
[多项式ln][多项式exp][多项式求幂][生成函数][DP][FNT] BZOJ 3684: 大朋友和多叉树
Vectorxj的博客
03-06 796
SolutionSolutionSolution 把DP写成生成函数的形式。f(x)=x+∑d∈Dfd(x)f(x)=x+∑d∈Dfd(x)f(x)=x+\sum_{d\in D}f^d(x)设g(f(x))=xg(f(x))=xg\big(f(x)\big)=x，有g(f(x))g(x)==f(x)−∑d∈Dfd(x)x−∑d∈Dxdg(f(x))=f(x)−∑d∈Dfd(x)g(x)=x−∑...
【Luogu4238】多项式取Ln模板题
დ♂Hany01's Blog Space♂ღ
07-28 241
Description 多项式取lnln\ln。 Solution lnf(x)=∫f′(x)f−1(x)ln⁡f(x)=∫f′(x)f−1(x)\ln f(x)=\int f'(x)f^{-1}(x) Code /************************************************ * Au: Hany01 * Date: Jul 2...
多项式对数函数（多项式求ln）小结
a_forever_dream的博客
05-19 1757
多项式对数函数 顾名思义，给出一个多项式 A(x)A(x)A(x)，求出一个多项式 G(x)G(x)G(x)，满足 G(x)=ln⁡A(x)(modxn)G(x)=\ln A(x)\pmod {x^n}G(x)=lnA(x)(modxn)。 设 F(x)=ln⁡(x)F(x)=\ln (x)F(x)=ln(x)，那么有 G(x)=F(A(x))G(x)=F(A(x))G(x)=F(A(x))，对等式两边同时求导，得到 G′(x)=F′(A(x))A′(x) G'(x)=F'(A(x))A'(x) G′(x)
python 求ln
最新发布
10-07
在Python中，可以使用math模块中的log函数来计算ln。log函数默认以e为底，因此可以直接使用log(x)来求ln(x)的值。下面是一个示例代码： import math x = 5 ln_x = math.log(x) print("ln(", x, ") = ", ln_x) 相关问题: 1. Python中如何计算其他对数函数，例如log10和log2？ 2. 如何在Python中使用对数函数来解决实
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值