多项式求ln

求

两边求导

所以对于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(ch)){
		ch=getchar();if(ch=='-')f=-1;
	}
	while(isdigit(ch)){
		cnt=cnt*10+ch-48;
		ch=getchar();
	}return cnt*f;
}
int a[2000003],b[2000003],c[2000003],r[2000003],d[2000003],e[2000003];
const int p=998244353;
int ksm(int x,int y){
	int sum=1;
	while(y){
		if(y&1)sum=sum*x%p;
		x=x*x%p;y>>=1;
	}return sum;
}
void NTT(int *A,int len,int inv){
	for(int i=0;i<len;i++){
		if(i<r[i])swap(A[i],A[r[i]]);
	}
	for(int mid=1;mid<len;mid<<=1){
		int wn=ksm(3,(p-1)/(mid<<1));
		for(int i=0;i<len;i+=mid*2){
			int omega=1;
			for(int j=0;j<mid;j++,omega=omega*wn%p){
				int x=A[i+j],y=A[i+j+mid]*omega%p;
				A[i+j]=(x+y)%p;A[i+j+mid]=((x-y)%p+p)%p;
			}
		}
	}
	if(inv==1)return;
	int gu=ksm(len,p-2);reverse(A+1,A+len);
	for(int i=0;i<len;i++){
		A[i]=A[i]*gu%p;
	}
}
void qiudao(int *A,int *B,int n){
	for(int i=1;i<n;i++){
		B[i-1]=A[i]*i%p;
	}B[n-1]=0;
}
void jifen(int *A,int *B,int n){
	for(int i=1;i<n;i++){
		B[i]=A[i-1]*ksm(i,p-2)%p;
	}B[0]=0;
}
void getinv(int *A,int *B,int deg){
	if(deg==1){
		B[0]=ksm(A[0],p-2);return;
	}
	getinv(A,B,(deg+1)>>1);
	int limit=1,l=0;
	while(limit<(deg<<1))limit<<=1,l++;
	for(int i=0;i<=limit;i++)r[i]=(r[i>>1]>>1)|((i&1)<<(l-1));
	for(int i=0;i<deg;i++)d[i]=a[i];for(int i=deg;i<limit;i++)d[i]=0;
	NTT(d,limit,1);NTT(B,limit,1);
	for(int i=0;i<limit;i++)B[i]=((2-d[i]*B[i]%p)+p)%p*B[i]%p;
	NTT(B,limit,-1);for(int i=deg;i<limit;i++)B[i]=0;
}
void solve(int n){
	qiudao(a,b,n);
	//for(int i=0;i<n;i++)cout<<b[i]<<" ";cout<<endl;//cin>>b[333];
	getinv(a,c,n);
	//for(int i=0;i<n;i++)cout<<c[i]<<" ";cout<<endl;
	int limit=1,l=0;while(limit<=n)limit<<=1,l++;
	for(int i=1;i<=n;i++)r[i]=(r[i>>1]>>1)|((i&1)<<(l-1));
	//for(int i=0;i<limit;i++)cout<<r[i]<<" ";cout<<endl;
	NTT(b,limit,1);NTT(c,limit,1);
	//for(int i=0;i<n;i++)cout<<b[i]<<" ";cout<<endl;
	//for(int i=0;i<n;i++)cout<<c[i]<<" ";cout<<endl;
	for(int i=0;i<limit;i++)b[i]=b[i]*c[i]%p;
	//for(int i=0;i<limit;i++)cout<<b[i]<<" ";cout<<endl;
	NTT(b,limit,-1);jifen(b,e,n);
	
}
signed main(){
	int n=in;
	for(int i=0;i<n;i++)a[i]=in;int m;
	for(m=1;m<=n;m<<=1);solve(m);
	for(int i=0;i<n;i++)printf("%lld ",e[i]);
	return 0;
}