SCOI2007排列perm

1072: [SCOI2007]排列perm

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 805  Solved: 497
[Submit][Status]

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,

Source

题解:

又一道恶心的状压DP，我对状压还不熟悉。。。

代码：

var f:array[0..10,0..1025,0..1025] of longint;
    i,j,nn,k,m,t,l,ans,d,n:longint;
    s,st:ansistring;
    a,calc,c:array[0..1025] of longint;
    b:array[0..10,0..1025] of longint;
procedure init;
 begin
   readln(s);
   st:=copy(s,1,pos(' ',s)-1);
   n:=length(st);
   for i:=1 to n do a[i]:=ord(st[i])-ord('0');
   delete(s,1,pos(' ',s));
   val(s,d);//writeln(st,' ',d);
   fillchar(b,sizeof(b),0);
   for i:=0 to 1<<n-1 do
     begin
      j:=i;calc[i]:=0;
      while j<>0 do
       begin
        inc(calc[i]);dec(j,j and (-j));
       end;
      j:=calc[i];
      inc(b[j,0]);b[j,b[j,0]]:=i;
     end;
 end;
procedure main;
 begin
   fillchar(f,sizeof(f),0);
   for i:=1 to n do f[1,1<<(i-1),a[i] mod d]:=1;
   for i:=2 to n do
    for l:=1 to b[i,0] do
     begin
      j:=b[i,l];
      for k:=1 to n do
       if j and (1<<(k-1))<>0 then
        for m:=0 to d-1 do
         inc(f[i,j,(10*m+a[k]) mod d],f[i-1,j-1<<(k-1),m]);
     end;
   ans:=f[n,1<<n-1,0];
   fillchar(c,sizeof(c),0);
   for i:=1 to n do inc(c[a[i]]);
   for i:=0 to 9 do
     if c[i]<=1 then continue
     else for j:=2 to c[i] do ans:=ans div j;
   writeln(ans);
 end;

begin
  assign(input,'input.txt');assign(output,'output.txt');
  reset(input);rewrite(output);
  readln(t);
  while t>0 do
   begin
    dec(t);
    init;
    main;
   end;
  close(input);close(output);
end.