方法:状压DP?
题解:这题其实没啥好写的,不算很难,推一推就能搞出来。
首先看到这个问题,对于被d整除这个条件,很容易就想到是取余数为0,所以想到可能状态中刚开始含有取余数。
先说我的第一个想法,f[i][j]表示选取i个数且此时的mod为j,这样的思想是第一下蹦出来的,当时想的就是在线来搞最终的答案。不过转瞬即发现,这TM不就是暴力吗魂淡!并没有什么卵用,于是开始想这个状态可不可以做什么优化。
显然第二维的j并不需要太大的优化,暂且先将其搁置一边,来考虑第一维的i怎么优化。
把滚动条拉到最下方发现这个i的范围最大是10,也就是说hash的空间为2^10-1,可以从状态的角度往下走,假设这个i代表hash状态,能不能表示他的转移方程呢?
首先按照一般的想法,对于每个i,可能从它的某一为1的位数变为零的状态转移过来,如果按照这个方向搞的话也许也能做出来,不过,我们更熟悉的操作是加上而并不是从原状态减去某个位,所以思路又转向到逆向更新。
也就是说我们可以初步写出来dp方程
f[(i|(1<<(k−1)))][(j∗10+a[k])%d]+=f[i][j]
其中更新的时候要保证i&(1<<(k-1))==0
然后我们做的只需要枚举i,枚举j,枚举k
嗯哼就是这样,然而这中途当然有很多自动剪掉的部分,所以复杂度不是很高。题目时间范围可承受。
但是并没有完,这时候你会发现如此的代码连第一个小样例都过不了,也就是说,如果有重复的数字的话,我们会算很多重复的情况,而怎么去重呢?只需要将答案除掉重复的数的个数的排列就好了,也就是阶乘!
后记:其实写这个题解的目的就是将自己解一道题的全过程写出来。
代码:
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 15
#define M 1100
using namespace std;
int tot,mod,l;
int a[N];
int cnt[N];
char s[N];
int f[1<<10][M];
int A[]={0,1,2,6,24,120,720,5040,40320,362880,3628800};
void dfs()
{
for(int i=0;i<(1<<l);i++)
{
for(int k=0;k<mod;k++)
{
if(f[i][k])
{
for(int j=1;j<=l;j++)
{
if((i&(1<<(j-1)))==0)
{
f[(i|(1<<(j-1)))][(k*10+a[j])%mod]+=f[i][k];
}
}
}
}
}
}
int main()
{
scanf("%d",&tot);
while(tot--)
{
scanf("%s",s);
l=strlen(s);
memset(cnt,0,sizeof(cnt));
for(int i=0;i<l;i++)
{
a[i+1]=s[i]-'0';
cnt[a[i+1]]++;
}
scanf("%d",&mod);
memset(f,0,sizeof(f));
f[0][0]=1;
dfs();
int ans=f[(1<<l)-1][0];
for(int i=0;i<=9;i++)
{
if(cnt[i])ans/=A[cnt[i]];
}
printf("%d\n",ans);
}
}