LightOJ1158 Anagram Division（状压DP）

题目问一个数字字符串的不重复全排列有几个能被d整除。

• dp[S][m]表示用字符集合S构成的%d为m的数字字符串个数
• dp[0][0]=0
• 我为人人转移，dp[S+{x}][(m*10+str[x]-'0')%d]+=dp[S][m]（x∉S）

最后的结果再除以各字符出现次数的阶乘就是答案了，即dp[2^strlen-1][0]/(t[0]!*t[1]!*t[2]!*t[3]!*t[4]!*t[5]!*t[6]!*t[7]!*t[8]!*t[9]!)。

#include<cstdio>
#include<cstring>
using namespace std;
int d[1<<10][1111],fact[11]={1};
int main(){
    for(int i=1; i<11; ++i) fact[i]=fact[i-1]*i;
    int t,m;
    char str[1111];
    scanf("%d",&t);
    for(int cse=1; cse<=t; ++cse){
        scanf("%s%d",str,&m);
        int n=strlen(str);
        memset(d,0,sizeof(d));
        d[0][0]=1;
        for(int i=0; i<(1<<n)-1; ++i){
            for(int j=0; j<m; ++j){
                if(d[i][j]==0) continue;
                for(int k=0; k<n; ++k){
                    if((i>>k)&1) continue;
                    d[i|(1<<k)][(j*10+str[k]-'0')%m]+=d[i][j];
                }
            }
        }
        int times[10]={0};
        for(int i=0; i<n; ++i) ++times[str[i]-'0'];
        for(int i=0; i<10; ++i){
            d[(1<<n)-1][0]/=fact[times[i]];
        }
        printf("Case %d: %d\n",cse,d[(1<<n)-1][0]);
    }
    return 0;
}