标签:状压DP
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模板题,f[i][j]表示状态为i时余数为j的个数
然后排列组合去重
Code
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define dep(i,a,b) for(int i=a;i>=b;i--)
#define LL long long
#define mem(x,num) memset(x,num,sizeof x)
using namespace std;
inline LL read()
{
LL f=1,x=0;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
const int maxn=20;
LL bin[20],T,d,len,a[maxn],v[maxn],tot[maxn],f[1<<11][1<<11];
char ch[maxn];
void dp()
{
rep(i,0,bin[len]-1)
rep(j,0,d-1)f[i][j]=0;
f[0][0]=1;
rep(i,0,bin[len]-1)
rep(j,0,d-1)
if(f[i][j])
rep(x,1,len)
if((bin[x-1]&i)==0)
f[i|bin[x-1]][(a[x]+j*10)%d]+=f[i][j];
}
int main()
{
T=read();
bin[0]=1;rep(i,1,19)bin[i]=bin[i-1]<<1;
while(T--){
scanf("%s",ch+1);d=read();
len=strlen(ch+1);
rep(i,0,9)v[i]=1,tot[i]=0;
rep(i,1,len){
a[i]=ch[i]-'0';
tot[a[i]]++;
v[a[i]]*=tot[a[i]];
}
dp();
rep(i,0,9)f[bin[len]-1][0]/=v[i];
printf("%lld\n",f[bin[len]-1][0]);
}
return 0;
}