這題是真的神仙啊...居然用的 stl 來卡常？

話說 998244353 真的可以一眼 NTT ？

noteskey

所以說只要推柿子就好了但是有的地方的推導根本就想不到...

我們令第 t 個答案為 \(ANS_t\over nm\) ，除去 nm 其實就是算期望時要除去的方案數

那麽有：

\[\begin{aligned}{}ANS_t=&\sum_{i=1}^{n} \sum_{j=1}^m (a_i+b_j)^t \\=&\sum_{i=1}^{n}\sum_{j=1}^{m}\sum_{k=0}^{t} \begin{pmatrix} t\\k \end{pmatrix} a_i^{k}b_{j}^{t-k} \\=&\sum_{k=0}^{t} \begin{pmatrix} t\\k \end{pmatrix} \sum_{i=1}^{n}a_i^{k} \sum_{j=1}^{m} b_{j}^{t-k} \\=& \sum_{k=0}^{t} \begin{pmatrix} t\\k \end{pmatrix} \Big( \sum_{i=1}^{n}a_i^{k} \Big) \Big(\sum_{j=1}^{m} b_{j}^{t-k} \Big) \\=&\sum_{k=0}^{t} t! \Big({ \sum_{i=1}^{n}a_i^{k}\over k! } \Big) \Big( {\sum_{j=1}^{m} b_{j}^{t-k} \over (t-k)!} \Big) \end{aligned}\]

那麽我們發現這個東西就是卷積的形式了，於是我們的任務就是快速算出後面兩個表達式所對應的的多項式

我們設：

\[\begin{aligned} f(k)=\sum_{i=1}^n a_i ^{k} \\ g(k)=\sum_{i=1}^{m}b_{i}^{k} \end{aligned}\]

那麽我們只要求出 \(f(1...t) ,g(1...t)\) 即可

但是求出來沒有這麽簡單...

我們考慮先求出：

\[f_j(x)=\sum_{i=1}^{k}a_j^{i}\]

然後我們令：

\[F(x)= \sum_{i=1}^{n} f_{i}(x)\]

那麽 \(F(x)\) 的第 k 項就表示了

\[\sum_{i=1}^{i} a_{i}^{k} x^k\]

即：

\[[x^{k}]F(x)=\sum_{i=1}^{i} a_{i}^{k}\]

那麽我們現在發現要求出 \(f_{j=1...n}(x)\) ，但是沒有什麽高效的辦法能解決這個問題

於是我們把 \(f_{j}(x)\) 中的上限 k 改掉，變成 \(\infty\) ，這樣我們就能用生成函數來解決了：

\[f_{j}(x)={1\over 1- a_i x} \]

然後就是把這些 \(f_j(x)\) 加起來了：

\[F(x)={1\over 1-a_1x} +{1\over 1-a_2x} + ...+{1\over 1-a_nx}\]

最後就是怎麽合並這些 \(f_j(x)\) 的問題了，我們發現：

