51nod1780 完美序列

　　首先我们需要排个序，如果排序之后相邻两项差超过1就是无解。统计一下每个数出现的次数num[i]。

　　然后我们很容易发现，一个完美序列，去掉所有权值大于某个值的数之后，还是完美的。

　　这样我们就考虑一路dp过去，将每一种数插进去。然后发现，当前数的决策之和上一个数有关。我们插的位置实际上只有两种可能：两边都是上一种数，一边是上一种数而且另一边没有数。这样状态比较好记录，只要记录上一种数有多少个连续的段和有有没有放在左右两边就好了。假如第i种数被分成了x段，考虑第一种情况，那就可以转移到第i+1种数被分成1..x-num[i]段，组合数算一下，应该比较好搞。有没有放在左右两边也可以讨论，非常简单。可以参考代码。

program j01;
const maxn=30086;mo=1000000007;
var a:array[0..maxn]of longint;
    c:array[0..200,0..200]of int64;
    num:array[0..maxn]of longint;
    cnt,i,j,k,n,sp:longint;
    ans:int64;
    f:array[0..maxn,0..100,0..3]of int64;

procedure sort(l,r:longint);
var i,j,x,y:longint;
begin
  i:=l;j:=r;x:=a[(i+j)div 2];
  repeat
    while a[i]<x do inc(i);
    while x<a[j] do dec(j);
    if i<=j then
    begin
      y:=a[i];a[i]:=a[j];a[j]:=y;
      inc(i);dec(j);
    end;
  until i>j;
  if i<r then sort(i,r);
  if l<j then sort(l,j);
end;

procedure pre;
var i:longint;
begin
  sort(1,n);
  cnt:=1;num[1]:=1;
  for i:=2 to n do
  begin
    if a[i]-a[i-1]>1 then
    begin
      writeln(0);halt;
    end;
    if a[i]<>a[i-1] then inc(cnt);
    inc(num[cnt]);
  end;
  c[0,0]:=1;
  for i:=1 to 200 do
  begin
    c[i,0]:=1;
    for j:=1 to i do c[i,j]:=(c[i-1,j]+c[i-1,j-1])mod mo;
  end;
end;

begin
  readln(n);
  for i:=1 to n do read(a[i]);
  pre;
  f[1,1,3]:=1;
  for i:=1 to cnt-1 do
  begin
    for j:=1 to num[i] do
    begin
      sp:=num[i]-j;
      if f[i,j,0]>0 then
        for k:=1 to sp do
          f[i+1,k,0]:=(f[i+1,k,0]+f[i,j,0]*c[sp,k]mod mo*c[num[i+1]-1,k-1]mod mo)mod mo;
      if f[i,j,1]>0 then
      begin
        for k:=1 to sp do
          f[i+1,k,0]:=(f[i+1,k,0]+f[i,j,1]*c[sp,k]mod mo*c[num[i+1]-1,k-1]mod mo)mod mo;
        for k:=0 to sp do
          f[i+1,k+1,1]:=(f[i+1,k+1,1]+f[i,j,1]*c[sp,k]mod mo*c[num[i+1]-1,k]mod mo)mod mo;
      end;
      if f[i,j,2]>0 then
      begin
        for k:=1 to sp do
          f[i+1,k,0]:=(f[i+1,k,0]+f[i,j,2]*c[sp,k]mod mo*c[num[i+1]-1,k-1]mod mo)mod mo;
        for k:=0 to sp do
          f[i+1,k+1,2]:=(f[i+1,k+1,2]+f[i,j,2]*c[sp,k]mod mo*c[num[i+1]-1,k]mod mo)mod mo;
      end;
      if f[i,j,3]>0 then
      begin
        for k:=1 to sp do
          f[i+1,k,0]:=(f[i+1,k,0]+f[i,j,3]*c[sp,k]mod mo*c[num[i+1]-1,k-1]mod mo)mod mo;
        for k:=0 to sp do
          f[i+1,k+1,1]:=(f[i+1,k+1,1]+f[i,j,3]*c[sp,k]mod mo*c[num[i+1]-1,k]mod mo)mod mo;
        for k:=0 to sp do
          f[i+1,k+1,2]:=(f[i+1,k+1,2]+f[i,j,3]*c[sp,k]mod mo*c[num[i+1]-1,k]mod mo)mod mo;
        for k:=2 to sp+2 do
          f[i+1,k,3]:=(f[i+1,k,3]+f[i,j,3]*c[sp,k-2]mod mo*c[num[i+1]-1,k-1]mod mo)mod mo;
      end;
    end;
  end;
  ans:=0;
  for i:=1 to num[cnt] do
    ans:=(ans+f[cnt,i,0]+f[cnt,i,1]+f[cnt,i,2]+f[cnt,i,3])mod mo;
  writeln(ans);
end.