Description
给定 G ( x ) G(x) G(x),求 F ( x ) ≡ e G ( x ) ( m o d x n ) F(x)\equiv e^{G(x)}\pmod {x^n} F(x)≡eG(x)(modxn)。
Solution
预备知识:牛顿迭代
求出 F ( x ) F(x) F(x)使得 G ( F ( x ) ) ≡ 0 ( m o d x n ) G(F(x))\equiv 0\pmod {x^n} G(F(x))≡0(modxn)。
假设我们已知 F n / 2 F_{n/2} Fn/2使得 G ( F n / 2 ( x ) ) = 0 G(F_{n/2}(x))=0 G(Fn/2(x))=0,我们在 F n / 2 F_{n/2} Fn/2泰勒展开,得:
G ( F n ( x ) ) ≡ ∑ i = 0 + ∞ G [ i ] ( F n / 2 ( x ) ) ( F n ( x ) − F n / 2 ( x ) ) i i ! ≡ 0 ( m o d x n ) G(F_n(x))\equiv \sum_{i=0}^{+\infty}\frac{G^{[i]}(F_{n/2}(x))(F_n(x)-F_{n/2}(x))^i}{i!}\equiv 0 \pmod {x^n} G(Fn(x))≡i=0∑+∞i!G[i](Fn/2(x))(Fn(x)−Fn/2(x))i≡0(modxn)
由于 i = 2 ⋯ + ∞ i=2\dots +\infty i=2⋯+∞时, ( F n ( x ) − F n / 2 ( x ) ) 2   m o d   x n (F_n(x)-F_{n/2}(x))^2 \bmod x^n (Fn(x)−Fn/2(x))2modxn一定为 0 0 0,我们解得:
F n ( x ) ≡ F n / 2 ( x ) − G ( F n / 2 ( x ) ) G ′ ( F n / 2 ( x ) ) F_n(x)\equiv F_{n/2}(x)-\frac{G(F_{n/2}(x))}{G'(F_{n/2}(x))} Fn(x)≡Fn/2(x)−G′(Fn/2(x))G(Fn/2(x))
于是我们用牛顿迭代就可以做多项式求逆、开方、exp等很多事情辣~
多项式exp
将式子化一下:
ln
(
F
(
x
)
)
≡
G
(
x
)
\ln (F(x))\equiv G(x)
ln(F(x))≡G(x)
代入牛顿迭代的式子,得:
F
n
(
x
)
≡
F
n
/
2
(
x
)
(
G
(
x
)
−
ln
(
F
n
/
2
(
x
)
)
+
1
)
F_n(x)\equiv F_{n/2}(x)(G(x)-\ln(F_{n/2}(x))+1)
Fn(x)≡Fn/2(x)(G(x)−ln(Fn/2(x))+1)
倍增计算即可。
Code
/************************************************
* Au: Hany01
* Date: Jul 28th, 2018
* Prob: Luogu4726 多项式exp
* Email: hany01@foxmail.com
* Inst: Yali High School
************************************************/
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
#define rep(i, j) for (register int i = 0, i##_end_ = (j); i < i##_end_; ++ i)
#define For(i, j, k) for (register int i = (j), i##_end_ = (k); i <= i##_end_; ++ i)
#define Fordown(i, j, k) for (register int i = (j), i##_end_ = (k); i >= i##_end_; -- i)
#define Set(a, b) memset(a, b, sizeof(a))
#define Cpy(a, b) memcpy(a, b, sizeof(a))
#define x first
#define y second
#define pb(a) push_back(a)
#define mp(a, b) make_pair(a, b)
#define SZ(a) ((int)(a).size())
#define INF (0x3f3f3f3f)
#define INF1 (2139062143)
#define debug(...) fprintf(stderr, __VA_ARGS__)
#define y1 wozenmezhemecaia
template <typename T> inline bool chkmax(T &a, T b) { return a < b ? a = b, 1 : 0; }
template <typename T> inline bool chkmin(T &a, T b) { return b < a ? a = b, 1 : 0; }
inline int read() {
static int _, __; static char c_;
for (_ = 0, __ = 1, c_ = getchar(); c_ < '0' || c_ > '9'; c_ = getchar()) if (c_ == '-') __ = -1;
for ( ; c_ >= '0' && c_ <= '9'; c_ = getchar()) _ = (_ << 1) + (_ << 3) + (c_ ^ 48);
return _ * __;
}
const int MOD = 998244353, g0 = 3, maxn = 1e5 + 5;
int rev[maxn << 2], powg[maxn << 2], ipowg[maxn << 2], ig0, invn;
inline int ad(int x, int y) { if ((x += y) >= MOD) return x - MOD; return x; }
inline LL Pow(LL a, LL b) {
static LL Ans;
for (Ans = 1; b; b >>= 1, (a *= a) %= MOD) if (b & 1) (Ans *= a) %= MOD;
return Ans;
}
inline void Init(int N) {
ig0 = Pow(g0, MOD - 2);
for (register int i = 1; i <= N; i <<= 1)
powg[i] = Pow(g0, (MOD - 1) / i), ipowg[i] = Pow(ig0, (MOD - 1) / i);
}
inline void pre(int n) {
static int N, cnt;
for (N = 1, cnt = 0; N < n; N <<= 1, ++ cnt);
rep(i, n) rev[i] = (rev[i >> 1] >> 1) | ((i & 1) << (cnt - 1));
invn = Pow(n, MOD - 2);
}
inline void NTT(int* a, int n, int ty) {
rep(i, n) if (i < rev[i]) swap(a[i], a[rev[i]]);
for (register int i = 2, p = 1; i <= n; p = i, i <<= 1) {
register int w0 = ty > 0 ? powg[i] : ipowg[i];
for (register int j = 0; j < n; j += i) {
register int w = 1;
rep(k, p) {
register int x = a[j + k], y = (LL)a[j + k + p] * w % MOD;
a[j + k] = ad(x, y), a[j + k + p] = ad(x, MOD - y);
w = (LL)w * w0 % MOD;
}
}
}
if (ty < 1) rep(i, n) a[i] = (LL)a[i] * invn % MOD;
}
inline void Mult(int* a, int* b, int* c, int n) {
static int A[maxn << 2], B[maxn << 2];
rep(i, n) A[i] = a[i], B[i] = b[i];
pre(n), NTT(A, n, 1), NTT(B, n, 1);
rep(i, n) A[i] = (LL)A[i] * B[i] % MOD;
NTT(A, n, -1);
rep(i, n) c[i] = A[i];
}
namespace Inv {
static int A[maxn << 2], B[maxn << 2], a[maxn << 2];
void Inv_(int* b, int n) {
if (n == 1) { b[0] = Pow(a[0], MOD - 2); return; }
Inv_(b, n >> 1);
rep(i, n) A[i] = a[i], B[i] = b[i];
For(i, n, n << 1) A[i] = B[i] = 0;
pre(n << 1), NTT(A, n << 1, 1), NTT(B, n << 1, 1);
rep(i, n << 1) A[i] = (LL)A[i] * B[i] % MOD * B[i] % MOD;
NTT(A, n << 1, -1);
rep(i, n) b[i] = ad(ad(b[i], b[i]), MOD - A[i]);
}
inline void Inv(int* x, int* y, int n) {
rep(i, n) a[i] = x[i], y[i] = 0;
Inv_(y, n);
}
}
inline void Int(int* a, int* b, int n) {
Fordown(i, n - 1, 0) b[i + 1] = Pow(i + 1, MOD - 2) * a[i] % MOD;//downto!!!
b[0] = 0;
}
inline void Der(int *a, int* b, int n) {
rep(i, n) b[i] = (LL)a[i + 1] * (i + 1) % MOD;
}
inline void Ln(int* a, int* b, int n) {
static int A[maxn << 2], B[maxn << 2];
rep(i, n) A[i] = B[i] = a[i];
For(i, n, n << 1) A[i] = B[i] = 0;
Inv:: Inv(A, A, n), Der(B, B, n << 1);
Mult(A, B, A, n << 1), Int(A, A, n << 1);
rep(i, n) b[i] = A[i];
}
namespace Exp {
static int a[maxn << 2], A[maxn << 2];
void Exp_(int* b, int n) {
if (n == 1) { b[0] = 1; return; }
Exp_(b, n >> 1), Ln(b, A, n);
rep(i, n) A[i] = ad(a[i], MOD - A[i]);
++ A[0], Mult(b, A, b, n << 1);
}
void Exp(int* x, int* y, int n) {
rep(i, n) a[i] = x[i], y[i] = 0;
Exp_(y, n);
}
}
int main()
{
#ifdef hany01
freopen("luogu4726.in", "r", stdin);
freopen("luogu4726.out", "w", stdout);
#endif
static int f[maxn << 2], n, N;
for (n = read(), N = 1; N < n; N <<= 1);
Init(N << 1);
rep(i, n) f[i] = read();
Exp:: Exp(f, f, N);
rep(i, n) printf("%d ", f[i]);
return 0;
}