【NTT模板】

没错，这就是我在网上找的一份NTT的板子…

使用方法如下：
多项式乘法模板

求多项式乘法，
A数组存第一个多项式的系数
B数组存第二个多项式的系数
也就是读入啦！
然后直接穿solve(A,B)，最后答案就会在A里面啦！
直接输出就好了。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 5e6+7;
const int mod = 998244353;
const int G = 3;
int n,m,L,R[maxn];
int A[maxn],B[maxn];
ll qpow(ll a,ll b)
{
    ll ans = 1;
    while(b>0)
    {
    	if(b&1) ans = ans*a%mod;
    	b>>=1;
    	a = a*a%mod;
	}
	return ans%mod;
}
void NTT(int *a,int f)
{
    for(int i = 0; i < n; i++) 
	{
		if(i < R[i]) swap(a[i],a[R[i]]);
	}
    for(int i = 1; i < n; i <<= 1)
	{
        ll gn = qpow(G,(mod - 1) / (i << 1));
        for(int j = 0; j < n; j += (i << 1))
		{
            ll g = 1;
            for(int k = 0; k < i; k++, g = g * gn % mod)
			{
                int x = a[j + k], y = g * a[j + k + i] % mod;
                a[j + k] = (x + y) % mod; 
				a[j + k + i] = (x - y + mod) % mod;
            }
        }
    }
    if(f==1) return;
    int inv = qpow(n,mod - 2); 
	reverse(a + 1,a + n);
    for(int i = 0; i < n; i++) a[i] = 1ll * a[i] * inv % mod;
}
void solve(int *A,int *B)
{
	m = n + m; 
	for(n = 1; n <= m; n <<= 1) L++;
    for(int i = 0; i < n; i++) R[i] = (R[i >> 1] >> 1) | ((i & 1) << (L - 1));
	NTT(A,1); 
	NTT(B,1);
    for(int i = 0; i < n; i++) A[i] = 1ll * A[i] * B[i] % mod;
    NTT(A,-1);
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i = 0; i <= n; i++) scanf("%d",A+i);
    for(int i = 0; i <= m; i++) scanf("%d",B+i);
    solve(A,B);
    for(int i = 0; i <= m; i++) printf("%d ",A[i]);
    return 0;
}