• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    初等數(shù)論在歷史年代學中的應(yīng)用

    2015-12-23 03:28:57羅見今
    關(guān)鍵詞:年代學干支甲子

    羅見今

    (內(nèi)蒙古師范大學科學技術(shù)史研究院,呼和浩特010022)

    初等數(shù)論在歷史年代學中的應(yīng)用

    羅見今

    (內(nèi)蒙古師范大學科學技術(shù)史研究院,呼和浩特010022)

    應(yīng)用初等數(shù)論中的同余式、連分數(shù)、漸近分數(shù)、不定方程等工具,分析考古年代學提出的若干問題,以敦煌漢簡、居延漢簡、《史記》等出土文物和史料為證據(jù),將《二十史朔閏表》、《三千五百年歷日天象》的記錄視為大時段日月運行的數(shù)據(jù)庫,證明或論述中國歷法的五個周期性質(zhì),從而闡明應(yīng)用數(shù)學進行年代考釋的基本原理和方法。

    初等數(shù)論;漢代簡牘;朔閏表;歷法的周期性質(zhì);考古年代學

    一、研究的背景和問題的提出

    當人們研究一個歷史問題或遇到一件出土文物的時候,首先須回答:它是什么朝代、什么時期的?有時還須進而了解:它是哪一年、哪一月的?這就是史學中的斷代問題,在整個歷史科學、特別是考古學中占有重要的地位。歷史、考古、文物工作者在長期的實踐和研究工作中具有豐富的經(jīng)驗,形成了年代學的基本理論,對大量年代學實際問題作出了科學而準確的判斷。人們現(xiàn)在所知關(guān)于歷史上的一切,幾乎都建立在年代學的框架之中。這個基本框架經(jīng)歷了千千萬萬史實的驗證,是不可動搖的。

    然而人們也不能忘記,史學家們所面對的是五千年悠久的歷史。由于自然和人的作用,關(guān)于古代紀年的信息所存無幾;隨著時間的流逝,“過去”離我們愈加遙遠、變得模糊不清。大量的問題不是“拂去塵封”就可以探知其真的。歷史是復雜的,它充滿了矛盾。人們有理由懷疑史書中個別年月日的記錄,因它是否后來所修,或有改寫的動機;也可以指謫歷譜中的錯誤,認為它主要靠推算而得,不見得是當時的應(yīng)用歷譜;或者把希望寄托在出土的文物之上:它的年號、月朔記錄應(yīng)當是真實的。然而,出土文物的記錄也可能出現(xiàn)牴牾,例如:1930年在今內(nèi)蒙古西部出土的“居延漢簡”中有一枝(編號137·3,224·18)為:[1]

    建平五年(前2)十二月丙寅朔乙亥誠北候長□充□言之官下詒詣(符號“□”表示原來木簡上的字難以辨認,“”表示以下部分斷缺。下略)。

    但是另有一枝(編號495·12A,506·20A)卻是:[2]

    兩簡同是同年十二月,但朔日(初一日)的干支卻相差25天!究竟哪一枝記錄可信呢?如果查一查史學界通用的、陳垣先生編《二十史朔閏表》,[3]得知建平五年(即元壽元年)十二月乙未朔,與兩簡所記又不相同!如果你繼續(xù)去查宋人劉羲叟的《長歷》、清人汪曰楨的《歷代長木輯要》[4]以及《漢書》[5](第3冊848頁記,該年十二月有庚子日),你會得出互相矛盾的結(jié)論而莫知所從。

    當然,這個例子較典型但并不普遍。由于分治、戰(zhàn)亂、改歷等因素的影響,有些歷史階段會出現(xiàn)各地使用歷譜不完全一致的情況;而竹木簡的記載也可能有筆誤、釋文錯誤等問題。絕大多數(shù)歷史記載、歷譜、文物的時間記錄基本是真實的,但又不是絕對正確的,需要進行具體的分析。

    另一方面,史學家對同一歷史事件年代的判斷可以出現(xiàn)重大分歧。最著名的一例是武王伐紂之年代,從西漢末的劉歆開始,到1997年為止,至少發(fā)表了106種論著,提出了多達44種不同的伐紂之年,其上下限為前1130年到前1018年。[6]如果史學家都爭論不已,一般人更是莫所適從了。所以國家撥巨款,從事“斷代工程”的年代學研究,以確定夏、商、周的年代。其成果當然對于中國歷史具有重要意義。

    事實上,除了這類重大的課題之外,應(yīng)用傳統(tǒng)的史學方法對于判定某些文物的年代、推算某些年代的歷譜等也可能出現(xiàn)錯誤。雖然無關(guān)重大歷史事件,但有悖于歷史的真實,也應(yīng)當予以糾正。在下文中將舉出具體例子說明。更多的情況是傳統(tǒng)的方法無法對文物精確斷代,只好劃出一個大致的時代上下限,使之處于模糊的狀態(tài)。

    隨著現(xiàn)代科學技術(shù)的進步,天文、數(shù)學、物理、化學等方法開始在年代學領(lǐng)域中取得巨大成果,成為斷代的有力工具。用天文學方法進行考古卓有成效,備受中外科學家的重視,已形成了“考古天文學”的分支。利用測定文物放射性的方法確定年代日臻完善,例如1947年利比創(chuàng)放射性同位素C14法測定年代已被大量使用。還有利用紅外線、電腦圖象識別等方法。其中,根據(jù)歷史、考古知識,利用歷法、數(shù)學方法,建立合適的數(shù)學模型,推導出一系列周期定理和判定定理,用來做考證年代的工具,簡便、適用、精確,與傳統(tǒng)方法結(jié)合緊密,值得系統(tǒng)討論。

    眾所周知,中國傳統(tǒng)歷法又稱為“干支歷”。相傳黃帝之臣大撓創(chuàng)甲子,乙丑,……癸亥,干支相配紀日、六十周而復始,是我國先人的偉大創(chuàng)造,其意義可與數(shù)字的發(fā)明相媲美。學術(shù)界一般認為傳統(tǒng)歷法自魯隱公元年(前722)以來,干支紀日至今未間斷,將近一百萬天的連續(xù)記錄,在世界歷法史上絕無僅有,并相信這一天的到來,更多的人們會體會到它的價值。我們把干支紀日的既不重復、又不間斷的性質(zhì)稱為“連續(xù)原理”,它是進行一切歷法、數(shù)學的計算、推理的基礎(chǔ)。

    干支紀日具有獨特之處?,F(xiàn)有一例:20世紀初出土的“敦煌漢簡”中有一枝歷譜簡(編號97),背面記有“閏(十二)月丙申朔大”[7]。因它是散簡,并不知是哪一年的。根據(jù)簡正面十二個月朔的記載,確認此閏月應(yīng)是閏十二月;經(jīng)考證知此歷譜簡屬西漢宣帝神爵三年(前59)??甲C有許多根據(jù);但是,僅根據(jù)“閏(十二)月丙申朔大”這一信息,就可以完全確定它為神爵三年簡,因為它是非常獨特的:在從公元前206年西漢始,至2000年的兩千二百多年中,僅有公元前59年是“閏十二月丙申朔大”。這種唯一性受到閏十二月周期性、丙申朔的周期性、大月的周期性的限定,往往被人忽視。同樣,大量無年號的朔閏干支記錄也都包含著各自的周期性和獨特性秘密,期待著人們?nèi)テ平?。相比之下,西歷如果記有“12月1日”,它可能屬于兩千多年中的任意一年,對于斷代則無所助益。由此凸顯出干支歷的特殊性,對于考釋年代,不啻為解決問題的一把鑰匙。

    中國歷法原被稱為“皇歷”,因制歷權(quán)為朝廷特有,并以皇帝名義頒行天下;后來又叫“農(nóng)歷”。在民間流行稱為“陰歷”(月亮歷),其實是不對的,應(yīng)是陰陽合歷(月亮—太陽歷),因它既含月象盈虧,又合寒暑時令。它的一年有354日(也有355日),閏年384日(也有383、385日);它的一個月大月有30日,小月29日。但實際上1回歸年有365.242 198 78日,而1朔望月有29.530 588 67日。中國歷代天文歷法家經(jīng)過千百年持續(xù)不斷的努力,精確地觀測天象,編制優(yōu)秀的歷法,使之與日月星辰的運行盡量相合。幾千年過去了,中國歷法作為一個整體,平均每年的日數(shù)、平均每月的日數(shù)分別與回歸年、朔望月逐漸相合,本文稱其為“趨同原理”。歷法中的許多周期性質(zhì)是由日、月的周期,干支周期共同決定的,它們建立在趨同原理的基礎(chǔ)之上。

    二、應(yīng)用初等數(shù)論的方法解決歷法問題

    為研究干支的周期性,須了解初等數(shù)論中的同余理論。對于整數(shù)a,b及正整數(shù)m有下列關(guān)系:

    a=mq1+r1(0≤r1≤m),b=mq2+r2,(0≤r2≤m)

    而r1=r2;則說:a,b對于模數(shù)的m為同余,也可以說:對于模數(shù)m,a,b同余,記作a≡b(mod m)上式叫做同余式。例如,42、50被8除,得到的余數(shù)都是2,于是說:42、50對于模數(shù)8為同余,記作42≡50(mod 8);又如-11、16被3除,得到的余數(shù)都是1,于是說-11、16對于模數(shù)3同余,記作-11≡16(mod 3)。由定義可證明a-b可被m整除。

    同余有諸多性質(zhì),包含深刻的內(nèi)容。用它研究干支,可以得到一些計數(shù)的結(jié)果。為此,先將甲子,乙丑,……,癸亥依次用01,02,…,60表示,記作A(或B,N等)每個A都由天干a(a=1,2,…,10)和地支b(b=1,2,…,12)構(gòu)成,于是有

    A=a(mod 10),A=b(mod 12),A=6a-5b(mod 60)

    如果從干支A(a1,b1)計數(shù)到B(a2,b2),得到N(≤60),那么N表示兩干支之間的距離,有

    任意兩干支日之間有多少天?這是一個復雜的問題,計算往往出錯。但實際問題非常之多。例如:

    建武六年七月戊戌朔乙卯甲渠障候敢言府書曰吏民毋得伐樹木……(EPT22.53A)[8]

    要問:乙卯是幾號?這是公元30年7月1日戊戌(5,11)到乙卯(2,4)日,上式表明可用兩種方法計算:乙卯B=52,戊戌A=35,據(jù)式可求出N=18;在不記得A,B的干支序號時,據(jù)式亦可求出N=18。

    通常這種求日期算法均用歷譜或60干支表,一個個口數(shù)心計;列出算式可制程序交給計算機去完成的。

    干支紀日具有順序意義、循環(huán)意義、數(shù)字意義。60甲子周期同回歸年和朔望月的日數(shù)(可以看作是無理數(shù))不可公度,也就是不能找出三者之間都恰好合適的公倍數(shù)。在初等數(shù)論中便借助于連分數(shù)、漸近分數(shù)的方法來解決。

    這種方法可以繼續(xù)進行下去,得到一個層層迭迭的分數(shù),叫做連分數(shù)。任何一個實數(shù)都可以表示成連分數(shù)的形式。上例用一種簡單的寫法可記作x=[1,1,1,1,…]。

    應(yīng)用迭代法,將試值x=1代入原式,得出x=1/2;將它再代入原式右邊,得出x=2/3;循此繼進,可求出一系列分數(shù):這同將試值x=1代入連分數(shù)的第1,2,…各層后展開的結(jié)果是相同的。其實,原式本來就是方程,求出的一個根是,即黃金分割,x=0.618 033 988 7。上列分數(shù)值逐漸向這個極限靠近,所以把它稱為漸近分數(shù)。

    我國西漢時已知求連分數(shù),計算歷法問題卓有成效。以下引用華羅庚[9]先生的一個例子:怎樣算農(nóng)歷的閏月?將回歸年和朔望月的日數(shù)相除:

    說明一年應(yīng)有12個朔望月,另外還余0.3683個月。應(yīng)用輾轉(zhuǎn)相除法將它展開成連分數(shù),得

    它的漸近分數(shù)為(用更精確的數(shù)據(jù)獲得)

    這表明:兩年一閏太多,三年一閏太少;八年三閏較多,十一年四閏較少;十九年七閏較多,334年123閏略少……古人常用十九年七閏;,積17個7閏之后,在下個11年中只作4閏,就更合于日、月的運行。以上對初等數(shù)論的同余和連分數(shù)作了簡介,對于了解以下內(nèi)容是必要的。

    三、兩列漸近分數(shù)

    我國古代天文歷法可分為觀象授時時期(上古至春秋中)、歷法準備時期(春秋中至戰(zhàn)國中)、歷法制定時期(戰(zhàn)國中至漢太初元年)。到《三統(tǒng)歷》臻于完備,載于《漢書·律歷志》。簡言之,漢代以來歷法、史書、文物的證據(jù)較先秦為多,在年代學中應(yīng)用數(shù)字方法可以找到歷史記錄作為證據(jù)。

    幾千年來中國傳統(tǒng)歷法是一個整體、一個系統(tǒng),雖然各個時期、每種歷法都有它的特點,但可作為分段演繹系統(tǒng)來考察。因此,有必要將陳垣《二十史朔閏表》等看作一個數(shù)據(jù)庫、把它作為數(shù)字研究的對象,認為歷譜中的記錄包含了日月運行的周期性、干支紀日的周期性等多方面的信息。本文將應(yīng)用天文年歷中精確的回歸年和朔望月數(shù)值,并以出土文物(主要是簡牘中的紀年)和《史記》、《漢書》等史書中的有關(guān)記載作為根據(jù)。

    今記回歸年α=365.24219878,朔望月β=29. 53058867,設(shè)γ=α/β-12=0.36826678。問題是,若經(jīng)過S年使得每年所余γ月盡可能的湊足整月數(shù)t,即Sγ≈t,則稱t為閏月數(shù)(也就是閏年數(shù))。即在S年內(nèi)應(yīng)有t個閏年。γ≈t/S,可將γ寫成連分數(shù),再展開為漸近分數(shù)ti/Si。承上,利用精確的α,β值,所得結(jié)果為:

    將這些理論計算的數(shù)值與朔閏表的實際結(jié)果相對照,人們會吃驚地發(fā)現(xiàn),兩者幾乎完全吻合!翻開陳垣《二十史朔閏表》,任取相差1021年或1708年的兩年,其間的閏年數(shù)恰為376年或629年!如果不合,往往是由于數(shù)錯造成。對于上列漸近分數(shù)中的最后一個,問題在于2729年中是否真的有1005閏?須查張培瑜的《三千五百年歷日天象》[10],如果從公元前721年(魯隱公元年)起算,1995年是第一個“千閏年”,2009年將迎來第1005閏,也可以看作是第六個千閏年。這充分說明我國歷代的天文歷法學家,以他們堅持不懈的共同努力,使歷法的精度漸次密合于天象;千閏年的出現(xiàn),是趨同原理的典型事件,具有深刻的天文—歷法意義和豐富的數(shù)學內(nèi)涵。[11]

    當然,不同歷法如何置閏可以小有區(qū)別,但目的都是使所制定的歷法與太陽運行相調(diào),稱為調(diào)日法。同時,不同歷法又必須安排好大小月的順序和兩者間的數(shù)量比,以使歷法與月亮圓缺相諧,這里稱為諧月法。因β-29.5=0.03058867,說明大月的總數(shù)應(yīng)當比小月多,設(shè)一個階段(比如q年)中大月數(shù)減小月數(shù)的差為諧月數(shù)φ(整數(shù)),那么諧月法就是要解決在q年中應(yīng)加入多少個諧月的問題。

    設(shè)m=30日為1諧月(它是1大月,也可能是蟬聯(lián)大月的后一個月),n=29.5日為1平月(大小月均值),ε=β-n=0.03058867日為1平月不足1朔望月之日數(shù),δ=m-β=0.46941133日為1諧月多于1朔望月之日數(shù)。顯然ε+δ=m-n=0.5。設(shè)比值μ為每過qi年應(yīng)增添的pi個諧(大)月,可以證明(此略):

    表示成漸近分數(shù),即有

    μ的意義不太直觀,最后一個分數(shù)表示1089年中應(yīng)有824個諧月,從月朔表中我們同樣獲得了有力的證據(jù)(關(guān)鍵是按定義計算諧月數(shù)要準確)。應(yīng)用這列分數(shù),可求出任一時段q內(nèi)諧大月數(shù)p,q大則精確。以《三統(tǒng)歷》(前104到公元85)行用的189年為例。從上列分數(shù)中選出28/37和3/4,存在以下等式表明在5個37年中增添5次28個諧大月之后,所余4年應(yīng)加入3個諧月(此時與的相對誤差約十萬分之六)。因此p=143對q=189是最佳的選擇。而《三統(tǒng)歷》的實際情況p值偏小,我們統(tǒng)計(從公元前104年的正月起算)為142諧月。這樣,就可以從一個方面對歷法精度作出估計。

    類似地,可將此法應(yīng)用到歷法研究的更多方面,還可用以分析木星(歲星)、金、火、水、土五星的運行,對于深入認識天文歷法,不失為一種有效的工具。歷法中常見的周期性質(zhì)也可以用連分數(shù)、漸近分數(shù)的方法來證明。

    人們有理由關(guān)注歷譜中不同性質(zhì)的干支出現(xiàn)的周期性,這些問題有:1.任一朔干支重現(xiàn)周期,2.任一月朔干支重現(xiàn)周期,3.任一閏月朔重現(xiàn)周期,4.某一年朔序(定義見下)重現(xiàn)周期,5.任一八節(jié)(定義見下)干支重現(xiàn)周期,等。以往的研究已經(jīng)提到“農(nóng)歷的年、月也有周期性?!讶崭芍У闹芷诤娃r(nóng)歷年月聯(lián)系起來看,卻找不到一個整數(shù)的公共周期。我們從實踐中,找到了大體上五年和三十一年的兩個周期?!保?3]這個結(jié)論不夠清晰、準確,也需要進行論證,但它指出了周期性的存在。

    以下列出上述五個周期的性質(zhì)定理,有的給出簡要論證,有的列舉它在考古年代學中的應(yīng)用。

    四、干支歷特有的與月朔、閏月朔相關(guān)的三個周期性質(zhì)

    性質(zhì)一、任一朔干支會在61、63或65個月后重現(xiàn)。

    設(shè)A為某月朔干支,在經(jīng)歷了p個大月和q個小月后重現(xiàn)(p,q>0)。30p+29q一定能被60整除,即p,q的最小值應(yīng)當滿足同余式30p≡-29q(mod 60)因29為素數(shù),必有q=30;此時p為奇數(shù),且q≠2k·30(k=1,2,…),故只取q=30。因大月數(shù)p>q,故可取p=31,33,35,…由歷法知在連續(xù)60多個月中如果小月僅有30個,則大月不會有37個。否則誤差就太大了。故p=31,33,35。由此性質(zhì)一得證。

    31×30+30×29=1800日=30甲子周;33× 30+30×29=1860日=31甲子;35×30+30× 29=1920日=32甲子。

    例如:公元前105年12月甲午朔歷61個月后到前100年12月甲午朔共經(jīng)1800日;從1991年正月丙辰朔到1995年12月丙辰朔共歷61個月。常見的是中間經(jīng)過63或65個月。上文所說“大體上五年”的周期應(yīng)被性質(zhì)一取代,它是進一步推理的基礎(chǔ)。

    性質(zhì)二、任一同月同朔干支可能在26、31或36年后重現(xiàn)??珊営洖?1n或31n±5年(n=1,2,…)。

    設(shè)a年u月(這里暫不考慮閏月)w朔,經(jīng)s年后為b年u月w朔,求最小s之值。顯然s=b-a為正整數(shù)。如果s年中含d個月,則據(jù)性質(zhì)一,當x, y,z為非負整數(shù)時,應(yīng)有65x+63y+61z=d;另一方面,設(shè)s年中閏t個月,則共有的月數(shù)d=12s+t,同時滿足d=平月數(shù)+諧月數(shù)??蓪⒉欢ǚ匠毯喕?,利用歷法中的條件淘汰不合用的解(過程較繁,從略),得出三個d值:

    d1=321月=12×26+9閏=300平月+21諧月=9480日=158甲子周

    d2=384月=12×31+12閏=360平月+24諧月=11340日=189甲子周

    d3=445月=12×36+13閏=420平月+25諧月=13140日=219甲子周

    統(tǒng)觀幾千年歷譜s1=26年和s3=36年很少出現(xiàn)。例如漢宣帝神爵二年(前60)正月丁未朔,經(jīng)26年后到漢元帝建昭五年(前34)正月丁未朔;又如漢章帝建初五年(80)正月辛亥朔,經(jīng)36年到漢安帝元初三年(116)正月辛亥朔等。較多出現(xiàn)的是31n年(n=1,2,…)的周期,對31年=384月而言,實際上由一組解確定,由性質(zhì)一來看,由4次65個月、1次63個月和1次61個月堆疊而成;這31年同時滿足閏12次、包含24個諧大月,用前文所述兩列漸近分數(shù)可解釋為:

    性質(zhì)二對于考證年代十分有用。在許多出土文物例如簡牘中,大部分未記年號,有一類僅記有月份和該月朔干支(或其月朔干支可以推算出來),叫做“月朔簡”[14],例如居延漢簡“六月癸巳朔”(148.37號),“五月廿八日戊戌”等,它們的年代大部分可以準確推出[15]。這是因為對居延發(fā)掘地各探方伴出紀年簡上下時限進行統(tǒng)計,其平均值約為18年;而性質(zhì)二能夠保證在31年(2,36年)內(nèi)同月同朔紀錄只有一個,這對考古工作實在太有用了,結(jié)合出土信息和簡上內(nèi)容。所定年代能經(jīng)得起多種考驗。在敦煌發(fā)掘地的情況也類似。事實上性質(zhì)二對干支歷都是適用的,因此對各歷史時期的文物年代鑒定能發(fā)揮它的作用。

    性質(zhì)三、閏月同朔干支經(jīng)常在51、52或57年后重現(xiàn)。

    設(shè)A年閏u月w朔,經(jīng)s年后為B年閏v月w朔,求常見的最小s之值。如果s年中含d個月,則據(jù)性質(zhì)一,當x,y,z為非負整數(shù)時,應(yīng)有65x+ 63y+61z=d另外,設(shè)s年中閏t個月,則共有的月數(shù)d=12s+t-u+v,同時滿足d=平月數(shù)+諧月數(shù)。

    這些都是必要條件。閏月朔的條件較復雜,用演繹的方法去求s會遇到很大困難,所以可用歸納的方法,對朔閏表中作為閏月朔的60干支分別進行統(tǒng)計。我們選取從公元前206年(西漢高祖元年)到公元815年(唐憲宗十年)共1021年,其間恰有376閏(與前列漸近分數(shù)376/1021相合,說明滿足趨同原理)。性質(zhì)三是統(tǒng)計的結(jié)果,它滿足已列出的必要條件,討論如下。

    d1=638月,表現(xiàn)為51或52年。在不定方程中x=7,y=0,z=3。51或52年中必有19閏,因

    故對51年,d1=638=51×12+19-u+v,即vu=7。說明凡過51年出現(xiàn)同朔之閏月,相差月數(shù)為7。對52年,d2=638=52×12+19-u+v,即v-u=5。說明凡過52年出現(xiàn)同朔之閏月,相差月數(shù)為5。

    在51或52年中應(yīng)有38個諧月,這時d1=638月=600平月+38諧月=18840日=314甲子周。

    另一結(jié)果是d2=705月=57年。在不定方程中x=8,y=1,z=2。因19×3=57故必有7×3=21閏,有

    d2=57×12+21-u+v=705,即v=u。說明凡過57年出現(xiàn)同朔之閏月,其閏月數(shù)應(yīng)相等。因此這是同月同朔干支的周期問題,它同時滿足性質(zhì)一和性質(zhì)二。事實上57年=26年+31年=321月+384月 =705月,是將兩者堆疊而成的。

    在57年中應(yīng)有45諧月,這時d2=705月=660平月+45諧月=20820日=347甲子周。

    以下舉出3例:

    晉武帝末年太熙元年(290)閏五月己丑朔,經(jīng)51年至晉成帝咸康七年(341)閏十二月己丑朔,歷638月,兩閏月份相差為7個月;

    唐玄宗天寶十三年(754)閏十一月壬戌朔,經(jīng)52年至唐憲宗元和元年(806)閏六月壬戌朔,歷638月,兩閏月份相差為5個月;

    漢宣帝地節(jié)元年(前69)閏正月己巳朔,經(jīng)57年至漢成帝延元年(前12)閏正月己巳朔,歷705月,兩閏月月份相同。

    當然,對于歷法系統(tǒng)而言,由于幾千年中存在的非歷法因素的影響,使得有些性質(zhì)會有例外。例如王莽建新朝(公元8),下令將當年十二月改為正月,第二年正月改為二月,等,直到莽亡(公元23)。所以在此期間的閏月數(shù)要比正常的多1。如莽歷閏五月實際是閏四月。

    在所考察的1021年中包括性質(zhì)三在內(nèi)全部閏月同朔干支相距的年代d均滿足57x+52y+51z=d,其中x,y,z為整數(shù),意即d是57,52和51的線性組合,我們已找出考察區(qū)間內(nèi)的全部d值,并求出每個d的(x,y,z)值。這里省略。

    例如前文所提到的漢宣帝神爵三年(前59)閏十二月丙申朔,要經(jīng)過d=533年后,在元徽三年(475)閏三月才出現(xiàn)第二次丙申朔,整個102年中閏月丙申朔才有三次,是閏十二月的只有一次。所以在2200多年中“閏十二月丙申朔”具有唯一性是不奇怪的了。d=533可分解為x=3,y=5,z=2。這種情況并非個別,例如漢成帝永始二年(前15)閏五月甲寅朔,要經(jīng)過d=635年后,在唐高祖武德四年(621)閏十月才出現(xiàn)第二次甲寅朔,整個1021年中閏月甲寅朔共有兩次??梢娭袊鴼v法朔閏干支紀日法具有相當?shù)奶厥庑?,以前沒有被揭示過。這里d=635年可分解為x=3,y=5,z=4。

    比57,52和51年小的d值有46年和47年,46年可分解為x=-1,y=1,z=1;47年為x=-1,y=2,z=0。它們出現(xiàn)的次數(shù)僅占全部d值出現(xiàn)總數(shù)的2.2%。

    對公元前后的年數(shù)相加時,因無公元零年,故須將公元前年數(shù)絕對值減1后再與公元年數(shù)相加。

    置閏在制歷中占有重要地位。太初元年(前104)之前,漢承秦制,以十月為一年之始,而將所有閏月排在歲尾,稱之謂“后九月”。以后歷法進步,按“無中(氣)置閏法”(不詳論),要據(jù)當時觀測、計算的結(jié)果來確定哪個月為閏月,后人所推出的歷譜就出現(xiàn)了不同的結(jié)果,如陳垣《二十史朔閏表》與劉曦叟的《長歷》就有諸多閏月記載不合;另外出土文物的記錄又出現(xiàn)于歷譜不合的情況,成為年代學研究中較突出的一個問題。所以有必要對閏月朔干支的周期性花費較多篇幅進行研究。上述必要條件還不是充分的,要作出置閏確否的判斷需要用綜合的方法,特別是要注重出土文物和史書的記載。

    例如諸歷表均記神爵元年(前61)閏四月壬午(19)朔。但以下證據(jù)證明此年應(yīng)閏三月壬子49朔:

    ①“神爵元年四月壬午朔己酉46佐解事敢言□”(居延漢簡109.7號)。按照當時記錄閏月的習慣,只寫“閏月”,不寫“閏幾月”。故此簡壬午決非閏月朔,它表明應(yīng)是閏三月壬子朔。

    ②陳夢家先生說:“居延漢簡10.27,5.10,332.26,10.33,10.30,10.32,10.29,10.31等八簡是相聯(lián)的一冊,記神爵元年二月至五月詔書之下達,其閏月丁巳54、庚申57皆屬于閏三月壬子49朔”。[18]

    ③“四月廿九日庚戌47寢兵五月大辛亥48一日壬子49夏至”(下略,居延漢簡179.10號)。此簡可據(jù)“五月大辛亥一日(朔)”按性質(zhì)二參照伴出紀年簡的上下時限考為神爵元年,也可據(jù)“壬子夏至”周期性得出相同的結(jié)論。它表明非閏四月。

    ④《漢書·王子候表第三》載:“新利候偃膠東戴王子神爵元年四月癸巳30封”(《漢書》2冊493頁),表明該年非閏四月。

    ⑤應(yīng)用數(shù)字方法,從集合論的角度看,用“相連兩歷法月同屬奇(或偶)集其間必有閏月”的判定法則(不詳論),可確知是年閏三月。綜上,前61年當閏三而非閏四,諸表皆誤。

    作此修改之后,經(jīng)過108年,到東漢光武帝廿四年(公元48)閏十月才出現(xiàn)壬子朔,d=108分解為x=1,y=0,z=1,說明這樣的修改滿足上述必要條件,從一方面證明它的可信度。

    史學中立論與駁論有嚴格要求,證“無三不立,無五不駁”,要駁倒成論須有大量可靠證據(jù)。

    五、干支歷特有的年朔序與分至八節(jié)的周期性質(zhì)

    性質(zhì)四、如果某兩年的年朔序相同,則這兩年相隔31n(n=3,4,6,7,10,13,14,17,20,21,24,27,…)年。其中n=4,7,10,14,較多,n=7最常見。

    所謂年朔序,指一年中朔日干支序列。例如據(jù)陳垣表漢高祖元年(前206)的年朔序為:乙卯乙酉甲寅甲申甲寅癸未癸丑壬午壬子辛巳辛亥庚辰,依次是前206年正月到十二月(此處暫不考慮十月歲首的規(guī)定)的朔日干支,將它譯為字符串,則成52 22 51 21 51 20 50 19 49 18 48 17。將全部歷譜的年朔序都譯為這樣的數(shù)串,要問哪一年將和它重復?

    我們回答:過961年(=31×31,這里n=31)后,在唐玄宗末年天寶十五年(756)的年朔序和它完全相同。初看起來令人難以置信:時隔近千年,其間出現(xiàn)各種歷法30余種,而且756年是著名天文學家一行改用《大衍歷》的末年。

    其實年朔序周期出現(xiàn)是較普遍的現(xiàn)象,它一直沒有引起重視。每次改歷,都要變更朔序,顯然對周期性有影響;但這種性質(zhì)突破各歷的差別,總會頑強地表現(xiàn)出來。不以大時段歷法的共同規(guī)律為研究對象則不能對它作出合理的解釋。本文是首次論述這些周期性質(zhì)。

    根據(jù)性質(zhì)二,同月同朔干支可能會在31n或31n5年后(n=1,2,…)重現(xiàn)。年朔序相同就是要求一年各月朔干支皆同,條件更嚴,31n年成了必要但非充分的條件,進而須對n作限定。其周期性十分復雜,可用對歷譜進行統(tǒng)計分析的方法。我們選取從公元前206年到公元805年的1010年作為考察區(qū)間,從三方面進行研究:①找出所有重復的年朔序;②對不同的再現(xiàn)周期進行分類統(tǒng)計;③建立一個模型解釋統(tǒng)計現(xiàn)象。

    結(jié)論是:在1010年中年朔序重復者有264例,涉及36年,占35.74%;同時,共有649個年號在千年之內(nèi)保持年朔序的獨特性、唯一性,占64.26%。將年朔序相同的公元年數(shù)寫在括號內(nèi),將重現(xiàn)周期(年數(shù))寫在括號前,按n值從小到大排到,每個n值只舉1例,有

    這個表并沒有寫完;還有個別特例。限于篇幅,不贅。

    把這些n值不同的年朔序作分類統(tǒng)計,可知它們出現(xiàn)的頻數(shù)各不相同,最常見的是n=7即相隔217年,可舉出98例,占264例的37.1%,其他依次為n=4(占12.1%),n=17(占9.85%),n=10(占9.47%)和n=14(占8.71%)。如果x,y為整數(shù),則值一般滿足3x+4y=n,|x-y|≤2。

    事實上除n為幾個常見值外,其他值都是據(jù)此推算出來。例如對x=5,y=3,有n=27;帶著“可能在837年(=31×27)后年朔序重現(xiàn)”的正確信息,才能發(fā)現(xiàn)從西漢惠帝二年(前193)過837年后到唐太宗貞觀十九年(645)年朔序重復出現(xiàn)。其他如當x=1,y=2時n=11;x=2,y=3時n=18等因出現(xiàn)頻數(shù)較小,或有其他個別例外。

    中國歷史上有多少不同的年朔序?這個問題還是未知的。當修改了某年某月朔干支,可能出現(xiàn)新的年朔序,但是獨特的年朔序占60%多,重現(xiàn)的年朔序占30%多,大體是不會改變的。了解這一基本知識可以使我們避免年代考證中的錯誤。包括一些著名的國學大師如羅振玉(1866-1940)、王國維(1877-1927)在內(nèi),因為忽略了年朔序的周期性,對敦煌漢簡中兩枚歷譜簡年代的確定,均推遲了217年(=31×7),這在年代學中當然是不能允許發(fā)生的。

    第1枚簡:疏勒298號[7](斯坦因編號)釋文為:“三日戊寅戊申戊寅丁未丁丑丙午丙子乙巳乙亥甲辰甲戌癸卯”。法國漢學家沙畹(E. Chavannes,1865-1918)和羅振玉[19]均考為永興元年(153),但實際上應(yīng)是元康元年(前65);這兩年的年朔序是相同的,相差217年。

    第2枚簡:疏勒697號[7](斯坦因編號)釋文為:“九月辛丑朔大(中略)辛巳立冬十月十一日乙未”,王國維、張鳳認為它是“建安十年(公元205)歷”,但應(yīng)當是永始四年(前13)。[20]

    在我們看來,當n=7時31×7(-65,153),(-13,205)只是周期為217年的兩例。其實從298號“三日”簡不難向前推兩天而求出“一日”簡的干支序列,即年朔序。它當然與公元153年朔序相同,但這僅是判斷年代的必要條件,其充分性并未證明,羅振玉在數(shù)學上犯了一個簡單的基本錯誤。我們學習高等數(shù)學,重要目的便是提交數(shù)學素養(yǎng)和邏輯思維能力。

    性質(zhì)五、分至八節(jié)干支的任一個都會在23年,80年或103年重現(xiàn)。

    中國傳統(tǒng)歷法是陰陽合歷。“分至”指春分、秋分、夏至、冬至,再加上立春、立夏、立秋、立冬,謂之“八節(jié)”,這是“陽歷”即太陽歷的內(nèi)容。重現(xiàn)24節(jié)氣,努力使它平均分布在一個回歸年(日)中,即節(jié)氣的確定應(yīng)符合太陽的運動,這樣耕種收獲就不誤農(nóng)時。

    《分至八節(jié)表》[10](前1500—公元2050年)載有3550年的2.84萬個干支數(shù)據(jù),平均每節(jié)的每個干支出現(xiàn)59.167次。由于規(guī)律性很明顯,不必去研究480種不同八節(jié)干支的全部,只須追蹤其中的幾個,就可以歸納出性質(zhì)五來;也不必列出所有(平均)59個年代,只須選擇其中一段,例如:

    壬戌春分:-190(23)-167(80)-87(23)-64(80)17(23)40(80)120(23)143(80)223…

    甲子秋分:-145(103)-42(103)62(103)62(103)165(103)268(103)371(23)394(80)474(23)497(80)577…

    癸巳冬至:-145(23)-122(80)-42(23)-19(80)62(23)85(80)165(23)188(80)268…

    其中負號表示公元前,括號內(nèi)的數(shù)字表示括號前后兩個年代的差。其規(guī)律性一目了然。

    很自然要問:為什么會出現(xiàn)性質(zhì)五的周期呢?23α=8400.5706日=140甲子周+0.5706日,8α=29219.3759日=487甲子周-0.6241日,103α=37619.9465日=627甲子周-0.0535日。

    在103年中只有1小時17分4.6秒的誤差,歷法當然要盡可能與之相合。不同歷法、出土文物等對某年八節(jié)干支的記錄也可能出現(xiàn)差別,但這不會影響性質(zhì)五的周期性。偶爾有57年的周期,因57α是347甲子周少1.1947日,誤差稍大,這就是它在歷法中較少出現(xiàn)的原因。

    事實上推算出上述幾萬個八節(jié)干支并不難,難的是一旦與文物、史書記載相忤時該如何解釋。由于規(guī)律性強,可將歷史年代劃分成若干時段(如秦漢440年可劃成三段),選定正確的初值(如前221年的八節(jié)干支),代入統(tǒng)一的八節(jié)干支算法公式,注意每段的修正值,則可推算出其間任一年的八節(jié)干支來。

    將性質(zhì)五應(yīng)用于簡牘等文物的年代考釋,以上述居延漢簡179.10號“壬子夏至”為例,據(jù)《分至八節(jié)表》:-164(103)-61(103)43(23)66(80)146…,只有公元前61年(神爵元年)符合此簡的年代。據(jù)《史記·匈奴傳》:太初三年(前102)“使強弩都尉路博德筑居延澤上”;[21]而建武八年(公元32)之后該地屯戍活動急劇減少,只出土了零星分散的紀年簡,所以神爵元年前后各103年皆不可能成為此簡的年代。

    又如上述敦煌漢簡697號“辛巳立冬”,有:-196(80)-116(23)-93(80)-13(23)11(80)91(23)114…,其中只有前13年(永始四年)符合此簡的年代,較近的公元11年(始建國三年)是王莽時代,與697號其他信息不符。

    六、結(jié)語

    中國傳統(tǒng)歷法用干支紀年月日時,形成幾千年的歷譜,包含了每個歷法階段的豐富信息,有的屬于天文、觀測,有的屬于數(shù)學、計算,有的屬于社會、人事,成為人類珍貴的科學遺產(chǎn);同時它又是歷史研究中必不可少的時間坐標,熟練而正確使用歷表是最基本的要求。因此要走進這座寶庫需要天文、歷法、數(shù)學的一些知識。

    顧炎武說:“三代之上人人皆知天文”?,F(xiàn)代分工和教學的內(nèi)容使大多數(shù)人缺乏起碼的天文知識,干支歷對青年人也日漸陌生。然而,中國歷法已有幾千年歷史,凡有華人的地方,年復一年隆重慶祝傳統(tǒng)的春節(jié)、元宵、端午、中秋等節(jié)日,將世代相傳。

    每當提到天文、歷法、數(shù)學時,學文科的學生不免立即聯(lián)想到非常復雜的公式、微積分方法等,從而產(chǎn)生畏懼心理。當然現(xiàn)代科學的確有許多領(lǐng)域不受專門訓練便不能得其門徑而入;但是古人也同樣不懂這些,他們卻靠堅持不懈觀測天象、不斷改進算法和進行主要是四則運算而制定了精密的歷法,沿用幾千年,這是發(fā)人深省的。

    本文所用的方法屬于初等數(shù)論的范疇,所提出的漸近分數(shù)和性質(zhì)定理具有實際應(yīng)用的價值。歷法研究一般是根據(jù)史書天文律歷志對天文常數(shù)、觀測數(shù)據(jù)和算法的記載去推算歷法,起止時間以某部歷法行用為始終,有一定局限性。本文則力圖以歷法系統(tǒng)為目標,去尋求適合于各種歷法的共同規(guī)律;[22]由于傳統(tǒng)歷法不是一個完整的演繹系統(tǒng),所以將歸納和演繹法相結(jié)合,研究統(tǒng)計的規(guī)律,推敲建立數(shù)學的模型,從而獲得一批周期性的結(jié)果。我們已據(jù)此研究了幾百枝簡牘的年代學問題,在史學和漢學的專業(yè)刊物如《文史》、《敦煌研究》、臺灣《漢學研究》上發(fā)表,將此性質(zhì)應(yīng)用于年代考釋初見成效。在人文社會科學中應(yīng)用數(shù)學方法所見多有,體現(xiàn)了數(shù)學向文科領(lǐng)域滲透的趨勢;而應(yīng)用初等數(shù)論研究歷史年代學,也將逐步完善,成為一種可以接受的輔助考釋方法。

    [1]謝桂華,等編.居延漢簡釋文合校上冊[M].北京:文物出版社,1987:228.

    [2]謝桂華,等編.居延漢簡釋文合校下冊[M].北京:文物出版社,1987:594.

    [3]陳垣.二十史朔閏表[M].北京:中華書局,1978.

    [4]汪曰楨.歷代長術(shù)輯要[M]//四部備要·子部.北京:中華書局,1936.

    [5]班固.漢書第3冊[M].北京:中華書局,1975:848.

    [6]北京師范大學國學研究所.武王克商之年研究[M].北京:北京師范大學出版社,1997.

    [7]林梅村,李均明.疏勒河流域出土漢簡[M].北京:文物出版社,1984.

    [8]甘肅省文物考古研究所等編.居延新簡——甲渠候官[M].北京:中華書局,1994.

    [9]華羅庚.高等數(shù)學引論.一卷一分冊[M].北京:科學出版社,1963:63.

    [10]張培瑜.三千五百年歷日天象[M].鄭州:河南教育出版社,1990.

    [11]羅見今.中國歷法中的“千閏年”[M]//黃留珠編.周秦漢唐文化研究.第四輯.西安:三秦出版社,2006:1-5.

    [12]呂子方.中國科學技術(shù)史論文集.上冊[M].成都:四川人民出版社,1983:35-57.

    [13]葛真.用日食、月相來研究西周的年代學[M].貴州工學院學報,1980(2):81-100.

    [14]羅見今.《居延新簡——甲渠候官》中月朔簡年代考釋[J].中國科技史料,1997(3):72-83.

    [15]羅見今.《居延漢簡》中月朔簡年代考釋[J].臺北:漢學研究,17,1999(2):169-185.

    [16]羅見今,關(guān)守義.敦煌漢簡中月朔簡年代考釋[J].敦煌研究,1998(1):158-165.

    [17]羅見今,關(guān)守義.敦煌、居延漢簡中與朔閏表不合諸簡考釋[M].文史,2000(1):57-72.

    [18]陳夢家.漢簡綴述[M].北京:中華書局,1980:233.

    [19]羅振玉,王國維.流沙墜簡[M].京都:京都東山學社, 1914:9.

    [20]羅見今,關(guān)守義.敦煌、居延若干歷簡年代考釋與質(zhì)疑[J].臺北:漢學研究,15,1997(2):37-50.

    [21]司馬遷.史記·匈奴列傳[M].北京:中華書局,1959:第9冊2916.

    [22]羅見今.中國歷法的五個周期性質(zhì)及其在考古年代學中的應(yīng)用[M]//黃留珠主編.周秦漢唐文化研究:第三輯.西安:三秦出版社,2004:6-18.

    (責任編輯:戢斗勇jidouyong@qq.com)

    The Application of Elementary Number Theory in the Historical Chronology

    LUOJian-jin
    (Institute for the HistoryofScience&amp;Technology,Inner Mongolia Normal University,Huhhot 010022,China)

    Based on some tools in elementary number theory,such as congruence,asymptotic fraction,the continued fraction Diophantine equation,etc.,this paper analyzes some questions presented in archaeological chronology.Taking Dunhuang bamboo slips,Ju Yan bamboo slips of Han Dynasty,the book“Historical Records”,unearthed relics and historical data as evidence,and the records of“Calendar of Twenty Dynasties in China”and“Calendar in Three Thousand Five Hundred Years”as a database of sun and moon movements, this paper proves and discusses five cycles properties in Chinese calendar so as to clarify the basic principles and methods of the archaeological chronology of mathematics.

    Elementary number theory;bamboo slips in Han Dynasty;Chinese Historical Calendar;the cyclical theorems of Chinese calendar;archaeological chronology

    K04

    A

    1008-018X(2015)06-0028-10

    2015-09-20

    羅見今(1942-),男,河南新野人,內(nèi)蒙古師范大學教授。

    猜你喜歡
    年代學干支甲子
    明 文震孟 行書甲子除夕詩卷
    中國書法(2023年7期)2023-09-14 09:59:45
    福安土白語匯的語言年代學考察
    干支梅
    草原歌聲(2020年3期)2021-01-18 06:51:58
    古代的“干支紀年法”
    鮮水河斷裂帶中南段同位素年代學及其地質(zhì)意義
    中國風
    北大漢簡《節(jié)》篇“冬夏至干支速算表”解讀
    海岱地區(qū)年代學新成果的思考與啟示
    東方考古(2018年0期)2018-08-28 10:05:04
    沈鵬《迎戊戌雙甲子》
    中華詩詞(2018年3期)2018-08-01 06:40:42
    《夢溪筆談》“甲子納音”構(gòu)造方法的數(shù)學分析
    ponron亚洲| 国产在线精品亚洲第一网站| 国产精品一区www在线观看| 男人和女人高潮做爰伦理| 久久久久国内视频| 99国产精品一区二区蜜桃av| 一级毛片电影观看 | 俄罗斯特黄特色一大片| 99久久精品热视频| 一进一出好大好爽视频| 久久精品夜夜夜夜夜久久蜜豆| 午夜影院日韩av| 精华霜和精华液先用哪个| 黄片wwwwww| 五月玫瑰六月丁香| 国产一区二区三区在线臀色熟女| 我要看日韩黄色一级片| 中文字幕熟女人妻在线| 美女xxoo啪啪120秒动态图| 成人午夜高清在线视频| 一本一本综合久久| 美女 人体艺术 gogo| 久久精品国产亚洲网站| videossex国产| 性插视频无遮挡在线免费观看| 毛片女人毛片| 日日摸夜夜添夜夜爱| 日本爱情动作片www.在线观看 | 国产单亲对白刺激| av黄色大香蕉| 18+在线观看网站| 搞女人的毛片| 成人特级黄色片久久久久久久| av国产免费在线观看| 天堂√8在线中文| 高清毛片免费观看视频网站| 欧美另类亚洲清纯唯美| 在线国产一区二区在线| 亚洲经典国产精华液单| 久久人人精品亚洲av| 波多野结衣高清作品| 亚洲国产精品成人综合色| 非洲黑人性xxxx精品又粗又长| av在线老鸭窝| 成人精品一区二区免费| 久久久久免费精品人妻一区二区| 婷婷亚洲欧美| 日本成人三级电影网站| 永久网站在线| 亚洲精品一区av在线观看| 成人国产麻豆网| 国产欧美日韩精品亚洲av| 激情 狠狠 欧美| 免费搜索国产男女视频| 日本欧美国产在线视频| 干丝袜人妻中文字幕| 国产aⅴ精品一区二区三区波| 日韩成人伦理影院| 成人特级黄色片久久久久久久| 久久精品国产清高在天天线| 精品久久久久久久末码| 亚洲在线自拍视频| 亚洲中文字幕一区二区三区有码在线看| 久久久久久久久大av| 亚洲精品久久国产高清桃花| 一级毛片aaaaaa免费看小| 国产女主播在线喷水免费视频网站 | 欧美性猛交黑人性爽| 一进一出抽搐动态| 搞女人的毛片| 欧美bdsm另类| 久久国内精品自在自线图片| 亚洲中文日韩欧美视频| 精华霜和精华液先用哪个| 精品免费久久久久久久清纯| 中文在线观看免费www的网站| av天堂中文字幕网| 日韩三级伦理在线观看| 美女 人体艺术 gogo| 在线播放无遮挡| av在线观看视频网站免费| 久久综合国产亚洲精品| 久久精品国产鲁丝片午夜精品| 欧美最新免费一区二区三区| 天堂动漫精品| 少妇熟女欧美另类| 黄色一级大片看看| 日韩制服骚丝袜av| 亚洲欧美清纯卡通| 小蜜桃在线观看免费完整版高清| 久久久午夜欧美精品| 听说在线观看完整版免费高清| 国产69精品久久久久777片| 欧美日本视频| 一区福利在线观看| 亚洲av免费高清在线观看| 啦啦啦韩国在线观看视频| 一进一出抽搐gif免费好疼| 国产一区二区在线av高清观看| 成人性生交大片免费视频hd| 欧美绝顶高潮抽搐喷水| 国产精品一及| 最近在线观看免费完整版| 99精品在免费线老司机午夜| 国产 一区精品| av视频在线观看入口| 精品99又大又爽又粗少妇毛片| 免费无遮挡裸体视频| 九九久久精品国产亚洲av麻豆| 一区二区三区免费毛片| 可以在线观看毛片的网站| 22中文网久久字幕| 天美传媒精品一区二区| 亚洲国产精品国产精品| 日韩一本色道免费dvd| 中文字幕熟女人妻在线| 亚洲人成网站在线播| 人妻久久中文字幕网| 亚洲国产精品合色在线| 国产男人的电影天堂91| 男人狂女人下面高潮的视频| 国产高清有码在线观看视频| 熟女电影av网| 国产视频一区二区在线看| 在线观看免费视频日本深夜| 12—13女人毛片做爰片一| 婷婷色综合大香蕉| 搡老岳熟女国产| 色哟哟·www| 亚洲精品久久国产高清桃花| 亚洲人成网站在线播| 欧洲精品卡2卡3卡4卡5卡区| 人人妻人人看人人澡| 在线天堂最新版资源| 国产精品无大码| 香蕉av资源在线| 亚洲高清免费不卡视频| 日本黄色片子视频| 国产在线男女| 亚洲中文字幕日韩| 成人三级黄色视频| 国产欧美日韩精品亚洲av| 一夜夜www| 日本欧美国产在线视频| 亚洲图色成人| 精品国内亚洲2022精品成人| 伦理电影大哥的女人| 免费在线观看影片大全网站| 成人国产麻豆网| a级毛片a级免费在线| 春色校园在线视频观看| 亚洲国产精品国产精品| 校园春色视频在线观看| 精品久久久久久久久久免费视频| 99视频精品全部免费 在线| 婷婷色综合大香蕉| 十八禁网站免费在线| av卡一久久| 午夜福利成人在线免费观看| 国产 一区 欧美 日韩| 成人综合一区亚洲| 色5月婷婷丁香| 亚洲欧美日韩高清在线视频| av专区在线播放| 色综合色国产| 日产精品乱码卡一卡2卡三| 在线观看美女被高潮喷水网站| 国产淫片久久久久久久久| 中文在线观看免费www的网站| 精品福利观看| 一级a爱片免费观看的视频| 最近视频中文字幕2019在线8| 国产精品久久久久久av不卡| 日韩高清综合在线| 99热这里只有是精品在线观看| 午夜精品一区二区三区免费看| 亚洲,欧美,日韩| 亚洲国产高清在线一区二区三| 日韩欧美精品v在线| 欧美又色又爽又黄视频| 色综合亚洲欧美另类图片| 亚洲无线观看免费| 美女xxoo啪啪120秒动态图| 少妇裸体淫交视频免费看高清| 欧美激情在线99| 99热这里只有精品一区| 内射极品少妇av片p| 亚洲av五月六月丁香网| 成人国产麻豆网| 俺也久久电影网| 一级毛片aaaaaa免费看小| 国产女主播在线喷水免费视频网站 | 一本精品99久久精品77| 欧美国产日韩亚洲一区| 一本精品99久久精品77| 国产色婷婷99| 精品久久久噜噜| 最近的中文字幕免费完整| 成人鲁丝片一二三区免费| 美女高潮的动态| 亚洲欧美精品综合久久99| 久久精品国产99精品国产亚洲性色| 成人一区二区视频在线观看| 国产精品野战在线观看| 人妻夜夜爽99麻豆av| 免费不卡的大黄色大毛片视频在线观看 | 精品久久久久久久久av| 国产一区亚洲一区在线观看| 日韩av在线大香蕉| 在线免费十八禁| 亚洲国产精品国产精品| 美女被艹到高潮喷水动态| 欧美+日韩+精品| 免费在线观看影片大全网站| 精品久久久久久久久亚洲| 久久人人爽人人爽人人片va| 亚洲精品影视一区二区三区av| 精品久久久久久成人av| 美女高潮的动态| 色视频www国产| 九九爱精品视频在线观看| av在线亚洲专区| 日韩欧美 国产精品| 国产午夜福利久久久久久| 日韩中字成人| 国产精品久久久久久av不卡| 亚洲精品乱码久久久v下载方式| 在线天堂最新版资源| 九色成人免费人妻av| 桃色一区二区三区在线观看| 白带黄色成豆腐渣| 国产视频一区二区在线看| 精品久久久久久久久亚洲| av视频在线观看入口| 亚洲最大成人手机在线| 51国产日韩欧美| 国产人妻一区二区三区在| 99久久成人亚洲精品观看| 真实男女啪啪啪动态图| 少妇熟女aⅴ在线视频| 国产精品一区二区三区四区免费观看 | 亚洲自偷自拍三级| 禁无遮挡网站| 内地一区二区视频在线| 日韩,欧美,国产一区二区三区 | 午夜视频国产福利| 亚洲第一区二区三区不卡| 香蕉av资源在线| 99久国产av精品国产电影| 老女人水多毛片| 日韩,欧美,国产一区二区三区 | 亚洲精华国产精华液的使用体验 | 精品国产三级普通话版| 亚洲av一区综合| 国产片特级美女逼逼视频| 国产中年淑女户外野战色| 欧美日本视频| 久久欧美精品欧美久久欧美| 欧美一区二区亚洲| 日本撒尿小便嘘嘘汇集6| 国产午夜精品论理片| 国产真实伦视频高清在线观看| 我的女老师完整版在线观看| 亚洲欧美日韩卡通动漫| 99热全是精品| 成人美女网站在线观看视频| 国产高清视频在线播放一区| 精品久久国产蜜桃| 秋霞在线观看毛片| 亚洲精品乱码久久久v下载方式| 国产精品国产高清国产av| 亚洲精品成人久久久久久| 长腿黑丝高跟| 国产精品美女特级片免费视频播放器| 色视频www国产| 久久精品国产亚洲av天美| 成人亚洲欧美一区二区av| av天堂在线播放| 国产精品三级大全| 成年免费大片在线观看| 蜜桃亚洲精品一区二区三区| 最近最新中文字幕大全电影3| av.在线天堂| 18禁黄网站禁片免费观看直播| 精品午夜福利在线看| 麻豆成人午夜福利视频| 国产成人a∨麻豆精品| 欧美成人a在线观看| 国产黄a三级三级三级人| 国产一区二区亚洲精品在线观看| 午夜福利视频1000在线观看| 久久草成人影院| 老熟妇乱子伦视频在线观看| 国产精品一区二区三区四区久久| 精品一区二区三区视频在线| 亚洲不卡免费看| 在线天堂最新版资源| 国产精品一区二区三区四区久久| 久99久视频精品免费| 免费人成在线观看视频色| 人人妻,人人澡人人爽秒播| 91在线精品国自产拍蜜月| 在线播放国产精品三级| 国产高潮美女av| 狂野欧美激情性xxxx在线观看| 国产 一区 欧美 日韩| 日本 av在线| 99久国产av精品| 国产精品久久久久久亚洲av鲁大| 在线免费观看不下载黄p国产| 人妻夜夜爽99麻豆av| 插阴视频在线观看视频| 高清午夜精品一区二区三区 | 欧美另类亚洲清纯唯美| 亚洲成人中文字幕在线播放| 亚洲欧美精品综合久久99| av中文乱码字幕在线| 午夜爱爱视频在线播放| 两个人的视频大全免费| 一本一本综合久久| 1000部很黄的大片| 日韩亚洲欧美综合| 国产高潮美女av| 国内精品宾馆在线| 三级国产精品欧美在线观看| 国内精品宾馆在线| 在线天堂最新版资源| 99九九线精品视频在线观看视频| 亚洲国产精品国产精品| 日韩精品青青久久久久久| videossex国产| 看片在线看免费视频| ponron亚洲| 亚洲高清免费不卡视频| 免费观看精品视频网站| 国产精品综合久久久久久久免费| 日韩欧美免费精品| 国产久久久一区二区三区| av专区在线播放| 国产男靠女视频免费网站| 欧美色视频一区免费| 成年版毛片免费区| 少妇的逼好多水| 成人美女网站在线观看视频| 久久人人爽人人爽人人片va| 一级毛片aaaaaa免费看小| 亚洲av成人av| 99久久精品一区二区三区| 波多野结衣高清作品| 老师上课跳d突然被开到最大视频| 国内揄拍国产精品人妻在线| 99热这里只有精品一区| 日韩制服骚丝袜av| 国产 一区 欧美 日韩| 亚洲av.av天堂| 日本三级黄在线观看| 亚洲成人中文字幕在线播放| 亚洲av电影不卡..在线观看| 日韩欧美精品v在线| 欧美成人一区二区免费高清观看| 可以在线观看毛片的网站| 看非洲黑人一级黄片| 国产人妻一区二区三区在| 亚洲av五月六月丁香网| 国产 一区 欧美 日韩| 不卡视频在线观看欧美| 色在线成人网| 欧美xxxx黑人xx丫x性爽| 国产在线男女| 丰满的人妻完整版| 精品无人区乱码1区二区| 插逼视频在线观看| a级毛片a级免费在线| 在线播放国产精品三级| 丰满人妻一区二区三区视频av| 久久精品国产鲁丝片午夜精品| 99久久精品热视频| 久久久a久久爽久久v久久| 91麻豆精品激情在线观看国产| 中文字幕人妻熟人妻熟丝袜美| 精品一区二区免费观看| 国产成人a区在线观看| 深夜精品福利| 啦啦啦韩国在线观看视频| 精品久久久久久久久av| 欧美最黄视频在线播放免费| 色噜噜av男人的天堂激情| 一个人观看的视频www高清免费观看| 俄罗斯特黄特色一大片| 可以在线观看的亚洲视频| 国产精品99久久久久久久久| 黑人高潮一二区| 夜夜爽天天搞| 乱码一卡2卡4卡精品| 欧美最新免费一区二区三区| 午夜a级毛片| 国模一区二区三区四区视频| 99热精品在线国产| 欧美三级亚洲精品| 91麻豆精品激情在线观看国产| 99久久无色码亚洲精品果冻| 免费看日本二区| 亚洲精品亚洲一区二区| 亚洲婷婷狠狠爱综合网| 国产色婷婷99| 2021天堂中文幕一二区在线观| 天天躁日日操中文字幕| 日韩欧美在线乱码| 免费不卡的大黄色大毛片视频在线观看 | 久久鲁丝午夜福利片| 女人十人毛片免费观看3o分钟| 国产精品一区二区三区四区免费观看 | 久久久久久久久中文| 国产精品爽爽va在线观看网站| 日韩成人伦理影院| 卡戴珊不雅视频在线播放| 免费av不卡在线播放| 欧美国产日韩亚洲一区| 91麻豆精品激情在线观看国产| 日本一本二区三区精品| 久久久精品94久久精品| 国产精品99久久久久久久久| 亚洲五月天丁香| 成人特级黄色片久久久久久久| 99在线视频只有这里精品首页| 91午夜精品亚洲一区二区三区| 国产精品一二三区在线看| 美女免费视频网站| 色av中文字幕| 男女啪啪激烈高潮av片| 亚洲va在线va天堂va国产| 日日撸夜夜添| 国产精品人妻久久久影院| 1024手机看黄色片| 亚洲无线观看免费| 国产精品av视频在线免费观看| a级毛片a级免费在线| 老师上课跳d突然被开到最大视频| 国产男靠女视频免费网站| 精品一区二区三区视频在线| 99久久九九国产精品国产免费| 国产高清有码在线观看视频| 女同久久另类99精品国产91| 亚洲成人久久爱视频| 亚洲中文字幕一区二区三区有码在线看| av卡一久久| 婷婷精品国产亚洲av| 久久久久国产网址| 一夜夜www| 亚洲精品一卡2卡三卡4卡5卡| 久久韩国三级中文字幕| 久久久久久久久久成人| 国产精品久久久久久精品电影| 97超视频在线观看视频| 欧美中文日本在线观看视频| h日本视频在线播放| 免费搜索国产男女视频| 久久久久国产精品人妻aⅴ院| 中国国产av一级| 啦啦啦韩国在线观看视频| 国产男人的电影天堂91| 黄色欧美视频在线观看| 有码 亚洲区| 午夜精品一区二区三区免费看| 美女黄网站色视频| 国产高清三级在线| 久久久久久伊人网av| 麻豆av噜噜一区二区三区| 国内精品美女久久久久久| 久久久久国内视频| 波野结衣二区三区在线| 亚洲五月天丁香| 波野结衣二区三区在线| 欧美激情在线99| 国产成人一区二区在线| 嫩草影院精品99| 国产午夜精品论理片| 日本免费a在线| 伊人久久精品亚洲午夜| 女人十人毛片免费观看3o分钟| 少妇的逼水好多| 国产精品,欧美在线| 国产精品久久久久久久久免| 久久久午夜欧美精品| 国产成人aa在线观看| 人妻夜夜爽99麻豆av| 久久九九热精品免费| 欧美三级亚洲精品| 欧美色视频一区免费| 色哟哟哟哟哟哟| 国产欧美日韩一区二区精品| 日韩欧美精品v在线| 久久久久精品国产欧美久久久| 男女边吃奶边做爰视频| 日韩精品中文字幕看吧| 国产片特级美女逼逼视频| 欧美日韩在线观看h| 又爽又黄a免费视频| 国产精品日韩av在线免费观看| 国产免费男女视频| 久久国内精品自在自线图片| 亚洲美女黄片视频| 乱码一卡2卡4卡精品| 露出奶头的视频| 网址你懂的国产日韩在线| 在线免费十八禁| 五月伊人婷婷丁香| 天堂网av新在线| 日韩欧美免费精品| 久久精品91蜜桃| 成年免费大片在线观看| 丰满人妻一区二区三区视频av| 在线播放无遮挡| 久久午夜福利片| ponron亚洲| 午夜精品国产一区二区电影 | 国产免费一级a男人的天堂| 国产精品久久视频播放| 黄色一级大片看看| 偷拍熟女少妇极品色| 亚洲中文日韩欧美视频| 亚洲国产欧美人成| 尾随美女入室| 国产亚洲91精品色在线| 精品久久久久久久久久久久久| 一a级毛片在线观看| 久久久久久久久久久丰满| 色吧在线观看| 亚洲在线自拍视频| 亚洲国产欧美人成| 尾随美女入室| 亚洲欧美日韩卡通动漫| 变态另类成人亚洲欧美熟女| 国产亚洲av嫩草精品影院| 国产高清视频在线播放一区| 在线免费观看不下载黄p国产| 亚洲乱码一区二区免费版| 精品日产1卡2卡| 久久99热这里只有精品18| 又黄又爽又刺激的免费视频.| 变态另类丝袜制服| 国产精品人妻久久久久久| 男女之事视频高清在线观看| 色综合色国产| 春色校园在线视频观看| 女同久久另类99精品国产91| 97热精品久久久久久| 熟妇人妻久久中文字幕3abv| 搡老妇女老女人老熟妇| 不卡视频在线观看欧美| 久久6这里有精品| av福利片在线观看| 最近手机中文字幕大全| 国产免费一级a男人的天堂| 最近2019中文字幕mv第一页| 最好的美女福利视频网| 特大巨黑吊av在线直播| 久久久成人免费电影| 国产探花在线观看一区二区| 欧美性感艳星| 伊人久久精品亚洲午夜| 久久天躁狠狠躁夜夜2o2o| 日韩av不卡免费在线播放| 亚洲国产精品国产精品| 成人漫画全彩无遮挡| 成年女人看的毛片在线观看| 成人二区视频| 国产一区二区在线av高清观看| 亚洲av二区三区四区| 日本撒尿小便嘘嘘汇集6| 免费不卡的大黄色大毛片视频在线观看 | 热99在线观看视频| 免费人成在线观看视频色| 99久久无色码亚洲精品果冻| 神马国产精品三级电影在线观看| av黄色大香蕉| 男女之事视频高清在线观看| ponron亚洲| 久久精品国产鲁丝片午夜精品| 欧美性感艳星| 两个人的视频大全免费| 日本在线视频免费播放| 深爱激情五月婷婷| 12—13女人毛片做爰片一| 免费黄网站久久成人精品| 嫩草影院精品99| 永久网站在线| 亚洲久久久久久中文字幕| 亚洲最大成人av| 国产午夜福利久久久久久| 精品99又大又爽又粗少妇毛片| 床上黄色一级片| 欧美日韩综合久久久久久| 干丝袜人妻中文字幕| 天堂√8在线中文| 国产白丝娇喘喷水9色精品| 免费在线观看成人毛片| 亚洲自偷自拍三级| 三级经典国产精品| 国产69精品久久久久777片| 天美传媒精品一区二区| 全区人妻精品视频| 99热全是精品| 精品国产三级普通话版| 中文字幕久久专区| 在线观看美女被高潮喷水网站| 中文亚洲av片在线观看爽| 久久精品国产亚洲网站| 18禁黄网站禁片免费观看直播| av黄色大香蕉| 成年女人永久免费观看视频| 欧美人与善性xxx| 天天躁日日操中文字幕| 国产女主播在线喷水免费视频网站 | 可以在线观看的亚洲视频|