[BZOJ1004]Cards

　　polya定理＋乘法逆元

#include<bits/stdc++.h>
using namespace std;
#define maxn 65
#define maxs 25
int sr,sb,sg,n,m,p;
int a[maxn][maxn],f[maxs][maxs][maxs],vis[maxn],cal[maxn];

int dp(int x){
    memset(vis,0,sizeof(vis));
    memset(f,0,sizeof(f));
    int cc=0;
    for(int i=1;i<=n;i++){
        if(!vis[i]){
            vis[i]=1;
            cal[++cc]=1;
            int ii=i;
            while(!vis[a[x][ii]]){
                vis[a[x][ii]]=1;
                cal[cc]++;
                ii=a[x][ii];
            }
        }
    }
    f[0][0][0]=1;
    for(int t=1;t<=cc;t++)
        for(int i=sr;i>=0;i--)
            for(int j=sb;j>=0;j--)
                for(int k=sg;k>=0;k--){
                    if(i>=cal[t])f[i][j][k]=(f[i][j][k]+f[i-cal[t]][j][k])%p;
                    if(j>=cal[t])f[i][j][k]=(f[i][j][k]+f[i][j-cal[t]][k])%p;
                    if(k>=cal[t])f[i][j][k]=(f[i][j][k]+f[i][j][k-cal[t]])%p;
                }
    return f[sr][sb][sg];
}
int main(){
    scanf("%d%d%d%d%d",&sr,&sb,&sg,&m,&p);
    n=sr+sb+sg;
    for(int i=1;i<=m;i++)
        for(int j=1;j<=n;j++)
            scanf("%d",&a[i][j]);
    m++;
    for(int i=1;i<=n;i++)a[m][i]=i;
    int ans=0;
    for(int i=1;i<=m;i++)
        ans=(ans+dp(i))%p;
    for(int i=1;i<=p-2;i++)
        ans=(ans*m)%p;
    printf("%d
",ans);
    return 0;
}