浙江大学学报(理学版), 2022, 49(6): 662-669 doi: 10.3785/j.issn.1008-9497.2022.06.003

数学与计算机科学

非线性边界条件下拟线性瞬态方程组的Phragmén-Lindelöf型二择一结果

李远飞,1, 肖胜中2, 曾鹏1, 欧阳柏平1

1.广州华商学院,广东 广州 511300

2.广东农工商职业技术学院, 广东 广州 510507

Phragmén-Lindelöf type alternative results for quasilinear transient equations with nonlinear boundary conditions

LI Yuanfei,1, XIAO Shengzhong2, ZENG Peng1, OUYANG Baiping1

1.Guangzhou Huashang College,Guangzhou 511300,China

2.Guangdong AIB College,Guangzhou 510507,China

收稿日期: 2020-04-13  

基金资助: 广东省普通高校创新团队项目.  2020WCXTD008

Received: 2020-04-13  

作者简介 About authors

李远飞(1982—),ORCID:https://orcid.org/0000-0002-9314-4104,男,博士,教授,主要从事偏微分方程研究. 。

摘要

考虑了一类定义在三维半无穷柱体上的拟线性方程组,其中假设方程的解在柱体的有限端和侧面满足非齐次条件。定义了“能量”表达式,通过限制非线性项,利用微分不等式技术, 推导了一阶微分不等式,解此不等式得到二择一结果,即证明了“能量”随与有限端距离的增大要么呈指数式(多项式)增加,要么呈指数式(多项式)衰减。同时,在衰减情形下得到了全能量的上界。

关键词: 拟线性方程 ; 非线性边界条件 ; Phragmén-Lindelöf型二择一 ; 空间衰减性

Abstract

In this paper, we consider a class of quasilinear equations defined on a three-dimensional semi-infinite cylinder, in which the solutions of the equations are assumed to satisfy the nonhomogeneous conditions on both the finite end and the side of the cylinder. An expression of "energy" is defined. By limiting the nonlinear terms and making full use of the differential inequality technique, we obtain a first order differential inequality. By solving this inequality, we prove the alternative results, i.e., the "energy" increases exponentially (polynomial) or decays exponentially (polynomial) with the distance from the finite end. Finally, in the case of decay, we get the upper bound of total energy.

Keywords: quasilinear equation ; nonlinear conditions ; Phragmén-Lindelöf type alternative ; spatial decay

PDF (485KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

李远飞, 肖胜中, 曾鹏, 欧阳柏平. 非线性边界条件下拟线性瞬态方程组的Phragmén-Lindelöf型二择一结果. 浙江大学学报(理学版)[J], 2022, 49(6): 662-669 doi:10.3785/j.issn.1008-9497.2022.06.003

LI Yuanfei, XIAO Shengzhong, ZENG Peng, OUYANG Baiping. Phragmén-Lindelöf type alternative results for quasilinear transient equations with nonlinear boundary conditions. Journal of Zhejiang University(Science Edition)[J], 2022, 49(6): 662-669 doi:10.3785/j.issn.1008-9497.2022.06.003

1856年SAINT-VENANT1提出的数学和力学方程,后被称为Saint-Venant原理,广受应用数学领域研究者的关注和深入研究。其含义为“能量”随与柱体有限端距离的增大逐渐衰减为零2-7。几乎所有的研究均需要方程的解在无穷远处衰减于零或趋近于瞬态层流的先验假设,且通常假设在柱体的侧面满足零边界条件。

20世纪90年代后,Phragmén-Lindelöf型二择一原理被提出,并逐步成为研究热点,不必再假设方程的解满足Saint-Venant原理,只要证明“能量”随与柱体有限端距离的增加要么无限增大要么衰减为零。Phragmén-Lindelöf型二择一原理在物理学、力学和生物学等领域具有巨大的应用前景8-14

通常的做法是先定义一个半无穷的柱体:

R=x1x2x3|x1x2Dx3>0

其中,Dx1Ox2平面上的一个有界区域且具有光滑的边界D。用Rz 表示R的一个子区域:

Rz={x1x2x3|x1x2Dx3z>0}

其中,zx3 轴上的一个动点。用Dz 表示Rx3=z处的横截面:

Dz=x1x2x3|x1x2D x3=z>0

假设方程的解在柱体的有限端满足非零边界条件,在柱体的侧面满足齐次Dirichlet边界条件或齐次Neumann边界条件。用微分不等式技术证明“能量”要么呈指数式(多项式)增加,要么呈指数式(多项式)衰减。如LIN13考虑了定义在R×(0,)上的稳态拟线性方程

(ρ(x,u,u)u,i),i=f(u)

其中,逗号表示求导,重复英文下标表示从1到3求和,如(ρ(x,u,u)u,i),i=i=13xiρ(x,u,u)uxi。在柱体的侧面,解满足齐次Dirichlet边界条件或齐次Neumann边界条件

u=0un=0 D×0

对非线性项做一定限制,可得到解的二择一定理。

考虑瞬态拟线性方程

ut=(ρ(x,u,u)u,i),i-f(u), (x,t)R×0

在柱体侧面施加非齐次Dirichlet边界条件

u+h(u)=0 D×0×0

其中,hu)为已知函数。在柱体的有限端,有

