数学建模集训知识大纲

杨利霞

数学建模集训知识大纲
2 J; z# q3 X" x" R6 @8 f数学建模集训知识大纲
# x4 S1 q, G& p# ?% Y0 u
文章目录

: [" X0 C. Z' n4 }: v: t数学建模集训知识大纲
9 S0 E' S! K) S  D. g2 P5 C评价算法
简单加权法
逼近于理想解的排序法(TOPSIS算法)
层次分析法
. _5 e6 N" f  y# q* K主成分分析法
模糊综合评价法
1 c1 ?- j0 Y3 i% {) c* ]* ~聚类分析法
. ]( ^5 _  `  p( }; G% G秩和比法
人工神经网络
熵权法
灰色关联度分析
$ o, r" z, }" V. g2 D预测算法
插值拟合
* }$ [: y5 \3 V回归模型预测
- d# Z% s5 z& x4 [8 e; {) R灰色预测GM
- g  i+ S: A; n( D$ Q时间序列
% ?1 Q; n$ X# o神经网络
. n- A0 p; T$ A7 i) \统计分析
: P0 ~# r' F% E3 }& e  W: o& n方差分析
  n$ s2 C. Z- J0 m2 ]. J" S. A回归分析
# ?* G) C$ F9 Q多元统计分析
分类问题
假设检验（非参数）
聚类
最优化方法
常用算法
" L; i# n6 i' X无约束优化
! @) d0 q: ?5 P) J有约束优化
: S& B% Q9 |, n二次规划
一般有约束非线性规划
计算机算法
" e1 @6 q' y+ q- q' {) i动态规划
回溯搜索
- \; u" a( w% A1 l- o. I1 d7 q2 U分治
贪心
: c; G. D' d# A- o6 J分枝定界
图论与网络
欧拉图
7 \- @1 Q/ z8 |" O二分图
网络流
$ I1 |' }8 O, u' V  B/ \/ V最短路径
TSP问题
& Y$ _7 t# D9 H" Q8 W8 v7 TP和NP问题
0 @* _, u' D2 e  g" a+ w' ^+ r# z背包问题
最小生成树
( n, b& ?( S6 j! |" M图的基本概念及其矩阵表示
0 ?; u4 K% b1 R$ }9 [1 L' x6 y常见的网络优化问题
/ K* l$ _0 N$ T: Q( q' P! u" O& O图像处理
& c. h) y! d+ t6 J9 [1 U. o$ q基本概念
数据类
图像类型
3 L6 v8 ?. ~% E8 K数据类之间转换
图像类型之间转换
' o' l' b' L2 J5 D) Z空间滤波器
频域变换
% @9 B$ E) W; f" `9 ^数字图像的水印防伪
' S0 G7 U3 W; e7 \  b4 n规划问题
. w$ e* f1 q4 v' \连续优化
离散优化
模拟和仿真
模拟
0 ?! p- Y% y) ~2 B7 C仿真
- q6 a+ Y% f/ p数值计算
常微分方程数值解
工具箱的使用
统计工具箱
. ^$ R* u( `+ d) \$ D' v( u拟合工具箱
) G8 l6 `! f+ I; G' x; z神经网络工具箱
4 y7 @, t7 d( @9 c9 W$ U小波工具箱
; S& ?9 H/ Z# l  l5 ~' b并行计算工具箱
  |# r/ T6 n1 h4 O优化工具箱
全局优化工具箱
0 F9 L/ h+ K1 Y2 D* O. S评价算法

  _% @6 V8 O, F  X+ ]简单加权法
: c* V$ q* m4 E8 {1 T
8 X( |+ h, I/ |线性加权
非线性加权
* W5 r2 z7 F5 A0 }4 ]$ k! G逼近于理想解的排序法(TOPSIS算法)

7 N: R  L& g3 Z层次分析法
" L( g3 A0 Y, o' `$ S& M7 x/ k
主成分分析法
0 I. Z! o, h# n
模糊综合评价法

- _- _! `5 P/ z: y聚类分析法
: l$ x$ W) A3 z' g9 p2 I1 w9 l& k5 `' _
秩和比法
% Y/ |1 F4 v2 v( T" F
人工神经网络
, O9 H, e: e& {6 Z% Z2 ^. Y" b
熵权法
0 ~0 F; w" z, r  ?1 @
灰色关联度分析
" n( b( i8 n! `
预测算法
7 X' i* E6 i* o) o) S" p% ^
0 q. T0 o5 X, U9 H; D/ q插值拟合

5 e8 j% X) U4 @) r+ Q小样本内部预测
回归模型预测
2 M) p/ I4 r  R
. i. }1 f8 P; s大样本内部预测
; h! ]; b3 R' P0 F灰色预测GM

