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