[CTS2019]随机立方体（容斥+组合数学）

这题七次方做法显然，但由于我太菜了，想了一会发现也就只会这么多，而且别的毫无头绪。发现直接做不行，那么，容斥！

f[i]为至少i个极值的方案，然后这里需要一些辅助变量，a[i]表示选出i个三维坐标均不相同的i个极大值的方案数，g[i]表示i个极大的数任意一个至少有一维坐标相同的点的个数，h[i]表示g[i]的极值可以同时存在的方案数，那么有f[i]=C(nml,g[i])a[i]h[i](nml-g[i])!。

a[i]很容易求得，就是(∏(n-j)(m-j)(l-j))/i!，其中j∈[0,i)，g[i]更好求，就是nml-(n-i)(m-i)(l-i)

然后要进行一些关于上升幂的运算，我这里打不出式子（因为太菜了不会LaTeX），所以就不打了。注意维护g[i]的前缀积，具体细节看code吧。

#include<bits/stdc++.h>

using namespace std;

const int N=5e6+,mod=;

int n,m,l,k,mn,ans,fac[N],inv[N],a[N],f[N],g[N],h[N],pre[N];

int calc(int x){return 1ll*(n-x)*(m-x)%mod*(l-x)%mod;}

int qpow(int a,int b)

{

    int ret=;

    while(b)

    {

        if(b&)ret=1ll*ret*a%mod;

        a=1ll*a*a%mod,b>>=;

    }

    return ret;

}

int C(int a,int b){return 1ll*fac[a]*inv[b]%mod*inv[a-b]%mod;}

int main()

{

    inv[]=inv[]=fac[]=;for(int i=;i<N;i++)inv[i]=1ll*(mod-mod/i)*inv[mod%i]%mod;

    for(int i=;i<N;i++)fac[i]=1ll*fac[i-]*i%mod,inv[i]=1ll*inv[i-]*inv[i]%mod;

    int T;scanf("%d",&T);

    while(T--)

    {

        scanf("%d%d%d%d",&n,&m,&l,&k);

        mn=min(n,min(m,l));

        a[]=h[]=pre[]=;

        for(int i=;i<=mn;i++)a[i]=1ll*a[i-]*calc(i-)%mod;

        for(int i=;i<=mn;i++)g[i]=(1ll*g[i-]+calc(i-)-calc(i)+mod)%mod;

        for(int i=;i<=mn;i++)pre[i]=1ll*pre[i-]*g[i]%mod;

        pre[mn]=qpow(pre[mn],mod-);

        for(int i=mn;i;i--)pre[i-]=1ll*pre[i]*g[i]%mod;

        for(int i=;i<=mn;i++)f[i]=1ll*a[i]*pre[i]%mod;

        ans=;

        for(int i=k;i<=mn;i++)

        if((i-k)&)ans=(ans-1ll*C(i,k)*f[i]%mod+mod)%mod;

        else ans=(ans+1ll*C(i,k)*f[i])%mod;

        printf("%d\n",ans);

    }

}