* _6 G" L3 w4 g' o$ J! j4 N. T' `小样本未来预测
时间序列
- e/ ~: A( q) t* S8 _
* E1 q% ^' W& S4 |; h) i3 j( T大样本的随机因素或周期特征的未来预测
* e9 ]* x% {6 t% O/ C5 n神经网络
3 A6 f! S. A& l3 ?0 Z
针对大样本的内部机理复杂的数据的未来预测
- c7 v1 M8 ?& S2 {# \统计分析

* k" w& P$ p% e/ g+ b8 V; N方差分析
0 o( L9 R$ d% E% }8 y
分析因变量的总误差中，除开随机误差以外，是否有类别变量（自变量）造成的处理误差，有多少误差是自变量造成的
分类
. z' B# A- z9 b/ ^2 D单自变量
& X1 S0 n* G' v# n  n3 [$ x- M- y单因素方差分析
两个自变量
无重复双因素分析：只考虑主效应，不考虑交互效应
可重复双因素分析：考虑主效应，也考虑交互效应
回归分析

多元统计分析

9 h2 @: U0 s- W1 c1 {分类
聚类分析
判别分析
综合评价
主成分分析
8 `( ^% y" C9 H! B3 S, C( t$ T从原来的坐标系转换到新的坐标系，第一个新坐标轴选择的是原始数据中方差最大的方向，第二个坐标轴选择的是和第一个新坐标轴正交且具有最大方差的方向。该过程一直重复，重复次数为原始数据中特征的数目。
+ ^) X9 x) m9 F8 ^$ K  a5 e因子分析
& m- [) D0 L+ c+ V. W1 k% A& M1 P隐变量和某些噪声的组合
典型相关分析
1 i5 D6 |% C( o3 o偏最小二乘回归
分类问题

8 w* d, S) N8 d. Z神经网络
$ f1 j3 U3 g, n, B% r2 N! H逻辑回归
& ?% r' F4 V. T, ?8 _判别分析
最邻近方法
朴素贝叶斯
; l" ]8 g3 i5 M  D4 B1 W3 q支持向量机
+ |; u; F7 s. ^4 K& m0 N% v) ^2 q) W决策树
" X6 C& {' G$ A8 ~) r集成学习
9 ^: j4 O3 D, i! K( M0 Z5 ]. ~3 X. k: {ROC曲线
0 z  r% N7 w) Y  Z假设检验（非参数）
' l9 t+ s4 C, ~7 s
分布拟合检验
秩和检验
. }# S% l' P( ]如果两个样本来自两个独立的但非正态或形态不清的两总体，要检验两样本之间的差异是否显著，不应运用参数检验中的t检验，而需要采用之和检验
配对检验
6 z! N% r+ r: w3 S* XK-S 检验
  e' R# P6 z9 A, NK-S检验不仅能够检验单个总体是否服从某一理论分布，还能够检验两总体分布是否存在显著差异。其原假设是：两组独立样本来自的两总体的分布无显著差异。
. Y8 R0 G" ?+ y  T: L( U( q8 x$ [9 EQ-Q图
用变量数据分布的分位数与所指定分布的分位数之间的关系曲线来进行检验的
P-P图
- g9 n5 m  e4 [5 Q& f; b( O! Q根据变量的累积比例与指定分布的累积比例之间的关系所绘制的图形。通过P-P图可以检验数据是否符合指定的分布。当数据符合指定分布时，P-P图中各点近似呈一条直线。
聚类
% @( a# O) U9 n; V- @
层次聚类
划分聚类
, ?. [7 e, a8 K( |  N! RK-means
  s7 W1 c% s6 @0 G0 v4 Q2 a" b6 [8 A7 ]" G网络聚类
聚类评价
共表型相关系数
相似性矩阵
共表矩阵
( v* V8 T( K' J! `; Q外部指标
+ b* c3 {+ ?3 y+ D) ^相仪表
调整兰德系数ARI
互信息MI
自助
2 Z3 Q( `6 o" ^7 U) C# M5 n2 `最优化方法
! v3 S7 C+ P8 c+ O- u8 ^; k
常用算法
' p2 M  Y  R' _4 |* m, B
6 c- U* Q; b/ @/ M) Z5 s模拟退火
神经网络
/ G8 a' \9 {+ H, V; ?/ ?遗传算法
% }1 d6 l  C% s  \无约束优化
  a" j4 b  g7 K0 J6 U  ?
9 G9 d" V- `; c6 B基本算法
共轭梯度法
在最优化方法中占有重要地位，最速下降法的优点是工作量小，缺点是收敛慢，适用于寻优过程中前期迭代或作为间插步骤，当接近极值点时，宜选用别种收敛快的算法
4 ^9 a2 |3 _, @+ T! T* ]牛顿法
如果f是对称正定矩阵A的二次函数，用牛顿法经过一次迭代就可到达最优点，如不是二次函数，则牛顿法不能一部达到极值点
- c3 l: w( S8 V' E1 z5 n牛顿法收敛速度虽然快，但要求Hessian矩阵可逆，要计算二阶导数和逆矩阵，就加大了计算机计算量和存储量。
拟牛顿法
修改了牛顿方向
/ F4 o+ R* k9 F7 r3 C4 R两种算法
( ~5 A% D/ h& PDFP
BFGS
有约束优化
6 A) P% M; c; V8 _& Q+ L
1 k' M6 I( x4 p# Z8 Y* B* e4 U1 x罚函数法
通过构造罚函数把约束问题转化为一系列无约束最优化问题，进而用无约束最优化方法求解，这类方法称为序列无约束最小化方法，简称SUMT法（分为SUMT外点法 其二为SUMT内点法）
近似规划法
' X0 e6 d  o! f/ X; j/ f3 r* n$ O! J二次规划

  N+ ]' L6 P1 E( M+ k5 aquadprog
一般有约束非线性规划
$ n8 k; p0 K! h9 E* c$ a
fmincon
9 ^3 B! [) Q" Q" l& C计算机算法
) i. s" O+ ~/ Q- X4 n
5 K9 y. }, D  _9 t动态规划
: Z4 g* @: |5 p/ M
8 k+ E/ n4 o  K9 A回溯搜索

