QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2623|回复: 0
打印 上一主题 下一主题

LDA主题模型简介及Python实现

[复制链接]
字体大小: 正常 放大
杨利霞        

5273

主题

82

听众

17万

积分

  • TA的每日心情
    开心
    2021-8-11 17:59
  • 签到天数: 17 天

    [LV.4]偶尔看看III

    网络挑战赛参赛者

    网络挑战赛参赛者

    自我介绍
    本人女,毕业于内蒙古科技大学,担任文职专业,毕业专业英语。

    群组2018美赛大象算法课程

    群组2018美赛护航培训课程

    群组2019年 数学中国站长建

    群组2019年数据分析师课程

    群组2018年大象老师国赛优

    跳转到指定楼层
    1#
    发表于 2022-9-7 15:38 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    LDA主题模型简介及Python实现( Z9 y7 a2 |/ ?; F

    8 g  ?! J+ ]0 R一、LDA主题模型简介7 @6 A# P8 }9 ?5 y/ Q. j9 P' B
            LDA主题模型主要用于推测文档的主题分布,可以将文档集中每篇文档的主题以概率分布的形式给出根据主题进行主题聚类或文本分类。, S2 j$ T0 U  w" J( u4 g  s8 ]% p2 v
    8 O. n  R4 c3 G- [% d
            LDA主题模型不关心文档中单词的顺序,通常使用词袋特征(bag-of-word feature)来代表文档。词袋模型介绍可以参考这篇文章:文本向量化表示——词袋模型 - 知乎, p2 l+ U3 ?& y3 M/ c
    5 J$ P' r5 h8 P+ Y7 ]1 }
            了解LDA模型,我们需要先了解LDA的生成模型,LDA认为一篇文章是怎么形成的呢?
    ! g. |7 ^7 ~5 B8 s$ e8 ^2 I& s: R' X  _8 o. m4 [& j) d
            LDA模型认为主题可以由一个词汇分布来表示,而文章可以由主题分布来表示。
    , r* B4 M6 P( T, W
    ) j! I; k3 j9 B( j3 N; T3 O        比如有两个主题,美食和美妆。LDA说两个主题可以由词汇分布表示,他们分别是:
    6 _+ Z! k! r' F# x& i( |, I
    0 E* y" a. i% h0 x{面包:0.4,火锅:0.5,眉笔:0.03,腮红:0.07}- d, c! p: C% X
    {眉笔:0.4,腮红:0.5,面包:0.03,火锅:0.07}7 U& }' B$ q  R
    $ Y1 G" X- M( |0 T, n! l
            同样,对于两篇文章,LDA认为文章可以由主题分布这么表示:8 _* v) |* ]" N1 l1 O4 |2 a# R

    : R# H3 q3 Q6 f  P《美妆日记》{美妆:0.8,美食:0.1,其他:0.1}( f% \; a- O* b$ L
    : c$ [$ D& C' R' p# K4 m
    《美食探索》{美食:0.8,美妆:0.1,其他:0.1}
      l' ?4 }* z9 `0 _7 F* v* i7 m. }8 R  @' l& V- b4 H: S9 s) D
            所以想要生成一篇文章,可以先以一定的概率选取上述某个主题,再以一定的概率选取那个主题下的某个单词,不断重复这两步就可以生成最终文章。
    3 t* D5 L0 M* {# x
    7 C( A9 |4 ^. \/ z2 x5 u        在LDA模型中,一篇文档生成的方式如下:
    , O: ^8 m5 E" B8 p2 \7 U& K9 _5 `% i: O9 ~8 c$ u

    % J4 _% i6 ]) D( x  J4 ~. G' a8 ^9 D) [7 L/ ~5 j
            其中,类似Beta分布是二项式分布的共轭先验概率分布,而狄利克雷分布(Dirichlet分布)是多项式分布的共轭先验概率分布。
    ) `  X7 v6 D% p" s9 r  T- D) ?4 I$ ?6 c5 s' s% r* D9 O" L

    . g/ p" }  b8 n( ^7 ^9 n* k6 n6 X0 i, _2 d
            如果我们要生成一篇文档,它里面的每个词语出现的概率为:1 E9 B+ ~* m- N! M6 ]
    # c) W: O1 a6 n3 {2 a

    6 l2 p) v% u( a+ n9 [+ R1 \% U. h# b# ~3 l6 w
             更详细的数学推导可以见:通俗理解LDA主题模型_结构之法 算法之道-CSDN博客_lda模型% V: }4 W0 T' u  j2 r( P
    ( ^, ]! K+ m- j6 x7 M$ q1 }
            看到文章推断其隐藏的主题分布,就是建模的目的。换言之,人类根据文档生成模型写成了各类文章,然后丢给了计算机,相当于计算机看到的是一篇篇已经写好的文章。现在计算机需要根据一篇篇文章中看到的一系列词归纳出当篇文章的主题,进而得出各个主题各自不同的出现概率:主题分布。/ S  `7 W/ ]8 z  x5 ^2 b, ~
    ) N' O( R& V3 q0 y
            至于LDA主题模型在计算机中具体是怎么实现的,我们也可以不必细究,现在有很多可以直接用来进行LDA主题分析的包,我们直接用就行。(没错,我就是调包侠)" ]" w$ y9 x. B$ N  a

    * v5 s9 w+ c5 J$ v1 V- M& s二、Python实现
    1 D% b' ^' l# ~        在用Python进行LDA主题模型分析之前,我先对文档进行了分词和去停用词处理(详情可以看我之前的文章:用python对单一微博文档进行分词——jieba分词(加保留词和停用词)_阿丢是丢心心的博客-CSDN博客_jieba 停用词)
    . [6 ~: h9 Z% }# n* C6 x8 S
    ) x9 i: v% _5 b# U8 I' z* o8 q        我下面的输入文件也是已经分好词的文件
    , O' E- v: I# Y' D- ~
    6 X. V2 z% }: t. C3 o( ^% X8 F1.导入算法包9 t! J" B5 ?# }1 L; b. o& w4 K
    import gensim# t" @$ T+ k! O* K; }- ]
    from gensim import corpora4 \6 T' y! Z2 p& x
    import matplotlib.pyplot as plt. r7 q; _9 J5 `% y" v$ D& U9 a% R
    import matplotlib. @# q7 v) \% }' z& p% S
    import numpy as np1 Y/ I5 Z& o# y; h
    import warnings
    " x# s$ ^' s7 ]& Fwarnings.filterwarnings('ignore')  # To ignore all warnings that arise here to enhance clarity
    7 t, _# t1 p: |5 o3 M, _7 @0 K$ W$ h  {: U
    from gensim.models.coherencemodel import CoherenceModel3 R0 d/ j& L- G/ U( J( ?$ R) e& ?
    from gensim.models.ldamodel import LdaModel
    5 a7 ^* U8 f1 o6 h" p5 j2.加载数据) I5 {3 ^1 p/ l5 N& V
            先将文档转化为一个二元列表,其中每个子列表代表一条微博:4 A9 `; j/ ?% I) d. q9 P- \) Z
    * b$ @$ H! y) V; w$ y
    PATH = "E:/data/output.csv"
    1 z: V; B% I5 ?  n/ x# _5 R$ ^
    & @: v( V, b; yfile_object2=open(PATH,encoding = 'utf-8',errors = 'ignore').read().split('\n')  #一行行的读取内容3 C: c8 k* u8 Q: [5 i
    data_set=[]  #建立存储分词的列表2 H2 F' c3 ~7 B+ W7 _7 `
    for i in range(len(file_object2)):/ Q1 E, p% R2 \" @
        result=[]
    8 i+ A; \+ u% H! r  R% u    seg_list = file_object2.split()
    ( a0 r: z! l% }" u    for w in seg_list :  #读取每一行分词
    - F* l6 Q4 S( W3 T+ q& X  X        result.append(w)
    2 ~2 I( }) [3 Q6 [    data_set.append(result)
    3 h- U3 [. o3 M5 l  ]+ Gprint(data_set)  P1 a+ C; P* B% u
            构建词典,语料向量化表示:
    " t( z9 Y# G+ @) ]+ z& l, p+ N3 d/ R6 ~. Z9 K5 g8 e
    dictionary = corpora.Dictionary(data_set)  # 构建词典0 x# h) p/ y. A
    corpus = [dictionary.doc2bow(text) for text in data_set]  #表示为第几个单词出现了几次. S% @. z& C$ j! o4 A5 p' p3 K
    3.构建LDA模型
    . a8 `! P1 B: I7 Y: |ldamodel = LdaModel(corpus, num_topics=10, id2word = dictionary, passes=30,random_state = 1)   #分为10个主题
    ( `* m' p# [# }! ^print(ldamodel.print_topics(num_topics=num_topics, num_words=15))  #每个主题输出15个单词$ j) d% U* ]) m; L+ j: p
            这是确定主题数时LDA模型的构建方法,一般我们可以用指标来评估模型好坏,也可以用这些指标来确定最优主题数。一般用来评价LDA主题模型的指标有困惑度(perplexity)和主题一致性(coherence),困惑度越低或者一致性越高说明模型越好。一些研究表明perplexity并不是一个好的指标,所以一般我用coherence来评价模型并选择最优主题,但下面代码两种方法我都用了。
    ( C0 D6 {: P: ]& Q0 Z4 i2 j3 W9 U
    #计算困惑度, l/ s, U; E' N' G9 G
    def perplexity(num_topics):
    5 C- f+ _4 N( j, b$ O: E    ldamodel = LdaModel(corpus, num_topics=num_topics, id2word = dictionary, passes=30)
    3 N0 Q& Z/ p) V! L    print(ldamodel.print_topics(num_topics=num_topics, num_words=15))
    7 j& E5 A9 W! s    print(ldamodel.log_perplexity(corpus))7 m6 z0 U. K/ Q% J, X
        return ldamodel.log_perplexity(corpus)
    ) B( `8 s3 e6 C: T8 p+ @#计算coherence
    ) R- ^8 d8 Z# _4 w+ ^+ D6 e" \def coherence(num_topics):* |1 ?1 l$ H; s+ m  O  R
        ldamodel = LdaModel(corpus, num_topics=num_topics, id2word = dictionary, passes=30,random_state = 1)
    * y' a/ E: [1 ~    print(ldamodel.print_topics(num_topics=num_topics, num_words=10))
    ( @: @8 ~6 e: ?9 q    ldacm = CoherenceModel(model=ldamodel, texts=data_set, dictionary=dictionary, coherence='c_v')
    6 _5 V* s( q( X2 W  [    print(ldacm.get_coherence())
    ) w& f- |! T2 V, ], `; R" R5 j    return ldacm.get_coherence()- W. N: @& _/ f& S9 d' a2 M
    4.绘制主题-coherence曲线,选择最佳主题数% d5 N+ ?" ~- X  I
    x = range(1,15)# e0 n/ p" F& l& C% s
    # z = [perplexity(i) for i in x]  #如果想用困惑度就选这个
    8 l/ N4 O. W" y9 `/ y+ G4 I0 my = [coherence(i) for i in x]) A' p8 L$ R# d- _0 S; K
    plt.plot(x, y): L0 E9 V9 T! y& u0 F+ i; \+ _
    plt.xlabel('主题数目'), u/ j6 }. w' t
    plt.ylabel('coherence大小')+ f! F% d6 b' R& Y% y
    plt.rcParams['font.sans-serif']=['SimHei']
    0 s2 R# ~* x" e  ?" ~+ ^8 hmatplotlib.rcParams['axes.unicode_minus']=False
    , |& F2 x/ Z; tplt.title('主题-coherence变化情况')
    3 s! ]* j. Q$ z) {) Y* w4 Hplt.show()9 o7 Q0 w' y7 c
            最终能得到各主题的词语分布和这样的图形:
    ' }# ?& S& S0 F: L  ^1 `
    * T0 o# t% |& ?0 @0 S* r
    9 c5 d: M1 N; a- ]. N. ]: |6 f
    " D6 Y8 j) f% w2 d# r 5.结果输出与可视化
    5 Q% d6 x. x* L        通过上述主题评估,我们发现可以选择5作为主题个数,接下来我们可以再跑一次模型,设定主题数为5,并输出每个文档最有可能对应的主题
    8 |* q1 \0 ?1 y* z0 e
    $ ^& _. m7 \6 c0 n. s' Wfrom gensim.models import LdaModel% v* j* A/ R5 f  m
    import pandas as pd5 v/ f) w2 c5 k% L5 [1 Z
    from gensim.corpora import Dictionary
    2 c  a" p: Q. S- p  g. G3 K# u' y* _from gensim import corpora, models
    + F# z" w. B$ K# c+ l* uimport csv' m( i! W" A2 F/ G! V4 T
    ! Q0 i* Z, a$ K' y  i3 U9 I
    # 准备数据
    : b$ @& c5 `. O, K: J& ^# @PATH = "E:/data/output1.csv"5 S- b9 |5 _$ w# F8 ^

    * i" \1 N; F- p" f1 t& j& I" qfile_object2=open(PATH,encoding = 'utf-8',errors = 'ignore').read().split('\n')  #一行行的读取内容
    0 A7 s( v/ U( r/ ]+ |+ p. P6 Ldata_set=[] #建立存储分词的列表
      _  C$ g  H+ j8 A2 L9 Ffor i in range(len(file_object2)):
    - n  i+ b( n1 ~3 S    result=[]
    " P. [; m( \# j  t, ~2 A: L" P    seg_list = file_object2.split()
    1 y  [) ]- J3 F" [- L" _4 o4 V    for w in seg_list :#读取每一行分词9 r1 I1 l3 G5 z' D3 `5 J' l! A" g
            result.append(w); ?+ X/ Y2 F! f; @4 U  d- P
        data_set.append(result)
    ; [& w9 V& t( p. A$ h" N
    7 t8 B# H  z, }2 U4 `# G# l; s7 Bdictionary = corpora.Dictionary(data_set)  # 构建词典( ]5 J; V$ N* l+ u/ `" V
    corpus = [dictionary.doc2bow(text) for text in data_set]
    7 h7 l; z) U4 y; x! i4 n$ X5 t2 H1 k5 q
    lda = LdaModel(corpus=corpus, id2word=dictionary, num_topics=5, passes = 30,random_state=1)9 \8 _8 m3 i: @0 @5 H
    topic_list=lda.print_topics()7 a6 C( X; r- R+ }9 z  b% e9 {
    print(topic_list)6 b. U* V" o: H: ]2 w

    . e3 `) ^, p% e0 zfor i in lda.get_document_topics(corpus)[:]:
    & L9 s8 `; H5 Q7 S3 ^. r    listj=[]) r- V& X+ r- j
        for j in i:
    3 n! O7 N" A5 t: D, ?        listj.append(j[1])9 V& L) d, R1 T% C
        bz=listj.index(max(listj)): P" x1 p2 T. |& w6 S
        print(i[bz][0])
    8 C0 }2 f* q7 I  e1 K
    4 ^& E3 q" o- N3 n        同时我们可以用pyLDAvis对LDA模型结果进行可视化:/ L4 d. x* I" K5 i: p3 Q) i

    & U, I: g% ~1 w$ a/ wimport pyLDAvis.gensim' s3 T- r4 ]+ r! Q: a+ F# s7 L+ }
    pyLDAvis.enable_notebook()
    & b) _( @+ K% g% r. Y! @5 [data = pyLDAvis.gensim.prepare(lda, corpus, dictionary)
    4 J3 m. E. k3 d# c+ M% a, opyLDAvis.save_html(data, 'E:/data/3topic.html')
    - y  H8 Y# m$ g8 u9 P( E0 ~$ n1 P        大概能得到这样的结果:
    ; \. I/ O4 Z4 E0 D5 F* F# o! D& \/ y  [" [+ Y

    ) d% `3 l9 @5 Q3 S0 P7 u. j, h# d, C- e  S  a' N9 }
    左侧圆圈表示主题,右侧表示各个词语对主题的贡献度。
    ; ^2 e% g( U; o! `
    + B5 E5 C$ O+ p# v$ {; r6 o所有代码如下:9 Q7 {% }% A% o
    import gensim7 j5 W9 `+ q. y3 j* W" V! z
    from gensim import corpora
    # b7 l0 W5 b4 a. X' ^- ?7 ~import matplotlib.pyplot as plt
    - H& l+ t9 B5 A; `import matplotlib
    : {+ z+ `7 Z3 S& T$ Gimport numpy as np8 j3 o2 R2 l- t8 `  n/ h9 F1 H3 F
    import warnings" p! p& X* D+ c8 X9 J
    warnings.filterwarnings('ignore')  # To ignore all warnings that arise here to enhance clarity
    ' v; D" |8 D/ p5 N% M
    ! t# A0 D6 n, U; N) M/ pfrom gensim.models.coherencemodel import CoherenceModel$ o: g, d  \8 }
    from gensim.models.ldamodel import LdaModel
      G! s& F, ]. j' [' p! C: ]0 o0 l/ t; A0 c+ r3 E, g
    : G. ?* \# B. t0 [2 r
    , j1 L. q5 }" a$ y( S
    # 准备数据; f, E! R# G7 c/ v! e0 H: C; U0 R) y
    PATH = "E:/data/output.csv"
    % Q/ V0 o9 k( M: Y- t
    $ I2 E0 B# \: S/ ~file_object2=open(PATH,encoding = 'utf-8',errors = 'ignore').read().split('\n')  #一行行的读取内容8 h- j0 N( r0 G
    data_set=[] #建立存储分词的列表% L0 D! z* q. Z! S2 t& u
    for i in range(len(file_object2)):
    2 c9 g8 U! F' C    result=[]
    ) y7 ^5 b& M2 s/ F2 K1 M    seg_list = file_object2.split()
    ( Y7 ], P" F, V# k) d8 y    for w in seg_list :#读取每一行分词
    ! J  N5 j8 O$ w, C! H6 o. \        result.append(w)
    5 l7 r; ~# f" V9 J. Y    data_set.append(result)
    5 r$ S6 f" Y$ E' t9 pprint(data_set)" l" ^6 c. W/ e2 _) [/ r) L
    , j) g' `; y& h& [' i
    # V7 [8 @- B8 {% c! I# o& s! _
    dictionary = corpora.Dictionary(data_set)  # 构建 document-term matrix
    1 P9 V6 W* R. G  x6 A; A5 |1 x% l  ?corpus = [dictionary.doc2bow(text) for text in data_set]
    ' L7 h& Y0 q& r$ U7 C+ I9 @#Lda = gensim.models.ldamodel.LdaModel  # 创建LDA对象
    2 s  @1 l. F/ d! D& U; J( c1 Y  Z% z5 K7 E$ |
    #计算困惑度
    ) {- y1 k+ M( l8 ~def perplexity(num_topics):* l- \: F6 W9 w. U
        ldamodel = LdaModel(corpus, num_topics=num_topics, id2word = dictionary, passes=30)
    : G" x" Y) w: A    print(ldamodel.print_topics(num_topics=num_topics, num_words=15))
    * F- v0 L! S$ W4 ], T    print(ldamodel.log_perplexity(corpus))
    1 l: w  x$ n, ~0 j7 S    return ldamodel.log_perplexity(corpus)8 o+ m! q' d/ z; v. f9 c

    2 x4 J+ k$ c6 }8 Y% ^) M#计算coherence
    7 M$ ~. c* {0 ?( t4 q4 `; _def coherence(num_topics):  `/ I: C* \7 ]2 A* H% r5 V. L& C
        ldamodel = LdaModel(corpus, num_topics=num_topics, id2word = dictionary, passes=30,random_state = 1)
    + w1 f1 R$ c+ p! D* B. I    print(ldamodel.print_topics(num_topics=num_topics, num_words=10))* F* q5 e) j# ?3 z
        ldacm = CoherenceModel(model=ldamodel, texts=data_set, dictionary=dictionary, coherence='c_v')1 S6 o7 R6 z) u6 R9 r( R
        print(ldacm.get_coherence())
    5 \( D  w9 ~  I    return ldacm.get_coherence()
    : `9 J: A' Q$ Z, |3 O5 }' i3 w$ E/ o9 ^" c& W# ^  T1 F
    # 绘制困惑度折线图; {5 [/ G5 `5 k0 g) r/ w
    x = range(1,15)+ z) M1 c% s- B# {
    # z = [perplexity(i) for i in x]
    ! y9 f. S8 J# H7 v+ R  p3 ry = [coherence(i) for i in x]
    - u2 S! B9 C2 N7 V$ Lplt.plot(x, y)- u, r/ s( x$ z2 {+ U- S1 p6 C+ b
    plt.xlabel('主题数目'); x' [8 i. L- c: V
    plt.ylabel('coherence大小')
    . R+ {) m. l) m, }+ D8 W2 uplt.rcParams['font.sans-serif']=['SimHei']
    4 f% C- t; u' ^7 T+ Rmatplotlib.rcParams['axes.unicode_minus']=False
    , t# N5 t: A. x- c4 p) eplt.title('主题-coherence变化情况')7 H7 I  ]* E- A3 A3 h
    plt.show()
    , S4 }2 f: K" B
    ( W% G6 @& ?3 l( t% P- Y- X1 L& o7 D4 \from gensim.models import LdaModel, A. l& O  Q! l; |  |' r
    import pandas as pd6 T; j. B3 i# _" o( X. y
    from gensim.corpora import Dictionary# ?+ w# ]* ~  C+ K3 Z' o5 K
    from gensim import corpora, models
    + y7 V1 r( k) Y& E+ [import csv
    2 x' G/ V4 G; ]" W0 M
    7 T- [9 f; O, A' i4 N) N. F3 r0 y) z# 准备数据7 X$ C) A& V1 J" w
    PATH = "E:/data/output1.csv": J8 S) ]# i) M# W8 h7 Z: }

    : P- d, h9 f  Ffile_object2=open(PATH,encoding = 'utf-8',errors = 'ignore').read().split('\n')  #一行行的读取内容
    & L8 f5 Z; z4 R) mdata_set=[] #建立存储分词的列表' K$ T  j* @: u, y+ ?
    for i in range(len(file_object2)):) w& w! H3 L' k/ }
        result=[]/ F0 _+ {3 ], Z. X+ e6 U" z* q
        seg_list = file_object2.split(): \6 \* z( O# o, u; z) O
        for w in seg_list :#读取每一行分词
    7 h& H8 D/ v% L! `" ]1 z8 _  w% U        result.append(w)0 I0 K+ v+ a7 g9 X! \$ I" E
        data_set.append(result)
    4 |' x) ~/ E- l) t- I1 S3 X4 r
    * Q9 v* w7 m8 U$ c) gdictionary = corpora.Dictionary(data_set)  # 构建 document-term matrix
    7 D' T0 H! V! n$ p% e3 u2 c9 ^$ vcorpus = [dictionary.doc2bow(text) for text in data_set]
    5 s  R- f: y  \) Z4 n4 j7 V" Y! B2 T/ l# R) k! g9 s4 d
    lda = LdaModel(corpus=corpus, id2word=dictionary, num_topics=5, passes = 30,random_state=1)
    6 }. w9 U; X0 {* ^topic_list=lda.print_topics()7 W" z/ \% i6 o1 }# v3 _" V: e! T
    print(topic_list)6 N/ a4 D0 d, N7 ~& l$ ~  U
    # Z" T  Z6 Y' I9 C! B+ N7 k
    result_list =[]% Y& h" Q( i' t! W" i" W% M: ^
    for i in lda.get_document_topics(corpus)[:]:/ s6 b! N! {* ^
        listj=[]
    1 Q8 Q2 w. Z- K0 q5 \    for j in i:" B2 [0 M) F' l- W
            listj.append(j[1]): `6 J8 [% _- l! E2 D
        bz=listj.index(max(listj))) U! h( f+ X) C! u; a; i+ F
        result_list.append(i[bz][0]), ]$ z2 u5 F8 _8 M8 A+ t0 _( O/ |
    print(result_list)& m, w3 S* q4 L* t1 V: d* X; f2 }

    2 q  `1 B6 e8 limport pyLDAvis.gensim
    & ]5 d2 k. A* T8 k6 ^  KpyLDAvis.enable_notebook(); D  Y- p- ~3 w/ j5 P0 {
    data = pyLDAvis.gensim.prepare(lda, corpus, dictionary)) p& r. A5 _3 |- e" R, O: r" x% g
    pyLDAvis.save_html(data, 'E:/data/topic.html')0 Y0 l$ p: f) {3 G
    有需要自取~
    & L( F. E; ?' Q, T- F————————————————, W! m5 P% s0 [8 P& O
    版权声明:本文为CSDN博主「阿丢是丢心心」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。+ j5 a7 N# Y: m
    原文链接:https://blog.csdn.net/weixin_41168304/article/details/122389948
    ' N8 c. B: l4 a4 J1 S! e
    8 E: {! m! D/ l7 }# y* k
    ) |) m2 D; d% }  u9 S' _! N
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏1 支持支持0 反对反对0 微信微信
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-5-30 04:03 , Processed in 0.454716 second(s), 50 queries .

    回顶部