bzoj 1072: [SCOI2007]排列perm

2016-06-21

数字 难道是数位DP 可惜不是， 哎n<=10 我竟然没想到是状压DP，

f[i][j]表示i这个状态%d==j的方案数有多少，转移的时候枚举哪一个做最后一位。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<cstdlib>
#include<map>
#define ll long long
#define M 2000009
using namespace std;
ll read()
{
    char ch=getchar();
    ll x=0,f=1;
    for(;ch<'0'||ch>'9';ch=getchar())
      if(ch=='-')
        f=-1;
    for(;ch>='0'&&ch<='9';ch=getchar())
      x=x*10+ch-'0';
    return x*f;
}
ll a[1024][1009],c[15],D[15];
int b[15],T;
int main()
{
    c[0]=1;
    for(int i=1;i<=10;i++)
      c[i]=c[i-1]*i;
    T=read();
    for(;T;T--)
      {
          memset(D,0,sizeof(D));
          char ch[15];
          scanf("%s",ch);
          int n=strlen(ch);
        for(int i=0;i<n;i++)
            {
              b[i]=ch[i]-'0';
            D[b[i]]++;
          } 
        int d=read();
        int S=(1<<n);
        memset(a,0,sizeof(a));
        a[0][0]=1;
        for(int i=1;i<S;i++)
          for(int j=0;j<n;j++)
              if(i&(1<<j))
                for(int k=0;k<d;k++)
                   a[i][(k*10+b[j])%d]=a[i-(1<<j)][k]+a[i][(k*10+b[j])%d];
        ll an=a[S-1][0];
        for(int i=0;i<=9;i++)
          an/=c[D[i]];
        printf("%lld
",an);
      }
    return 0;
}