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的倍数。
【限制】
100%的数据满足：s的长度不超过10, 1<=d<=1000, 1<=T<=15
分析：
觉得还是很像数位dp
但是实际上是一道状压dp
基本上不会状压dp，所以这道就算是状压first吧
把每个数选与不选表示成01串
最终状态就是2^strlen(s)-1
一看这个d<=1000
嗯很好，我们可以枚举这个余数了
设计状态：f[i][j]表示i状态，余数是j的个数
i状态就是已经选了哪些数
然而
直接递推会出现重复计算，
利如数字串为001，可以组合出001，010，100三种不重复的数字，
但是直接计算的话3！=6，
因为有两个0，所以他们对答案的贡献是2！，
那么我们只要在计算完之后，除去每个数字个数的阶乘就可以了。
代码详解：
num表示每个数出现的状态
jc表示每个数出现次数的阶乘
a表示原序列中这一位上具体是什么数
预处理完成之后，就是程序的主体部分：dp
我们已经明确了最终状态是
tt=(1<<len)-1；一个状态（例如1000010）
只能转移到比ta含有的1的个数多一个的状态（自己领悟一下）
所以我们只要循环完所有的状态，就可以完成pd了
第二维循环余数
第三维循环是枚举这一次我们要选哪一个数，
如何判断这一位选没选过：
没选过
((1<<(k-1))&i)==0
选过
((1<<(k-1))&i)==1
tip
d的前导零要处理一下
这里写代码片
#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
int f[1200][1005];
char s[11];
int zt[2000],len,sum[10],a[15],d,jc[15];
int main()
{
int T;
scanf("%d",&T);
while (T--)
{
scanf("%s",&s);
len=strlen(s);
memset(f,0,sizeof(f));
for (int i=0;i<=9;i++)
jc[i]=1,sum[i]=0,a[i]=0;
for (int i=0;i<len;i++)
{
a[i+1]=s[i]-'0';
sum[a[i+1]]++;
jc[a[i+1]]*=sum[a[i+1]];
}
scanf("%s",&s);
int i=0;
while (s[i]=='0') i++;
d=0;
for (int j=i;j<strlen(s);j++) d=d*10+s[j]-'0';
f[0][0]=1;
int tt=(1<<len)-1;
for (int i=0;i<=tt;i++)
for (int j=0;j<d;j++)
if (f[i][j])
{
for (int k=1;k<=len;k++)
if (((1<<(k-1))&i)==0)
f[i|(1<<(k-1))][(j*10+a[k])%d]+=f[i][j];
}
for (int i=0;i<=9;i++) f[tt][0]/=jc[i];
printf("%d
",f[tt][0]);
}
return 0;
}