数学建模集训知识大纲

杨利霞

数学建模集训知识大纲
数学建模集训知识大纲
8 a+ \5 W5 [" ]: I- g. @
文章目录
' p% o" z+ b& d
4 r5 G/ [9 F6 h" g5 d数学建模集训知识大纲
评价算法
! \. V" U- r6 C4 [5 ?! o简单加权法
逼近于理想解的排序法(TOPSIS算法)
" s! N& q$ ~; M! P层次分析法
1 O4 R; S, b$ [) t9 k' q0 E主成分分析法
8 |  v9 r- r. k( L  b1 y4 J模糊综合评价法
聚类分析法
6 b" `4 v  v- G. Z( ^秩和比法
5 r7 K. j: W( K7 r8 \7 g) j人工神经网络
熵权法
灰色关联度分析
- F( Q  e0 |. J! c) }预测算法
插值拟合
回归模型预测
2 p. W, O# H3 J: C* ?5 e灰色预测GM
时间序列
' D5 V, f! P& Z神经网络
7 d' L7 }; W0 _( |; K! L4 P统计分析
方差分析
回归分析
. y* V9 T, N- D0 B0 q" P& }- G多元统计分析
分类问题
假设检验（非参数）
聚类
6 n& L! f. N& r: C最优化方法
) ~. J! O; w5 \! r' D常用算法
' L( L. f% o" i无约束优化
有约束优化
二次规划
一般有约束非线性规划
* s( b, _) g  z9 Q计算机算法
# L: S3 ~$ c9 ]( x" G$ u* R% ^动态规划
回溯搜索
分治
5 F; q! {0 n: W5 b7 q贪心
分枝定界
4 C' c+ K) u0 u" M2 [3 R图论与网络
欧拉图
二分图
网络流
* ^7 S0 V% }  I- I最短路径
TSP问题
P和NP问题
背包问题
/ I3 f3 s7 `3 g) {最小生成树
图的基本概念及其矩阵表示
常见的网络优化问题
图像处理
  Y$ L; x" A0 j2 r7 B基本概念
, Y5 Y! c  F; ]0 e# r  T数据类
图像类型
数据类之间转换
4 X% k/ @" Y1 r1 ]图像类型之间转换
3 p) u+ E. S' [: r空间滤波器
频域变换
% w2 |6 \8 W  N5 l3 p  c数字图像的水印防伪
规划问题
, m. ]' [! O4 e5 ?, @  z! T  `连续优化
离散优化
8 R9 c+ f+ t9 Y7 a* J& m模拟和仿真
; O8 g6 H0 E3 u" c9 L0 j  C( N$ i模拟
: T1 K+ M: {; G# ^/ S仿真
* P2 t" g$ n# b0 e) W2 {数值计算
常微分方程数值解
工具箱的使用
2 p. `9 [% o8 F$ G: b/ m- R* T1 @统计工具箱
# u( s4 i. O& h. p0 u拟合工具箱
& b) t- n( `  ?, Z神经网络工具箱
5 f) C) F* V, |. d% [( y7 ]& a小波工具箱
并行计算工具箱
优化工具箱
全局优化工具箱
评价算法

简单加权法
( H8 g+ k  W! k- A$ O
% I; g. I3 H! P' y- o线性加权
( b9 p4 b, X: a& {0 D* T0 F非线性加权
逼近于理想解的排序法(TOPSIS算法)

" \! h3 X% Z" k0 @$ V9 \层次分析法

主成分分析法

模糊综合评价法

聚类分析法
: i. a& X: ?# y
秩和比法

( Z* C+ h% y, F8 W, P  m人工神经网络

$ u1 _4 s1 C# j3 n) \熵权法

灰色关联度分析

预测算法

插值拟合

  x4 i. }+ U4 c: Y7 u$ K小样本内部预测
0 Y; i8 |! C( o' Y回归模型预测

大样本内部预测
9 |( Y3 i7 Z6 |$ J7 B灰色预测GM
! x; X( S% f/ c9 T: [/ Y8 j) `
# }" n, x* \( H小样本未来预测
6 R6 C, B" {8 S6 l8 M. L时间序列

大样本的随机因素或周期特征的未来预测
; Q6 [, W9 |3 p神经网络
5 _6 n$ j9 W) l' ^
$ W  I7 A0 C" S9 E$ L5 D) C针对大样本的内部机理复杂的数据的未来预测
2 d+ V9 q) ^6 ]2 f统计分析
( }, \, O4 i% \* D6 \2 I
方差分析
: l& U8 ]1 M% I/ j. @6 R
( h$ f  v3 w2 Z分析因变量的总误差中，除开随机误差以外，是否有类别变量（自变量）造成的处理误差，有多少误差是自变量造成的
分类
8 [' j- A: X7 Z  C: d9 A3 ~单自变量
, J8 B3 L0 H( A" V* \/ F: k单因素方差分析
6 N3 s9 Z1 ?2 p4 j两个自变量
1 n& B& l6 |) i无重复双因素分析：只考虑主效应，不考虑交互效应
" N) S' _- S8 [1 h" W' k" |7 r' c可重复双因素分析：考虑主效应，也考虑交互效应
7 H7 A9 t- O. H) }回归分析
1 S. p( n/ V* o* X4 {
1 z5 A- [. N+ m1 \( @* a; ^多元统计分析
8 c: L7 ]& ?( P7 C5 O2 o
$ f1 o2 u3 |; X% h, y7 S( m; X分类
聚类分析
% X( p  Y/ I. |6 e  o判别分析
8 k7 D7 |4 K. J5 M  h综合评价
主成分分析
) i& ]( C$ U/ d+ I从原来的坐标系转换到新的坐标系，第一个新坐标轴选择的是原始数据中方差最大的方向，第二个坐标轴选择的是和第一个新坐标轴正交且具有最大方差的方向。该过程一直重复，重复次数为原始数据中特征的数目。
因子分析
隐变量和某些噪声的组合
! k8 @  k$ P  B3 V1 C$ J典型相关分析
偏最小二乘回归
分类问题

神经网络
逻辑回归
/ w4 u4 K# V6 Q9 @* k# |判别分析
最邻近方法
朴素贝叶斯
支持向量机
% U. y$ g( `3 P/ @& ^& u# `决策树
$ [' W# W' I/ X集成学习
5 u# `+ }5 R+ Y$ p# U. d" K" bROC曲线
假设检验（非参数）

: A1 ^) D8 H8 J/ F8 @4 q分布拟合检验
秩和检验
如果两个样本来自两个独立的但非正态或形态不清的两总体，要检验两样本之间的差异是否显著，不应运用参数检验中的t检验，而需要采用之和检验
配对检验
$ \- I. k2 |* F4 V* fK-S 检验
& X* v5 v$ U! |) J4 cK-S检验不仅能够检验单个总体是否服从某一理论分布，还能够检验两总体分布是否存在显著差异。其原假设是：两组独立样本来自的两总体的分布无显著差异。
1 I0 K8 t9 b$ `, oQ-Q图
用变量数据分布的分位数与所指定分布的分位数之间的关系曲线来进行检验的
- K& ^1 z! i. ?, \& k, jP-P图
根据变量的累积比例与指定分布的累积比例之间的关系所绘制的图形。通过P-P图可以检验数据是否符合指定的分布。当数据符合指定分布时，P-P图中各点近似呈一条直线。
, \" z6 L" X; r, U3 S. L聚类

1 K5 M. d- q6 Y' ^. U4 @层次聚类
: ^6 l$ x5 F  u2 S$ p. k3 k划分聚类
K-means
- W2 U( n' Q. g7 C" U1 B网络聚类
聚类评价
" q" O! d2 i2 ]0 p, R共表型相关系数
相似性矩阵
3 S! s/ C( x0 t  V+ t* c4 ?7 X- t$ H共表矩阵
外部指标
3 ^% M5 v  o2 C& f相仪表
2 O0 D" B" t$ r' D& S  O9 I" R3 k0 S/ F调整兰德系数ARI
互信息MI
自助
8 Z, T: f  v- p6 I最优化方法

' k0 h4 [) J8 i) I9 D( B常用算法

# _. n. _) J5 X1 Z7 J模拟退火
神经网络
遗传算法
无约束优化

2 c) A2 x4 b1 [" o& `/ `基本算法
共轭梯度法
在最优化方法中占有重要地位，最速下降法的优点是工作量小，缺点是收敛慢，适用于寻优过程中前期迭代或作为间插步骤，当接近极值点时，宜选用别种收敛快的算法
牛顿法
如果f是对称正定矩阵A的二次函数，用牛顿法经过一次迭代就可到达最优点，如不是二次函数，则牛顿法不能一部达到极值点
6 c+ q# B# }5 o6 l牛顿法收敛速度虽然快，但要求Hessian矩阵可逆，要计算二阶导数和逆矩阵，就加大了计算机计算量和存储量。
4 P3 R/ Z6 K6 C5 x拟牛顿法
9 @% j4 ]4 r7 j$ c( D6 V( ^修改了牛顿方向
# ^3 m/ P7 H2 w6 `! F. i/ v( f% f两种算法
DFP
. V3 Q5 q! m+ x1 A  E: BBFGS
/ g, n/ e; h' }% j- q) P# j/ O有约束优化
6 S5 \0 o6 ?" J& d$ c, T
* G+ M! M, w3 ?罚函数法
通过构造罚函数把约束问题转化为一系列无约束最优化问题，进而用无约束最优化方法求解，这类方法称为序列无约束最小化方法，简称SUMT法（分为SUMT外点法 其二为SUMT内点法）
7 B# q% V) v8 G5 _' H! k近似规划法
$ d7 @6 j' j* n8 t二次规划

% V6 x  j: R) ]1 V4 Iquadprog
- r. f6 D) e) Y, a4 e+ w' d# U$ a一般有约束非线性规划
' T. d( t( ]" z3 @( l4 q- m
fmincon
2 K, ?! }* l: s% J9 _( U计算机算法

3 _. x7 N$ m* j& _7 p$ T动态规划

回溯搜索

; H6 Y3 ~' _$ V: u$ ?6 D分治

* Y" x+ `8 A# U6 x  a+ ?: [$ P贪心

( t7 k, r4 [. k9 \) Y9 ]& y* `; k分枝定界
& @& b$ [9 B7 `. ]" C
+ f6 |& J. u1 b$ p9 Q: h4 Y3 G图论与网络
7 D9 Y, {" D- Z6 L5 r( E8 |
欧拉图
9 U6 c% x2 l4 e5 {, J) B% U8 @
' v- l1 Z+ S) z  {: t( r二分图
- @; [0 m7 C* v) W
- r" ]1 n) T5 [9 g& ?1 C5 ~# n, I匈牙利算法
网络流
# D: _, Z5 {, |
最小费用最大流
+ ?, I0 i$ f: h" S% M% @最短路径
+ _, D1 y+ [5 `. ^/ i, |" W
TSP问题

, \9 b8 e3 q) b' f+ k! p- V近似算法和启发式算法
二边逐次修正法
P和NP问题
9 S6 H. y- @: P0 X
; F" i5 ?3 ^% J9 [+ U7 f: |. Y背包问题

最小生成树

Prime算法
Kruskal算法
图的基本概念及其矩阵表示
8 p1 x$ M2 ^& c( v" n
# T9 Q2 f# V% y, l' P无向图
有向图
! u: ^3 U5 d- w3 v4 b7 M9 |( s0 b9 a完全图
/ _- x/ P( V7 B% A; V# Z二分图
: N2 |+ \" y3 u  K0 x4 C) _/ R图与网络的数据结构表示方法
邻接矩阵表示法
3 U; D) k4 W% f% {关联矩阵表示法
. Y" K, C4 G: T- E弧表表示法
. w8 \) T4 r5 b/ {! K+ ~邻接表表示法
星形表示法
! C8 @" ~8 V- Q7 J; W  ^* F常见的网络优化问题

3 A# p: x) [- z! h9 [7 K最短路问题
$ \7 h0 k& O# Y6 t7 c0 \2 C公路连接问题
5 h7 L5 b' q& n) V+ N1 \+ c中国邮递员问题（欧拉图）
旅行商问题TSP（哈密顿图）
2 h4 o- V2 D# C运输问题
图像处理
: E9 }, f' H6 ^
' F$ K' q" d& G# U. ^& @: {基本概念
) z- q! J" t% T
图像分类
离散图像
# b& o+ F  G/ z: c* W用一个数字序列表示的图像，0和1
, r  b" N' a6 I- i+ v* \0 e  T. d+ K连续图像
二维坐标系中具有连续变化的图像
$ ^2 {& o  A% U8 R/ w/ D矩阵中的元素称为像素，以256灰色等级的数字图像为例，一般由8位，即一个字节表示灰度值
! j. J4 f+ q/ m! r9 U+ u3 t+ o灰度值量化为对应灰度等级
2 U6 C. G) a# ^; L6 V等间隔量化（一般采用这个）
非等间隔量化
2 u* T) P0 o8 T. X7 {7 V& A  A' ]数据类

数值数据类
double
uint8
/ Z9 t4 L( \6 D& Y2 ?7 Buint16
  |  m8 H4 @4 ?1 G, ?7 Euint32
int8
int16
( s/ k5 r& R* T/ @5 Oint32
: P6 ?5 O# |  x  y1 f" ^- asingle
6 n( g( g4 q" P' Fchar
( B9 o5 k$ k% _2 t% ^* c字符类（逻辑数据类）
logical
图像类型
; [- B6 k! b5 a- `* U/ W1 b( H/ Y; b
二值图像
8 u$ \" G* \0 S# h6 j二维矩阵由０、１构成，０为黑色，１为白色
5 W1 }$ C, s# c3 x通常用于OCR
& J- U% M5 z2 K! ^一般二值图像是逻辑数组，只有０和１的uint8类数组，并不会认为是二值图像，需要使用logical函数 B= logical(A);
- T( s/ I6 u7 J+ x5 m- W8 v- {灰度图像
, I' M6 g0 n! ?, _. ~; Z即人们常说的256灰度图像，0表示纯黑色，255表示纯白色
) K. S. ?. b# W/ p5 O二值图像可以看成是灰度图像的一个特例
0 O0 m8 j/ |# j3 l8 ?索引图像
# h/ j. n  V* E+ K+ y' @* |有两个分量。即数据矩阵X和彩色映射矩阵map，矩阵map是一个大小为m*3且由范围在[0,1]之间的浮点值构成的double数组，map数组的长度同它所定义的颜色数目相等，map数组的每一行都定义单色的红绿蓝三个分量，X则是索引矩阵
真彩色RGB图像
3 x6 ^: I  m! X( M7 ?是彩色像素的一个mn3数组，其中每一个彩色像素点，都是在特定空间位置的彩色图像相对应的红、绿、蓝三个分量
. v+ A# V" O, c+ m0 M6 s# k6 l. Z数据类之间转换

im2uint8
) G* S8 C) w2 u, B" B, G( L, qim2uint16
mat2gray：将输入转换为double，范围为[0,1]
im2double
im2bw
) Y: G0 b# ^+ h7 L图像类型之间转换
9 l/ a! b) P0 Z: K# S% `
2 V9 |+ ?+ [# w" k% qind2fray
0 x6 f' H  w& _3 G$ hgray2ind
rgb2ind
2 R- B7 t$ R$ Y8 {. |( B* G: oind2rgb
/ ]  g9 Q, l/ ^4 Q1 Mntec2rgb
rgb2ntsc
4 A* k# h8 G/ K$ z. ]使用imtool命令查看一个图像文件的信息
3 L$ _% U, M! k/ H1 P空间滤波器
6 k  p+ |  y! M# x6 z
2 ~7 _; F! m* t* n' w( @线性滤波器
' A+ R; d0 H% x+ p- U使用拉普拉斯滤波器增强图像
+ I1 _4 @1 h3 z, k( w使用fspecial生成过滤器以及imfilter的使用
非线性滤波器
7 Q8 ^; k5 e7 C9 f2 E/ ]一个工具是ordfilt2函数，可以生成统计排序滤波器 g=ordfilt2(f,order,domain)
8 D+ h& k* Q+ g, i" Q2 {ordfilt2函数生成图像g的方式：使用邻域的一组排序元素中的第order个元素来代替f中的每个元素，该淋雨则由domain中的非零元素指定
' W& b4 A( ?& d* ]数字图像处理中最著名的统计排序滤波器是中值滤波器，对应第50个百分位，使用g=ordfilt2(f,median(1:m*n),ones(m,n));创建中值滤波器
工具箱提供了二维中值滤波函数g=medfilt2(f,[m,n])
& X. z' {1 N; j0 C0 _9 V频域变换

1 \" H- Y& W% T/ p& q8 V) \/ u为了有效的对图像进行处理和分析，需要将原定义的图像空间的图像以某种形式转换到频域空间，利用频域空间的特有性质方便的进行一定的加工，最后转换回图像空间
傅里叶变换
  \2 |' u* D4 X# v& t/ @将图像从空域变换到频域
二维连续傅里叶变换
4 o: e/ w8 H3 ^+ r) P: ?( P; d) ^7 A. K3 \二维离散傅里叶变换（DFT）
" D( X# s, ?# t) C2 P5 g基于离散傅里叶变换的频域滤波
离散余弦变换DCT
$ F, @4 r; k2 X图像处理中常用的变换算法，通过DCT变换，可以将图像空间域上的信息变换到频率域上
  f6 w3 ?* R+ h0 r8 m两种实现
9 [- D4 |1 W( Z) g; H3 L基于快速傅里叶变换FFT的算法，通过工具箱提供的dct2
另一种是DCT变换矩阵，工具箱提供了dctmtx函数来计算变换矩阵
图像保真和质量
图像处理中为了增加压缩率有时会放弃图像细节或者其他不太重要的内容，为了衡量解码图像相对于原始图像的偏离程度，这些测度一般称为保真度准则
. c% S' k5 e; B主要准则
客观保真度准则
当所损失的信息量可用编码输入图与解码输出图的函数表示时，可以认为是基于客观保真度准则的
均方根误差
均方信噪比（SNR）
均方根误差越小，峰值信噪比越大，处理的图像质量越好
主观保真度准则
& ~, o: h' i+ a8 c) z2 A5 S用分数代表主观评价{很差，较差，相同，稍好，较好，很好}
( Z! t+ Y2 l7 c: J9 S1 O数字图像的水印防伪

1 B5 Z( q$ w4 t+ T+ N9 q" i  g: q内嵌水印的特点
6 ~/ U3 k/ b* \/ \+ I透明性
' X- U6 M0 ^' Y( v鲁棒性
2 C4 H' K; z' v2 R1 {能够承受施加于图像的变换操作，不会因变换处理而丢失
安全性
1 x" k4 G" O! s" X& @# U9 I空间域水印
将水印信息嵌入到载体图像的空间域特性上
频率域水印
将水印信息嵌入到载体图像的变换域系数等特征上
8 U* I2 Q2 P3 z( B; P% h基于矩阵奇异值分解的数字水印算法
2 r& k7 R& o7 o& Y: M$ T6 \奇异值分解定理
3 z  b( p5 G  \0 xWeyl定理
在原矩阵上加一个小小的扰动，看是否矩阵奇异值的变化是否会超过扰动矩阵的最大奇异值，从而判断稳定性
* W* B8 Q. {) G4 k* p. S水印嵌入
水印提取
# c$ C! Y; f" u% z2 y- Z! V嵌入的逆过程
基于DCT变换的水印算法
DCT变换是实数域变换，对实系数处理更加方便，不会使相位信息发生变化，另外，DCT变换是有损图像压缩JPEG的核心，基于DCT变换的图像水印将兼容JPEG图像压缩
1 i+ ]& r; B* m% w" n水印嵌入算法
  F2 ^5 @" {2 r5 e- Y" c% L3 |) A图像加密
. Q$ Q2 X% G2 z2 ]图像隐藏
规划问题
1 s" Z  `. U1 H6 S! C9 C% I# T! C$ `
连续优化

线性规划 LP
非线性规划 NLP
二次规划 QP
离散优化

整数线性规划
% b  o$ M: W& n, P' G9 _整数非线性规划
+ E9 J( e5 ?# [- x纯整数规划
混合整数规划
1 @& k, P! t( B" [一般整数规划
0-1整数规划
2 `8 ?$ ~9 g- Q$ X8 y' O整数规划的特殊情形，要求线性规划模型中的决策变量只能取值为0和1
模拟和仿真

模拟

( [3 I" H. `7 B7 o0 S拟合
) S' B9 W4 z- J, c/ ~5 F
! M! l) r  z' m% I( L9 H非线性最小二乘法
Isqcurvefit
1 c' c9 ^( f6 J7 [) w+ S7 ?Isqnonlin
拟合与统计回归
; w8 d0 x$ [0 W0 u8 e- l+ S线性回归
非线性回归
; L1 z- k! o: L( O  C7 P6 R( N, j[beta,t,J]=nlinfit(x,y,‘model’,beta0)
( b9 j4 q. |' @常用解法：线性最小二乘法
2 [! ^) g3 \& f- ja=polyfit(x,y,m)
插值

2 i4 ?) m5 z& k3 W一维插值
9 S* {4 ]4 Q% a% p3 U  T拉格朗日插值
3 h( x+ X5 C* O# [( y
" ]9 d& k# C$ [+ J3 syy=lagrange(x,y,xx)
分段插值
5 C2 u- ^7 ^- J5 A
yy=interp1(x , y ,xx , ‘method’)
, a" s% t/ E& I  L三次样条插值

yy=spline(x,y,xx)
二维插值

最邻近插值
分片线性插值
6 X0 Z: q" D! g  ]. {. v3 k双线性插值
仿真
! D; F8 Z5 Y7 X+ X, f
动态仿真
, q- [' X# a3 j) l连续系统仿真
1 e) A/ I) |! u* r+ O4 o时间步长法
离散系统仿真
事件步长法
静态仿真
蒙特卡罗方法
随机数产生
! H7 y9 c. i9 b9 b均匀随机数
4 L  W4 F4 [! F% S- h( [其他方法
逆变换法
舍选法
近似抽样法
指数分布
泊松分布
标准正态分布
1 A+ n! f$ [* w# G2 c4 D% e正态分布
1 C9 w0 B2 r6 I7 o" ]. Q" p- h二项分布
排队问题
; A" `3 U$ o% X- A可靠性问题
- t4 h: o% Z& e9 E% g$ y. }$ ~蒙特卡罗方法
元胞自动机
数值计算

常微分方程数值解

$ D4 \) @4 X- l. ]* n' u$ V1 E# M微分方程解析解
8 Z; \9 j$ T# M0 c. d8 S# j& Gu=dsolve(‘Du=1+u^2’,‘t’)
; R8 ]' G7 G8 Z0 X) `" O: \! T欧拉方法
向前欧拉公式
  r  o, _5 |5 S/ A0 M! q: K向后欧拉公式
( a2 a: e# U$ S6 w5 m, [" I) \梯形公式
( d) J9 @1 |1 ?- X' _8 z/ k龙格—库塔方法
matlab求解
高阶微分方程转一阶常微分方程
3 E* r3 ?' F9 H& h% S3 ^刚性常微分方程与非刚性常微分方程
; ~/ N' ]5 [& r6 b9 [工具箱的使用
# g6 C+ I8 \8 t- P# d# Y1 ]
统计工具箱
4 d( k! c  g9 M1 C1 n
拟合工具箱
! x, f3 u/ }8 E2 S+ V6 q
, C- w; a0 [2 M神经网络工具箱
( F6 e+ L% g% K+ {
小波工具箱

并行计算工具箱
3 Z# f  H1 ^0 ~
% j8 @# J# E/ F% [0 N1 k" e7 t优化工具箱

, B! C1 F, `; ?5 v3 @0 Z4 H, U一元函数极小值
1 [0 M1 x# ^: M+ U! ifminbnd
# ?, L; b" g0 I& q* a8 l& Q4 [1 k无约束极小
fminunc
为无约束优化提供了大型优化和中型优化算法
为中型优化算法的搜索方向提供了4种算法，由options中的HessUpdate控制，默认为bfgs（拟牛顿法的DFP公式） 另外，还有dfp（拟牛顿法的DFP公式），还有steepdesc（最速下降法）
为中型优化算法的步长一维搜索提供了两种算法，由options中的参数LineSearchType控制，默认是quadcubic （缺省值，混合的二次和三次插值）另外还有cubicpoly（三次多项式插值）
2 j$ B, J+ ]% y- G7 }fminsearch
1 \9 Y! ^9 [+ I& p线性规划
; B& {( {# w! m* K% Jlinprog
二次规划
. Y, i5 V, L. g: t$ t& A( qquadprog
约束极小
, G1 N: `, Q4 ~0 {# s0 xfmincon
达到目标问题
fgoalattain
- t1 p: n' v2 z! X$ B* O极小极大问题
+ {" }+ ]5 L3 I6 mfminimax
/ J5 y" w$ j; Y8 \6 r' G全局优化工具箱
$ A. _0 u+ l: b9 B. }1 i( ]————————————————
原文链接：https://blog.csdn.net/SinclairWang/article/details/100891757
; j7 y- Y/ r0 w; \1 z2 `+ c5 w
; ]0 f9 @7 w8 T4 J1 }
1 @0 J/ z: k: X4 A* v

柠檬草lll

发表回复很全面

柠檬草lll

发表回复谢谢分享