BZOJ1072

1072: [SCOI2007]排列perm

Time Limit: 10 Sec   Memory Limit: 162 MB
Submit: 758   Solved: 471
[ Submit][ Status]

Description

给一个数字串s和正整数d, 统计s有多少种不同的排列能被d整除（可以有前导0）。例如123434有90种排列能被2整除，其中末位为2的有30种，末位为4的有60种。

Input

输入第一行是一个整数T，表示测试数据的个数，以下每行一组s和d，中间用空格隔开。s保证只包含数字0, 1, 2, 3, 4, 5, 6, 7, 8, 9.

Output

每个数据仅一行，表示能被d整除的排列的个数。

Sample Input

7
000 1
001 1
1234567890 1
123434 2
1234 7
12345 17
12345678 29

Sample Output

1
3
3628800
90
3
6
1398

HINT

在前三个例子中，排列分别有1, 3, 3628800种，它们都是1的倍数。 【限制】 20%的数据满足：s的长度不超过5, 1<=T<=5 50%的数据满足：s的长度不超过8 100%的数据满足：s的长度不超过10, 1<=d<=1000, 1<=T<=15,

【题解】可以用状态压缩DP,f[ i ][ j ][ k ],i 我们用二进制来理解，因为s 至多只有10为，所以在二进制中为1，表示选了这个数，0表示没有选。 j 表示当前状态下，最后一位是那个数。k 表示当前的余数。

   转移方程为，f[ i ][ j ][(p*10+a[j])%d]=f[ i^(1<<(j-1))][k][p];是从减去第 j 个数的状态转移过来的。详见代码。对了其中有个需要注意的地方，就是两个重复数字的情况。我用了两个数组判掉了。

#include<cstdio>
#include<cstring>

int f[1024][11][1000],a[11],w[10],twice[10];
 
int main()
{
	int t,s,d;scanf("%d",&t);
	while(t--){
		char ch;int l=0;
		while(ch=getchar(),ch=='\n' || ch==' ' || ch=='\r');
		while(ch!=' ')a[++l]=(int)ch-48,ch=getchar();
		scanf("%d",&d);
		memset(f,0,sizeof(f));
		for(int i=1;i<=(1<<l)-1;i++)
		{
		  memset(w,0,sizeof(w));
		  for(int j=1;j<=l && (1<<(j-1))<=i;j++)
		   if((1<<(j-1))&i && (!w[a[j]])){
			 w[a[j]]=1;
			 if(!(i^(1<<(j-1))))f[i][j][a[j]%d]=1;
			  else{
					memset(twice,0,sizeof(twice));
					for(int k=1;k<=l && (1<<(k-1))<=i;k++)
			  		if((j-k)&&((1<<(k-1))&i)&&(!twice[a[k]]))
			  		{
						twice[a[k]]=1;
						for(int p=0;p<=d-1;p++)
					    f[i][j][(p*10+a[j])%d]+=f[i^(1<<(j-1))][k][p];
					}
			     }
		    }
		}
		int ans=0,t=(1<<l)-1;
		for(int i=1;i<=l;i++)ans+=f[t][i][0];
		printf("%d\n",ans);
	}
	return 0;
}