分治

2 W4 T* h% K; c贪心

. A) A9 c2 }: O0 f! Y分枝定界

; N% @2 G! f8 x/ m9 A图论与网络
8 B+ Y5 F  G- r* |. C5 N
: `* i0 c& a' R4 `. d3 U* K欧拉图
/ ]9 P$ J" E6 F0 p6 X, W
( w+ g% d# _$ X  h) G6 d& m. y二分图

匈牙利算法
6 l7 s  H+ W# w  Z网络流

8 |; ~5 y  O; Q+ t  l最小费用最大流
最短路径

$ D0 a, k* S$ O+ U2 J% gTSP问题
% a* L8 f3 I7 s, T! X
1 Z$ W. |9 f6 R4 r  S近似算法和启发式算法
# h6 t+ t- b- Z3 m二边逐次修正法
P和NP问题
% Z! ?: B+ I2 z: B8 G
背包问题
  Q) e% S+ s& ^
# f7 ^' Z. l6 R, h: S最小生成树

Prime算法
Kruskal算法
图的基本概念及其矩阵表示
' q4 ~: W9 E$ S  K
无向图
% P2 r9 L- H9 ]) X0 \9 `有向图
; R7 f" Y- P& F, e; R% c4 p完全图
二分图
0 U' I: h& \( W1 D9 l* U图与网络的数据结构表示方法
邻接矩阵表示法
. a" S5 e3 c8 L5 E* b( U- W关联矩阵表示法
弧表表示法
邻接表表示法
星形表示法
常见的网络优化问题
  e% i% s' J8 c3 }9 |1 j2 _
6 D# S" {- r- z9 x- e1 }1 M/ M1 R% |最短路问题
& ~( D6 c7 S8 {0 u公路连接问题
5 v; R% B1 D$ X! I中国邮递员问题（欧拉图）
旅行商问题TSP（哈密顿图）
运输问题
9 Y9 ~& p: C" u9 w; f图像处理
& s6 A- R+ V( K2 P2 P
+ t( v/ {) Y3 M0 S6 }6 Q基本概念

6 b* L: D- E) i+ W图像分类
离散图像
* P' x$ O" A$ X3 Y: U8 V: L) Z3 d用一个数字序列表示的图像，0和1
连续图像
% n" R9 J/ P, U& ?, Y二维坐标系中具有连续变化的图像
) E7 z3 L2 ?: }4 o矩阵中的元素称为像素，以256灰色等级的数字图像为例，一般由8位，即一个字节表示灰度值
! u5 I2 e, V  W# g; F8 R灰度值量化为对应灰度等级
% l0 S% b( {, K! @3 i3 \% U3 X等间隔量化（一般采用这个）
非等间隔量化
数据类

/ S' b5 q3 l1 b( I2 J+ g$ Y7 J数值数据类
3 Z: ]5 x3 ]5 g5 b# Ldouble
uint8
uint16
uint32
int8
int16
int32
7 x6 A- m+ E5 [8 Gsingle
9 N8 [$ A$ h, _" V6 F3 F3 n% d9 G* xchar
字符类（逻辑数据类）
* R: I& `2 G/ f8 Zlogical
图像类型

