BZOJ 1072 SCOI2007 排列perm 狀壓DP

標籤：bzoj   bzoj1072   狀壓dp   

題目大意：給定n個數字，求這些數位全排列中有多少數能被d整除

令f[i][j]為狀態為i，餘數為j的方案數

枚舉最高位轉移

小心爆int

#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>using namespace std;int n,d,ans,f[1<<10][1<<10],digit[1<<10],tens[10],cnt[10],factorial[11];char s[20];int Get_Digit(int x){int re=0;while(x)re+=x&1,x>>=1;return re;}void State_Compressed_DP(int x){int i,j;if(x==1023)++x,--x;for(i=0;i<n;i++)if(1<<i&x){int temp=tens[digit[x]-1]%d*(s[i]-'0')%d;for(j=0;j<d;j++)f[x][(j+temp)%d]+=f[1<<i^x][j];}}int main(){int T,i;for(i=0;i<1<<10;i++)digit[i]=Get_Digit(i);tens[0]=1;for(i=1;i<10;i++)tens[i]=tens[i-1]*10;factorial[0]=1;for(i=1;i<=10;i++)factorial[i]=factorial[i-1]*i;for(cin>>T;T;T--){memset(f,0,sizeof f);scanf("%s%d",s,&d);n=strlen(s);f[0][0]=1;for(i=1;i<1<<n;i++)State_Compressed_DP(i);ans=f[(1<<n)-1][0];memset(cnt,0,sizeof cnt);for(i=0;i<n;i++)cnt[s[i]-'0']++;for(i=0;i<=9;i++)ans/=factorial[cnt[i]];cout<<ans<<endl;}}

BZOJ 1072 SCOI2007 排列perm 狀壓DP