bzoj 1072: [SCOI2007]排列perm【状压dp】

先写了个next_permutation结果T了，于是开始写状压

设f[s][i]为选取状态为s，选的数模d为i的方案数，去重的话直接除以每个数字的出现次数的阶乘即可

#include<iostream>

#include<cstdio>

#include<cstring>

using namespace std;

const int N=20;

int T,d,n,a[N],fac[N],c[N],f[2005][2005];

char s[N];

int main()

{

	scanf("%d",&T);

	while(T--)

	{

		scanf("%s%d",s+1,&d);

		n=strlen(s+1);

		for(int i=0;i<=9;i++)

			fac[i]=1,c[i]=0;

		for(int i=1;i<=n;i++)

		{

			a[i]=s[i]-'0';

			c[a[i]]++;

			fac[a[i]]*=c[a[i]];

		}

		for(int i=0;i<(1<<n);i++)

			for(int k=0;k<d;k++)

				f[i][k]=0;

		f[0][0]=1;

		for(int i=0;i<(1<<n);i++)

			for(int j=0;j<d;j++)

				for(int k=1;k<=n;k++)

					if(!((1<<(k-1))&i))

						f[i|(1<<(k-1))][(a[k]+j*10)%d]+=f[i][j];

		for(int i=0;i<=9;i++)

			f[(1<<n)-1][0]/=fac[i];

		printf("%d\n",f[(1<<n)-1][0]);

	}

	return 0;

}

巴特西

bzoj 1072: [SCOI2007]排列perm【状压dp】

最新文章

热门文章