* l8 y9 U. y; H) U9 z' G2 y5 c二值图像
- l5 F. S! V: b4 ~6 M0 p二维矩阵由０、１构成，０为黑色，１为白色
通常用于OCR
一般二值图像是逻辑数组，只有０和１的uint8类数组，并不会认为是二值图像，需要使用logical函数 B= logical(A);
. M- M" D* L. A' d: ?灰度图像
即人们常说的256灰度图像，0表示纯黑色，255表示纯白色
6 p% Y  Q" a; c$ A5 n! J二值图像可以看成是灰度图像的一个特例
$ Q# b1 ~0 u  Q' e( X索引图像
有两个分量。即数据矩阵X和彩色映射矩阵map，矩阵map是一个大小为m*3且由范围在[0,1]之间的浮点值构成的double数组，map数组的长度同它所定义的颜色数目相等，map数组的每一行都定义单色的红绿蓝三个分量，X则是索引矩阵
真彩色RGB图像
是彩色像素的一个mn3数组，其中每一个彩色像素点，都是在特定空间位置的彩色图像相对应的红、绿、蓝三个分量
. L. N' {+ Y0 V) n4 ?6 W数据类之间转换

im2uint8
im2uint16
8 Z) K! q4 ?' Q8 tmat2gray：将输入转换为double，范围为[0,1]
im2double
8 f2 S* r5 K; `" n7 A9 nim2bw
图像类型之间转换
5 T( Z9 S- _6 G/ s/ j; y
2 {1 Y6 X& S" \4 y( |/ g7 V2 |ind2fray
6 H; H$ S5 k+ z/ q! p9 E, bgray2ind
) E8 G- V7 D2 E( w" `0 H, I; b" s* y' |rgb2ind
ind2rgb
ntec2rgb
rgb2ntsc
使用imtool命令查看一个图像文件的信息
空间滤波器