u=g(x1,x2,t)    (x1,x2,t)D×0

初始条件为

u(x,0)=0,  xR

式(1)中,=(x1,x2,x2)表示梯度算子,f(u)为大于零的非线性函数,满足

uf(u)k1u2qq-1,    q>1,k1>0

且假设ρ满足

0<c1ρM1+M2(ρ|u|2)1q

其中,c1M1M2均为大于零的常数。

式(2)中,hu)满足

uh(u)k2|u|2p,    p>12,k2>0

式(3)中,g为大于零的给定函数且满足兼容性条件g=0,在D0×0gx1x2, 0)=0。

实际情况是研究的物理模型很难满足齐次边界条件,因此本研究更有意义。据知,目前除文献[811]外,很少见关注非齐次边界条件下二择一抛物问题的文献。文献[8]研究的是三维柱体上调和方程的二择一拋物问题,文献[11]则假设调和方程在柱体侧面满足不同的非零边界条件,得到解的二择一定理。本文的模型更复杂,无法直接应用文献[811]中的方法。

1 基本不等式

在研究二择一问题时,大多文献通常用wL2积分控制wL2积分,即若w|D=0,则

λDw2dADw,αw,αdA

其中,λ为二维固定膜问题的最小特征值,w,αw,α表示α=12w,αw,α。显然此不等式在本文条件下不成立。

为证明本文的主要结果, 结合文献[58]的方法, 推导常用的微分不等式。

引理1 存在一个依赖于区域D的大于零的常数C1,使得

D|w|2dAC1D|w|dl2+Dw,αw,αdA

证明PD内一点。令P1P2分别表示过点P平行于x1坐标轴的直线与D的交点, 令Q1Q2分别为过点P平行于坐标轴x2的直线与D的交点。首先,注意到

w(P)=w(P1)+P1Pwx1dx1=w(P2)-PP2wx1dx1

所以

|w(P)|=12[w(P1)+w(P2)]+P1P2wx1dx1

类似地

|w(Q)|=12[w(Q1)+w(Q2)]+Q1Q2wx1dx1

式(8)和式(9),有

D|w|2dA12D|w|dl+Dw,αdA2

利用不等式

(a+b)2a2+b2,    a,b>0

及Ho¨lder不等式, 可得

D|w|2dA12D|w|dl2+|D|Dw,αw,αdA

其中,|D|表示区域D的面积。现取C1=12max{1,|D|},引理1得证。

接下来,定义能量表达式,然后利用微分不等式技术推导关于此能量表达式的一阶微分不等式,得到解的二择一结果。

式(1),可得恒等式

0tz0zDξu[uη-(ρ(x,u,u)u,i),i+f(u)]dAdξdη=0

其中,z0(0)为坐标轴x3上的点,分部积分,可得

12z0zDξu2dAdξ|η=t+0tz0zDξρuh(u)dldξdη+0tz0zDξρu,iu,idAdξdη-0tDzρuu,3dAdη+0tDz0ρuu,3dAdη+0tz0zDξuf(u)dAdξdη=0

E(z, t)=0tDzρuu,3dAdη

式(10),可得

E(z, t)-E(z0, t)=12z0zDξu2dAdξ|η=t+0tz0zDξρu,iu,idAdξdη+0tz0zDξuf(u)dAdξdη+0tz0zDξρuh(u)dldξdη

