石業(yè)嬌, 孟憲濤

(1. 大連海洋大學(xué) 應(yīng)用技術(shù)學(xué)院 遼寧 大連 116300; 2. 沈陽師范大學(xué) 數(shù)學(xué)與系統(tǒng)科學(xué)學(xué)院, 沈陽 110034)

理論與應(yīng)用研究

用代數(shù)方法探討四圓相交區(qū)域填充數(shù)字問題

石業(yè)嬌1, 孟憲濤2

(1. 大連海洋大學(xué) 應(yīng)用技術(shù)學(xué)院 遼寧 大連 116300; 2. 沈陽師范大學(xué) 數(shù)學(xué)與系統(tǒng)科學(xué)學(xué)院, 沈陽 110034)

受在三圓相交區(qū)域中填充數(shù)字問題以及用代數(shù)方法求解幻方問題的啟發(fā),設(shè)計了利用線性代數(shù)方法在四圓相交區(qū)域中填充數(shù)字問題。首先,建立了填充問題的約束方程組,根據(jù)需要將約束方程組變形為5種形式,即所謂5個約束條件;然后,對約束條件進(jìn)行討論,得出四圓重疊區(qū)域中心位置的數(shù)與兩圓重疊區(qū)域的4個數(shù)字之和的奇偶性,以及三圓重疊區(qū)域的4個數(shù)字之和與只屬于一個圓區(qū)域的4個數(shù)字之和的奇偶性,以約束條件為基礎(chǔ),兼顧數(shù)字的對稱性與互補性,采用試驗的方法,考慮3種情況下的不同取值,得到相應(yīng)問題的15個解;最后,給出了相對于每一個解,每一個圓中所包含的7個數(shù)字之和的上限與下限,給出相應(yīng)的證明。

四圓相交區(qū)域; 填充數(shù)字; 約束方程組; 求解

在幾個圓構(gòu)成的相交區(qū)域中填充數(shù)字問題與幻方問題頗為類似,幻方為中國人首創(chuàng),這一點在漢朝的《數(shù)術(shù)記遺》中有明確的記載。我國宋朝數(shù)學(xué)家楊輝對幻方的研究頗有建樹,他的研究往往是針對于某階幻方直接給出構(gòu)造方法,即直接給出幻方的解,雖然巧妙致極,但總有一點不知其所以然的感覺。筆者曾用線性代數(shù)的方法研究了三階幻方的所有解。受幻方問題的啟發(fā),本文則嘗試?yán)镁€性代數(shù)方法探討在四個圓相交區(qū)域里填充數(shù)字問題。

1 問題描述

四圓相交,圍成13個區(qū)域如圖1?，F(xiàn)將1,2,…,13這13個數(shù)字填在圖1中,每一個區(qū)域里一個數(shù)字,要求各個圓中所包含的7個數(shù)字之和相等。此問題雖然與三圓相交區(qū)域填充數(shù)字類似,卻比解決三圓相交區(qū)域填充數(shù)字問題困難許多。因為圖中的變數(shù)為13個,根據(jù)問題描述要求可能確定的約束方程的個數(shù)卻遠(yuǎn)少于變量的個數(shù)。因此問題必須在所建立的約束條件基礎(chǔ)上結(jié)合試驗方法加以解決。

圖1 四圓相交填充Fig.1 Filling method of four circles intersection

2 約束條件

記S={a1,a2,…,a13}={1,2,…,13}。設(shè)圖1的填法滿足問題要求,設(shè)每個圓中所包含的7個數(shù)字之和為M,即

(1)

得到第一個約束條件

(2)

從式(2)推得

(3)

知a1+a3+a10+a12+a13與a2+a4+a9+a11+a13分別為奇數(shù)。由a1,a3,a10,a12,a13與a2,a4,a9,a11,a13分別為圖1中過a13的直線上的5個數(shù),得第2個約束條件為這5個數(shù)之和為奇數(shù)。

由式(1)可得

(4)

代入式(2)中第1式便得到

(5)

