【BZOJ1004】Cards

polya 계수 + DP의 입문 문제, 인터넷에서 거의 모든 polya 계수를 소개하는 웹 페이지를 직시할 수 없습니다...그룹 수학 교재나 국가 훈련팀의 강연 원고를 사도록 추천합니다...
제목 입력에서 "데이터를 입력하면 임의로 여러 번 카드를 씻을 수 있다. 모두 이 m종의 카드를 씻는 방법 중의 하나로 대체할 수 있고 모든 카드를 씻는 방법에 대해 원래의 상태로 돌아갈 수 있는 방법이 존재한다."모든 방안에 변하지 않는 교환을 더하면 하나의 교환군이다. 이런 중요한 조건들은 왜 이렇게 함축적인가...
사실 카드를 씻는 방식은 저장할 필요가 없어요. 하나의DP에 기록해도 돼요...처음 써봤는데 그렇게 생각 안 했어...이분 쾌속 멱을 쓰는 것이 귀찮아졌고, 폭력을 다 쓴 후에야 길이의 차이가 많지 않다는 것을 발견하였다.

#include<cstdlib>
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<queue>
using namespace std;
int SA,SB,SC,m,p,n;
struct data
{
	int a[65];
	data()
	{
		memset(a,0,sizeof(a));
	}
}A[65];

void Init()
{
	scanf("%d%d%d%d%d",&SA,&SB,&SC,&m,&p);
	n=SA+SB+SC;
	for(int i=1;i<=m;i++)
	{
		for(int j=1;j<=n;j++)
		{
			scanf("%d",&A[i].a[j]);
		}
	}
	m++;
	for(int j=1;j<=n;j++)
	{
		A[m].a[j]=j;
	}
	return ;
}
int dp[25][25][25];
int qkpower(int a,int x,int mo)
{
	int ans=1;
	for(int i=1;i<=x;i++)
	{
		ans=(ans*a)%mo;
	}
	return ans;
}
void work()
{
	bool v[65];
	int cc[65],cnt,sum=0;
	for(int i=1;i<=m;i++)
	{
		for(int j=1;j<=n;j++)
		{
			v[j]=0; cc[j]=0;
		}
		cnt=0;
		int t;
		for(int j=1;j<=n;j++)if(!v[j])
		{
			cnt++; 
			t=A[i].a[j];
			while(t!=j)
			{
				cc[cnt]++; v[t]=1; t=A[i].a[t];
			}
			cc[cnt]++; v[t]=1; 
		}
		for(int u=0;u<=SA;u++)
		{
			for(int v=0;v<=SB;v++)
			{
				for(int w=0;w<=SC;w++)
				{
					dp[u][v][w]=0;
				}
			}
		}
		dp[0][0][0]=1;
		for(int k=1;k<=cnt;k++)
		{
			for(int v=SB;v>=0;v--)
			{
				for(int w=SC;w>=0;w--)
				{
					for(int u=SA;u>=0;u--)
					{
						if(cc[k]+u<=SA)
						{
							dp[u+cc[k]][v][w]=(dp[u+cc[k]][v][w]+dp[u][v][w])%p;
						}
						if(cc[k]+v<=SB)
						{
							dp[u][v+cc[k]][w]=(dp[u][v+cc[k]][w]+dp[u][v][w])%p;
						}
						if(cc[k]+w<=SC)
						{
							dp[u][v][w+cc[k]]=(dp[u][v][w+cc[k]]+dp[u][v][w])%p;
						}
					}
				}
			}
		}
		sum=(sum+dp[SA][SB][SC])%p;
	}
	sum=(sum*qkpower(m,p-2,p))%p;
	printf("%d
",sum);
	return ;
}
int main()
{
	freopen("in.txt","r",stdin);
	Init();
	work();
	return 0;
}