求导,可得

zE(z, t)=12Dzu2dA|η=t+0tDzρu,iu,idAdη+0tDzuf(u)dAdη+0tDzρuh(u)dldη

利用Ho¨lder不等式和式(6),由式(11),可得

|E(z, t)|0tDzρu,32dAdη120tDzρu2dAdη12M10tDzρu,32dAdη120tDzu2dAdη12+M20tDzρu,32dAdη12×0tDz(ρ|u|2)1qu2dAdη12=I1+I2

由引理1和Schwarz不等式,可得

|I1|M1C10tDz|u|dldη2+0tDzu,αu,αdAdη120tDzρu,32dAdη12M1C10tDz|u|dldη0tDzρu,32dAdη12 +M1C10tDzu,αu,αdAdη120tDzρu,32dAdη12M1C1(Lt)2p-12p0tDz|u|2pdldη12p×0tDzρu,32dAdη12 +M1C1c1×0tDzρu,αu,αdAdη120tDzρu,32dAdη12M1C1(Lt)2p-12pp+1k2-12p×1c1p+12p0tDzρuh(u)dldηp+12p+p0tDzρu,32dAdηp+12p +M1C12c10tDzρu,iu,idAdη,
|I2|12M20tDz(ρ|u|2)1qu2dAdη+12M20tDzρu,32dAdη12M20tDzρu,32dAdη+0tDzρ|u|2dAdη1q0tDz|u|2qq-1dAdηq-1q12M20tDzρu,32dAdη+12qM20tDzρ|u|2dAdη+q-12k1qM20tDzuf(u)dAdη,

其中,LD的周长。将式(15)和式(16)代入式(14),再由式(13),可得

|E(z, t)|m1zE(z, t)+m2zE(z, t)p+12p

其中,

m1=max12M2,q-12k1qM2,M1C12c1 m2=M1C1(Lt)2p-12pp+1k2-12pmax1c1p+12p,p

2 二择一结果

分以下3种情形分析式(17)。

(i) 12<p<1

(A) 假设存在一个z0,使得Ez0t)>0。由于zE(z, t)0,所以对所有的zz0>0,有Ez,t)>0。此时对式(17)右边第1项,利用Young不等式,可得

m1zE(z, t)m1zE(z, t)p+12p×3p-1p+1×zE(z, t)p+14p×2-2pp+1m13p-1p+1zE(z, t)p+12p+m12-2pp+1zE(z, t)p+14p

所以

E(z, t) m3zE(z, t)p+14p+m4zE(z, t)p+12p

其中, m3=m13p-1p+1, m4=m2+m13p-1p+1

于是由式(19), 可得

E(z, t)+m324m4m4zE(z, t)p+14p+m32m42

化简,可得

Ez1m4E(z,t)+m324m42-m32m44pp+1

将右边的项移到左边并从z0z积分,有

m411m4E(z,t)+m324m42-m32m44pp+1×d1m4E(z,t)+m324m42z-z0

所以

2m4z0z1m4E(z,t)+m324m42-m32m41-3pp+1×d1m4E(z,t)+m324m42+m3z0z1m4E(z,t)+m324m42-m32m4-4pp+1×d1m4E(z,t)+m324m42z-z0

积分并舍弃左边一些非正项,可得

p+11-pm4E(z,t)m4+m324m42-m32m42-2pp+1(z-z0)-p+13p-1×m3E(z0,t)m4+m324m42-m32m41-3pp+1

利用不等式

a+ba+b,    a,b>0

式(22)可改写为

E(z,t)m5z-z0-p+13p-1×m3E(z0,t)m4+m324m42-m32m41-3pp+11+p1-p
m5=1-p1+p1+p1-pm4-2p1-p
limz[z-1+p1-pE(z,t)]m5

表明如果E(z,t)在某点z0处大于零,那么对任意的z>0 E(z,t)无界,且其增长速度等于或快于z-1+p1-p。将式(24)代入式(20),可得

limzz-2p1-pEzm5m42p1-p

(B) 若对任意的z>0,有E(z,t)<0,则式(17)可改写为

-E(z,t)m1Ez+m2Ezp+12p

对右边第2项,应用Young不等式,可得

