数学建模集训知识大纲

杨利霞

数学建模集训知识大纲
* H7 J' w) d  \数学建模集训知识大纲

文章目录
/ x  Q; Q9 I% s) u- }5 C8 ^
数学建模集训知识大纲
评价算法
简单加权法
逼近于理想解的排序法(TOPSIS算法)
层次分析法
主成分分析法
5 \3 x, q' v) {: K+ W0 _, ^! E: W( c模糊综合评价法
4 X# h! J7 ~& W; m2 X6 ^聚类分析法
' L) |) C& b) l' W$ \秩和比法
7 F  Q7 m/ |+ ~1 H- Q9 ]4 P% b人工神经网络
& O$ J6 n! |- x) o! k熵权法
灰色关联度分析
预测算法
9 W" }/ E2 K7 f$ E- A插值拟合
% V* E6 }( T. `" t" ~( L) \回归模型预测
灰色预测GM
! c1 d+ m7 A/ m" j时间序列
神经网络
* B0 ?2 y! P* d6 Z; J- Q统计分析
方差分析
1 }, Z3 p5 M# S  i回归分析
多元统计分析
分类问题
假设检验（非参数）
聚类
最优化方法
6 `% s. [( U# J% V( ?7 H* G常用算法
无约束优化
有约束优化
二次规划
一般有约束非线性规划
计算机算法
动态规划
6 G$ K' N! D! i回溯搜索
分治
1 J) h+ c; j0 ]. M贪心
分枝定界
' U6 e. K* _& a: O! D+ {# S( `图论与网络
# A, I6 S8 C! l8 t* _$ Y欧拉图
二分图
4 U: Z( J. l" R网络流
最短路径
8 P) v2 n8 [3 s; g+ \* ~  J( F4 oTSP问题
& S- _9 i" f4 n- N4 UP和NP问题
背包问题
. E5 G5 y# Q: q1 x+ w& w最小生成树
; O8 K1 p7 k1 F$ D# J图的基本概念及其矩阵表示
; `( M/ s7 m" n; e5 P- ^常见的网络优化问题
图像处理
基本概念
数据类
图像类型
7 @* Q: g9 h( h3 E  Z数据类之间转换
, S1 M2 B0 Z9 T, L, w' U6 u( z5 Q$ W图像类型之间转换
2 E! B! |5 g7 S& T空间滤波器
频域变换
数字图像的水印防伪
( t5 ~) ?7 V( S0 |& b规划问题
; y. d. C  |) y. g6 D连续优化
离散优化
模拟和仿真
0 _* F0 H  a& R4 R模拟
" B1 G* u/ Q$ Y# r; D仿真
数值计算
/ W9 N* ?% L: D7 I常微分方程数值解
工具箱的使用
4 u4 B1 q, Z2 X4 }2 N8 Y0 L统计工具箱
. W) I  Y7 O# ~: e  R/ j2 n6 B拟合工具箱
神经网络工具箱
6 Z& T- F7 R# P# ?$ Y' a小波工具箱
! [3 k) M: k; A- Z并行计算工具箱
1 r8 f$ Z# `% k8 M5 @6 w优化工具箱
# J4 H+ ~7 a" f! d* o全局优化工具箱
评价算法

( d! b& y4 l9 [# A! E% h: \简单加权法
- M( \' i( d" d6 W* \/ P7 y8 o
" L6 U, Y) B! H5 F2 t- [' e线性加权
非线性加权
逼近于理想解的排序法(TOPSIS算法)
( X7 l/ r. g, x6 Z6 Q5 ]8 H  a/ |
层次分析法

5 X- A3 }( L. p" _" X$ F+ W3 Y1 z主成分分析法
9 g$ z8 i* e4 l# [
模糊综合评价法

聚类分析法
6 {0 k5 C9 `4 o* X7 K* \
% x7 s5 A. j8 [7 O- S& G0 X: v秩和比法
/ h- n/ _( m+ O% o: Q
- v7 b+ V' ^$ i2 q& y9 p人工神经网络
0 ~+ R# ~: m3 v7 Z( E$ y2 D
1 s! S  f6 ?5 A& o4 Y' i5 _/ k4 S熵权法
1 ?1 g5 M2 H" s4 p6 K
灰色关联度分析

预测算法

/ i( V; c- `/ q0 M插值拟合
5 W4 F' _1 S  l( D
' i, t$ Y( L% p小样本内部预测
回归模型预测

大样本内部预测
灰色预测GM
/ [0 Y& R, p/ F1 T" `; R5 C
小样本未来预测
+ X/ i& g+ c. C3 h4 y" u9 o0 _% G3 N2 P时间序列
1 e& s" e6 o0 y0 i4 ?
大样本的随机因素或周期特征的未来预测
神经网络

针对大样本的内部机理复杂的数据的未来预测
统计分析

' [2 j; W6 t+ z1 K  q方差分析

分析因变量的总误差中，除开随机误差以外，是否有类别变量（自变量）造成的处理误差，有多少误差是自变量造成的
  L( X# N& L( e! x6 ]分类
. x; P. ~' s. ^0 z5 q) ]% |单自变量
单因素方差分析
3 q% K7 ~8 t9 s- O" V8 U% ]两个自变量
无重复双因素分析：只考虑主效应，不考虑交互效应
可重复双因素分析：考虑主效应，也考虑交互效应
回归分析
7 Y- T1 k$ v% {
多元统计分析

6 [0 D: d( H( Q+ y& G6 p分类
聚类分析
判别分析
综合评价
2 _4 f! r9 v# ]+ \1 W, x# c主成分分析
从原来的坐标系转换到新的坐标系，第一个新坐标轴选择的是原始数据中方差最大的方向，第二个坐标轴选择的是和第一个新坐标轴正交且具有最大方差的方向。该过程一直重复，重复次数为原始数据中特征的数目。
因子分析
' }: t2 n9 q6 `& h, Q2 s7 ~* P隐变量和某些噪声的组合
典型相关分析
; Q$ h2 E- D6 ^* S( o偏最小二乘回归
分类问题

神经网络
. \& @* f2 X6 ]0 v- ?逻辑回归
2 `: k* N+ ^9 v: e- j1 C4 Z判别分析
最邻近方法
( G) [3 O/ A. E3 u2 c- h0 E& Q朴素贝叶斯
支持向量机
8 j7 x. [" g% I5 f2 e决策树
集成学习
! o* ?. q+ A4 k; |ROC曲线
- G$ U: A! H- I" d假设检验（非参数）
8 k1 v- F$ W2 _! U0 b( ^1 H  X# e
分布拟合检验
秩和检验
如果两个样本来自两个独立的但非正态或形态不清的两总体，要检验两样本之间的差异是否显著，不应运用参数检验中的t检验，而需要采用之和检验
% o( m- ~9 r9 R配对检验
K-S 检验
K-S检验不仅能够检验单个总体是否服从某一理论分布，还能够检验两总体分布是否存在显著差异。其原假设是：两组独立样本来自的两总体的分布无显著差异。
Q-Q图
用变量数据分布的分位数与所指定分布的分位数之间的关系曲线来进行检验的
: _; e2 P# S- Q7 {# ^# g! @P-P图
" M; S, x! ]* J- L根据变量的累积比例与指定分布的累积比例之间的关系所绘制的图形。通过P-P图可以检验数据是否符合指定的分布。当数据符合指定分布时，P-P图中各点近似呈一条直线。
聚类
2 N2 V" {2 T5 r5 a
' v% s3 x) J: V' o1 e# y: S层次聚类
划分聚类
: V0 y) \1 a- i+ WK-means
网络聚类
聚类评价
1 S$ t$ A% X" \( \3 h共表型相关系数
! \* U+ [3 T/ s. K/ \相似性矩阵
) @5 V( Q: q2 F* e& B. W0 G共表矩阵
" n6 H& p$ T! _2 u/ q& n+ F外部指标
# a$ g3 N2 j# S, x6 n6 Q相仪表
" j9 R" a( A4 w( ~* Z7 U& E* R' h调整兰德系数ARI
: Q1 F% ]1 J9 Y! x互信息MI
自助
8 d% I/ b2 _1 |# J0 l) }! \# p% ]最优化方法
1 a8 \* Y1 ^, X3 s0 N) t+ d4 S
; x' z4 V5 F- w4 ~- Q常用算法

; k- f( {1 n0 w* z' l0 ]模拟退火
神经网络
. `8 ~( \" ]8 d8 ?+ e# I遗传算法
2 y# K( K5 ]/ L. Y) e无约束优化

基本算法
共轭梯度法
8 @5 u9 h6 Q) D在最优化方法中占有重要地位，最速下降法的优点是工作量小，缺点是收敛慢，适用于寻优过程中前期迭代或作为间插步骤，当接近极值点时，宜选用别种收敛快的算法
4 R* W5 |% e& m- A. X; l牛顿法
, @/ c7 B; _) o5 A) R- W如果f是对称正定矩阵A的二次函数，用牛顿法经过一次迭代就可到达最优点，如不是二次函数，则牛顿法不能一部达到极值点
牛顿法收敛速度虽然快，但要求Hessian矩阵可逆，要计算二阶导数和逆矩阵，就加大了计算机计算量和存储量。
拟牛顿法
. R) _4 i# ?0 [, @: \修改了牛顿方向
两种算法
DFP
/ m( V; x2 e" o' K. m  r( J3 ^- s3 JBFGS
8 d" s" O+ S/ v. W! k7 M3 t有约束优化
( }9 p& e7 i5 R9 ^, ]' B4 v" f
! f& f4 S' [6 X2 _罚函数法
' X5 k9 f' e  T- ?6 E: g通过构造罚函数把约束问题转化为一系列无约束最优化问题，进而用无约束最优化方法求解，这类方法称为序列无约束最小化方法，简称SUMT法（分为SUMT外点法 其二为SUMT内点法）
近似规划法
6 y6 N1 c' A' a二次规划
2 n. I  D4 {2 i0 @" K  [; K, L7 o; n
quadprog
一般有约束非线性规划
/ @' a6 ]4 |7 T" Z2 l) ]
fmincon
计算机算法

& |+ u1 ~, a  Q$ ~; X" O* a动态规划

回溯搜索

+ m! S: b/ w. a7 S, H6 E7 T$ [, U分治

! N# V3 C7 s5 |/ O, M- s  N  v贪心
# X( o3 r$ A: a
  L# ]) b5 j# P& `! `2 t分枝定界

图论与网络

欧拉图
: w+ H" e, u" m) R+ T7 m6 [
二分图
* g$ a* G9 d9 F' B, j
2 V9 P3 F% z% ~+ |, A; F! |  o匈牙利算法
网络流

最小费用最大流
最短路径
4 v$ y( i( l* D+ r
8 l* ^0 Q# P6 v$ [8 i& i1 q3 dTSP问题
) a( h! I) e) g/ e/ R8 j% b
, W" i) Q6 I0 g2 u$ N) z近似算法和启发式算法
二边逐次修正法
  f6 M% k$ \3 A# ^( Z$ aP和NP问题
9 p3 o0 f# O* L# G. R) o% J
1 s# i4 V# a4 @/ N2 i背包问题

. F$ @3 O! N# J4 a# {0 a5 O" Y! M% V% o最小生成树
; L+ Q! ]& l1 T* o& i/ K# P' p
Prime算法
Kruskal算法
. w0 h' Y4 M* V) i. X图的基本概念及其矩阵表示

; M* {9 m+ c7 K2 |- Q无向图
2 u! S5 d% t; I5 m有向图
完全图
  Z; @! ~  P/ g( j, w1 `二分图
图与网络的数据结构表示方法
邻接矩阵表示法
; r1 G) m" ?% \4 w& d0 }关联矩阵表示法
# ]9 N- V. g6 T9 G弧表表示法
! B8 S+ j) e! k! B: d& a+ J% e7 I2 H邻接表表示法
1 z! e( |. o" B! Z& y5 @星形表示法
常见的网络优化问题
1 r1 U% m! _0 B
' t; T% J- y6 p6 K9 T7 _最短路问题
, }- j$ n' [4 w2 v; v公路连接问题
中国邮递员问题（欧拉图）
) {% x4 @% ^, R+ a旅行商问题TSP（哈密顿图）
# Y: f: }1 ^: n! v3 q运输问题
* Z/ W9 S, ~' p8 o$ t! O( S图像处理
9 j, b5 K/ \2 `/ H& _
基本概念
( L7 v8 P4 g" u, L0 l0 ^1 F
2 R; ]: Q5 a: R/ M) |图像分类
离散图像
5 U2 l9 q( [3 M/ G- w用一个数字序列表示的图像，0和1
- T" v& }/ P3 N- c) H连续图像
4 E& {3 i$ d0 C1 j二维坐标系中具有连续变化的图像
矩阵中的元素称为像素，以256灰色等级的数字图像为例，一般由8位，即一个字节表示灰度值
灰度值量化为对应灰度等级
, g. g0 Q; z$ i7 A等间隔量化（一般采用这个）
9 d' K" d# S/ A  l$ E非等间隔量化
- w# U" S$ ~3 C数据类

+ M3 a2 M6 ], r: y数值数据类
double
uint8
& B  O  t) t( r0 Y& Z+ [( t) l& ruint16
uint32
int8
; Y* z" K. s* x4 j' C' c" Eint16
int32
single
char
, _) j7 ]4 \+ ]7 u+ S3 ^字符类（逻辑数据类）
logical
图像类型

二值图像
二维矩阵由０、１构成，０为黑色，１为白色
通常用于OCR
一般二值图像是逻辑数组，只有０和１的uint8类数组，并不会认为是二值图像，需要使用logical函数 B= logical(A);
灰度图像
- _$ D, M  W. F  e即人们常说的256灰度图像，0表示纯黑色，255表示纯白色
* f, Q- b% P9 A0 q- ^4 ?二值图像可以看成是灰度图像的一个特例
索引图像
有两个分量。即数据矩阵X和彩色映射矩阵map，矩阵map是一个大小为m*3且由范围在[0,1]之间的浮点值构成的double数组，map数组的长度同它所定义的颜色数目相等，map数组的每一行都定义单色的红绿蓝三个分量，X则是索引矩阵
. a. ~& r2 j5 o& ?真彩色RGB图像
是彩色像素的一个mn3数组，其中每一个彩色像素点，都是在特定空间位置的彩色图像相对应的红、绿、蓝三个分量
* {5 I3 m) }% p$ I) [数据类之间转换
0 K1 c, E( b& b: }( A4 }7 Y# Q
im2uint8
im2uint16
mat2gray：将输入转换为double，范围为[0,1]
1 {1 Z. }* h$ p" d4 Him2double
im2bw
; b( r3 _2 J; M图像类型之间转换
8 `0 v! o' t/ ^5 D$ e+ y0 a
4 A& j0 Z4 Q0 _. y2 P5 C& Nind2fray
$ M- [+ T: r- z* f  _gray2ind
% e, r" G  J& n4 B# ]5 qrgb2ind
ind2rgb
ntec2rgb
5 v+ E8 \! C; M! `# R# srgb2ntsc
: N) N3 ]. ^6 J0 O' A( v使用imtool命令查看一个图像文件的信息
" W  O" a3 k; B0 X4 f+ {空间滤波器

: V2 d* f4 o% m8 R# q( V, i线性滤波器
使用拉普拉斯滤波器增强图像
) i! T2 t, [: r/ e% X. r使用fspecial生成过滤器以及imfilter的使用
非线性滤波器
  g% w( }8 E( @( M一个工具是ordfilt2函数，可以生成统计排序滤波器 g=ordfilt2(f,order,domain)
6 y8 [) E* |7 n) vordfilt2函数生成图像g的方式：使用邻域的一组排序元素中的第order个元素来代替f中的每个元素，该淋雨则由domain中的非零元素指定
' {1 L* |: r( P数字图像处理中最著名的统计排序滤波器是中值滤波器，对应第50个百分位，使用g=ordfilt2(f,median(1:m*n),ones(m,n));创建中值滤波器
工具箱提供了二维中值滤波函数g=medfilt2(f,[m,n])
$ c3 y! g1 D7 R' y) i: F% c  k. N2 {频域变换
* `+ r2 G# H3 O3 D
- s: N- i/ @7 N9 N  k为了有效的对图像进行处理和分析，需要将原定义的图像空间的图像以某种形式转换到频域空间，利用频域空间的特有性质方便的进行一定的加工，最后转换回图像空间
傅里叶变换
将图像从空域变换到频域
$ \) n( O: e# L$ B' W, Y$ x二维连续傅里叶变换
6 ]( `# U8 H4 B9 {) q二维离散傅里叶变换（DFT）
基于离散傅里叶变换的频域滤波
; w- y* U3 Q$ A0 f- O3 M离散余弦变换DCT
图像处理中常用的变换算法，通过DCT变换，可以将图像空间域上的信息变换到频率域上
* l( h+ @# V2 i$ J6 [两种实现
5 ^& `1 z" d  w0 v6 D基于快速傅里叶变换FFT的算法，通过工具箱提供的dct2
9 `; T* W- Q; w. o+ `另一种是DCT变换矩阵，工具箱提供了dctmtx函数来计算变换矩阵
( k. g& K& B8 i9 N2 X1 {2 t. v' W; m图像保真和质量
- V2 ?, `: |# x4 S  X7 q, j图像处理中为了增加压缩率有时会放弃图像细节或者其他不太重要的内容，为了衡量解码图像相对于原始图像的偏离程度，这些测度一般称为保真度准则
. B/ t9 V" R* P% d5 v9 h! n主要准则
客观保真度准则
6 ~0 R: u/ v+ b0 [4 }5 O当所损失的信息量可用编码输入图与解码输出图的函数表示时，可以认为是基于客观保真度准则的
均方根误差
$ q9 b3 L! x+ {. E- }" h; S* W' y' Z均方信噪比（SNR）
均方根误差越小，峰值信噪比越大，处理的图像质量越好
9 w# L3 Y( R% L  f) o0 a+ N# `- Y( J主观保真度准则
; [' o3 D$ F4 Y! V  N7 V用分数代表主观评价{很差，较差，相同，稍好，较好，很好}
数字图像的水印防伪
9 t. u6 |+ j# B& F6 r3 t2 T
9 p5 Z& E$ z- ~0 f* K内嵌水印的特点
透明性
) M0 P2 {1 H$ C; i& G% `1 `鲁棒性
- o6 a2 K& v( B% ^能够承受施加于图像的变换操作，不会因变换处理而丢失
安全性
$ _3 W* d4 ^4 p" P' R  E' Y空间域水印
' k/ I% o$ i# ?  E将水印信息嵌入到载体图像的空间域特性上
频率域水印
将水印信息嵌入到载体图像的变换域系数等特征上
4 g. t& A+ V/ F' g& [+ E  @) f基于矩阵奇异值分解的数字水印算法
( c/ p1 f* L2 I* S5 C奇异值分解定理
Weyl定理
在原矩阵上加一个小小的扰动，看是否矩阵奇异值的变化是否会超过扰动矩阵的最大奇异值，从而判断稳定性
水印嵌入
水印提取
, M' V* F* q5 V  @# Z1 ]嵌入的逆过程
基于DCT变换的水印算法
DCT变换是实数域变换，对实系数处理更加方便，不会使相位信息发生变化，另外，DCT变换是有损图像压缩JPEG的核心，基于DCT变换的图像水印将兼容JPEG图像压缩
水印嵌入算法
图像加密
图像隐藏
5 t9 i7 X4 K/ T/ i) w) T$ @: R$ r7 B规划问题

连续优化
3 ~  b- j8 Z6 J+ j0 Y3 ~0 B" U7 q
线性规划 LP
非线性规划 NLP
, {3 q3 D: Y! a5 T- W9 d二次规划 QP
离散优化
( W- }* K& y. ?* v3 h( r
$ o& e, O' g( e; g) p3 z& k# }整数线性规划
整数非线性规划
5 y& P9 X/ a* b% [纯整数规划
混合整数规划
- u0 D  x' x! |一般整数规划
+ r/ h. u; j, Q' @+ P; x0-1整数规划
整数规划的特殊情形，要求线性规划模型中的决策变量只能取值为0和1
模拟和仿真

模拟

4 ~4 Q6 W+ K0 p) d) z" W拟合
+ M2 ?* ?9 a  s( t
非线性最小二乘法
Isqcurvefit
6 g' M# M1 Y, F* ~8 X2 OIsqnonlin
拟合与统计回归
线性回归
非线性回归
[beta,t,J]=nlinfit(x,y,‘model’,beta0)
: ]7 O* W  l4 Q/ z2 a常用解法：线性最小二乘法
* w7 ]6 f, \+ G- T$ l; A( }$ na=polyfit(x,y,m)
1 b7 I1 b$ J& K4 A/ z' U9 @% X插值
/ g! ?  l, V& H% J; @
一维插值
: a1 z: ~1 i& U& K$ f3 x拉格朗日插值

yy=lagrange(x,y,xx)
分段插值

2 w# ?. L. T4 ~: S9 o( }+ Wyy=interp1(x , y ,xx , ‘method’)
! x. _) p. D7 f" m  N) X/ s三次样条插值
8 w: A! t* D" f( f  r  z$ L% G% A
# ]5 A! H5 W2 |2 v+ j9 g7 j* jyy=spline(x,y,xx)
二维插值

- a4 B8 T: ^" ]/ @# {最邻近插值
分片线性插值
; }2 y0 \5 x; M9 I1 Q双线性插值
- f! {; B/ n# A/ v! p, y; ^仿真

动态仿真
' a9 {7 c7 R7 }0 V8 d! r' ~0 ^) l连续系统仿真
7 A2 i% s* |8 j) {2 ]& F时间步长法
7 N3 G& g; B' h$ w# H2 c离散系统仿真
事件步长法
0 K$ D  u( T  h1 b3 z静态仿真
蒙特卡罗方法
随机数产生
均匀随机数
其他方法
8 L9 l, W7 o# Q逆变换法
, t. A% n1 x2 K- l. r- F舍选法
- l  W5 r4 F2 a+ r5 |- v  |近似抽样法
& y6 r/ c& `+ F; e% @9 c指数分布
泊松分布
标准正态分布
正态分布
3 ?7 v; O& a+ c2 J0 W# D/ ]0 B二项分布
排队问题
可靠性问题
8 g& m6 I1 s8 a# K蒙特卡罗方法
$ s* d. H( ]$ x+ C4 {' W3 u: S元胞自动机
数值计算

9 b4 b- M# c( g, X7 v常微分方程数值解
9 W$ y" G( B4 X
: n. k$ t7 C* b$ o2 |- c微分方程解析解
( m8 q8 Y" C) O5 F# w% _u=dsolve(‘Du=1+u^2’,‘t’)
1 ?0 O) C0 z9 L% H, \. I欧拉方法
6 r. P9 j7 i1 e8 L2 e向前欧拉公式
向后欧拉公式
梯形公式
龙格—库塔方法
matlab求解
9 F- L' p0 w  }% d; c1 I高阶微分方程转一阶常微分方程
, H. I- S7 Y& v9 ]$ E刚性常微分方程与非刚性常微分方程
3 c/ {, Q1 L" j( V) o工具箱的使用

* D& j  v  B$ _0 y/ N6 t) n统计工具箱
3 ?* w& `% R6 L5 g6 u, n
# A( Z" K: p- \6 _9 }5 H% {! {) i拟合工具箱
1 Q' p4 `0 F- X+ j) v+ o& P
神经网络工具箱

小波工具箱
& D; y% u, ?/ d2 N6 I
4 U3 I/ f# }- h3 I0 J' v2 J4 V并行计算工具箱

* t4 i) ?# K% N! w2 N8 l- {优化工具箱

一元函数极小值
fminbnd
+ {) W$ ?' D0 N# O无约束极小
! V* F% a! r' N* K" X9 tfminunc
为无约束优化提供了大型优化和中型优化算法
: G! O" m$ N/ J* ~% m6 L为中型优化算法的搜索方向提供了4种算法，由options中的HessUpdate控制，默认为bfgs（拟牛顿法的DFP公式） 另外，还有dfp（拟牛顿法的DFP公式），还有steepdesc（最速下降法）
& @9 u, w. y, k/ ^4 O为中型优化算法的步长一维搜索提供了两种算法，由options中的参数LineSearchType控制，默认是quadcubic （缺省值，混合的二次和三次插值）另外还有cubicpoly（三次多项式插值）
fminsearch
8 F; Q6 \9 R% S& ]' M5 I3 e线性规划
, T! J- f7 Q- X* O: M+ `( Vlinprog
二次规划
! l( p% C2 \; y+ B+ b; Gquadprog
& q  v& T6 ~" |9 T/ t( E" T约束极小
fmincon
达到目标问题
fgoalattain
, h: j, f. D4 d, A1 i- A0 M极小极大问题
fminimax
全局优化工具箱
————————————————
) J. G/ M# {9 r% J  L  n原文链接：https://blog.csdn.net/SinclairWang/article/details/100891757

$ P" P9 v6 A: V8 `  U6 i' s5 o

柠檬草lll

发表回复很全面

柠檬草lll

发表回复谢谢分享
2 q! }( H5 F% N, L0 T