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

    逐次Taylor替換與一類冪指函數(shù)不等式的機(jī)器證明

    2019-11-26 05:36:06陳世平
    關(guān)鍵詞:展開式代數(shù)區(qū)間

    陳世平,陳 果

    (1.中國民航飛行學(xué)院德陽校區(qū)-四川省商貿(mào)學(xué)校,四川 德陽 618000;2.西南科技大學(xué)信息工程學(xué)院,四川 綿陽 612000)

    0 引 言

    Nishizawa 證明了[10]:

    Branko Male?evic 等證明了文獻(xiàn)[10]提出的公開問題[11]:

    文獻(xiàn)[12-17]討論了形如f(x,trans1(x),…,transn(x))>0 的超越多項(xiàng)式不等式的機(jī)器證明,其中f(x,x1,…,xn)是n+1 元多項(xiàng)式,trans1(x)可以是初等超越函數(shù),如sin x、cos x、tan x、arcsin x、arccos x、arctan x、ex、ln x 等,也可以是形如trans1(trans2(x))的復(fù)合函數(shù),利用一次或多次Taylor 替換,目標(biāo)不等式的證明轉(zhuǎn)化為一系列單變量的多項(xiàng)式不等式的驗(yàn)證,然后由代數(shù)不等式證明軟件BOTTEMA 完成.算法對于常見的超越多項(xiàng)式不等式非常有效,而且是“可讀的”,可以輸出能夠理解的證明過程.這種方案被稱為逐次Taylor 替換法.

    本文的目的是討論定理2 中不等式的機(jī)器證明,這是一個(gè)比文獻(xiàn)[12-17]中引用的例子更復(fù)雜、超越多項(xiàng)式不等式范圍外形如sin(x)/x>u(x)v(x)的冪指函數(shù)不等式問題,運(yùn)用逐次Taylor 替換并結(jié)合人工證明和巧妙的變量代換,實(shí)現(xiàn)了該不等式的機(jī)器證明.通過與文獻(xiàn)[11]不同的途徑,也證明了文獻(xiàn)[10]提出的公開問題,盡管結(jié)論是已知的結(jié)果,但方法本身對于同類不等式具有示范性.

    此外,算法輸出的證明過程中除了有的需要判定其正定性的一元有理多項(xiàng)式次數(shù)較高從而驗(yàn)證計(jì)算量較大外,其余內(nèi)容都容易“手工”驗(yàn)算.

    1 逐次Taylor替換

    不等式的自動證明問題,一直以來都是數(shù)學(xué)機(jī)械化和自動推理領(lǐng)域研究熱點(diǎn)和難點(diǎn)問題.近年來,代數(shù)不等式的自動證明吸引了大批學(xué)者研究,也取得了豐富的成果,已有專著《不等式機(jī)器證明與自動發(fā)現(xiàn)》問世,特別是楊路教授提出了降維算法,據(jù)之編制的通用程序BOTTEMA 能夠成批驗(yàn)證不等式,但BOTTEMA 只能處理代數(shù)不等式或先把幾何不等式轉(zhuǎn)化為代數(shù)不等式,尚不能解決一般的超越不等式自動證明問題[18-22].當(dāng)然關(guān)于超越不等式已有了大量的研究成果和證明方法,但這些成果和方法尚不能適應(yīng)數(shù)學(xué)機(jī)械化和推理自動化的需要,故有必要對此類不等式的自動證明進(jìn)行專門的研究,探索新的算法.

    1.1 Taylor替換

    給定(n+1)元實(shí)系數(shù)多項(xiàng)式環(huán)R[x,x1,…,xn],定義該環(huán)上一個(gè)映射:hom:f(x,x1,…,xn)→F(x),將變量xi用某個(gè)超越函數(shù)transi(x)代替.

    定義1.1 給定(n+1)元實(shí)系數(shù)多項(xiàng)式f(x,x1,…,xn),它在映射hom 下的像F(x)=hom(f)=f(x,trans1(x),…,transn(x))稱為超越函數(shù)多項(xiàng)式,稱transi(x)為其超越因子.

    如x+sin(x)、x+cos(sqrt(x))、x+esin(x)、x+arctan(x)-arcsin(x)都是超越函數(shù)多項(xiàng)式,sin(x)、cos(sqrt(x))、esin(x)、arctan(x)、arcsin(x)就是函數(shù)的超越因子.

    定義 1.2 對于超越函數(shù) F(x),在某個(gè)區(qū)間 I,如果有代數(shù)函數(shù)列{T_min(n,F(xiàn))}和{T_max(n,F(xiàn))},存在自然數(shù) n0,當(dāng) n≥n0時(shí),滿足:

    1)T_min(n+1,F(xiàn))>T_min(n,F(xiàn)),且當(dāng) n→∞ 時(shí),T_min(n,F(xiàn))→F(x);

    2)T_max(n+1,F(xiàn))<T_max(n,F(xiàn)),且當(dāng) n→∞ 時(shí),T_max(n,F(xiàn))→F(x).

    我們稱{T_min(n,F(xiàn))}和{T_max(n,F(xiàn))}為 F(x)在區(qū)間 I 的下限多項(xiàng)式列和上限多項(xiàng)式列,T_min(n,F(xiàn))為 F(x)的下限多項(xiàng)式,T_max(n,F(xiàn))為 F(x)的上限多項(xiàng)式,n0稱為臨界值.

    顯然,F(xiàn)(x)的下限多項(xiàng)式和上限多項(xiàng)式滿足:

    T_min(n0,F(xiàn))(x)<T_min(n0+1,F(xiàn))(x)<T_min(n0+2,F(xiàn))(x)<…<F(x)<…<T_max(n0+2,F(xiàn))(x)<T_max(n0+1,F(xiàn))(x)<T_max(n0,F(xiàn))(x),也就是說若 F(x)存在有下限多項(xiàng)式列和上限多項(xiàng)式列,則有一個(gè)多項(xiàng)式套來逼近 F(x):(T_min(n0,F(xiàn))(x),T_max(n0,F(xiàn))(x))?(T_min(n0+1,F(xiàn))(x),T_max(n0+1,F(xiàn))(x))?(T_min(n0+2,F(xiàn))(x),T_max(n0+2,F(xiàn))(x))?……?{F(x)}.

    下面我們討論如何計(jì)算或構(gòu)造超越函數(shù)多項(xiàng)式F(x)=f(x,trans(x))的下限多項(xiàng)式和上限多項(xiàng)式(其中f(x,y)為二元多項(xiàng)式).

    為方便,我們將函數(shù)f(x)在0 點(diǎn)的Taylor 展開式中的前n 項(xiàng)之和記為taylor(f,n).

    定義1.3 如果超越函數(shù)f(x)在某個(gè)區(qū)間(0,T]滿足以下條件:(1)大于0;(2)在0 點(diǎn)的Taylor 展開式收斂于f(x);(3)taylor(f,n)可以表示為交錯(cuò)級數(shù)f1(x)-f2(x)+f3(x)+…+(-1)n-1fn(x),其中fn(x)=An*x^mn,0<An≤1,mn-1<mn;(4)存在常數(shù)n0,當(dāng)n≥n0時(shí),taylor(f,n)>0,fn(x)>fn+1(x)>0.我們稱f(x)在區(qū)間(0,T]內(nèi)可以規(guī)范展開,常數(shù)n0稱為其臨界值.若還滿足An≤1/(n-1)!,則稱f(x)在區(qū)間(0,T]內(nèi)可以強(qiáng)規(guī)范展開.

    常見的基本初等超越函數(shù)大都可以在相應(yīng)的區(qū)間內(nèi)規(guī)范展開:

    arctan(x)=x-x3/3+x5/5-x7/7+…+(-1)n-1x2n-1/(2n-1)+……在區(qū)間(0,1]可以規(guī)范展開,臨界值n0=1;

    e-x=1-x+x2/2!-x3/3!+…+(-x)n-1/(n-1)?。?,在區(qū)間(0,T](T 為任意正數(shù),下同)可以規(guī)范展開,臨界值n0=min{n∈N,當(dāng)0<x≤T 時(shí)taylor(e-x,n)>0,taylor(e-x,n+1)>0,n>T};

    ln(1+x)=x-x2/2+x3/3-……+在區(qū)間(0,1]可以規(guī)范展開,臨界值n0=1;

    sin(x)=x-x3/3!+x5/5!-……+在區(qū)間(0,T]可以規(guī)范展開,其中臨界值n0=min{n∈N,當(dāng)0<x≤T 時(shí)taylor(sin(x),n)>0,taylor(sin(x),n+1)>0 且2n(2n+1)>T2},特別地,當(dāng)時(shí),n0=1;

    cos(x)=1-x2/2!+x4/4!-……+在區(qū)間(0,T]可以規(guī)范展開,其中臨界值n0=min{n∈N,當(dāng)0<x≤T 時(shí)taylor(cos(x),n)>0,taylor(cos(x),n+1)>0 且 2n(2n-1)>T2}.我們注意到,當(dāng)時(shí),taylor(cos(x),2n)<cos(x)=0,所以 cos(x)在內(nèi)不能規(guī)范展開;

    (1+x)a=1+ax+a(a-1)x2/2!+a(a-1)(a-2)x3/3!+C(a,n)xn+…,在區(qū)間(0,1)可以規(guī)范展開(其中 C(a,n)=a(a-1)(a-2)…(a-n+1)).當(dāng) 0<a<1 時(shí),n0=1(展開式中f1(x)=1+ax).

    其中 e-x、sin(x)、cos(x)在相應(yīng)區(qū)間內(nèi)還可強(qiáng)規(guī)范展開.

    也有一些基本初等超越函數(shù)不能規(guī)范展開,如:

    ln(1-x)=-x-x2/2-x3/3-……-xk/k-……(0<x<1);

    tan(x)=ΣB2n(-4)n(1-4n)x(2n-1)/(2n)!=x+(1/3)x3+(2/15)x5+(17/315)x7+(62/2835)x9+……

    引理1.1 如果f(x)區(qū)間I 可以規(guī)范展開,臨界值為n0,則當(dāng)n>n0時(shí),在區(qū)間I 上

    1)taylor(f,2n-2)<taylor(f,2n)<f(x);

    2)taylor(f,2n-1)>taylor(f,2n+1)>f(x);

    3)當(dāng)n→∞時(shí),taylor(f,n)→f(x).

    為方便,引入如下記號:

    acrtan_max(x,n)=taylor(arctan(x),2n+1),arctan_min(x,n)=taylor(arctan(x),2n),則當(dāng) x∈(0,1],對任意 n∈N,arctan_max(x,n)>arctan(x)>arctan_min(x,n);

    ln_max(x,n)=taylor(ln(1+x),2n+1),ln_min(x,n)=taylor(ln(1+x),2n),則當(dāng)x∈(0,1),對任意 n∈N,ln_max(x,n)>ln(1+x)>ln_min(x,n);

    sin_max(x,n)=taylor(sin(x),2n+1),sin_min(x,n)=taylor(sin(x),2n),則當(dāng)x∈對任意 n∈N,sin_max(x,n)>sin(x)>sin_min(x,n);

    cos_max(x,n)=taylor(cos(x),2n+1),cos_min(x,n)=taylor(cos(x),2n),則當(dāng)對任意 n∈N,cos_max(x,n)>cos(x)>cos_min(x,n);

    用f+和f-分別表示多項(xiàng)式f 展開后的正項(xiàng)和負(fù)項(xiàng)之和,顯然f=f++f-,并且下面的引理成立:

    引理1.2 假設(shè)實(shí)函數(shù)T1(y)>0,T2(y)>0 且T1(y)<x<T2(y),則

    f+(T1(y),y)+f-(T2(y),y)<f(x,y)<f+(T2(y),y)+f-(T1(y),y).

    定理1.1 F(x)=f(x,trans(x)),trans(x)在區(qū)間I 可以規(guī)范展開,則對任意n∈N,

    T_max(n,F(xiàn))=f+(x,taylor(trans(x),2n-1))+f-(x,taylor(trans(x),2n))為 F(x)的上限多項(xiàng)式;

    T_min(n,F(xiàn))=f+(x,taylor(trans(x),2n))+f-(x,taylor(trans(x),2n-1))為 F(x)的下限多項(xiàng)式.

    定理 1.2 設(shè){T_min(n,F(xiàn))}和{T_max(n,F(xiàn))}分別是超越函數(shù)多項(xiàng)式 F(x)在區(qū)間 I 的下限多項(xiàng)式列和上限多項(xiàng)式列,則在相應(yīng)區(qū)間內(nèi),

    1)若存在 n0使得 T_min(n0,F(xiàn))≥0 成立,則 F(x)>0 成立;

    2)若存在 n0使得 T_max(n0,F(xiàn))≤0 成立,則 F(x)<0 成立;

    3)若存在 n0使得 T_max(n0,F(xiàn))>0 和 T_min(n0,F(xiàn))<0 都不成立,則 F(x)在區(qū)間 I的符號不固定.

    上述方案運(yùn)用Taylor 展開式構(gòu)造超越函數(shù)多項(xiàng)式的上下限多項(xiàng)式列,建立一個(gè)逼近目標(biāo)函數(shù)的一元多項(xiàng)式套,從而將目標(biāo)不等式的證明轉(zhuǎn)化為一系列的一元多項(xiàng)式不等式的驗(yàn)證,我們將這種方案稱為Taylor 替換法.

    1.2 逐次Taylor替換

    一般的超越函數(shù)多項(xiàng)式可能包含多個(gè)超越因子,超越因子可能還是由超越函數(shù)或代數(shù)函數(shù)復(fù)合而成,經(jīng)過一次Taylor 替換后的函數(shù)可能還是包含超越因子,還需要再次甚至多次的Taylor 替換.我們以超越函數(shù)多項(xiàng)式F(x)=f(x,trans1(x),trans2(x))和F(x)=f(x,trans1(trans2(x)))為例來描述這個(gè)過程:

    1、F(x)=f(x,trans1(x),trans2(x)),其中超越函數(shù)trans1(x)、trans2(x)在區(qū)間I 內(nèi)可以規(guī)范展開,臨界值分別為n1、n2.

    令f1(x,trans2(x),m)=f+(x,taylor(trans1(x),2*m),trans2(x))+f-(x,taylor(trans1(x),2*m-1),trans2(x));

    令f2(x,m,n)=f1+(x,taylor(trans2(x),2*n),m)+f1-(x,taylor(trans2(x),2*n-1),m);

    顯然當(dāng)m>n1時(shí),f1(x,trans2(x),m)<f(x,trans1(x),trans2(x)),當(dāng)m→∞時(shí),f1(x,trans2(x),m)→f(x,trans1(x),trans2(x));

    當(dāng)n>n2時(shí),對任意m>n1,f2(x,m,n)<f1(x,trans2(x),m),當(dāng)n→∞時(shí),f2(x,m,n)→f1(x,trans2(x),m);

    即當(dāng)m>n1,n>n2時(shí),f2(x,m,n)<f(x,trans1(x),trans2(x)),當(dāng)m→∞,n→∞時(shí),f2(x,m,n)→F(x).

    特別地,T_min(n,F(xiàn))=f2(x,n,n)就是F(x)的下限多項(xiàng)式,{T_min(F,n)}是其下限多項(xiàng)式列,其臨界值n0=max{n1,n2}.

    類似的方法可以求得F(x)的上限多項(xiàng)式列.

    同樣的方法我們可以計(jì)算含更多“超越因子”形如f(x,trans1(x),trans2(x),trans3(x),…,transm(x))的超越函數(shù)的上下限多項(xiàng)式列.

    2、F(x)=f(x,trans1(trans2(x))),其中trans2(x)在區(qū)間I 可以規(guī)范展開,其臨界值為n2,trans2(x)在區(qū)間I 上的值域?yàn)閰^(qū)間J,trans1(x)在J 內(nèi)可以規(guī)范展開,其臨界值為n1.

    令trans2(x)=u,得到F1(x,u)=f(x,trans1(u)),

    令f1(x,m,u)=f+(x,taylor(trans1(u),2m))+f-(x,taylor(trans1(u),2m-1)),

    顯然,當(dāng)m>n1時(shí)f1(x,m,u)<F1(x,u)=f(x,trans1(u)),且當(dāng)m→∞時(shí),f1(x,m,u)→F1(x,u);

    若trans2(x)是代數(shù)函數(shù),則f1(x,m,trans2(x))就是F(x)的下限多項(xiàng)式,否則

    令f2(x,m,n)=f1+(x,m,taylor(trans2(x),2n))+f1-(x,m,taylor(trans2(x),2n-1)),

    當(dāng)n>n2時(shí),f2(x,m,n)<f1(x,m,trans2(x))<F(x)=f(x,trans1(trans2(x))),且當(dāng)n→∞時(shí),f2(x,m,n)→f1(x,m,trans2(x));

    從而當(dāng)m→∞,n→∞時(shí),f2(x,m,n)→F(x)=f(x,trans1(trans2(x))).

    特別地,T_min(F,n)=f2(x,n,n)是F(x)的下限多項(xiàng)式,{T_min(F,n)}是其下限多項(xiàng)式列.其臨界值n0=max{n1,n2}.

    類似的方法可以求得F(x)的上限多項(xiàng)式列.

    同樣的方法我們可以計(jì)算含形如trans1(trans2(trans3(x)))經(jīng)過更多次復(fù)合而成的超越因子的超越函數(shù)的上下限多項(xiàng)式列.

    我們把上述方案稱為逐次Taylor 替換,一個(gè)可以計(jì)算一般超越函數(shù)多項(xiàng)式f(x,trans1(x),…,transn(x))的上下限多項(xiàng)式列的逐次Taylor 替換過程可以遞歸地定義如下:

    算法 1.1 STS(Successive Taylor Substitution)

    輸入①超越函數(shù)多項(xiàng)式F(x);②n∈N;

    輸出F(x)的下(上)限多項(xiàng)式.

    1)若 F(x)為代數(shù)函數(shù)

    return[F(x),F(xiàn)(x)];

    算法結(jié)束;

    2)設(shè)trans(x)是F(x)的一個(gè)超越因子且可以表示為F(x)=f(x,trans(x))

    2.1)若trans(x)在相應(yīng)區(qū)間可以規(guī)范展開,則

    f1(x,n)←f+(x,taylor(trans(x),2n))+f-(x,taylor(trans(x),2n-1));

    T_min(n,F(xiàn))←STS(f1(x,n))[1];

    f2(x,n)←f+(x,taylor(trans(x),2n-1))+f-(x,taylor(trans(x),2n));

    T_max(n,F(xiàn))←STS(f2(x,n))[2];

    return[T_min(n,F(xiàn)),T_max(n,F(xiàn))];

    2.2)若trans(x)具有形如trans1(trans2(x))的復(fù)合形式,令F(x)=f(x,trans1(u)),其中u=trans2(x),假設(shè)trans1(u)在相應(yīng)區(qū)間可以規(guī)范展開.

    f1(x,n,u)←f+(x,taylor(trans1(u),2n))+f-(x,taylor(trans1(u),2n-1));

    T_min(n,F(xiàn))←STS(f1(x,n,trans2(x)))[1];

    f2(x,n,u)←f+(x,taylor(trans1(u),2n-1))+f-(x,taylor(trans1(u),2n));

    T_max(n,F(xiàn))←STS(f2(x,n,trans2(x)))[2];

    return[T_min(n,F(xiàn)),T_max(n,F(xiàn))];

    3)END.

    算法STS 還不能直接處理類似tan(x)、ln(1-x)、ex這類不能規(guī)范展開的超越因子,在證明過程中需要結(jié)合人工證明和做一些的變量代換.

    超越函數(shù)多項(xiàng)式的“超越因子”還可能就是超越常數(shù),也需要有理化.本文的處理方法是將超越數(shù)視為一個(gè)超越函數(shù),構(gòu)造一個(gè)逼近該超越數(shù)的有理數(shù)端點(diǎn)區(qū)間套.比如對于超越數(shù)t,我們構(gòu)造滿足如下性質(zhì)的區(qū)間列{[t1(n),t2(n)]}:

    1)t1(1)<t1(2)<…<t1(n)<…<t<…t2(n)<…t2(2)<t2(1);

    2)當(dāng)n→∞時(shí),t1(n)→t,t2(n)→t;

    3)端點(diǎn)t1(n)和t2(n)均為有理數(shù).

    若f(x,t)是系數(shù)含超越數(shù)t 的一元多項(xiàng)式,則

    T_max(n,f)=f+(x,t1(n))+f-(x,t2(n))為f(x,t)的上限多項(xiàng)式;

    T_min(n,f)=f+(x,t2(n))+f-(x,t1(n))為f(x,t)的下限多項(xiàng)式.

    輸入:n∈N;

    輸出:有理數(shù)端點(diǎn)區(qū)間[t1(n),t2(n)],滿足:對任意p≤n,t1(1)<t1(2)<…<t1(p)<<t2(p)<…<t2(2)<t2(1);當(dāng)n→∞ 時(shí),t1(n)→,t2(n)→.

    1)t1←16*subs(x=1/5,arctan_min(x,n))-4*subs(x=1/239,acrtan_max(x,n));

    2)t2←16*subs(x=1/5,acrtan_max(x,n))-4*subs(x=1/239,tarctan_min(x,n));

    3)return[t1,t2];算法結(jié)束.

    2 定理2的證明

    令f3(x)=numer(f2(x))*denom(f2(x))(其中numer(f(x))和denom(f(x))分別取分式的分子和分母,下同),顯然sgn(f2(x))=sgn(f3(x)).

    令g1(x,n)=subs(ln(x/sin(x))=B(x,n),f3+)+subs(ln(x/sin(x))=A(x,n),f3-),則對任意n∈N,g1(x,n)<f3(x);

    令g2(x,n)=numer(g1(x))*denom(g1(x));

    令g3(x,n)=subs(sin(x)=sin_min(x,n),g2+)+subs(sin(x)=sin_max(x,n),g2-),則對任意n∈N,g3(x,n)<g2(x,n);

    令g4(x,n)=subs(cos(x)=cos_min(x,n),g3+)+subs(cos(x)=cos_max(x,n),g3-),則對任意n∈N,g4(x,n)<g3(x,n),此時(shí)g4(x,n)還含有超越數(shù)作為系數(shù);

    令g5(x,n)=subs(=_to_ra(tn)[1],g4+)+subs(=_to_ra(tn)[2],g4-),則對任意n∈N,g5(x,n)是關(guān)于x 的一元有理多項(xiàng)式,且g5(y,n)<g4(x,n)<g3(x,n)<g2(x,n);

    使用BOTTEMA 的xprove 指令:xprove(g5(x,n)>0,[x<_to_rat(n)[2]/3]), 當(dāng)n=5 時(shí),不等式g5(x,n)>0 成立,即當(dāng)x∈(0,x<_to_ra(t5)[2]/3)(0/3],g2(x,5)>0,從而f3(x)>g1(x,5)>0,繼而f2(x)>0,得到f1(x)在(0/3]內(nèi)單調(diào)上升.

    證明f1(x)=-ln(u(x))-1/v(x)*ln(x/sin(x)),由于u(0)=1,v(0)=且當(dāng)x→0+,ln(x/sin(x))→0,所以當(dāng)x→0+,f1(x)→0.而由引理3.1,當(dāng)x∈(0,/3],f1(x)關(guān)于x單調(diào)上升,從而f1(x)>0.

    證明 令w3(x)=nume(rw1(x))*denom(w1(x)),得到

    令w4((x,n)=subs(=_to_ra(tn)[2],w3+)+subs(=_to_ra(tn)[1],w3-),則w4(x,n)是關(guān)于x 的有理多項(xiàng)式,且對任意n∈N,w4(x,n)>w3(x);

    使用BOTTEMA 的xprove 指令:xprove(w4(x,n)<0,[x<_to_ra(tn)[2]/2]),當(dāng)n=2 時(shí),w4(x,n)<0 成立,即在(0/2]內(nèi)0>w4(x,2)>w3(x),從而w1(x)<0.

    證明 令f4(x)=-(ln(x/sin(x))+w2(x)/w1(x)),由引理2.2 知,在(0,/2]內(nèi)w1(x)<0,所以在(/3,/2]內(nèi)sgn(f4(x))=sgn(f2(x)).

    令f5(x)=f4(`x),則f5(x)為包含x、sin(x)、cos(x)、的分式.由于 cos(x)在(/3,/2)內(nèi)不能規(guī)范展開,作變量替換x=/2-y,則f5(x)<0 在(/3,/2)內(nèi)成立等價(jià)于在(0/6)內(nèi)g(y)=f5(/2-y)<0 成立.

    令g1(y)=numer(g(y))*denom(g(y));

    令g2(y,n)=subs(sin(y)=sin_max(y,n),g1+)+subs(sin(y)=sin_min(y,n),g1-),則對任意n∈N,g2(y,n)>g1(y);

    令g3(y,n)=subs(cos(y)=cos_max(y,n),g2+)+subs(cos(y)=cos_min(y,n),g2-),則對任意n∈N,g3(y,n)>g2(y,n);

    令g4(y,n)=sub(s=_to_ra(tn)[2],g3+)+sub(s=_to_ra(tn)[1],g3-),則g4(y,n)是關(guān)于y 的一元有理多項(xiàng)式,且g4(y,n)>g3(y,n)>g2(y,n)>g1(y);

    使用BOTTEMA 的xprove 指令:xprove(g4(y,n)<0,[y<_to_rat(n)[2]/6]), 當(dāng)n=8 時(shí),不等式g4(y,n)<0 成立,即在(0/6)內(nèi)g1(y)<0 成立,g(y)<0 成立.從而在(/3/2)內(nèi)f5(x)<0 成立.由此得到f4(x)在(/3/2)內(nèi)單調(diào)下降.

    下面討論f2(/2)的符號.令

    h1=nume(rf2(/2))*denom(f2(/2))=其中和ln(/2)需要有理化,令t=/2-1,則 t∈(0,1),ln(1+t) 可以規(guī)范展開,即對任意 n∈N,ln_max(t,n)>ln(1+t)=ln(/2)>ln_min(t,n),令 E(n)=subs(t=/2-1,ln_max(t,n)),F(xiàn)(n)=subs(t=/2-1,ln_min(t,n)),則 E(n)>ln(/2)>F(n).

    令h2(n)=subs(ln(/2)=E(n),h1+)+subs(ln(/2)=F(n),h1-);

    令h3(n)=subs(=_to_ra(tn)[2],h2(n)+)+subs(=_to_ra(tn)[1],h2(n)-),則對任意n∈N,h3(n)>h2(n)>h1.h3(n)是含n 的有理式且容易得到h3(5)<0,所以h1<0,從而f2(/2)<0.

    f2(/2)<0 就得到f4(/2)<0 成立.又由引理2.1 知f2(/3)>0,所以f4(/3)>0,從而在(/3,/2)內(nèi)存在唯一x0使得f4(x0)=0,并且在(/3,x0)內(nèi)f4(x)>0,從而f2(x)>0,在(x0,/2)內(nèi)f4(x)<0,從而f2(x)<0.

    而f2(x)=f1(`x),所以f1(x)在(/3,x0)內(nèi)單調(diào)上升,在(x0/2)內(nèi)單調(diào)下降.由推論2.1 知f1(/3)>0,從而f1(x)>0 在(/3,x0]內(nèi)成立;而f1(/2)=0,所以f1(x)>0 在(x0/2)內(nèi)成立.所以f1(x)>0 在(/3,/2)內(nèi)成立.

    由推論2.1 和引理2.3 知f1(x)>0 在(0,/2)內(nèi)成立,從而 (fx)>0 在(0,/2)內(nèi)成立,即定理2 成立.

    3 結(jié)束語

    逐次Taylor 替換是處理超越函數(shù)多項(xiàng)式不等式自動證明的有效方案,其主要思路是借助Taylor 展開式建立一個(gè)逼近目標(biāo)函數(shù)的多項(xiàng)式套,從而將證明轉(zhuǎn)化為一系列的一元多項(xiàng)式不等式的驗(yàn)證,然后借助BOTTEMA 等代數(shù)不等式證明工具完成最后的工作.算法可以處理含多個(gè)不同超越因子和復(fù)合超越因子的超越函數(shù)多項(xiàng)式不等式,在證明過程中若結(jié)合人工證明和巧妙的代換方法,更能發(fā)揮逐次Taylor 替換算法的效能.

    本文解決了一個(gè)一類形如sin(x)/x>u(x)v(x)的冪指函數(shù)不等式的機(jī)器證明問題,是逐次Taylor 替換成功運(yùn)用的案例.雖然結(jié)論是已知結(jié)果,但方法本身對同類不等式具有示范性,能夠?yàn)槲墨I(xiàn)[1-11]中其它冪指函數(shù)不等式(或冪函數(shù)不等式)的機(jī)器提供證明方案.

    猜你喜歡
    展開式代數(shù)區(qū)間
    解兩類含參數(shù)的復(fù)合不等式有解與恒成立問題
    你學(xué)會“區(qū)間測速”了嗎
    兩個(gè)有趣的無窮長代數(shù)不等式鏈
    泰勒展開式在函數(shù)中的應(yīng)用
    Hopf代數(shù)的二重Ore擴(kuò)張
    什么是代數(shù)幾何
    科學(xué)(2020年1期)2020-08-24 08:08:06
    函數(shù)Riemann和式的類Taylor級數(shù)展開式
    區(qū)間對象族的可鎮(zhèn)定性分析
    一個(gè)非平凡的Calabi-Yau DG代數(shù)
    對一道冪級數(shù)展開式例題的思考
    国产午夜福利久久久久久| 舔av片在线| 中文字幕人妻丝袜一区二区| 香蕉av资源在线| 欧美日韩黄片免| 免费av观看视频| 亚洲不卡免费看| 亚洲av免费高清在线观看| 色综合站精品国产| 国产不卡一卡二| 1000部很黄的大片| 小蜜桃在线观看免费完整版高清| 97超级碰碰碰精品色视频在线观看| 最近最新中文字幕大全免费视频| 热99在线观看视频| 亚洲专区中文字幕在线| 亚洲最大成人中文| 午夜精品一区二区三区免费看| 99热这里只有精品一区| 好男人在线观看高清免费视频| 少妇裸体淫交视频免费看高清| 国产91精品成人一区二区三区| 免费看a级黄色片| 又紧又爽又黄一区二区| 99精品在免费线老司机午夜| 国产v大片淫在线免费观看| 青草久久国产| 日本爱情动作片www.在线观看| 日韩伦理黄色片| 国产一级毛片七仙女欲春2| 尾随美女入室| 国产精品麻豆人妻色哟哟久久 | 国产亚洲av片在线观看秒播厂 | 亚洲一区高清亚洲精品| 国产精品女同一区二区软件| 亚州av有码| 久久精品人妻少妇| 少妇熟女aⅴ在线视频| 最近最新中文字幕大全电影3| 日韩 亚洲 欧美在线| 国产高清有码在线观看视频| 午夜激情久久久久久久| 美女主播在线视频| 国内少妇人妻偷人精品xxx网站| 欧美不卡视频在线免费观看| 午夜视频国产福利| 精品人妻一区二区三区麻豆| 亚洲人与动物交配视频| 卡戴珊不雅视频在线播放| 国产国拍精品亚洲av在线观看| 久久精品久久久久久噜噜老黄| 日韩精品青青久久久久久| 精品国产一区二区三区久久久樱花 | 亚洲伊人久久精品综合| 神马国产精品三级电影在线观看| 看十八女毛片水多多多| 国产激情偷乱视频一区二区| 久久久色成人| 精品久久久久久久末码| 91久久精品国产一区二区三区| 99热网站在线观看| 欧美日韩精品成人综合77777| 久久久亚洲精品成人影院| 夫妻性生交免费视频一级片| 色5月婷婷丁香| 久久久精品免费免费高清| 高清在线视频一区二区三区| 晚上一个人看的免费电影| 一级av片app| 国产精品久久久久久av不卡| 99热网站在线观看| 欧美日韩在线观看h| 男女啪啪激烈高潮av片| 欧美三级亚洲精品| 欧美最新免费一区二区三区| 91狼人影院| 99久久精品一区二区三区| 日本三级黄在线观看| 嘟嘟电影网在线观看| 欧美xxxx性猛交bbbb| 男插女下体视频免费在线播放| 一本久久精品| 久久6这里有精品| 欧美一区二区亚洲| 精品久久久久久久久av| 国产精品无大码| 大话2 男鬼变身卡| 三级男女做爰猛烈吃奶摸视频| 欧美+日韩+精品| 日日摸夜夜添夜夜添av毛片| 久久精品熟女亚洲av麻豆精品 | 成人一区二区视频在线观看| 免费黄网站久久成人精品| 观看美女的网站| 亚洲av中文av极速乱| 中文字幕av在线有码专区| 性插视频无遮挡在线免费观看| 人人妻人人澡欧美一区二区| 亚洲av中文字字幕乱码综合| 日本欧美国产在线视频| 免费播放大片免费观看视频在线观看| 搞女人的毛片| a级毛色黄片| 激情五月婷婷亚洲| 国产精品一区二区性色av| 久久久久久久国产电影| 美女xxoo啪啪120秒动态图| 午夜福利网站1000一区二区三区| 特大巨黑吊av在线直播| 在线免费观看不下载黄p国产| 国产高清有码在线观看视频| 国产精品女同一区二区软件| 91久久精品国产一区二区三区| 亚洲av成人av| 99久久中文字幕三级久久日本| 最新中文字幕久久久久| 精华霜和精华液先用哪个| 自拍偷自拍亚洲精品老妇| 永久免费av网站大全| 22中文网久久字幕| 国产免费福利视频在线观看| 一区二区三区高清视频在线| 熟妇人妻久久中文字幕3abv| 啦啦啦韩国在线观看视频| av黄色大香蕉| 内地一区二区视频在线| 免费看日本二区| 91精品国产九色| 丝瓜视频免费看黄片| 丝袜喷水一区| 少妇的逼水好多| 十八禁网站网址无遮挡 | 亚洲国产精品国产精品| 最近最新中文字幕大全电影3| av专区在线播放| 成人欧美大片| 国产亚洲精品av在线| 免费大片黄手机在线观看| 九九在线视频观看精品| 午夜福利高清视频| 日韩人妻高清精品专区| 亚洲精品亚洲一区二区| 91午夜精品亚洲一区二区三区| 汤姆久久久久久久影院中文字幕 | 免费av毛片视频| av免费在线看不卡| 国产精品人妻久久久影院| 国产精品一区二区三区四区久久| 一个人免费在线观看电影| 免费大片18禁| 成人午夜高清在线视频| 大香蕉97超碰在线| 简卡轻食公司| 久久精品久久久久久久性| 亚洲精品一区蜜桃| 赤兔流量卡办理| 欧美日韩综合久久久久久| 成人毛片a级毛片在线播放| 97人妻精品一区二区三区麻豆| 亚洲怡红院男人天堂| 国产精品国产三级国产av玫瑰| 色哟哟·www| 国产精品精品国产色婷婷| av卡一久久| 久久久久九九精品影院| 国产精品国产三级国产专区5o| 美女高潮的动态| 亚洲国产精品国产精品| 欧美成人精品欧美一级黄| 插阴视频在线观看视频| av免费观看日本| 久久精品国产鲁丝片午夜精品| 成年免费大片在线观看| 成人无遮挡网站| 亚洲高清免费不卡视频| 亚洲精品成人av观看孕妇| 老师上课跳d突然被开到最大视频| 麻豆成人av视频| 寂寞人妻少妇视频99o| 日日撸夜夜添| 久久99热6这里只有精品| 亚洲经典国产精华液单| 久久久精品欧美日韩精品| 中文字幕免费在线视频6| 亚洲欧美精品自产自拍| 一区二区三区四区激情视频| 一级爰片在线观看| 美女主播在线视频| 亚洲欧美中文字幕日韩二区| 内射极品少妇av片p| 一级毛片 在线播放| 最近的中文字幕免费完整| 欧美激情国产日韩精品一区| 美女黄网站色视频| 免费看日本二区| 99视频精品全部免费 在线| 床上黄色一级片| 国产单亲对白刺激| av网站免费在线观看视频 | 精品人妻一区二区三区麻豆| 亚洲av中文av极速乱| 91精品伊人久久大香线蕉| 国产 一区 欧美 日韩| 男女边摸边吃奶| 最近中文字幕高清免费大全6| 中文在线观看免费www的网站| 亚洲精品视频女| 免费观看性生交大片5| 99热这里只有是精品50| 少妇丰满av| 在线 av 中文字幕| .国产精品久久| 欧美另类一区| 国产亚洲av片在线观看秒播厂 | 精品欧美国产一区二区三| av在线亚洲专区| 亚洲av男天堂| 女的被弄到高潮叫床怎么办| 国产免费福利视频在线观看| 少妇裸体淫交视频免费看高清| 最近最新中文字幕大全电影3| 国产精品国产三级专区第一集| 亚洲婷婷狠狠爱综合网| 国产成人一区二区在线| 麻豆国产97在线/欧美| 国产v大片淫在线免费观看| 欧美日韩一区二区视频在线观看视频在线 | 3wmmmm亚洲av在线观看| 一级黄片播放器| 成人亚洲欧美一区二区av| 淫秽高清视频在线观看| 十八禁国产超污无遮挡网站| 中文欧美无线码| 久久国内精品自在自线图片| 日韩,欧美,国产一区二区三区| 欧美精品国产亚洲| 国产久久久一区二区三区| 中文字幕av成人在线电影| 天天躁日日操中文字幕| 久久精品国产亚洲av涩爱| 亚洲一级一片aⅴ在线观看| 建设人人有责人人尽责人人享有的 | 好男人视频免费观看在线| 亚洲人成网站在线观看播放| 日韩制服骚丝袜av| 国产精品久久久久久久久免| 大话2 男鬼变身卡| 久久韩国三级中文字幕| 中文字幕av在线有码专区| 亚洲综合色惰| 精品99又大又爽又粗少妇毛片| av网站免费在线观看视频 | 伦精品一区二区三区| 一级毛片 在线播放| 麻豆av噜噜一区二区三区| 日本免费在线观看一区| 日韩av在线大香蕉| 欧美激情国产日韩精品一区| 国产成人精品婷婷| 久久久a久久爽久久v久久| 亚洲精品日韩av片在线观看| 亚洲国产精品成人久久小说| 女人久久www免费人成看片| 91久久精品国产一区二区三区| 国产乱人偷精品视频| 国产黄频视频在线观看| 97精品久久久久久久久久精品| 久久久久久久久久人人人人人人| 欧美性感艳星| 日韩电影二区| 日韩制服骚丝袜av| 欧美极品一区二区三区四区| av在线观看视频网站免费| 日本欧美国产在线视频| 精品国产一区二区三区久久久樱花 | 久久久久国产网址| 男女边吃奶边做爰视频| 夜夜爽夜夜爽视频| 97超视频在线观看视频| 日日啪夜夜爽| 久久国产乱子免费精品| 成年免费大片在线观看| 日日啪夜夜撸| videossex国产| 欧美变态另类bdsm刘玥| 国产亚洲午夜精品一区二区久久 | 精品国产一区二区三区久久久樱花 | 国产精品福利在线免费观看| 丰满乱子伦码专区| a级毛色黄片| 日韩,欧美,国产一区二区三区| 精品一区二区免费观看| 久久久久网色| 国产亚洲精品av在线| 久久精品久久久久久噜噜老黄| 好男人在线观看高清免费视频| 欧美高清成人免费视频www| 亚洲熟女精品中文字幕| 美女内射精品一级片tv| 久久久a久久爽久久v久久| 特大巨黑吊av在线直播| 精品酒店卫生间| 日韩大片免费观看网站| 亚洲成人久久爱视频| 欧美激情国产日韩精品一区| 亚洲av男天堂| 国产麻豆成人av免费视频| 亚洲性久久影院| 99热网站在线观看| 亚洲成人久久爱视频| 熟女电影av网| 插阴视频在线观看视频| 国产精品.久久久| 少妇裸体淫交视频免费看高清| 天堂av国产一区二区熟女人妻| 高清日韩中文字幕在线| 别揉我奶头 嗯啊视频| 91av网一区二区| 久久人人爽人人片av| 午夜福利网站1000一区二区三区| 亚洲高清免费不卡视频| 波野结衣二区三区在线| 观看美女的网站| 97超碰精品成人国产| 亚洲熟妇中文字幕五十中出| 乱码一卡2卡4卡精品| 亚洲欧美成人综合另类久久久| 天堂俺去俺来也www色官网 | 亚洲高清免费不卡视频| 免费av不卡在线播放| 在线观看美女被高潮喷水网站| 午夜免费观看性视频| 国产精品一区二区三区四区久久| av在线老鸭窝| 99久久中文字幕三级久久日本| 内地一区二区视频在线| 草草在线视频免费看| 中文字幕制服av| 亚洲国产精品sss在线观看| 青春草国产在线视频| 精品人妻偷拍中文字幕| 国产麻豆成人av免费视频| 我的老师免费观看完整版| 美女cb高潮喷水在线观看| 一个人看的www免费观看视频| 日日干狠狠操夜夜爽| 欧美人与善性xxx| 欧美性猛交╳xxx乱大交人| 秋霞伦理黄片| 国产精品三级大全| 日韩一区二区视频免费看| 一级片'在线观看视频| 国产综合懂色| 亚洲国产精品专区欧美| 最近最新中文字幕大全电影3| 久久久成人免费电影| 午夜精品国产一区二区电影 | 美女黄网站色视频| 国产一区有黄有色的免费视频 | 久热久热在线精品观看| 国产一区二区在线观看日韩| 欧美成人精品欧美一级黄| 深夜a级毛片| 免费看av在线观看网站| 亚洲精品成人久久久久久| 亚洲精品中文字幕在线视频 | 毛片女人毛片| 人妻夜夜爽99麻豆av| 国产精品熟女久久久久浪| 国产熟女欧美一区二区| 男人和女人高潮做爰伦理| 欧美xxⅹ黑人| 成年人午夜在线观看视频 | 性插视频无遮挡在线免费观看| 精品酒店卫生间| 亚洲aⅴ乱码一区二区在线播放| 日本熟妇午夜| 欧美激情在线99| 大香蕉久久网| av卡一久久| 国产精品一区二区在线观看99 | 熟妇人妻久久中文字幕3abv| 欧美一区二区亚洲| 成人一区二区视频在线观看| 免费播放大片免费观看视频在线观看| 高清av免费在线| 欧美日韩国产mv在线观看视频 | 国产黄a三级三级三级人| 女的被弄到高潮叫床怎么办| 我的老师免费观看完整版| av在线播放精品| 久久精品熟女亚洲av麻豆精品 | 嫩草影院入口| 国产午夜精品一二区理论片| 久久这里只有精品中国| 最近2019中文字幕mv第一页| 老司机影院毛片| 狂野欧美激情性xxxx在线观看| 精品久久久久久久久久久久久| 男人舔奶头视频| 国语对白做爰xxxⅹ性视频网站| 亚洲无线观看免费| av又黄又爽大尺度在线免费看| 亚洲人成网站高清观看| 麻豆国产97在线/欧美| 午夜视频国产福利| 青青草视频在线视频观看| 精品少妇黑人巨大在线播放| 亚洲欧美日韩卡通动漫| 亚洲成人av在线免费| 欧美日韩国产mv在线观看视频 | 国产伦在线观看视频一区| 又粗又硬又长又爽又黄的视频| 亚洲人成网站在线播| 在线观看一区二区三区| videos熟女内射| 国产精品国产三级国产av玫瑰| 国产免费又黄又爽又色| 精品一区在线观看国产| 不卡视频在线观看欧美| videossex国产| 免费观看av网站的网址| 国产美女午夜福利| 黄片无遮挡物在线观看| 蜜桃久久精品国产亚洲av| 免费av毛片视频| 国产成人精品一,二区| 亚洲精品色激情综合| 欧美97在线视频| 久久久久久国产a免费观看| 国产亚洲精品av在线| 综合色丁香网| 网址你懂的国产日韩在线| 一区二区三区乱码不卡18| 永久网站在线| 秋霞在线观看毛片| 国产真实伦视频高清在线观看| 久久久久久久大尺度免费视频| 日韩制服骚丝袜av| 大又大粗又爽又黄少妇毛片口| 色综合站精品国产| 国产成人freesex在线| 六月丁香七月| 日本-黄色视频高清免费观看| 欧美3d第一页| 五月伊人婷婷丁香| videos熟女内射| 国产国拍精品亚洲av在线观看| 人人妻人人看人人澡| 日本av手机在线免费观看| 免费无遮挡裸体视频| 春色校园在线视频观看| 男的添女的下面高潮视频| 国产精品国产三级国产专区5o| 五月伊人婷婷丁香| 国产 一区 欧美 日韩| 97人妻精品一区二区三区麻豆| 亚洲av电影不卡..在线观看| 亚洲国产精品专区欧美| 美女被艹到高潮喷水动态| 丰满乱子伦码专区| 黄色配什么色好看| 国产亚洲精品久久久com| 日韩一本色道免费dvd| 又爽又黄无遮挡网站| 亚洲怡红院男人天堂| 精品久久久噜噜| 噜噜噜噜噜久久久久久91| 欧美激情国产日韩精品一区| 国模一区二区三区四区视频| 午夜福利高清视频| 如何舔出高潮| 国产精品一及| 午夜免费男女啪啪视频观看| 女人十人毛片免费观看3o分钟| 国产熟女欧美一区二区| 亚洲精品第二区| 黄片wwwwww| 国产 一区 欧美 日韩| 国产一区二区亚洲精品在线观看| 日本熟妇午夜| 大香蕉97超碰在线| 2021天堂中文幕一二区在线观| 中文字幕亚洲精品专区| 日韩三级伦理在线观看| 国产伦理片在线播放av一区| 大片免费播放器 马上看| 麻豆久久精品国产亚洲av| 99热这里只有是精品在线观看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲在线观看片| 狂野欧美激情性xxxx在线观看| 少妇高潮的动态图| 国产精品熟女久久久久浪| 日产精品乱码卡一卡2卡三| 精品久久久久久成人av| 久久久国产一区二区| 欧美潮喷喷水| 国产成人精品婷婷| 好男人在线观看高清免费视频| 全区人妻精品视频| 日本av手机在线免费观看| 十八禁网站网址无遮挡 | 亚洲美女视频黄频| 非洲黑人性xxxx精品又粗又长| 国内精品美女久久久久久| 91久久精品国产一区二区成人| 中文在线观看免费www的网站| 国产精品无大码| 日韩一区二区三区影片| 国产精品一区二区在线观看99 | 少妇人妻精品综合一区二区| 国产高潮美女av| 亚洲国产色片| 国产精品久久久久久久电影| 精品久久久久久久末码| 看十八女毛片水多多多| 观看美女的网站| 亚洲性久久影院| 狂野欧美激情性xxxx在线观看| 波多野结衣巨乳人妻| 少妇的逼好多水| 日韩一区二区视频免费看| 美女黄网站色视频| 99热这里只有是精品50| 蜜臀久久99精品久久宅男| 免费看av在线观看网站| 人人妻人人澡欧美一区二区| 亚洲精品第二区| av在线观看视频网站免费| 青春草国产在线视频| 日产精品乱码卡一卡2卡三| 狂野欧美白嫩少妇大欣赏| 两个人的视频大全免费| 性色avwww在线观看| 高清午夜精品一区二区三区| 久久久久九九精品影院| 一级毛片黄色毛片免费观看视频| 婷婷色综合www| 日本黄大片高清| 男人爽女人下面视频在线观看| av在线蜜桃| 午夜爱爱视频在线播放| 久久久亚洲精品成人影院| 久久精品夜色国产| 免费高清在线观看视频在线观看| 日本黄色片子视频| 国产真实伦视频高清在线观看| 久热久热在线精品观看| 日韩伦理黄色片| 亚洲精品一区蜜桃| 一级毛片 在线播放| 国产成人a区在线观看| 精品一区在线观看国产| 精品不卡国产一区二区三区| 青春草国产在线视频| 最近最新中文字幕免费大全7| 日韩成人伦理影院| 99久久精品国产国产毛片| 深夜a级毛片| 夜夜爽夜夜爽视频| 日日摸夜夜添夜夜爱| 久久热精品热| 精品亚洲乱码少妇综合久久| 国产视频首页在线观看| 国产老妇女一区| 美女高潮的动态| 精品午夜福利在线看| 国产伦在线观看视频一区| 2021少妇久久久久久久久久久| 亚洲精品日韩av片在线观看| 成人一区二区视频在线观看| 性插视频无遮挡在线免费观看| 国产黄色免费在线视频| 极品少妇高潮喷水抽搐| 欧美日韩视频高清一区二区三区二| 国产精品久久久久久久电影| 亚洲精品456在线播放app| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产又色又爽无遮挡免| 性色avwww在线观看| 欧美日韩在线观看h| 精品一区二区三卡| 色播亚洲综合网| 婷婷色麻豆天堂久久| 中文在线观看免费www的网站| 少妇裸体淫交视频免费看高清| 我的老师免费观看完整版| 熟妇人妻不卡中文字幕| 狂野欧美激情性xxxx在线观看| 一级毛片黄色毛片免费观看视频| 国产一级毛片在线| 色哟哟·www| 精品少妇黑人巨大在线播放| 嫩草影院入口| 国产精品综合久久久久久久免费| 久久精品久久久久久久性| 国产精品av视频在线免费观看| 丝袜喷水一区| 精品久久国产蜜桃| 国产午夜精品一二区理论片| 日本免费a在线| .国产精品久久| 男人爽女人下面视频在线观看| 国产av在哪里看| 熟妇人妻不卡中文字幕| av又黄又爽大尺度在线免费看| 久久久欧美国产精品| 亚洲欧美日韩卡通动漫| 波野结衣二区三区在线| 亚洲国产色片| 国产精品无大码| 极品少妇高潮喷水抽搐| 国产一区二区亚洲精品在线观看| 日韩国内少妇激情av|