4 F& G, A& Q8 N5 x  a7 f$ k线性滤波器
使用拉普拉斯滤波器增强图像
$ F. v4 a7 A! j  X+ K使用fspecial生成过滤器以及imfilter的使用
+ f9 q/ E& X& n- C& Q" W非线性滤波器
* l! G$ w: A' ]% j- C' ?, `一个工具是ordfilt2函数，可以生成统计排序滤波器 g=ordfilt2(f,order,domain)
7 ^1 `  T& `' D* s, M3 r! g3 T) Zordfilt2函数生成图像g的方式：使用邻域的一组排序元素中的第order个元素来代替f中的每个元素，该淋雨则由domain中的非零元素指定
数字图像处理中最著名的统计排序滤波器是中值滤波器，对应第50个百分位，使用g=ordfilt2(f,median(1:m*n),ones(m,n));创建中值滤波器
工具箱提供了二维中值滤波函数g=medfilt2(f,[m,n])
频域变换

6 h1 f+ I& X) b/ }为了有效的对图像进行处理和分析，需要将原定义的图像空间的图像以某种形式转换到频域空间，利用频域空间的特有性质方便的进行一定的加工，最后转换回图像空间
傅里叶变换
将图像从空域变换到频域
) Q, p! j% g8 S! |0 A( z; Q# S二维连续傅里叶变换
二维离散傅里叶变换（DFT）
6 ]' R. X+ M  @( L0 {5 k/ l3 ~基于离散傅里叶变换的频域滤波
& X6 l# M6 h( _1 y2 G' M0 s0 k离散余弦变换DCT
# m" b# ^4 [+ m8 `8 F2 z1 m* I图像处理中常用的变换算法，通过DCT变换，可以将图像空间域上的信息变换到频率域上
* u, _1 I5 b9 P两种实现
基于快速傅里叶变换FFT的算法，通过工具箱提供的dct2
另一种是DCT变换矩阵，工具箱提供了dctmtx函数来计算变换矩阵
图像保真和质量
图像处理中为了增加压缩率有时会放弃图像细节或者其他不太重要的内容，为了衡量解码图像相对于原始图像的偏离程度，这些测度一般称为保真度准则
- u3 h5 W7 o  T  ]1 c# h1 i主要准则
5 S8 B4 n3 |2 B) q客观保真度准则
当所损失的信息量可用编码输入图与解码输出图的函数表示时，可以认为是基于客观保真度准则的
% [5 F+ t1 w+ U4 h0 f均方根误差
) G1 q5 u, I) E5 m% _均方信噪比（SNR）
均方根误差越小，峰值信噪比越大，处理的图像质量越好
/ {+ y+ J6 w% @) U( J主观保真度准则
3 M* L, r- J$ V( i+ r1 ~* Y用分数代表主观评价{很差，较差，相同，稍好，较好，很好}
% T4 G, l) R5 L2 f数字图像的水印防伪

; D4 Q' Q8 l2 i; [  n' T' T! O内嵌水印的特点
' t: Q+ `2 s$ b8 p5 h' H: ]: J5 {透明性
* M, G4 t# Q+ E2 z鲁棒性
7 W5 N9 \( h5 U6 o能够承受施加于图像的变换操作，不会因变换处理而丢失
: d: h8 \& g2 w5 S' d% n, t" R安全性
空间域水印
8 ]4 M% v4 U" O* }. W. r/ _将水印信息嵌入到载体图像的空间域特性上
频率域水印
将水印信息嵌入到载体图像的变换域系数等特征上
基于矩阵奇异值分解的数字水印算法
* C) e) c$ |# d# L' O奇异值分解定理
  F. H+ ~: s. J$ D: eWeyl定理
7 R( T$ [/ S  i' a& m在原矩阵上加一个小小的扰动，看是否矩阵奇异值的变化是否会超过扰动矩阵的最大奇异值，从而判断稳定性
水印嵌入
# J  _; R, L9 ^7 S) n* J水印提取
3 o) ]+ x2 E0 W5 Q) z5 \3 `嵌入的逆过程
基于DCT变换的水印算法
, t5 _' q+ C1 fDCT变换是实数域变换，对实系数处理更加方便，不会使相位信息发生变化，另外，DCT变换是有损图像压缩JPEG的核心，基于DCT变换的图像水印将兼容JPEG图像压缩
水印嵌入算法
' H3 F, E2 S1 O! U0 e) U. ~图像加密
" @, e: A, k' E+ K1 i# N! `3 a$ f图像隐藏
规划问题

连续优化

+ I- l# @) ^8 R  Z2 d线性规划 LP
$ D3 T! h! N7 Y非线性规划 NLP
二次规划 QP
离散优化

整数线性规划
整数非线性规划
纯整数规划
* e  a" w% u3 Y* @' U混合整数规划
一般整数规划
0-1整数规划
整数规划的特殊情形，要求线性规划模型中的决策变量只能取值为0和1
模拟和仿真