\[ln(1-a_ix)'={-a_i\over 1-a_ix}\]

發現這個東西和我們要求的東西蠻像的，那麽這裏又有：

\[\begin{aligned}G(x)=&\sum_{i=1}^{n} ln(1-a_ix)' \\ F(x)=&-x·G(x)+n \end{aligned}\]

那麽我們的任務就已經是求出 G(x) 了

\[G(x)=\sum_{i=1}^{n} ln(1-a_ix)'= ln(\prod_{i=1}^n (1-a_ix))'\]

然後我們就要在優秀的時間內求出\(\prod_{i=1}^n (1-a_ix)\) ，然後取個 Ln 再求個 Inv ，右移一位讓最低項為 n ， \(F(x)\) 就構造好了

至於對 b 數組也是同樣的求法，兩個要求的多項式求出來後，我們還要讓兩個多項式除去對應的階乘，然後再去卷積

卷完積後也和最開始的式子一樣，第 t 項乘上 t! ，然後別忘了我們還要除去 nm 才是最終的答案...

最後就是順序輸出就好了

code

這沒救了的壓行【霧

//by Judge
#include<bits/stdc++.h>
#define Rg register
#define fp(i,a,b) for(Rg int i=(a),I=(b)+1;i<I;++i)
#define fd(i,a,b) for(Rg int i=(a),I=(b)-1;i>I;--i)
#define ll long long
using namespace std;
const int mod=998244353;
const int iG=332748118;
const int M=524288+3;
typedef int arr[M];
#ifndef Judge
#define getchar() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
#endif
char buf[1<<21],*p1=buf,*p2=buf;
inline int inc(int x,int y){return (x+y)%mod;}
inline int mul(int x,int y){return 1ll*x*y%mod;}
inline int dec(int x,int y){return (x-y+mod)%mod;}
inline int read(){ int x=0,f=1; char c=getchar();
    for(;!isdigit(c);c=getchar()) if(c=='-') f=-1;
    for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;
} char sr[1<<21],z[20];int CCF=-1,Z;
inline void Ot(){fwrite(sr,1,CCF+1,stdout),CCF=-1;}
inline void print(int x,char chr='\n'){
    if(CCF>1<<20)Ot();if(x<0)sr[++CCF]=45,x=-x;
    while(z[++Z]=x%10+48,x/=10);
    while(sr[++CCF]=z[Z],--Z);sr[++CCF]=chr;
} int n,m,t,qq,rev,d,limit;
arr a,b,f,g,fac,ifac,lg,r[20],D[21],G[2];
inline int qpow(Rg int x,Rg int p=mod-2,Rg int s=1){
    for(;p;p>>=1,x=mul(x,x)) if(p&1) s=mul(s,x); return s;
}
inline void prep(int n=262143){ for(limit=1;limit<=n;limit<<=1);
    fac[0]=fac[1]=1; fp(i,2,limit) fac[i]=mul(fac[i-1],i);
    fp(d,1,19){ lg[1<<d]=d; fp(i,0,(1<<d)-1)
        r[d][i]=(r[d][i>>1]>>1)|((i&1)<<(d-1));
    }
    for(Rg int t=(mod-1)>>1,i=1,x,y;i<limit;i<<=1,t>>=1){
        x=qpow(3,t),y=qpow(iG,t),G[0][i]=G[1][i]=1;
        fp(k,1,i-1) G[1][i+k]=mul(G[1][i+k-1],x);
        fp(k,1,i-1) G[0][i+k]=mul(G[0][i+k-1],y);
    }
}
inline void init(int n){
    for(limit=1;limit<n;limit<<=1); d=lg[limit];
}
inline void NTT(int* a,int tp){
    fp(i,0,limit-1) if(i<r[d][i]) swap(a[i],a[r[d][i]]);
    for(Rg int mid=1,I=2;mid<limit;mid<<=1,I<<=1)
        for(Rg int j=0,x,y;j<limit;j+=I) fp(k,0,mid-1)
            x=a[j+k],y=mul(a[j+k+mid],G[tp][mid+k]),
            a[j+k]=inc(x,y),a[j+k+mid]=dec(x,y);
    if(tp) return ; int inv=qpow(limit);
    fp(i,0,limit-1) a[i]=mul(a[i],inv);
}
inline void Mul(int* a,int* b,int* c,int n){
    static arr A,B; fp(i,0,n-1) A[i]=a[i],B[i]=b[i];
    fill(A+n,A+limit,0),fill(B+n,B+limit,0),NTT(A,1),NTT(B,1);
    fp(i,0,limit-1) A[i]=mul(A[i],B[i]); NTT(A,0),copy(A,A+limit,c);
}
void Inv(int* a,int* b,int n){
    if(n==1) return b[0]=qpow(a[0]),void();
    Inv(a,b,n>>1),init(n<<1); static arr c,d;
    fp(i,0,n-1) c[i]=a[i],d[i]=b[i];
    fp(i,n,limit-1) c[i]=d[i]=0; NTT(c,1),NTT(d,1);
    fp(i,0,limit-1) c[i]=mul(c[i],mul(d[i],d[i])); NTT(c,0);
    fp(i,0,n-1) b[i]=dec(inc(b[i],b[i]),c[i]);
}
inline void Derv(int* a,int* b,int n){
    fp(i,1,n-1) b[i-1]=mul(a[i],i); b[n-1]=0;
}
inline void Ln(int* a,int* b,int n){ static arr c,d;
    memset(c,0,(n+1)<<2),memset(d,0,(n+1)<<2);
    Derv(a,c,n),Inv(a,d,n),init(n<<1);
    Mul(c,d,b,n),fill(b+n,b+limit,0);
}
inline void solv(int* a,int d,int l,int r){
    if(l==r) return D[d][0]=1,D[d][1]=mod-a[l],void(); static arr A,B;
    int mid=(l+r)>>1; solv(a,d,l,mid),solv(a,d+1,mid+1,r),init(r-l+2);
    fp(i,0,mid-l+1) A[i]=D[d][i]; fp(i,0,r-mid) B[i]=D[d+1][i];
    fill(A+mid-l+2,A+limit,0),fill(B+r-mid+1,B+limit,0);
    Mul(A,B,D[d],r-l+1),fill(D[d]+r-l+2,D[d]+limit,0);
}
int main(){ n=read(),m=read(),rev=qpow(mul(n,m)),prep();
    fp(i,1,n) a[i]=read(); fp(i,1,m) b[i]=read();
    solv(a,1,1,n),copy(D[1],D[1]+n+1,f);
    solv(b,1,1,m),copy(D[1],D[1]+m+1,g);
    t=read(),qq=max(max(n,m),t); init(qq+1),qq=limit;
    Ln(f,f+1,qq),Ln(g,g+1,qq),f[0]=n,g[0]=m;
    fp(i,1,qq-1) f[i]=dec(mod,f[i]),g[i]=dec(mod,g[i]);
    fac[0]=ifac[0]=ifac[1]=1; fp(i,1,qq-1) fac[i]=mul(fac[i-1],i);
    fp(i,2,qq-1) ifac[i]=mul(mod-mod/i,ifac[mod%i]);
    fp(i,1,qq-1) ifac[i]=mul(ifac[i],ifac[i-1]);
    fp(i,0,qq-1) f[i]=mul(f[i],ifac[i]),g[i]=mul(g[i],ifac[i]);
    fp(i,qq,qq<<1) f[i]=g[i]=0; init(qq<<1),NTT(f,1),NTT(g,1);
    fp(i,0,limit-1) f[i]=mul(f[i],g[i]); NTT(f,0);
    fp(i,1,t) print(mul(f[i],mul(fac[i],rev))); return Ot(),0;
}

題解 P4705 【玩遊戲】