[atARC061F]Card Game for Three

记录每一次操作的玩家为操作序列（去掉第一次），需要满足：$a$的个数为$n$且以$a$为结尾，$b$和$c$的个数分别不超过$m$和$k$

其所对应的概率：每一个字符恰好确定一张卡牌，因此即$3^{n+m+k-|s|}$

暴力枚举$b$和$c$的个数，即$sum_{i=0}^{m}sum_{j=0}^{k}{n+i+j-1choose i}{n+j-1choose j}3^{m+k-i-j}$

将组合数用阶乘展开，即$sum_{i=0}^{m}sum_{j=0}^{k}frac{(n+i+j-1)!3^{m+k-i-j}}{(n-1)!i!j!}$

枚举$i+j=s$，即$frac{sum_{s=0}^{m+k}(n+s-1)!3^{m+k-s}sum_{0le ile m,0le s-ile k}frac{1}{i!(s-i)!}}{(n-1)!}$

关于后者显然是一个多项式乘法的形式，但由于模数是$10^{9}+7$，需要写一个拆系数fft/三模数ntt，复杂度为$o(nlog n)$且常数较大

#include<bits/stdc++.h>
using namespace std;
#define N 300005
#define M (1<<20)
#define mod 1000000007
#define K 40000
#define cd complex<double>
#define PI acos(-1.0)
int n,m,k,ans,fac[M],inv[N],mi[M],rev[M];
cd a1[M],a2[M],b1[M],b2[M],c1[M],c2[M],c3[M],w[M];
void fft(cd *a,int p){
    for(int i=0;i<M;i++)
        if (i<rev[i])swap(a[i],a[rev[i]]);
    for(int i=2;i<=M;i<<=1)
        for(int j=0;j<M;j+=i)
            for(int k=0;k<(i>>1);k++){
                cd s=w[(i>>1)+k];
                if (p)s=conj(s);
                cd x=a[j+k],y=a[j+k+(i>>1)]*s;
                a[j+k]=x+y;
                a[j+k+(i>>1)]=x-y;
            }
    if (p){
        for(int i=0;i<M;i++)a[i]/=M;
    }
}
int main(){
    fac[0]=inv[0]=inv[1]=mi[0]=1;
    for(int i=1;i<M;i++)fac[i]=1LL*fac[i-1]*i%mod;
    for(int i=2;i<N;i++)inv[i]=1LL*(mod-mod/i)*inv[mod%i]%mod;
    for(int i=2;i<N;i++)inv[i]=1LL*inv[i-1]*inv[i]%mod;
    for(int i=1;i<M;i++)mi[i]=3LL*mi[i-1]%mod;
    scanf("%d%d%d",&n,&m,&k);
    for(int i=0;i<M;i++)rev[i]=(rev[i>>1]>>1)+((i&1)*(M>>1));
    for(int i=1;i<M;i<<=1)
        for(int j=0;j<i;j++)w[i+j]=cd(cos(j*PI/i),sin(j*PI/i));
    for(int i=0;i<=m;i++)a1[i]=cd(inv[i]/K,0);
    for(int i=0;i<=m;i++)a2[i]=cd(inv[i]%K,0);
    for(int i=0;i<=k;i++)b1[i]=cd(inv[i]/K,0);
    for(int i=0;i<=k;i++)b2[i]=cd(inv[i]%K,0);
    fft(a1,0);
    fft(a2,0);
    fft(b1,0);
    fft(b2,0);
    for(int i=0;i<M;i++){
        c1[i]=a1[i]*b1[i];
        c2[i]=a1[i]*b2[i]+a2[i]*b1[i];
        c3[i]=a2[i]*b2[i];
    }
    fft(c1,1);
    fft(c2,1);
    fft(c3,1);
    for(int i=0;i<=m+k;i++){
        int s1=llround(c1[i].real())%mod;
        int s2=llround(c2[i].real())%mod;
        int s3=llround(c3[i].real())%mod;
        int s=(1LL*K*K%mod*s1+1LL*K*s2+s3)%mod;
        ans=(ans+1LL*fac[n+i-1]*mi[m+k-i]%mod*s)%mod;
    }
    ans=1LL*ans*inv[n-1]%mod;
    printf("%d",ans);
}

事实上还有更好的做法，由于$frac{1}{i!(s-i)!}=frac{{schoose i}}{s!}$，代入即$frac{sum_{s=0}^{m+k}s!(n+s-1)!3^{m+k-s}sum_{0le ile m,0le s-ile k}{schoose i}}{(n-1)!}$

考虑后者，也就是$sum_{i=l}^{r}{schoose i}$，定义为$f(s,l,r)$，我们将这一行的每一个数乘以2（即有2个），在将非首尾的数利用杨辉三角合并，即$2f(s,l,r)=f(s,l+1,r)+{schoose l}+{schoose r}$

利用这个式子，注意到每一次$s$增加区间两端点变化为为$o(1)$，就可以线性的完成，复杂度为$o(n)$

#include<bits/stdc++.h>
using namespace std;
#define N 1000005
#define mod 1000000007
int n,m,k,ans,fac[N],inv[N],mi[N];
int c(int n,int m){
    if (n<m)return 0;
    return 1LL*fac[n]*inv[m]%mod*inv[n-m]%mod;
}
int main(){
    fac[0]=inv[0]=inv[1]=mi[0]=1;
    for(int i=1;i<N;i++)fac[i]=1LL*fac[i-1]*i%mod;
    for(int i=2;i<N;i++)inv[i]=1LL*(mod-mod/i)*inv[mod%i]%mod;
    for(int i=2;i<N;i++)inv[i]=1LL*inv[i-1]*inv[i]%mod;
    for(int i=1;i<N;i++)mi[i]=3LL*mi[i-1]%mod;
    scanf("%d%d%d",&n,&m,&k);
    //l=max(s-k,0),r=min(m,s)
    int s=1,l=0,r=0;
    for(int i=0;i<=m+k;i++){
        if (i){
            s=((2*s%mod+mod-c(i-1,l))%mod+mod-c(i-1,r))%mod;
            l++;
            while (max(i-k,0)<l)s=(s+c(i,--l))%mod;
            while (l<max(i-k,0))s=(s+c(i,l++))%mod;
            while (r<min(m,i))s=(s+c(i,++r))%mod;
            while (min(m,i)<r)s=(s+c(i,r--))%mod;
        }
        ans=(ans+1LL*fac[n+i-1]*mi[m+k-i]%mod*s%mod*inv[i])%mod;
    }
    ans=1LL*ans*inv[n-1]%mod;
    printf("%d",ans);
}