https://codeforces.com/problemset/problem/401/D
思路:看到这个18开始往数位dp搞了。但是是个状压。
定义dp[s][j]:当前取数状态为s,且modm=j的方案数。
对于排列,我们考虑其枚举s中的一个数,然后从之前状态转移,且枚举到的这个数放在最后一个位置加进去。这样第二部分的mod就可以顺推,前一个状态为k,现在加入某个数的状态可以算出。
之前的数为X,X=km+j;选择了w[i]之后,X=10km+10j+w[i]; 余数就变成了:(10j+w[i])%m。
但是这样会重复转移。
举例:n=221 111会从101 转移一次,还会从011转移一次。然而这两个组合其实都是2、1,所以会算重。
因此我们对这次状态转移,对于取过的数放最后一位进行标记不再转移即可。(这部分也可以最后用组合数学的同色球除掉)
最后注意前导零的状态去除。
#include<iostream>
#include<vector>
#include<queue>
#include<cstring>
#include<cmath>
#include<map>
#include<set>
#include<cstdio>
#include<algorithm>
#define debug(a) cout<<#a<<"="<<a<<endl;
using namespace std;
const int maxn=19;
typedef long long LL;
inline LL read(){LL x=0,f=1;char ch=getchar(); while (!isdigit(ch)){if (ch=='-') f=-1;ch=getchar();}while (isdigit(ch)){x=x*10+ch-48;ch=getchar();}
return x*f;}
LL dp[1<<maxn][102];
int num[30],tot=0;
bool vis[maxn];
int main(void){
LL n,m;n=read();m=read();
vector<LL>v;
LL n1=n;
while(n1){
v.push_back(n1%10);
n1/=10;
}
reverse(v.begin(),v.end());
for(auto i:v) num[++tot]=i;
dp[0][0]=1;
n=tot;
for(LL s=0;s<(1<<n);s++){
memset(vis,0,sizeof(vis));
for(LL j=1;j<=n;j++){
if(s==(1LL<<(j-1))&&!num[j]) break;///去掉前导零
if(!(s&(1LL<<(j-1)))||vis[num[j]]) continue;
vis[num[j]]=1;///重复转移,也可以组合方法除掉
for(LL k=0;k<m;k++){
dp[s][(k*10+num[j])%m]+=dp[s^(1LL<<(j-1))][k];
}
}
}
printf("%lld\n",dp[(1<<n)-1][0]);
return 0;
}