m2Ezp+12p=m2δEz3p-12pEz2δ-3p-11-p1-p2p3p-11-pm2δEz+1-p2pm2δ-3p-11-pEz2

其中,δ为大于零的任意常数,代入式(26),可得

-E(z,t)m6Ez+m7Ez2

其中,m6=m1+3p-12pm2δm7=1-p2pm2δ-3p-11-p,解得

Ez-E(z,t)m7+m624m72-m62m7

积分,可得

2m7- E(z,t)m7 + m624m72-m6+m6ln- E(z,t)m7 + m624m72 - m62m7z0 -z

所以

-E(z,t)m7Φ2(0,t)e-2zm6+m6Φ(0,t)e-zm6

其中,

Φ(0,t) =-E(0,t)m7+m624m72-m62m7e2m7m6-E(0,t)m7+m624m72

显然limz-E(z,t)=0。对式(13)从0到积分,易得

-E(z,t)=12Rzu2dAdξη=t+0tRzρu,iu,idAdξdη+0tRzuf(u)dAdξdη+0tzDξuh(u)dldξdη

将二择一结果总结为以下定理。

定理1 假设u是式(1)~式(4)的解,且12<p<1,则要么

limzz-2p1-pEzm5m42p1-p

要么-E(z,t)m7Φ2(0,t)e-2zm6+m6Φ(0,t)e-zm6

Ez的定义见式(13),-E(z,t)的定义见式(31),mi(i=4,5,6,7)为大于零的常数。

注1m6的定义,知m6和常数δ有关, 只要选择δ足够小,衰减率1m6可足够大。

(ii) p=1

式(17)中,取p=1,并记m8=m1+m2,有

E(z,t)m8Ez

(A) 假设存在一个z0,使得E(z0,t)>0由于E(z,t)z0,所以对任意的zz0,有E(z,t)>0所以

E(z,t)m8Ez

z0z积分,可得

E(z,t)E(z0,t)ez-z0m8

(B) 若对任意的z>0,E(z,t)<0,则

-E(z,t)m8Ez,

从0到z积分,可得

-E(z,t)-E(0,t)e-zm8

式(13)和式(31),可得到以下定理。

定理2 假设u是式(1)~式(4)的解,且p=1,则当z沿x3轴趋近于时,u要么呈指数式增长,要么呈指数式衰减,即要么式(33)成立,要么式(35)成立。

(iii) p>1

(A) 假设存在一个z0,使得E(z0,t)>0由于E(z,t)z0,所以对所有的zz0,有E(z,t)>0式(17)右边第1项,应用Young不等式,可得

m2Ezp+12p=m2Ezp-12pEz1pp-1pm2Ez12+1pm2Ez

所以

E(z,t)m9Ez12+m10Ez

其中,m9=p-1pm2,m10=m1+1pm2

采用与式(18)~式(25)类似的方法,可得

E(z,t)m11e2z

其中,

m11=m10m9m10E(z0,t)m10+m924m102-m92m10×em92m10-E(z0,t)m10+m924m102-m92m10-2z0

(B) 若对任意的zz0,有E(z,t)<0,则式(17)可改写为

-E(z,t)m1Ez+m2Ezp+12p

应用Young不等式, 可得

m1Ezp+12p=m1Ezp+1p×p-1p+1Ezp+12p×2p+1p-1p+1m1Ezp+1p+2p+1m1Ezp+12p

于是式(39)可写为

-E(z,t)m12Ezp+12p+m13Ezp+1p

其中,m12=m2+2p+1m1,m13=p-1p+1m1

采用与情形(i)和情形(ii)类似的方法,可得

-E(z,t)m13Q0+p-1m12(p+1)z2(1+p)1-p+m12Q0+p-1m12(p+1)z1+p1-p

不难发现

limzz-2p1-p[-E(z,t)]p-1(1+p)m121+p1-p=m14

定理3 假设u是式(1)~式(4)的解,且p>1,则当z沿x3轴趋近于时,u要么呈指数式增长,要么呈指数式衰减,即要么式(33)成立,要么式(35)成立。

显然,在衰减情形下,Φ(0)Q0均与-E(0,t)有关。因此,要使定理1~定理3都有意义,须证明-E(0,t)有界。

3 全能量估计

为在衰减情形下,推导全能量E(0,t)的上界,首先假设

