bzoj 1072: [SCOI2007]排列perm 状压dp

code:

#include <bits/stdc++.h>
#define N 1005
using namespace std;
void setIO(string s)
{
string in=s+".in";
string out=s+".out";
freopen(in.c_str(),"r",stdin);
// freopen(out.c_str(),"W",stdout);
}
char str[N];
int num[N],f[N<<1][N],poww[N],tax[11],fac[N],vis[11];
void solve()
{
memset(f,0,sizeof(f));
memset(tax,0,sizeof(tax));
memset(vis,0,sizeof(vis));
int mod,i,j,n,k;
scanf("%s%d",str,&mod);
n=strlen(str);
fac[0]=poww[0]=1;
for(i=1;i<=10;++i) fac[i]=fac[i-1]*i;
for(i=1;i<=n;++i) poww[i]=poww[i-1]*10%mod;
for(i=0;i<n;++i) num[i]=str[i]-'0', ++tax[num[i]];
f[0][0]=1;
for(i=0;i<(1<<n);++i)
{
for(j=0;j<mod;++j)
{
int sz=0;
for(k=0;k<n;++k) if(i&(1<<k)) ++sz;
for(k=0;k<n;++k)
{
if(i&(1<<k)) continue;
else
{
int nx=i|(1<<k);
int cur=num[k];
f[nx][(j+cur*poww[sz]%mod)%mod]+=f[i][j];
}
}
}
}
int ans=f[(1<<n)-1][0];
for(i=0;i<n;++i) if(!vis[num[i]]) vis[num[i]]=1, ans/=fac[tax[num[i]]];
printf("%d\n",ans);
}
int main()
{
// setIO("input");
int T;
scanf("%d",&T);
while(T--)
{
solve();
}
return 0;
}

  

上一篇:[置顶] Android访问控制系统测试与评估


下一篇:Mysql Cluster配置基本篇