怪盗基德之最后の战 vijos1614

vijos1614

 

我们来考虑这样一个方阵：

在左下角放上1，然后每个格子的左边即为这个格子中数的4倍，上面为5倍，则整个方阵就成为：

【下面是4*4的例子】

125	500	2000	8000
25	100	400	1600
5	20	80	320
1	4	16	64

于是我们就发现，剩下的牌都被完美地填进了这个方阵之中！

而实验证明，当k=12时，方阵已全部包含2^22内数字。因此方阵最大只需12*12。

接下去我们要解决的就是在一个方阵中取若干个格子，任意两个格子都没有公共边。

显而易见，第i行怎么种，对前i-1行是没有影响的，因而符合无后效性原则。当然啦，这道题是没有所谓的最优解的。因此我们想到了什么？那就是递推~~

状态非常好写：f[i,j]表示第i行是j状态。j是二进制数，即一串01串，第i位是1表示这个格子取。那么f[i,j]=sigma(f[i-1,k])其中j和k能够并存。

那么怎么判断j和k能够并存呢？j和k能够并存说明它们没有同一位上是1,即j and k=0。

当然还要判断j状态是否可行。比如说n=300，则第二行的状态1010是不可取的状态。也就是说，某列可以用能否取化为一个串0000…011…1，设它为g[i]【g[i]的01串中0为能取，1为不能取】。则必须满足g[i] and j=0。可以预处理a[i,j]表示第i行j状态是否可行。最后只需把第n行能取的状态的总数全加起来即可。

【注意】还要判断，某行是否有两个连起来的格子被取，这个可以这样判断：

对状态j，只要j shr 1 and j=0即可

由于相邻的两个格子不能同时取，每行的状态被限制在了521个。同时还有不少格子本来就不能取，使得状态数进一步地减少了，因而最终时间复杂度最差也是O（512^2*12）完全可以承受了。

const nn=100000000; var map:array[0..13,0..13]of boolean; c:array[1..13]of longint; g:array[0..13,0..133]of int64; f:array[0..13,0..4096]of int64; bo:array[0..13,0..4096]of boolean; i,j,ans,n:longint; procedure dfs(i,j,now,new:longint); begin if j>12 then begin bo[i+1,new]:=true; f[i+1,new]:=f[i+1,new]+f[i,now]; if f[i+1,new]>nn then f[i+1,new]:=f[i+1,new] mod nn; exit; end; dfs(i,j+1,now,new); if (map[i+1,j])and((now and c[j])=0) then dfs(i,j+2,now,new+c[j]); end; begin g[1,1]:=1; for j:=2 to 12 do g[1,j]:=g[1,j-1]*4; for i:=2 to 12 do for j:=1 to 12 do g[i,j]:=g[i-1,j]*5; readln(n); for i:=1 to 12 do for j:=1 to 12 do if g[i,j]<=n then map[i,j]:=true; c[1]:=1; for i:=2 to 13 do c[i]:=c[i-1]*2; bo[0,0]:=true; f[0,0]:=1; for i:=0 to 11 do for j:=0 to c[13]-1 do if bo[i,j] then begin dfs(i,1,j,0); end; for j:=0 to c[13]-1 do if bo[12,j] then begin ans:=ans+f[12,j]; if ans>nn then ans:=ans mod nn; end; writeln(ans mod nn); end.