数学建模社区-数学中国

标题: 残缺棋盘 [打印本页]

作者: 韩冰    时间: 2004-10-4 05:16
标题: 残缺棋盘
<b>残缺棋盘</b>. h! p; y3 \  q' l  t
<>残缺棋盘(defective chessboard)是一个有2k×2k 个方格的棋盘,其中恰有一个方格残缺。图2 - 3给出k≤2时各种可能的残缺棋盘,其中残缺的方格用阴影表示。注意当k= 0时,仅存在一种可能的残缺棋盘(如图1 4 - 3 a所示)。事实上,对于任意k,恰好存在22k 种不同的残缺棋盘。
/ x" ~  [  A8 E9 p6 v
1 ^. Y0 l3 v# M4 n1 V5 ^残缺棋盘的问题要求用三格板(t r i o m i n o e s)覆盖残缺棋盘(如图1 4 - 4所示)。在此覆盖中,两个三格板不能重叠,三格板不能覆盖残缺方格,但必须覆盖其他所有的方格。在这种限制条件下,所需要的三格板总数为( 22k -1 ) / 3。可以验证( 22k -1 ) / 3是一个整数。k 为0的残缺棋盘很容易被覆盖,因为它没有非残缺的方格,用于覆盖的三格板的数目为0。当k= 1时,正好存在3个非残缺的方格,并且这三个方格可用图1 4 - 4中的某一方向的三格板来覆盖。
/ s& y2 o) |2 B+ _, k& A' \  p
3 k! c* i' i8 D0 C) O用分而治之方法可以很好地解决残缺棋盘问题。这一方法可将覆盖2k×2k 残缺棋盘的问题转化为覆盖较小残缺棋盘的问题。2k×2k 棋盘一个很自然的划分方法就是将它划分为如图1 4 - 5 a所示的4个2k - 1×2k - 1 棋盘。注意到当完成这种划分后, 4个小棋盘中仅仅有一个棋盘存在残缺方格(因为原来的2k×2k 棋盘仅仅有一个残缺方格)。首先覆盖其中包含残缺方格的2k - 1×2k - 1 残缺棋盘,然后把剩下的3个小棋盘转变为残缺棋盘,为此将一个三格板放在由这3个小棋盘形成的角上,如图14-5b 所示,其中原2k×2k 棋盘中的残缺方格落入左上角的2k - 1×2k - 1 棋盘。可以采用这种分割技术递归地覆盖2k×2k 残缺棋盘。当棋盘的大小减为1×1时,递归过程终止。此时1×1的棋盘中仅仅包含一个方格且此方格残缺,所以无需放置三格板。- [4 z# P7 F! x7 d
& {$ G- ]5 D' X) \6 o
可以将上述分而治之算法编写成一个递归的C++ 函数Ti l e B o a r d (见程序1 4 - 2 )。该函数定义了一个全局的二维整数数组变量B o a r d来表示棋盘。B o a r d [ 0 ] [ 0 ]表示棋盘中左上角的方格。该函数还定义了一个全局整数变量t i l e,其初始值为0。函数的输入参数如下:
, \7 ^5 u# H, s+ P
& }' i( q9 G/ b? tr 棋盘中左上角方格所在行。
' M/ v$ W' s9 A4 a: |5 H- ?) F8 I; D: L6 e& J) L. M
? tc 棋盘中左上角方格所在列。8 D. l$ @: R+ q) m% d7 z, p4 H
( A$ G$ H& F+ o2 C
? dr 残缺方块所在行。# H: V. N1 x/ `; e) _
1 s( K; A9 d! V1 d3 M4 Q
? dl 残缺方块所在列。5 c, R. l8 f$ f0 ^8 i

" [4 n4 U$ K1 S? size 棋盘的行数或列数。9 T- g4 G+ m: N7 m: F* E

+ H) Y4 Q, X' p+ E" STi l e B o a r d函数的调用格式为Ti l e B o a r d(0,0, dr, dc,size),其中s i z e = 2k。覆盖残缺棋盘所需要的三格板数目为( s i z e2 -1 ) / 3。函数TileBoard 用整数1到( s i z e2-1 ) / 3来表示这些三格板,并用三格板的标号来标记被该三格板覆盖的非残缺方格。/ h) n4 ?+ o% V* N/ @
7 {- A/ E% e; G, V7 y  u: X
令t (k) 为函数Ti l e B o a r d覆盖一个2k×2k 残缺棋盘所需要的时间。当k= 0时,s i z e等于1,覆盖它将花费常数时间d。当k &gt; 0时,将进行4次递归的函数调用,这些调用需花费的时间为4t (k-1 )。除了这些时间外, if 条件测试和覆盖3个非残缺方格也需要时间,假设用常数c 表示这些额外时间。可以得到以下递归表达式:/ f9 o; \, l9 B6 V4 u4 i& P) A
& o/ J: D( N. _' }# k
程序14-2 覆盖残缺棋盘
( l; L, X2 ^* b2 D0 K- n7 X# B( j! d! S! ]8 K- ]
void TileBoard(int tr, int tc, int dr, int dc, int size)# ?, ~5 v- u* p" z1 a

0 t6 ?% e& L/ b- K{// 覆盖残缺棋盘) Q) m4 C# {7 I9 C& n* W0 A
) N6 X, k% H" c+ c! O* c4 z/ C
if (size == 1) return;3 k% I: }4 B2 l0 ]. B3 G, B
& D% h8 k+ w  x- V0 ~' A
int t = tile++, // 所使用的三格板的数目
& E$ C" t  C: J- w$ `& f) y4 @' U5 Z. @1 D
s = size/2; // 象限大小
( ~- ~& X" d# O/ t
& ?3 p* \" E1 e/ /覆盖左上象限
) \. V) M5 X/ e9 o0 z
. \) Z! d, Q3 ]% U1 s0 sif (dr &lt; tr + s &amp;&amp; dc &lt; tc + s)
4 O. u. h" \  C/ }! R
$ B% a) C4 a( @. m4 R# |' D' H  H// 残缺方格位于本象限9 B( @& c8 J! _0 {# i

; p; p2 N: |1 _7 C) BTi l e B o a r d ( t r, tc, dr, dc, s);, F% P) p: n. n# V

