Codeforces 401D Roman and Numbers 状压DP

题意:给出数n,m.问有多少个x 满足x%m==0 x和n包含相同的数字(无前导0,顺序可以不同)?
n<=1e18,m<=100.4s时限

最多18位数字,用18位二进制表示当前状态,第i位为1则表示n的第i位数字被使用,答案为dp[(1<<18)-1][0]]
按最后一个添加进集合的数字转移状态,注意(a[1],a[3]) (a[3],a[1])被重复计算 答案在除以每个数字出现次数的排列.

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1<<19;
ll n,m,dp[N][101],cnt[20],f[20];
ll a[20],pos;
void init()
{
	pos=0;
	memset(cnt,0,sizeof(cnt));
	while(n)
	{
		a[pos++]=n%10;
		cnt[n%10]++;
		n/=10;
	}
}
int main()
{
	f[0]=1;
	for(int i=1;i<20;i++)
		f[i]=f[i-1]*(ll)i;
	while(cin>>n>>m)
	{
		memset(dp,0,sizeof(dp));
		dp[0][0]=1;
		init();
		int s=1<<pos;
		for(int i=0;i<s;i++)
		{
			for(int j=0;j<m;j++)
			{
				for(int k=0;k<pos;k++)
				{
					if((i>>k)&1)
						continue;
					if(i==0&&a[k]==0)
						continue;
					
					dp[i|1<<k][(j*10+a[k])%m]+=dp[i][j];
				}
			}
		}
		ll ans=dp[s-1][0];
		for(int i=0;i<=9;i++)
			ans/=f[cnt[i]];
		cout<<ans<<endl;
	}
	return 0;
}