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

设f[i][j]表示当前状态为i余数为j的方案数，其实这个还是蛮好想得，就是一开始脑残，非要加上一维表示现在填到了第k位数字，结果时间复杂度死活过不了，真是智障qaq。

转移f[i|(1<<j)][(k*10+s[j]-'0')%d]+=f[i][k]

考虑去重，其实就是用总的排列数除以每一种重复数字的个数的排列就好了。

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
int f[1024][1010],x[10],cnt[10],d;
int T;char s[10];

int main(){
	scanf("%d",&T);
	while(T--){
		memset(f,0,sizeof(f));f[0][0]=1;
		memset(cnt,0,sizeof(cnt));
		scanf("%s%d",s,&d);int len=strlen(s);
		for(int i=0;i<len;i++)cnt[s[i]-'0']++;
		for(int i=0;i<1<<len;i++){
			for(int k=0;k<d;k++){
				for(int j=0;j<len;j++)if(!(1<<j&i))
					f[i|(1<<j)][(k*10+s[j]-'0')%d]+=f[i][k];
			}
		}
		int ans=f[(1<<len)-1][0];
		for(int i=0;i<10;i++)if(cnt[i])
			for(int j=1;j<=cnt[i];j++)ans/=j;
		cout<<ans<<endl;
	}
	return 0;
}