: m" U; P. ?; C0 A  ^* c9 C/ {else {// 本象限中没有残缺方格7 L# M# O: l- L- b3 W- D
8 E0 C& h5 Q! e
// 把三格板t 放在右下角2 Q% o7 J! i: f0 w
% k+ D; f" y  B: i1 A
Board[tr + s - 1][tc + s - 1] = t;
* y* D+ P* G/ Z7 Z1 X/ n- d0 [3 R' X0 ?5 R
// 覆盖其余部分
6 N  A3 H3 e; Z' F5 O. M* W/ \4 ~
$ T% e" a2 G/ F0 Z: D' E. iTi l e B o a r d ( t r, tc, tr+s-1, tc+s-1, s);}% C8 }+ q& V; C3 b' O/ d
2 s4 @4 W+ r# ~
/ /覆盖右上象限
" c, U2 \! Z. o8 c0 K- Y" W5 Q' N0 d! }* j
if (dr &lt; tr + s &amp;&amp; dc &gt;= tc + s)
9 Z, Q) D) I4 s
  E4 M" {5 L- E- f( p5 j// 残缺方格位于本象限
0 X/ B" C1 `; l; e# x4 A+ b' V9 ]3 ~
Ti l e B o a r d ( t r, tc+s, dr, dc, s);
9 W" a6 t! K/ `6 O. b
9 n' Y  {# J$ J0 Zelse {// 本象限中没有残缺方格
6 x; U2 e7 o  D* _( M0 L" j
, [3 r) j$ J& N: R// 把三格板t 放在左下角
3 l8 C& M) p/ m+ v
. }: ?5 A2 @9 bBoard[tr + s - 1][tc + s] = t;, {/ R) u- y! P, z
! b" W8 c! }; G4 D
// 覆盖其余部分
9 r) N0 R8 |) p, n& n2 @: x% S9 E
0 z+ J" t# K  S8 H- j* e4 g8 vTi l e B o a r d ( t r, tc+s, tr+s-1, tc+s, s);}
: Y3 a7 u: \( n. h$ Z- z  C" q& Z1 {7 d9 E0 R6 p
/ /覆盖左下象限9 m9 H" T  w  I7 c7 E) E& d: }

& e( w8 ^2 W3 @: R$ Qif (dr &gt;= tr + s &amp;&amp; dc &lt; tc + s)7 ^- X& C; q- X0 U1 S2 S( k7 B
% x, p1 k% u( Z% @5 `
// 残缺方格位于本象限
7 r' }0 V2 M4 |/ [& a
9 \7 V+ s: Y  ]9 W- dTileBoard(tr+s, tc, dr, dc, s);7 d. ]) A2 {+ K, F8 Q: k3 H
; Y3 a, n+ Z! [2 R/ @; [
else {// 把三格板t 放在右上角3 \% }# G; ?. t1 v3 f* A- c5 \

  K9 g( f5 e5 C& H. I6 _" y' bBoard[tr + s][tc + s - 1] = t;
5 W% J9 k& M6 t5 L  A5 D4 B
2 k  w( Q2 e% ~( A# A/ ?// 覆盖其余部分
0 t* C% Z6 i0 L! o: h7 ~4 d4 v7 c# T* G1 [5 L
TileBoard(tr+s, tc, tr+s, tc+s-1, s);}
1 I6 n  o1 Z9 {# O9 Z+ ?% G$ @3 c  I  ?
// 覆盖右下象限6 l4 u& e  J3 b* S& K  y+ l

: ~0 y+ c2 M( Y7 N# a+ Sif (dr &gt;= tr + s &amp;&amp; dc &gt;= tc + s)  r+ ?+ j3 Y; N1 ~9 B# T
( M: ]" ?; }& ?$ X0 X- b
// 残缺方格位于本象限( W* Y& p( [# [# q1 v6 L9 q

  _# E1 |9 c. l3 a/ tTileBoard(tr+s, tc+s, dr, dc, s);* t. z6 R$ V2 p
6 s4 o5 N5 k) R5 I4 i& C1 ^  i6 d
else {// 把三格板t 放在左上角% [  N+ `# G1 d  ?9 i; d8 s' {

: ^# y' R2 f2 ~1 V* x1 cBoard[tr + s][tc + s] = t;
3 l. |, L: f7 _- w, I
, Z) C& w% {% D# D. f" a// 覆盖其余部分0 b* I7 S2 V7 u# C
/ S5 e5 b+ |) n" ?
TileBoard(tr+s, tc+s, tr+s, tc+s, s);}
* Z6 T* q" \' ?0 p1 _) Z6 x3 d: M0 n( m8 K3 W+ M; p: {
}
6 k# ~; ]% j2 O& v2 T  {8 _* H2 g
( c; H) z; e! s; B* K) cvoid OutputBoard(int size), A( p6 b. |! {6 X$ a: M0 y
6 ~* j2 t8 C8 t  D* K! q
{- u( Z3 w5 V& C# d! q' Q! W

" D( S5 i  [# c$ @& h% \for (int i = 0; i &lt; size; i++) {
! G+ D5 e7 L/ V! Y9 @. h
( A7 \5 S: ^( u+ I8 w0 d! ifor (int j = 0; j &lt; size; j++)- S( q2 C  e$ I. E- \: q2 m
* t, h) q* K" ?5 P% ]5 t
cout &lt;&lt; setw (5) &lt;&lt; Board[j];
; M) c) d4 Y4 _7 E* f2 P6 e7 V  M  M0 U. x! d7 [
cout &lt;&lt; endl;2 c6 n! |3 D+ k5 o7 M
0 o+ R! Q. C& L! s. U
}
8 J8 _0 r; R( p
/ j9 ~" y9 g2 i) u, w}2 C5 U) s0 ?9 a) }' o! K2 O0 r
$ r0 y. R; y4 u6 P6 U1 |: C
可以用迭代的方法来计算这个表达式(见例2 - 2 0),可得t (k )= ( 4k )= (所需的三格板的数目)。由于必须花费至少( 1 )的时间来放置每一块三格表,因此不可能得到一个比分而治之算法更快的算法。</P>
作者: yammay    时间: 2005-4-3 17:22
<>太感谢啦。找了好久才找到这个程序。真是真是太感谢啦。有点激动</P>
作者: cupidvenus    时间: 2005-9-27 13:13
没看明白,图呢?
作者: wendy28    时间: 2005-9-29 00:46
我们学过了,很经典的一道题!




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5