觀察式(5)可知,當(dāng)a13為奇數(shù)時,3a13+91為偶數(shù),于是a5+a6+a7+a8為偶數(shù);當(dāng)a13為偶數(shù)時,3a13+91為奇數(shù),于是a5+a6+a7+a8為奇數(shù)。注意到a5+a6+a7+a8所處的位置是圖1中兩圓重疊的區(qū)域,便得到第3個約束條件為“圖1中四圓重疊區(qū)域即中心位置的數(shù)a13與圖1中處在兩圓重疊區(qū)域的4個數(shù)字之和a5+a6+a7+a8的奇偶性相反”。

由式(4)有

(6)

可見當(dāng)a1+a2+a3+a4為奇數(shù)時,a9+a10+a11+a12為奇數(shù);當(dāng)a1+a2+a3+a4為偶數(shù)時,a9+a10+a11+a12為偶數(shù),有第4個約束條件為“圖1中三圓重疊區(qū)域的四數(shù)字之和與只屬于1個圓區(qū)域的4個數(shù)字之和具有相同的奇偶性”。

對方程組(2)的增廣矩陣A進(jìn)行初等變換,得式(2)的同解方程組

(7)

方程組(7)即是第5個約束條件。

應(yīng)該注意的是,以上得到的5個約束條件并不是相互“獨立”的,如式(2)與式(7)是同解方程組,于是約束條件1與約束條件5實際上是同一個約束條件。同理約束條件2、3與約束條件4也屬于約束條件一變形。因而這5個約束條件實際就是約束條件1(即式(2))的5種不同的表現(xiàn)形態(tài)。之所以如此,是為了在以下的對圖1的解法探討中便于從不同側(cè)面加以約束,更方便求解。

3 求解舉例

滿足方程組(7)的解為問題的解。求問題的解,必須在約束條件的基礎(chǔ)上輔之以試驗方法。試驗取值要充分關(guān)注1,2,…,13這13數(shù)的分布特點,注意數(shù)字之間的對稱性與互補性。

如果取a13=1,考慮了數(shù)字間的對稱與互補,取a9=4,a10=5,a11=2,a12=3填入圖1中,再考慮補償關(guān)系及相關(guān)約束條件,取a5=6,a6=7,a7=8,a8=9填入圖1中,把這些數(shù)代入式(7),求得

圖2 滿足方程組解的第1種填充方法Fig.2 The first filling method of satisfying equations

將求得的a1,a2,a3,a4的值填入圖1,得圖2。

檢驗圖2中各圖所包含的7個數(shù)字之和,得M=38。于是a1=13,a2=12,a3=11,a4=10,a5=6,a6=7,a7=8,a8=9,a9=4,a10=5,a11=2,a12=3,a13=1為方程組(2)的一個解,即圖2為一種填法。

取a9=5,a10=4,a11=3,a12=3。考慮約束條件,取a5=10,a6=11,a7=12,a8=13,代入到式(7)中,得a1=8,a2=9,a3=6,a4=7,填入四圓相交區(qū)域得圖3。經(jīng)檢驗知,圖3為滿足問題要求的一種填法,每個圓內(nèi)所含有的7個數(shù)字之和均為42。

取a9=8,a10=9,a11=6,a12=7,a5=2,a6=3,a7=4,a8=5。代入式(7)中,得a1=13,a2=12,a3=11,a4=10,填入四圓相交區(qū)域得圖4。檢驗可知M=42,即圖4為問題的一個解。

圖3 滿足方程組解的第2種填充方法

圖4 滿足方程組解的第3種填充方法

取a9=9,a10=8,a11=7,a12=6,a5=10,a6=11,a7=12,a8=13。由式(7)求得a1=4,a2=5,a3=2,a4=3,填入四圓相交區(qū)域得圖5。檢驗之,M=50。于是圖5為符合問題要求的填法。

如果取a13=13,考慮數(shù)字間對稱與互補,取a9=2,a10=3,a11=4,a12=1,a5=6,a6=5,a7=8,a8=7。由式(7)求得a1=11,a2=10,a3=9,a4=12,得圖6填法。經(jīng)檢驗,M=44,因此圖6填法是問題的一個解。

圖5 滿足方程組解的第4種填充方法

圖6 滿足方程組解的第5種填充方法

如果取a13=7時,可給出3種填法,對應(yīng)的M值分別為40,44,58,……。

以此類推,可得出四圓相交區(qū)域填充數(shù)字問題的15種解法,這里略述。

4 結(jié) 語

以上給出四圓相交區(qū)域填充數(shù)字問題的幾種解法,發(fā)現(xiàn)每一種解法中M值均滿足38≤M≤60。事實上,由式(5)可推得