; z' c7 m; e- ^) I& @6 I0 x. k模拟
" f- i- U* M8 z6 o& |7 Z* N' |
拟合
/ H$ E1 ^4 d. a$ L8 L6 D
0 F* \6 v6 Q1 B4 K4 _, d- }非线性最小二乘法
Isqcurvefit
5 }9 T+ ]& f6 P( F. I$ KIsqnonlin
" x7 X9 m/ j6 J& Y3 \) _拟合与统计回归
4 \9 A- X( ^- I" L0 [* m! A线性回归
/ w1 u( f( {) B非线性回归
[beta,t,J]=nlinfit(x,y,‘model’,beta0)
常用解法：线性最小二乘法
a=polyfit(x,y,m)
插值
1 T/ G+ g; A& D( v+ y
7 n& t3 k  c2 H/ g; P1 i& N一维插值
拉格朗日插值

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

yy=interp1(x , y ,xx , ‘method’)
三次样条插值

yy=spline(x,y,xx)
& f: _" i4 C* ]3 w8 b二维插值
9 a# U* m9 g( y1 o/ O! l5 n$ z) O
: X; S8 G5 b" R( ~最邻近插值
分片线性插值
! {. {/ R. e( K9 u) {双线性插值
仿真

/ Y6 Z  L7 g$ x动态仿真
9 h  L- P2 c" D, r* v9 T连续系统仿真
5 {0 T$ b: }6 U8 u+ E时间步长法
离散系统仿真
& N  h8 `% H. M0 y/ v9 [. G" ^事件步长法
" n* |  Q$ C; d: E1 A静态仿真
) _6 Y: Z% W' U9 z' K; o. r8 l蒙特卡罗方法
随机数产生
4 P2 D  p6 h0 u0 z' {9 t& j; ?均匀随机数
其他方法
逆变换法
舍选法
5 }/ M: R" S4 x! y, q# H5 t近似抽样法
" T* ]0 N' l/ b7 ?4 S指数分布
泊松分布
标准正态分布
正态分布
二项分布
排队问题
可靠性问题
蒙特卡罗方法
3 {% d% P8 b+ w& W. d' ]元胞自动机
数值计算

常微分方程数值解

$ a* S# u$ v! |0 s微分方程解析解
u=dsolve(‘Du=1+u^2’,‘t’)
欧拉方法
, N: L" ?# q! Q+ P$ [4 e' d0 H0 y( U向前欧拉公式
向后欧拉公式
梯形公式
龙格—库塔方法
- e, F4 @6 V! F- xmatlab求解
高阶微分方程转一阶常微分方程
* R! ?5 D2 j! W, p' q刚性常微分方程与非刚性常微分方程
工具箱的使用
% {- o1 V) D8 ?% q$ N4 f
- \# a0 e" y8 A0 `* Y- S. w# B统计工具箱

; ]0 i2 m6 A! h# g拟合工具箱
% }* C" I  m/ M9 i- f5 o2 \! J3 V' c
7 s1 t8 p& V# v) f" S1 d神经网络工具箱
0 b2 W- n9 p( c2 Y5 h- a$ r
( ?+ J6 |# I4 m小波工具箱
6 ^' ?. V  W' a# e0 c
并行计算工具箱
$ a4 F, L9 g/ W$ ~- e
. W2 T5 S0 `7 Z1 Y优化工具箱
% n' K) ]1 ?$ f6 |+ Y. b4 K
一元函数极小值
5 L) t& A; r9 E7 Z3 J5 a2 mfminbnd
无约束极小
fminunc
% E! i; c% k/ V; F) T+ t为无约束优化提供了大型优化和中型优化算法
) _4 ]( F; x5 D6 r) a为中型优化算法的搜索方向提供了4种算法，由options中的HessUpdate控制，默认为bfgs（拟牛顿法的DFP公式） 另外，还有dfp（拟牛顿法的DFP公式），还有steepdesc（最速下降法）
+ e/ b$ @, Q" ~! ?为中型优化算法的步长一维搜索提供了两种算法，由options中的参数LineSearchType控制，默认是quadcubic （缺省值，混合的二次和三次插值）另外还有cubicpoly（三次多项式插值）
fminsearch
线性规划
& k1 r* d  h4 W( k! }9 O4 `linprog
二次规划
* k5 p  U$ Y# ~+ u( H2 nquadprog
- W* B& p  x8 A约束极小
fmincon
+ h& r: _* ^0 ^8 M* {4 d达到目标问题
fgoalattain
' k3 c' y  T1 ~6 k' \7 o极小极大问题
fminimax
) G' }/ |1 J* q  X全局优化工具箱
————————————————
) B. @- E: e0 J1 B: ~原文链接：https://blog.csdn.net/SinclairWang/article/details/100891757

2 b+ T; w; o' k  R
9 `7 s! s/ \! x

柠檬草lll

发表回复很全面

柠檬草lll

发表回复谢谢分享