uf(u)αf(u)q1,    uh(u)βh(u)p1,    α,β>0

只要取q1=2pp+1>1p1=2p2p-1>1以及α,β足够小,此条件与式(5)和式(7)并不矛盾,在式(11)中,取z=0,可得

-E(0,t)=-0tD0ρu,3udAdη

类似地,在式(31)中取z=0,可得

-E(0,t)=12Ru2dAdξη=t+0tRρu,iu,idAdξdη+0tRuf(u)dAdξdη+0t0Dξuh(u)dldξdη

引入辅助函数

S(x1,x2,x3,t)=g(x1,x2,t)e-σx3

其中,σ为大于零的任意常数。对式(45)分部积分,由式(1),有

-E(0,t)=-0tD0ρu,3SdAdη=0t0Dξ(ρu,3S),idAdξdη=0tR(ρu,i),iSdAdξdη+0t0DξρSh(u)dldξdη+0tRρu,iS,idAdξdη=RuSdAdξη=t+0tRuStdAdξdη+0t0DξρSh(u)dldξdη+0tR f(u)SdAdξdη+0tRρu,iS,idAdξ

利用Ho¨der不等式和算术几何平均不等式,可得

RuSdAdξη=t14Ru2dAdξη=t+RS2dAdξη=t
0tRρu,iS,idAdξ120tRρu,iu,idAdξ+120tRρS,iS,idAdξ

再由式(5),可得

0tRuStdAdξdη0t Ru2qq-1dAdξdηq-12q×0tRSt2qq+1dAdξdηq+12qk1-q-12q0tRuf(u)dAdξdηq-12q×0tRSt2qq+1dAdξdηq+12qq-12qk1-q-12qε10tRuf(u)dAdξdη+q+12qk1-q-12qε1-q-1q+10tRSt2qq+1dAdξdη

其中,ε1为大于零的任意常数。再由式(44),可得

0t0DξρSh(u)dldξdη0t0Dξh(u)p1dldξdη1p1×0t0Dξ(ρS)p1p1-1dldξdηp1-1p1β-1p10t0Dξuh(u)dldξdη1p1×0t0Dξ(ρS)p1p1-1dldξdηp1-1p1β-1p11p1ε20t0Dξuf(u)dldξdη+β-1p1p1-1p1ε2-1p1-10t0Dξ(ρS)p1p1-1dldξdη

其中,ε2为大于零的任意常数。

类似地,可得

0tRf(u)SdAdξdηα-1q1ε3q10tRuf(u)dAdξdη+α-1q1q1-1q1ε3-1q1-10tRSq1q1-1dAdξdη

其中,ε3为大于零的任意常数。取适当的ε1ε2ε3使得

q-12qk1-q-12qε1+α-1q11q1ε3=12,    β-1p11p1ε2=12

将式(48)~式(52)代入式(47),再由式(46),可得

-E(0,t)-12E(0,t)+F(t)

其中,

F(t)=RS2dAdξη=t+120tRρS,iS,idAdξdη+q+12qk1-q-12qε1-q-1q+10tRSt2qq+1dAdξdη+β-1p1p1-1p1ε2-1p1-10t0Dξ(ρS)p1p1-1dldξdη+α-1q1q1-1q1ε3-1q1-10tRSq1q1-1dAdξdη

式(53),可得

-E(0,t)2F(t)

再由Φ0Q0的定义,易得其上界。

4 总 结

首先推导了引理1,假设式(1)~式(4)在柱体侧面满足非齐次Dirichlet条件,利用引理1并运用微分不等式技术,得到了解的二择一结果。事实上,经必要的修正,本文的结论也适合非齐次Neumann边界条件。显然,本文结论是对文献[13]的推广,研究方法和结果可为进一步研究复杂模型提供借鉴。

今后,可继续进行更深层次的研究,如研究文献[16]中的模型

(h(u))t=[(up+1)u,i]+k(t)b(t)f(t)

以及文献[17]中变系数热量方程

ut=Δum-V(x)um+up,m,    p>1

等。

http://dx.doi.org/10.3785/j.issn.1008-9497.2022.06.003

参考文献

A J C B de SAINT-VENANT.

Mémoiresur la flexion des prismes

[J]. Journal de Mathematiques Pures et Appliquees, 185612):89-189.

[本文引用: 1]

