bzoj1801

题目就是每行每列最多放两个炮的意思；

首先不难想到状态压缩dp，但是当n,m<=100的时候显然会跪掉；

考虑每行最多就2个点，状压dp浪费了大量的空间

由于每行最多两个点，我们可以直接用f[i,j,k]表示状态到第i行放1个炮有j列，2个炮有k列

方程具体见程序

 const re=;

 var f:array[..,..,..] of int64;

     i,j,n,p,k,m:longint;

     ans:int64;

 function c2(x:int64):int64;

   begin

     c2:=x*(x-) div  mod re;

   end;

 begin

   readln(n,m);

   f[,,]:=;

   ans:=;

   p:=;

   for i:= to n do

   begin

     p:=-p;

     for j:= to m do

       for k:= to m-j do

       begin

         f[p,j,k]:=f[-p,j,k];

         if j> then

           f[p,j,k]:=(f[p,j,k]+f[-p,j-,k]*(m-j-k+) mod re) mod re;

         if (k>) and (j<m) then

           f[p,j,k]:=(f[p,j,k]+f[-p,j+,k-]*(j+) mod re) mod re;

         if (j>) then

           f[p,j,k]:=(f[p,j,k]+f[-p,j-,k]*c2(m-j-k+) mod re) mod re;

         if (j>) and (k>) then

           f[p,j,k]:=(f[p,j,k]+f[-p,j,k-]*(m-j-k+)*j mod re) mod re;

         if (k>) and (j<m-) then

           f[p,j,k]:=(f[p,j,k]+f[-p,j+,k-]*c2(j+) mod re) mod re;

       end;

   end;

   for i:= to m do

     for j:= to m-i do

       ans:=(ans+f[p,i,j]) mod re;

   writeln(ans);

 end.

 const re=;

 var f:array[..,..,..] of int64;
     i,j,n,p,k,m:longint;
     ans:int64;
 function c2(x:int64):int64;
   begin
     c2:=x*(x-) div  mod re;
   end;

 begin
   readln(n,m);
   f[,,]:=;
   ans:=;
   p:=;
   for i:= to n do
   begin
     p:=-p;
     for j:= to m do
       for k:= to m-j do
       begin
         f[p,j,k]:=f[-p,j,k];
         if j> then
           f[p,j,k]:=(f[p,j,k]+f[-p,j-,k]*(m-j-k+) mod re) mod re;
         if (k>) and (j<m) then
           f[p,j,k]:=(f[p,j,k]+f[-p,j+,k-]*(j+) mod re) mod re;
         if (j>) then
           f[p,j,k]:=(f[p,j,k]+f[-p,j-,k]*c2(m-j-k+) mod re) mod re;
         if (j>) and (k>) then
           f[p,j,k]:=(f[p,j,k]+f[-p,j,k-]*(m-j-k+)*j mod re) mod re;
         if (k>) and (j<m-) then
           f[p,j,k]:=(f[p,j,k]+f[-p,j+,k-]*c2(j+) mod re) mod re;
       end;
   end;
   for i:= to m do
     for j:= to m-i do
       ans:=(ans+f[p,i,j]) mod re;
   writeln(ans);
 end.