Tyvj4879（dp+倍增+NTT）

题面
题意：有n组数，每个数为1到6，每组有m个，每组内不计顺序。问总共至少有x个6的方案数。模998244353。n，m≤400。

看到这个熟悉的数字，大概就是个NTT了，然后考虑哪里有卷积。设g[i]为一组内有i个6的方案数，f[i]为两组内有i个6的方案数，有

$$f[i]=\sum_{j=0}^ig[j]*g[i-j]$$
就是个卷积了。

对于有n组，我们考虑倍增，先求出n/2组的g数组，然后ntt合并即可。

考虑怎么求g[i]。由于每组内无关顺序，根据套路，我们可以强行设定单调。
设h[i][j]为第1~i个数非严格单调，第i个数为j的方案数。字面意思转移即可。

我的直觉认为，在倍增途中，若每一次卷积的次数界都动态变化，则总复杂度可以省一个log。

#include <iostream>
#include <fstream>
#include <algorithm>
#include <cmath>
#include <ctime>
#include <cstdio>
#include <cstdlib>
#include <cstring>

using namespace std;
#define mmst(a, b) memset(a, b, sizeof(a))
#define mmcp(a, b) memcpy(a, b, sizeof(b))

typedef long long LL;

const int N=800400;
const LL p=998244353,g=3;

int n,rev[N];

LL cheng(LL a,LL b)
{
    LL res=1ll;
    for(;b;b>>=1,a=a*a%p)
    if(b&1)
    res=res*a%p;
    return res;
}

void init(int lim)
{
    int k=-1;
    n=1;
    while(n<lim)
    n<<=1,k++;
    for(int i=0;i<n;i++)
    rev[i]=(rev[i>>1]>>1) | ((i&1)<<k);
}

void ntt(LL *a,int ops)
{
    for(int i=0;i<n;i++)
    if(i<rev[i])
    swap(a[i],a[rev[i]]);
    for(int l=2;l<=n;l<<=1)
    {
        int m=(l>>1);
        LL wn;
        if(ops)
        wn=cheng(g,(p-1)/l);
        else
        wn=cheng(g,p-1-(p-1)/l);
        for(int i=0;i<n;i+=l)
        {
            LL w=1ll;
            for(int k=0;k<m;k++)
            {
                LL t=a[i+k+m]*w%p;
                a[i+k+m]=(a[i+k]-t+p)%p;
                a[i+k]=(a[i+k]+t)%p;
                w=w*wn%p;
            }
        }
    }
    if(!ops)
    {
        LL Inv=cheng(n,p-2);
        for(int i=0;i<n;i++)
        a[i]=a[i]*Inv%p;
    }
}

int nn,mm,x,y;
LL f[401][6],gg[N],ans[N],by[N];

void work(int b)
{
    int csj=mm+1;
    for(;b;b>>=1)
    {
        init(csj*2);
        if(b&1)
        {

            ntt(ans,1);

            for(int i=0;i<n;i++)
            by[i]=gg[i];
            ntt(by,1);

            for(int i=0;i<n;i++)
            ans[i]=(ans[i]*by[i])%p;
            ntt(ans,0);
        }

        ntt(gg,1);
        for(int i=0;i<n;i++)
        gg[i]=gg[i]*gg[i]%p;
        ntt(gg,0);
        csj*=2;
    }
}

int main()
{
    cin>>nn>>mm>>x>>y;
    for(int i=1;i<=mm;i++)
    {
        int tu;
        scanf("%d",&tu);
        if(tu==y)
        x--;
    }

    f[0][1]=1ll;
    for(int i=1;i<=mm;i++)
    for(int j=1;j<=5;j++)
    for(int k=j;k<=5;k++)
    f[i][k]=(f[i][k]+f[i-1][j])%p;

    for(int i=0;i<=mm;i++)
    {
        LL hy=0;
        for(int j=1;j<=5;j++)
        hy=(hy+f[i][j])%p;
        gg[mm-i]=hy;
    }

    ans[0]=1ll;

    work(nn);

    LL biu=0;
    for(int i=x;i<n;i++)
    biu=(biu+ans[i])%p;

    cout<<biu<<endl;

    return 0;
}