BZOJ_P1072 [SCOI2007]排列perm(状态压缩动态规划)

来源：互联网发布：linux stopped 编辑：程序博客网时间：2024/05/26 12:05

BZOJ传送门
Time Limit: 10 Sec Memory Limit: 162 MB
Submit: 1484 Solved: 931
[Submit][Status][Discuss]
Description
给一个数字串s和正整数d, 统计s有多少种不同的排列能被d整除（可以有前导0）。例如123434有90种排列能被2整除，其中末位为2的有30种，末位为4的有60种。

Input
输入第一行是一个整数T，表示测试数据的个数，以下每行一组s和d，中间用空格隔开。s保证只包含数字0, 1, 2, 3, 4, 5, 6, 7, 8, 9.

Output
每个数据仅一行，表示能被d整除的排列的个数。

Sample Input
7
000 1
001 1
1234567890 1
123434 2
1234 7
12345 17
12345678 29

Sample Output
1
3
3628800
90
3
6
1398

HINT
在前三个例子中，排列分别有1, 3, 3628800种，它们都是1的倍数。
【限制】
100%的数据满足：s的长度不超过10, 1<=d<=1000, 1<=T<=15

Source

因为s长度不超过10，将数字的使用变成状态，长为1024，f[i][j]中i表示状态，j表示余数；
每次枚举S中1的位数，然后状态转移，将第i位删去变成倒数第tot位，然后就简单咯
f[S][(j+s[i]*pow(10,tot-1)]+=f[S^(1<

#include<cstdio>#include<cstring>#include<algorithm>#include<iostream>using namespace std;#define N 1025int f[N][N];int d;inline int pow(int a,int b,int res=1){while(b){if(b&1) res*=a%d;a*=a%d;b>>=1;}return res;}int main(){    int t,s[15],cnt[10],v[10],l,lim;scanf("%d",&t);char ch=getchar();    while(t--){        for(int i=0;i<10;i++) cnt[i]=0,v[i]=1;l=0,d=0;        while(ch>'9'||ch<'0') ch=getchar();        while(ch>='0'&&ch<='9') s[++l]=ch-'0',cnt[s[l]]++,v[s[l]]*=cnt[s[l]],ch=getchar();        while(ch>'9'||ch<'0') ch=getchar();        while(ch>='0'&&ch<='9') d=d*10+ch-'0',ch=getchar();        memset(f,0,sizeof(f));lim=(1<<l);f[0][0]=1;        for(int S=1,tot=0;S<lim;S++,tot=0){            for(int i=0;i<l;i++) if(S&(1<<i)) tot++;tot=pow(10,tot-1);            for(int i=0;i<l;i++)                if(S&(1<<i)){                    for(int j=0,w=(1<<i);j<d;j++){                        f[S][(j+s[i+1]*tot)%d]+=f[S^w][j];                    }                }        }        for(int i=0;i<10;i++) f[lim-1][0]/=v[i];        printf("%d\n",f[lim-1][0]);    }    return 0;}

0 0