(5+6+7+8)+2(9+10+11+12)+3×13+91=240

以及

(9+8+7+6)+2(5+4+3+2)+3×1+91=152

即有

152≤4M≤240

從而

38≤M≤60

以上考慮3種情況下的不同取值,得到相應(yīng)問題的15個解,還可以考慮其他情況下的相應(yīng)問題的解。

[ 1 ]吳鶴齡. 幻方與素數(shù)----娛樂數(shù)學(xué)兩大經(jīng)典命題[M]. 北京:科學(xué)出版社, 2008.

[ 2 ]許莼舫. 古算趣味[M]. 北京:中國青年出版社, 1955.

[ 3 ]錢寶琮. 中國數(shù)學(xué)史[M]. 北京:科學(xué)出版社, 1964.

[ 4 ]姜偉,劉彥佩. 幾類4-正則平面圖的最小折數(shù)縱橫擴張[J]. 沈陽師范大學(xué)學(xué)報(自然科學(xué)版), 2007,25(2):129-134.

[ 5 ]李冠林,顧大全. 幻方的實現(xiàn)方法研究[J]. 微型電腦應(yīng)用, 2010,26(1):17-18.

[ 6 ]盛中平,王曉輝,朱本喜. 有理插值的擴展方程組與約束方程組[J]. 高等學(xué)校計算數(shù)學(xué)學(xué)報, 2005(1):85-96.

[ 7 ]陳劍南. 素數(shù)階均衡完美幻方若干問題初探[J]. 計算機工程與應(yīng)用, 2009,45(21):179-182.

[ 8 ]歐陽錄. 最佳拉丁方與高級原幻方[J]. 數(shù)學(xué)理論與應(yīng)用, 2001,21(3):22-28.

[ 9 ]楊富鋒. 構(gòu)造奇次同心幻方的一種方法[J]. 數(shù)學(xué)的實踐與認(rèn)識, 2006,36(5):192-199.

[10]祝寶滿,龔和林. 非素數(shù)階幻方的構(gòu)造[J]. 數(shù)學(xué)的實踐與認(rèn)識, 2008(15):207-214.

[11]李紅林,劉昆. 基于圓的對稱等分點畫線算法的圓內(nèi)區(qū)域填充[J]. 云南民族大學(xué)學(xué)報(自然科學(xué)版), 2013(4):292-294,298.

[12]李向利. 幾類帶界約束方程組的算法研究[D]. 西安:西安電子科技大學(xué),2011.

[13]段云嶺. 非線性方程組的解法:局部弧長法[J]. 力學(xué)學(xué)報, 1997(1):117-123.

[14]尹釗,賈尚暉. Moore-Penrose廣義逆矩陣與線性方程組的解[J]. 數(shù)學(xué)的實踐與認(rèn)識, 2009(9):239-244.

Algebraic method for filling numerical problems of intersection region of four circles

SHIYejiao1,MENGXiantao2

(1. Applied Technology College, Dalian Ocean University, Dalian 116300, China; 2. School of Mathematics and System Science, Shenyang Normal University, Shenyang 110034, China)

Inspired by the filling numbers of three circles in the intersectional region and algebraic method to solve the problem of magic square, the problem of filling numbers in intersectional regions of four circles is designed based on linear algebra methods. Firstly, the constraint equations of the filling problem are established, and they are deformed into five forms, namely, the five constraints. Then, upon discussion of the constraints, odevity of the number in intersectional region of four circles and summation of the four numbers in intersectional region of two circles as well as the odevity of summation of the four numbers in intersectional region of three circles and the four numbers in only one circle are obtained. On the basis of constraints, considering the symmetry and complementation of the numbers and different values in three cases, the 15 solutions of the problem are obtained by the method of experiment. Finally, the upper and lower bounds of the sum of 7 numbers contained in each circle are given to each solution. Proofs are provided respectively.

intersectional region of four circles; fill numbers; constraint equations; solve

2017-04-27。

國家自然科學(xué)基金資助項目(11201313)。

石業(yè)嬌(1970-),女,遼寧大連人,大連海洋大學(xué)副教授。

1673-5862(2017)03-0335-04

O151.26

A

10.3969/ j.issn.1673-5862.2017.03.014