LIU YLI Y FLIN Y Wet al.

Spatial decay bounds for the channel flow of the Boussinesq equations

[J]. Journal of Mathematical Analysis and Applications, 20113811): 87-109. DOI:10.1016/j.jmaa.2011.02.066

[本文引用: 1]

LIU YQIU HLIN C H.

Spatial decay bounds of solutions to the Navier-Stokes equations for transient compressible viscous flow

[J]. Journal of the Korean Mathematical Society, 2011486): 1153-1170. DOI:10.4134/JKMS.2011.48.6.1153

LI Y FLIN C H.

Spatial decay for solutions to 2-D Boussinesq system with variable thermal diffusivity

[J]. Acta Applicandae Mathematicae, 20181541): 111-130. DOI:10.1007/s10440-017-0136-z

AMES K APAYNE L ESCHAEFER P W.

Spatial decay estimates in time-dependent Stokes flow

[J]. SIAM Journal on Mathematical Analysis,1993246): 1395-1413. DOI:10.1137/0524081

[本文引用: 1]

KNOPS R JQUINTANILLA R.

Spatial decay in transient heat conduction for general elongated regions

[J].Quarterly of Applied Mathematics,201776611-625. DOI:10.1090/qam/1497

QUINTANILLA R.

Some remarks on the fast spatial growth/decay in exterior regions

[J]. Zeitschrift für angewandte Mathematik und Physik, 20197083. DOI: 10.1007/s00033-019-1127-x

[本文引用: 1]

HORGAN C OPAYNE L E.

Phragmén-Lindelöf type results for harmonic functions with nonlinear boundary conditions

[J]. Archive for Rational Mechanics & Analysis, 1993122123-144. DOI:10.1007/BF00378164

[本文引用: 5]

LIN CPAYNE L E.

Phragmén-Lindelöf alternative for a class of quasilinear second order parabolic problems

[J]. Differential and Integral Equations,199583): 539-551.

LIN CPAYNE L E.

A Phragmén-Lindelöf type results for second order quasilinear parabolic equation in R 2

[J]. Zeitschrift für angewandte Mathematik und Physik, 199445294-311. DOI:10.1007/BF00943507

LESEDUARTE M CQUINTANILLA R.

Phragmén-Lindelöf alternative for the Laplace equation with dynamic boundary conditions

[J]. Journal of Applied Analysis and Computation, 201774):1323-1335. DOI: 10.11948/2017081

[本文引用: 3]

PAYNE L ESCHAEFER P W.

Some Phragmén-Lindelöf type results for the biharmonic equation

[J]. Zeitschrift für angewandte Mathematik und Physik, 199445414-432. DOI:10.1007/BF00945929

LIN C H.

A Phragmén-Lindelöf alternative for a class of second order quasilinear equations in R 3

[J]. Acta Mathematica Scientia, 1996162): 181-191. DOI:10.1016/S0252-9602(17)30793-2

[本文引用: 2]

LIU YLIN C H.

Phragmén-Lindelöf type alternative results for the stokes flow equation

[J]. Mathematical Inequalities & Applications, 200694): 671-694. doi:10.7153/mia-09-60

[本文引用: 1]

李远飞.

海洋动力学中二维黏性原始方程组解对热源的收敛性

[J].应用数学和力学,2020413): 339-352. DOI:10.21656/1000-0887.400176

LI Y F.

Convergence results on heat source for 2D viscous primitive equations of ocean dynamics

[J]. Applied Mathematics and Mechanics, 2020413): 339-352. DOI:10.21656/1000-0887.400176

李远飞.

Robin边界条件下更一般化的非线性抛物问题全局解的存在性

[J]. 应用数学学报, 2018412): 257-267. DOI:10.12387/C2018020

[本文引用: 1]

LI Y F.

Blow-up and global existence of the solution to some more general nonlinear parabolic problems with Robin boundary conditions

[J]. Acta Mathematicae Applicatae Sinica, 2018412): 257-267. DOI:10.12387/C2018020

[本文引用: 1]

YANG C XZHAO L ZZHENG S N.

The critical Fujita exponent for the fast diffusion equation with potential

[J]. Journal of Mathematical Analysis and Applications, 20133982): 879-885. DOI:10. 1016/j.jmaa.2012.09.050

[本文引用: 1]

/