第一道 状态压缩dp

f[mask | 1 << i][(k*10+a[i]) % m] += f[mask][k];

using namespace std;
#define long long long
long f[1 << 18][100];
int a[20];
int main(){
long n, m;
int N = 0;
cin >> n >> m;
for( ; n>=0; n /=10) a[N++] = n%10;
sort(a,a+N);
int mask, k;
f[0][0] = 1;
for( mask = 0; mask < (1 << N); mask ++)
for( k = 0; k < m; k ++) if( f[mask][k] ){                //f[mask][k] 即能否产生这样的情况
long cnt = f[mask][k], pre = -1;                      //前面已经对a[] 进行排序，这里加入pre后可以去除重复
for(int i = N-1; i >= 0; i --) if(~mask >> i & 1 ){   //mask 中的第i位 即第i个数是否用过
if(a[i] == pre) continue;                         //去重复
if(a[i] == 0 && mask == 0) continue;              //不是很理解
f[ mask | 1 << i][(k*10 + a[i]) % m] += cnt;
pre = a[i];
}
}
cout << f[(1 << N)-1][0] << endl;
return 0;
}