【以前的空间】bzoj 1072 [SCOI2007]排列perm

又颓废了一个下午，最近撸mc撸到丧失意识了，玩的有点恶心，于是找水题做，瞧不起颓废的自己啊。

another水题。

这题题意很明显啦，就是找数字排列后组成的数去mod d=0后有多少种。

普通的搜索的话，是会tle的（应该是o（n！）没错？）。注意到长度n还是比较小的，于是想到状压dp。

状态就是每个数取和不取组成的结果（就是00110表示第3,4个数取了啦，学过状压都知道）。

然后转移就是f[i,j,k]表示现在取到第i个数状态为i余数为j有多少种情况，

那么f[i,j,(k*10+a[i])mod d]=Σf[i-1，j',k]（也就是同余的东西啦，在123后加一个数字4，那么1234 mod d=（（123 mod d *10）+4 ）mod d ）。k的范围就是0-（d-1）啦。其实就是一个01背包，然后i是可以去掉的。

最后就是处理重复的问题了，之前好像也做过类似的，不过反过来的……重复的话可以这么想，比如有4个2重复，那么对于第4个2来说，加上它后重复的情况就是之前的情况*4，对于第三个2来说，加上它后重复的情况就是之前的情况*3，对于第二个2就是之前的情况*2（其实就是个排列组合……四个2有序号要比没有序号多A（4,4)=4！）。

2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
var
b:array[0..10,0..2000]of longint;
c,a:array[0..20]of longint;
f:array[0..2000,0..1500]of longint;
i,j,k,l,n,m,t,tot,top:longint;

procedure into;
var
i,j,k,len:longint;
s:string;
begin
readln(s);
len:=length(s);
i:=pos(' ',s);
val(copy(s,i+1,len-i),tot);
delete(s,i,len-i+1);
n:=length(s);
fillchar(c,sizeof(c),0);
for i:=1 to n do begin
a[i]:=ord(s[i])-ord('0');
inc(c[a[i]]);
end;
top:=1<<n-1;
fillchar(b,sizeof(b),0);
for i:=0 to top do begin
j:=i;
k:=0;
while j<>0 do begin
inc(k);
dec(j,j and (-j));
end;
j:=k;
inc(b[j,0]);
b[j,b[j,0]]:=i;
end;
end;

procedure work;
var
i,j,k,l,m,ans:longint;
begin
fillchar(f,sizeof(f),0);
for i:=1 to n do
f[1<<(i-1),a[i] mod tot]:=1;
for i:=2 to n do
for j:=1 to b[i,0] do begin
k:=b[i,j];
for l:=1 to n do
if k and (1<<(l-1)) <> 0 then
for m:=0 to tot-1 do
inc(f[k,(10*m+a[l]) mod tot],f[k-1<<(l-1),m]);
end;
ans:=f[top,0];
for i:=0 to 9 do
if c[i]>1 then
for j:=2 to c[i] do ans:=ans div j;
writeln(ans);
end;

begin
readln(t);
while t>0 do begin
dec(t);
into;
work;
end;
end.

View Code