51nod 1514 美妙的序列

题意：

某个1~n的排列如果满足：
在1~n-1 这些位置后面将序列断开，使得总可以从右边找到一个数，并且该数不大于左边的所有数，则称该序列为“美妙的”。
给出n，求长度为n的“美妙的序列”的数量。

题解：

首先第一步，读懂题！换句话说，就是当左边所有数都比右边的最小值小时不合法。
那么就按照题意dp，考虑容斥，枚举从哪个点开始不合法，得到$f[n]=n!-\sum_{i=1}^{n-1}f(i)*(n-i)!$
注意到这是个卷积形式，直接分治fft/多项式求逆即可。
code：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#define LL long long
using namespace std;
const LL p=998244353,yg=3;
LL bin[400000],a[400000],b[400000],c[400000],n;
LL pow(LL a,LL b)
{
    LL ans=1;
    while(b)
    {
        if(b&1) ans=ans*a%p;
        a=a*a%p;b>>=1;
    }
    return ans;
}
void ntt(LL *a,LL n,LL op)
{
    for(LL i=0;i<n;i++) bin[i]=(bin[i>>1]>>1)|((i&1)*(n>>1));
    for(LL i=0;i<n;i++) if(i<bin[i]) swap(a[i],a[bin[i]]);
    for(LL i=1;i<n;i<<=1)
    {
        LL wn=pow(yg,op==1?(p-1)/(i<<1):p-1-(p-1)/(i<<1)),t,w;
        for(LL j=0;j<n;j+=i<<1)
        {
            w=1;
            for(LL k=0;k<i;k++)
            {
                t=a[j+i+k]*w%p;w=w*wn%p;
                a[i+j+k]=(a[j+k]-t+p)%p;a[j+k]=(a[j+k]+t)%p;
            }
        }
    }
    if(op==-1)
    {
        LL Inv=pow(n,p-2);
        for(LL i=0;i<n;i++) a[i]=a[i]*Inv%p;
    }
}
void pol_inv(LL deg,LL *a,LL *b)
{
    if(deg==1) {b[0]=pow(a[0],p-2);return;}
    pol_inv((deg+1)>>1,a,b);
    LL N=1;while(N<deg<<1) N<<=1;
    copy(a,a+deg,c);fill(c+deg,c+N,0);
    ntt(c,N,1);ntt(b,N,1);
    for(LL i=0;i<N;i++)
    {
        b[i]=(2-b[i]*c[i]%p)*b[i]%p;
        b[i]+=b[i]<0?p:0;
    }
    ntt(b,N,-1);fill(b+deg,b+N,0);
}
LL fac[100010],inv[100010];
void pre()
{
    fac[0]=fac[1]=inv[0]=inv[1]=1;
    for(LL i=2;i<=100000;i++) fac[i]=fac[i-1]*i%p,inv[i]=(p-p/i)*inv[p%i]%p;
    for(LL i=2;i<=100000;i++) inv[i]=inv[i-1]*inv[i]%p;
}
int main()
{
    n=100000;
    pre();for(LL i=1;i<=n;i++) a[i]=fac[i];
    pol_inv(n+1,fac,b);
    LL N=1;while(N<(n+1)<<1) N<<=1;
    ntt(a,N,1);ntt(b,N,1);
    for(LL i=0;i<N;i++) a[i]=a[i]*b[i]%p;
    ntt(a,N,-1);
    LL T;scanf("%lld",&T);
    while(T--)
    {
        LL x;scanf("%lld",&x);
        printf("%lld\n",a[x]);
    }
}