[Bzoj1072][SCOI2007]排列perm(状压dp)

题目链接：https://www.lydsy.com/JudgeOnline/problem.php?id=1072

状压dp，dp[i][j]表示状态为i，对d取余为j的方案数，则有dp[i|(1<<k)][(j*10+a[k])%d]+=dp[i][j],其中k为不在状态i中的数。

同时，因为有重复的数字，所以最后还要除掉重复数字。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 3e5 + 10;
const ll mod = 1e9 + 7;
const ll INF = 1e18;
int dp[1050][1010];
int main() {
    int t;
    scanf("%d", &t);
    while (t--) {
        int d;
        char s[15];
        int vis[15];
        int a[15];
        scanf("%s%d", s, &d);
        int len = strlen(s);
        memset(dp, 0, sizeof(dp));
        memset(vis, 0, sizeof(vis));
        for (int i = 0; i < len; i++)
            a[i] = s[i] - '0', vis[a[i]]++;
        dp[0][0] = 1;
        for (int i = 0; i < (1 << len); i++) {
            for (int j = 0; j < d; j++) {
                for (int k = 0; k < len; k++) {
                    if ((i&(1 << k)) == 0)
                        dp[i | (1 << k)][(j * 10 + a[k]) % d] += dp[i][j];
                }
            }
        }
        int ans = dp[(1 << len) - 1][0];
        for (int i = 0; i <= 9; i++)
            for (int j = 1; j <= vis[i]; j++)
                ans /= j;
        printf("%d\n", ans);
    }
    return 0;
}