1072: [SCOI2007]排列perm 状压DP

我太弱了。。据说是状压裸题。。我还是不会做。。。T T

考虑用 fi,j 表示当前状态为 i ,且 mod d=j 的方案数。
若存在 k 使得 (1<<(k1)) and i==0 就可以转移。
最后全排列去重

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int d,len;
int a[15],v[15],tot[15],f[1025][1025];
char s[15];
inline int read()
{
    int a=0,f=1; char c=getchar();
    while (c<'0'||c>'9') {if (c=='-') f=-1; c=getchar();}
    while (c>='0'&&c<='9') {a=a*10+c-'0'; c=getchar();}
    return a*f;
}
inline void dp()
{
    memset(f,0,sizeof(f));
    f[0][0]=1;
    for (int i=0;i<(1<<len);i++)
        for (int j=0;j<d;j++)
            if (f[i][j])
                for (int k=1;k<=len;k++)
                    if (((1<<(k-1))&i)==0)
                        f[i|(1<<(k-1))][(a[k]+j*10)%d]+=f[i][j];
}
int main()
{
    int testcase=read();
    while (testcase--)
    {
        scanf("%s",s+1); d=read();
        len=strlen(s+1);
        for (int i=0;i<=9;i++) v[i]=1,tot[i]=0;
        for (int i=1;i<=len;i++)
        {
            a[i]=s[i]-'0';
            tot[a[i]]++;
            v[a[i]]*=tot[a[i]];
        }
        dp();
        for (int i=0;i<=9;i++) f[(1<<len)-1][0]/=v[i];
        printf("%d\n",f[(1<<len)-1][0]);
    }
    return 0;
}

你可能感兴趣的:(1072: [SCOI2007]排列perm 状压DP)