codeforces 856C Eleventh Birthday

$Description$

$$\begin{gathered} 给定n个数，求将这n个数按拼在一起形成一个大数，这个大数是 \\ 是11的倍数的方案数。 \end{gathered}$$

$Solution$

$对于这题，我们首先得想到一个数学性质$
$10\equiv -1(mod11)$
$$\begin{gathered} 对于这个性质，比如a,b,c三个数拼接成abc形式， \\ abc=a\ast 10^{len(b+c)}+b\ast 10^{len(c)}+c \\ 若满足a\ast 10^{len(b+c)}+b\ast 10^{len(c)}+c\equiv 0(mod11) \\ 则可化为a\ast (-1{)}^{len(b+c)}+b\ast (-1{)}^{len(c)}+c\equiv 0(mod11) \\ 即是a,b,c本身\ast 各自所在的奇偶位是11的倍数时，成立 \\ (举例,a=9,b=2,c=4,组成的数9\ast 1+2\ast -1+4\equiv 0，故924是11的倍数) \end{gathered}$$

$$\begin{gathered} 我们就若设a\ast 1，a在奇数位，a\ast -1,a在偶数位,那么我们现在只需要去找 \\ 哪些数在偶数位，哪些数在奇数位的情况，在数列中，我们插入一个len为 \\ 偶数的数字，不会改变后续的奇偶位情况情况，而插入一个len为奇数的则会全部 \\ 反过来 \end{gathered}$$
$$\begin{gathered} 我们把所有len为奇数的数统计出来，设置这些数在奇数位或偶数位的时候，和为k时 \\ 有多少方案，那么剩下的len为偶数的数，直接插入就好了，不会影响数字本身所在的 \\ 奇偶位,设len为奇数的数有n1个,命名为a_i,len为偶数有n2个，命名为b_i \end{gathered}$$
$设f[i][j][k]为len为奇数中前i个数,有j个数在偶数位时,和为k的方案有多少$
$设g[i][j][k]为len为偶数中前i个数,有j个数在偶数位时,和为k的方案有多少$
$len为奇数的最终状态一定时f[n1][n1/2][0..10],len为偶数可以任意$
状态转移方程为
$f[i][j][k]=f[i-1][j][k-a[i]]+f[i-1][j-1][k+a[i]]$
$g[i][j][k]=g[i-1][j][k-b[i]]+g[i-1][j-1][k+b[i]]$
$$\begin{gathered} 而f[n1][n1/2][0..10]对于每一个只是表示n1/2个偶数下,但对于每个具体位置,n1/2个位于偶数位的数和n1-n1/2个位于奇数位的数是可以互相换位的,故还需乘以 \\ (n1/2)!\ast (n1-n1/2)! \\ 再思考，对于len为偶数中的情况，g[n2][j][k]中，有j个数在偶数位，那么这偶数位 \\ 得位于n1-n1/2个数之后，变成了一个组合数问题,j个球放在n1-n1/2个盒子里，一个 \\ 盒子可放任意多个球，盒子可以为空，即是多重集组合数，j个球可任意交换， \\ 则再乘j!,对于n2-j个放在奇数为得数同理，只是奇数位多了一个第一位可以放, \\ 则是n2-j个球放在n1/2+1个盒子 \\ OK,完全solved!! \end{gathered}$$

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
//#pragma GCC optimize(3, "Ofast", "inline")
#define fir(i, a, b) for (int i = a; i <= b; i++)
#define mov(i) (1ll << i)
#define MAX_INF 0x3f3f3f3f
ll qpow(ll a, ll b, ll p)
{
    ll ans = 1;
    while (b)
    {
        if (b & 1)
            ans = ans * a % p;
        a = a * a % p;
        b >>= 1;
    }
    return ans;
}
const int maxn = 2e3 + 10;
int a[maxn],b[maxn];
int mod = 998244353;
int f[2][maxn][12],g[2][maxn][12];//直接开数组MLE，状态压缩以下
inline int add(ll a,ll b)
{
    return (a+b+mod)%mod;
}
ll fac[maxn], invfac[maxn], t[maxn];
void init(int n, ll MOD)
{
    fac[0] = 1;
    mod = MOD;
    for (int i = 1; i <= n; i++)
        fac[i] = fac[i - 1] * i % mod;
    invfac[n] = qpow(fac[n], mod - 2, mod);
    for (int i = n; i >= 1; i--)
        invfac[i - 1] = invfac[i] * i % mod;
}
ll C(ll n, ll m)//求组合数
{
    return n >= m ? fac[n] * invfac[n - m] % mod * invfac[m] % mod : 0;
}
//多重集排列数
ll cal(ll n,ll m)
{
    if(!m) return (n==0);//n,m可能均为0
    return fac[n]*C(n+m-1,m-1)%mod;
}
int main()
{
  //  input();
    int T;
    init(maxn-1,mod);
    fac[0]=1;
    cin >> T; 
    while(T--)
    {
        int n;
        cin>>n;
        int cnta=0,cntb=0;
        fir(i,1,n)
        {
            int x;
            scanf("%d",&x);
            if(to_string(x).size()&1)
            a[++cnta]=x%11;
            else 
            b[++cntb]=x%11;
        }
        memset(f,0,sizeof(f));
        memset(g,0,sizeof(g));
        f[0][0][0]=1;g[0][0][0]=1;
        for(int i=1;i<=cnta;i++)
        {
            int curi=i&1,lasti=!(i&1);
            for(int j=0;j<=i;j++)
            {
                fir(k,0,10)
                {
                    f[curi][j][k]=add(f[curi][j][k],f[lasti][j][(k-a[i]+11)%11]);
                    if(j)
                    f[curi][j][k]=add(f[curi][j][k],f[lasti][j-1][(k+a[i])%11]);
                }
            }
            memset(f[lasti],0,sizeof(f[lasti]));
        }
        
        for(int i=1;i<=cntb;i++)
        {
            int curi=i&1,lasti=!(i&1);
            for(int j=0;j<=i;j++)
            {
                fir(k,0,10)
                {
                    g[curi][j][k]=add(g[curi][j][k],g[lasti][j][(k-b[i]+11)%11]);
                    if(j)
                    g[curi][j][k]=add(g[curi][j][k],g[lasti][j-1][(k+b[i])%11]);
                }
            }
            memset(g[lasti],0,sizeof(g[lasti]));
        }
        ll ans=0;
        ll temp=fac[cnta/2]*fac[cnta-cnta/2]%mod;
        int m=cnta/2+1;
       // cout<<"temp:"<<temp<<endl;
        for(int i=0;i<=cntb;i++)
        {
            for(int j=0;j<11;j++)
            {
                int p=cntb-i;
                ll val=(ll)g[cntb&1][i][j]*f[cnta&1][cnta/2][(11-j)%11]%mod*temp%mod*cal(p,m)%mod*cal(i,cnta-cnta/2)%mod;
                ans=(ans+val)%mod;
            }
        }
        printf("%lld\n",ans);
    }
}