CTSC2019

（upd：随机立方体AC）

太菜了只会部分分。以后慢慢补坑吧……

随机立方体：

30分：

正常人都能想到的的转移状态（我的确是弱智），从大往小填数，记录有多少个极大值点和三个方向上各占了多少。转移可以计算。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;

#define upd(x,y) x+=y,x%=P

typedef long long i64;

const int N=;
const int P=;

int f[][N][N][N][N];

int pwr(int x,int a){
    int s=;
    while(a){
        if(a&)s=(i64)s*x%P;
        x=(i64)x*x%P;
        a>>=;
    }
    return s;
}

int main()
{
    int T;
    scanf("%d",&T);
    while(T--){
        int r,c,l,n,m;
        scanf("%d%d%d%d",&r,&c,&l,&m);
        n=r*c*l;
        if(m>r&&m>c&&m>l){
            printf("0\n");
            continue;
        }
        memset(f,,sizeof(f));
        int w=,p=;
        f[w][][][][]=;
        int i,j,x,y,z;
        for(i=n;i;i--){
            swap(w,p);
            memset(f[w],,sizeof(f[w]));
            for(j=;j<=n-i+&&j<=m;j++)
                for(x=;x<=n-i+&&x<=r;x++)
                    for(y=;y<=n-i+&&y<=c;y++)
                        for(z=;z<=n-i+&&z<=l;z++){
                            upd(f[w][j][x][y][z],(i64)f[p][j][x][y][z]*max(,x*y*z-(n-i))%P);
                            upd(f[w][j][x][y][z],(i64)f[p][j][x-][y][z]*(r-x+)*y*z%P);
                            upd(f[w][j][x][y][z],(i64)f[p][j][x][y-][z]*x*(c-y+)*z%P);
                            upd(f[w][j][x][y][z],(i64)f[p][j][x][y][z-]*x*y*(l-z+)%P);
                            upd(f[w][j][x][y][z],(i64)f[p][j][x][y-][z-]*x*(c-y+)*(l-z+)%P);
                            upd(f[w][j][x][y][z],(i64)f[p][j][x-][y][z-]*(r-x+)*y*(l-z+)%P);
                            upd(f[w][j][x][y][z],(i64)f[p][j][x-][y-][z]*(r-x+)*(c-y+)*z%P);
                            upd(f[w][j][x][y][z],(i64)f[p][j-][x-][y-][z-]*(r-x+)*(c-y+)*(l-z+)%P);
                        }
        }
        p=;
        for(i=;i<=n;i++)
            p=(i64)p*i%P;
        printf("%d\n",(i64)f[w][m][r][c][l]*pwr(p,P-)%P);
    }
    return ;
}

AC：

容斥，钦定几个点为极大值点，其他不管，然后用组合数作为系数同时算本质相同的方案。

#include<iostream>
#include<cstdio>
using namespace std;
typedef long long i64;
const int N=5e6+;
const int P=;

int pwr(int x,int a){
    int s=;
    for(;a;a>>=,x=1ll*x*x%P)
        if(a&)s=1ll*s*x%P;
    return s;
}

int inv[N],iinv[N];
int sinv[N],siinv[N];
int v[N],f[N];

int C(int n,int m){
    return 1ll*inv[n]*iinv[m]%P*iinv[n-m]%P;
}

int main()
{
    int i;
    inv[]=;
    for(i=;i<N;i++)
        inv[i]=1ll*inv[i-]*i%P;
    iinv[N-]=pwr(inv[N-],P-);
    for(i=N-;i;i--)
        iinv[i-]=1ll*iinv[i]*i%P;
    int T;
    scanf("%d",&T);
    while(T--){
        int c,r,h,m;
        scanf("%d%d%d%d",&c,&r,&h,&m);
        int n=min(min(c,r),h);
        int sv=1ll*c*r%P*h%P;
        for(i=;i<=n;i++)
            v[i]=(P+sv-1ll*(c-i)*(r-i)%P*(h-i)%P)%P;
        for(i=;i<=n;i++)
            f[i]=1ll*C(r,i)*C(c,i)%P*C(h,i)%P*pwr(inv[i],)%P;
        sinv[]=;
        for(i=;i<=n;i++)
            sinv[i]=1ll*sinv[i-]*v[i]%P;
        siinv[n]=pwr(sinv[n],P-);
        for(i=n;i;i--)
            siinv[i-]=1ll*siinv[i]*v[i]%P;
        int ans=,d;
        for(i=m,d=;i<=n;i++,d=P-d)
            (ans+=1ll*d*C(i,m)%P*f[i]%P*siinv[i]%P)%=P;
        printf("%d\n",ans);
    }
    return ;
}

珍珠：

O(nd)：

只跟奇偶性有关。

#include<iostream>
#include<cstdio>
using namespace std;
typedef long long i64;
const int P=;
const int S=5e7+;
const int N1=4e3+;

int f[N1][N1];

int main()
{
    int d,n,m,i,j;
    scanf("%d%d%d",&d,&n,&m);
    if((i64)d*n<S){
        f[][]=;
        for(i=;i<=n;i++)
            for(j=;j<=d;j++){
                if(j)f[i][j]+=(i64)f[i-][j-]*(d-j+)%P,
                f[i][j]%=P;
                if(j!=d)f[i][j]+=(i64)f[i-][j+]*(j+)%P,
                f[i][j]%=P;
            }
        int ans=;
        for(i=;i<=n-*m&&i<=d;i++)
            ans+=f[n][i],ans%=P;
        printf("%d",ans);
        return ;
    }
    return ;
}