记录编号	392577	评测结果	AAAAAAAAAA
题目名称	[HNOI 2008]Cards	最终得分	100
用户昵称	sxysxy	是否通过	通过
代码语言	C++	运行时间	0.125 s
提交时间	2017-04-08 10:05:00	内存使用	0.47 MiB

#include <cstdio>
#include <cstring>
#include <cstdarg>
#include <string>
#include <algorithm>
#include <vector>
#include <list>
#include <queue>
#include <deque>
#include <cctype>
using namespace std;
const int MAXN = 110;
int pow_mod(int a, int x, int p){
  int ans = 1;
  for(a %= p; x; a = a*a%p, x >>= 1)if(x&1)ans = ans*a%p;
  return ans;
}
inline int inv(int x, int p){
  return pow_mod(x, p-2, p); //求逆元
}
int f[MAXN/3][MAXN/3][MAXN/3], g[MAXN];
int st[MAXN], size;
bool b[MAXN];
int n, MOD;
int sr, sg, sb;
int calc(){
  int size = 0; memset(b, false, sizeof(b)); memset(st, 0, sizeof(st));
  for(int i = 1; i <= n; i++){
    if(!b[i]){
      size++;
      for(int j = i; !b[j]; j = g[j])b[j] = true, st[size]++;
          //把在同一个循环节里面的物品看作同一个。
    }
  }
  memset(f, 0, sizeof(f));
  f[0][0][0] = 1;
  for(int i = 1; i <= size; i++) //背包方案dp
    for(int r = sr; ~r; r--)for(int g = sg; ~g; g--)for(int b = sb; ~b; b--){
      if(r >= st[i])f[r][g][b] = (f[r][g][b]+f[r-st[i]][g][b])%MOD;
      if(g >= st[i])f[r][g][b] = (f[r][g][b]+f[r][g-st[i]][b])%MOD;
      if(b >= st[i])f[r][g][b] = (f[r][g][b]+f[r][g][b-st[i]])%MOD;
    }
  return f[sr][sg][sb];
}
int main(){
  freopen("hnoi_cards.in", "r", stdin); 
  freopen("hnoi_cards.out", "w", stdout);
  int m;
  scanf("%d %d %d %d %d", &sr, &sb, &sg, &m, &MOD);
  n = sr+sb+sg; //这个不要忘，迷之wa了半小时...
  int ans = 0;
  for(int i = 1; i <= m; i++){
    for(int i = 1; i <= n; i++)scanf("%d", g+i); //一个置换
    ans = (ans+calc())%MOD; //在这个置换下的方案。满足加法原理。
  }
  for(int i = 1; i <= n; i++)g[i] = i; 
  //还有恒等变换(如果之前给出的有恒等变换，这里再做一次也是不影响结果的)
  ans = (ans+calc())%MOD;
  ans = ans*inv(m+1, MOD)%MOD;
  printf("%d\n", ans);
  return 0;
}