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[s][j]表示状态s下【s表示已经选择了哪些数】余数为j的方案数，那么f[s | (1<<i-1)][(j * 10 + a[i])%d] += f[s][j]

很明显，状态s下可以通过在末尾添加一个不在状态中的i号数来转移到s|(1<<i-1)这个状态

设立一个状态（突然不知道怎么用语言描述这个状态）DP(I,J)，IDP(I,J)，I表示当前枚举的状态，JJ表示状态I对应的数字的数值，这样设立状态后很容易得出下面的状态转移方程：

其中判断语句的含义是在II状态中，KK号没有选择出来。

状态转移结束后，由于数串中可能存在重复的数字（样例已经给出来了），这个时候我们就会有许多重复的计算。这个问题很好解决，我们根据排列的知识将最后的Ans/=Cnt[I] Cnt[I] 记录数字I在数串中出现的次数）就可以了。

C++代码/暴力

#include<bits/stdc++.h>

using namespace std;

int a[],b[];

int n ,d;

int ans ;

void dfs(int i,long long x){

    if(i > n){

        if(x % d == ) ans ++;

        return;

    }

    for(int j = ;j <= ;j ++){

        if(b[j]) {

            --b[j];

        dfs(i + ,x *  + j);

        ++b[j];

    }

    }

}

int main(){

    int t;

    cin >> t;

    while(t--){

        string str;

        cin >> str;

        ans = ;

        n = str.size();

        memset(a,,sizeof a);

        memset(b,,sizeof b);

        for(int i = ; i < str.size() ; i++){

            a[i] = str[i] - '';

            b[a[i]]++;

        }

        cin >> d;

        dfs(,);

        cout << ans << endl;

    }

}

C++代码/状压

#include <cmath>

#include <cstdio>

#include <cstring>

#include <iostream>

#include <algorithm>

#define DB double

#define SG string

#define LL long long

#define DP(A,B) DP[A][B]

#define Fp(A,B,C,D) for(A=B;A<=C;A+=D)

#define Fm(A,B,C,D) for(A=B;A>=C;A-=D)

#define Clear(A) memset(A,0,sizeof(A))

using namespace std;

const LL Mod=1e9+;

const LL Max=2e3+;

const LL Inf=1e18;

LL T,M,Num[Max],Cnt[Max],DP[Max][Max];

inline LL Read(){

    LL X=;char CH=getchar();bool F=;

    while(CH>''||CH<''){if(CH=='-')F=;CH=getchar();}

    while(CH>=''&&CH<=''){X=(X<<)+(X<<)+CH-'';CH=getchar();}

    return F?-X:X;

}

inline void Write(LL X){

    if(X<)X=-X,putchar('-');

    if(X>)Write(X/);

    putchar(X%+);

}

int main(){

    LL I,J,K,L;

    T=Read();

    while(T--){

        Clear(Cnt);Clear(DP);

        char CH[];scanf("%s",CH+);

        LL Length=strlen(CH+);M=Read();

        Fp(I,,Length,){

            Num[I]=CH[I]-'';

            Cnt[Num[I]]++;

        }DP(,)=;

        K=(<<Length)-;

        Fp(I,,K,){

            Fp(J,,M-,){

                Fp(L,,Length,){

                    if((I&(<<(L-)))==){

                        DP(I|(<<(L-)),((J<<)+(J<<)+Num[L])%M)+=DP(I,J);

                    }

                }

            }

        }LL Ans=DP(K,);

        Fp(I,,,){

            Fp(J,,Cnt[I],){

                Ans/=J;

            }

        }Write(Ans);putchar('\n');

    }

    return ;

}

巴特西

排列perm HYSBZ - 1072（状压dp/暴力）