1072: [SCOI2007]排列perm - BZOJ

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的倍数。 【限制】 20%的数据满足：s的长度不超过5, 1<=T<=5 50%的数据满足：s的长度不超过8 100%的数据满足：s的长度不超过10, 1<=d<=1000, 1<=T<=15,

状压dp

f[s,i]表示状态为s（用二进制表示选的数字的情况），余数为i的方案数

最后把重复的去掉，即除以每一种数字的个数的阶乘

 var
         flag,f:array[..,..]of longint;
         x,y:array[..]of longint;
         t,time,d,n:longint;
         s:array[..]of longint;

 procedure init;
 var
         c:char;
 begin
         read(c);
         n:=;
         while c<>' ' do
           begin
             inc(n);
             s[n]:=ord(c)-ord('');
             read(c);
           end;
         readln(d);
 end;

 procedure work;
 var
         head,tail,i,j,k,nx,ny:longint;
 begin
         flag[,]:=time;
         f[,]:=;
         head:=;
         tail:=;
         x[]:=;
         y[]:=;
         while head<=tail do
           begin
             for i:= to n do
               if x[head] and (<<(i-))= then
               begin
                 nx:=x[head]+<<(i-);
                 ny:=(y[head]*+s[i])mod d;
                 if flag[nx,ny]<>time then
                 begin
                   flag[nx,ny]:=time;
                   f[nx,ny]:=;
                   inc(tail);
                   x[tail]:=nx;
                   y[tail]:=ny;
                 end;
                 inc(f[nx,ny],f[x[head],y[head]]);
               end;
             inc(head);
           end;
         for i:= to n do
           begin
             k:=;
             for j:= to i do
               if s[j]=s[i] then inc(k);
             f[<<n-,]:=f[<<n-,] div k;
           end;
         if flag[<<n-,]=time then writeln(f[<<n-,])
         else writeln();
 end;

 begin
         readln(t);
         while t> do
           begin
             dec(t);
             inc(time);
             init;
             work;
           end;
 end.