前置知识:求导,积分,多项式求逆
模板
题目大意:给出
n
−
1
n−1
n−1 次多项式
A
(
x
)
A(x)
A(x),求一个
m
o
d
x
n
\bmod x^n
mod xn 下的多项式
B
(
x
)
B(x)
B(x),满足
B
(
x
)
≡
ln
A
(
x
)
(
m
o
d
x
n
)
B(x)\equiv\lnA(x)\pmod{x^n}
B(x)≡lnA(x)(modxn).
系数对
998244353
998244353
998244353 取模
这题不难,关键是要想到思路。
多项式求对数怎么是一个多项式呢?所以要想办法把 ln \ln ln 消去。
两边同时求导得
B
′
(
x
)
≡
A
′
(
x
)
A
(
x
)
(
m
o
d
x
n
)
B'(x)\equiv\dfrac{A'(x)}{A(x)}\pmod{x^n}
B′(x)≡A(x)A′(x)(modxn)
此时发现 A ′ ( x ) A'(x) A′(x) 和 A ( x ) A(x) A(x) 都是已知的。所以可以使用多项式求逆求出 A ′ ( x ) A ( x ) \dfrac{A'(x)}{A(x)} A(x)A′(x)。
这时再两边同时积分,得到
B
(
x
)
≡
∫
A
′
(
x
)
A
(
x
)
(
m
o
d
x
n
)
B(x)\equiv\int\dfrac{A'(x)}{A(x)}\pmod{x^n}
B(x)≡∫A(x)A′(x)(modxn)
这样对 A ′ ( x ) A ( x ) \dfrac{A'(x)}{A(x)} A(x)A′(x) 积分就能得到 B ( x ) B(x) B(x) 了
注意常数项为 0 0 0
代码
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=(1<<18)+1;
const ll mod=998244353,g=3;
int len=1,n;
ll a1[N],w,wn,a[N],inva[N],da[N],ln[N];
ll ksm(ll a,ll b)
{
ll ans=1;
while(b){
if(b&1) ans=ans*a%mod;
a=a*a%mod;
b>>=1;
}
return ans;
}
void change(ll num[])
{
for(int i=1,j=len/2;i<len-1;i++){
if(i<j) swap(num[i],num[j]);
int k=len/2;
while(j>=k) j-=k,k>>=1;
if(j<k) j+=k;
}
}
void ntt(ll num[],int fl)
{
for(int i=2;i<=len;i<<=1){
if(fl==1) wn=ksm(g,(mod-1)/i);
else wn=ksm(g,mod-1-(mod-1)/i);
for(int j=0;j<len;j+=i){
w=1;
for(int k=j;k<j+i/2;k++){
ll u=w*num[k+i/2]%mod,t=num[k];
num[k]=(t+u)%mod;
num[k+i/2]=(t-u+mod)%mod;
w=w*wn%mod;
}
}
}
if(fl==-1){
ll inv=ksm(len,mod-2);
for(int i=0;i<len;i++) num[i]=num[i]*inv%mod;
}
}
void getinv(int n,ll a[],ll ans[])
{
if(n==1){ans[0]=ksm(a[0],mod-2);return;}
getinv((n+1)/2,a,ans);
len=1;
while(len<2*n) len*=2;
for(int i=0;i<n;i++) a1[i]=a[i];
for(int i=n;i<len;i++) a1[i]=0;
change(a1),change(ans);
ntt(a1,1),ntt(ans,1);
for(int i=0;i<len;i++) ans[i]=ans[i]*(2-ans[i]*a1[i]%mod+mod)%mod;
change(ans),ntt(ans,-1);
for(int i=n;i<len;i++) ans[i]=0;
}
void getln(int n,ll a[],ll ln[])
{
for(int i=1;i<n;i++) da[i-1]=a[i]*i;
getinv(n,a,inva);
len=1;
while(len<2*n) len*=2;
change(da),change(inva);
ntt(da,1),ntt(inva,1);
for(int i=0;i<len;i++) ln[i]=da[i]*inva[i]%mod;
change(ln),ntt(ln,-1);
for(int i=len-1;i>=0;i--) ln[i+1]=ksm(i+1,mod-2)*ln[i]%mod;
ln[0]=0;
}
int main()
{
scanf("%d",&n);
for(int i=0;i<n;i++) scanf("%lld",&a[i]),a[i]%=mod;
getln(n,a,ln);
for(int i=0;i<n;i++) printf("%lld ",ln[i]);
}