数学建模戒烟问题

.z.第2页，共10页2014高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了"全国大学生数学建模竞赛章程"和"全国大学生数学建模竞赛参赛规则"〔以下简称为“竞赛章程和参赛规则〞，可从全国大学生数学建模竞赛下载〕。我们完全明白，在竞赛开场后参赛队员不能以任何方式〔包括、电子、网上咨询等〕与队外的任何人〔包括指导教师〕研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料〔包括网上查到的资料〕，必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进展公开展示〔包括进展网上公示，在书籍、期刊和其他媒体进展正式或非正式发表等〕。我们参赛选择的题号是〔从A/B/C/D中选择一项填写〕：B题 我们的参赛报名号为〔如果赛区设置报名号的话〕：79所属学校〔请填写完整的全名〕：理工大学参赛队员(打印并签名)：1.2.3.指导教师或指导教师组负责人(打印并签名)： 〔论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上容请仔细核对，提交后将不再允许做任何修改。如填写错误，论文可能被取消评奖资格。〕日期：2014年9月3日赛区评阅编号〔由赛区组委会评阅前进展编号〕：-.z.2014高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号〔由赛区组委会评阅前进展编号〕：赛区评阅记录〔可供赛区评阅时使用〕：评阅人评分备注全国统一编号〔由赛区组委会送交全国前编号〕：全国评阅编号〔由全国组委会评阅前进展编号〕：-.z.摘要本文针对戒烟和各种相关因素关系，解决影响戒烟成功的因素有哪些的问题，我们利用E*cel，spss，Matlab软件对相关数据进展分析，利用主成分分析法、层次分析法建立模型，并进展可靠性检验，得到影响戒烟成功的重要因素，对有志于戒烟的人士提供戒烟对策和建议。针对问题一，首先观察数据发现数据残缺，运用spss中缺失值替换的线性插值法将缺失数据补充完整。然后求出在不同年龄段、不同性别、不同调整CO浓度下、不同每日抽烟数的条件下的发病累加率，并运用E*cel，Matlab，Spss的等软件做出图像，结合图像对发病累加率的分布进展分析，得到累加发病率的根本趋势为随着每日抽烟数和调整CO浓度的增加而提高，男性比女性的累加发病率略低。针对问题二，首先用spss做其它各变量与戒烟天数的相关性分析，得到显著性〔双侧〕值，CO浓度和戒烟天数0.01水平显著相关，距离最后一支烟的分钟数和调整的CO浓度与戒烟天数0.05水平显著相关；然后，根据题意，CO浓度和距离最后一支烟的分钟数共同影响调整的CO浓度；可以得到每日抽烟数和调整的CO浓度会影响戒烟时间〔天数〕长短。最后，对相关变量的数据进展分组，在不同区间对戒烟天数求平均值，然后用spss回归分析中曲线估计对数据进展拟合，再用MATLAB拟合求回归方程中系数的置信区间来分析拟合效果。得到每日抽烟和戒烟天数的回归方程为，=0.948，F=12.038，Sig=0.078，可得拟合效果较好；调整的CO浓度和戒烟天数的回归方程为，=0.785，F=12.799，Sig=0.005，拟合的较好。针对问题三，我们认为CO浓度和距抽最后一支烟的分钟数是来控制调整CO浓度的，这两个因素可以用调整CO浓度说明，因此在考虑戒烟成功主要因素时没有针对这两个因素讨论，则可以假定戒烟成功受年龄，性别，每日抽烟数以及调整CO浓度影响。就这四个因素建立层次分析，探究每个因素对戒烟成功的影响。先用主成分分析法对影响因素探究，再运用层次分析法，根据主成分分析得到的累计奉献率来建立判断矩阵，并算出各成分的权向量。结果显示影响戒烟成功的因素主要为每天抽烟数，调整CO浓度和年龄。其比重分别为37.15%，25.00%，23.22%。通过一致性检验，判断矩阵具有满意一致性，可以为模型提供可靠分析。针对问题四，我们根据前三问所得结果向有志于戒烟的人士提供相应的戒烟对策及建议，撰写的报告具体见5.4问题四的求解。关键字：线性插值法相关性分析回归分析主成分分析层次分析法-.z.一、问题的重述众所周知，吸烟不仅危害自身安康，而且由此引起的被动吸烟更是危害公众身心安康的主要原因。为此，如何帮助相关人士摆脱烟瘾的困扰也就成为一个重要的研究课题。本文研究数据涉及234人，他们都自愿表示戒烟但还未戒烟。在他们戒烟的这一天，测量了每个人的CO〔一氧化碳〕水平并记下他们抽最后一支烟到CO测定时间.。CO的水平提供了一个他们先前抽烟数量的客观指标，但其值也受到抽最后一支烟的时间的影响,因此抽最后一支烟的时间可以用来调整CO的水平。记录下研究对象的性别、年龄及自述每日抽烟支数。这个调查跟踪1年,考察他们一直保持戒烟的天数,由此估计这些人中再次吸烟的累加发病率,也就是原吸烟者戒烟一段时间后又再吸烟的比例.其中假设原烟民戒烟的可信度是很低的〔更恰当地说多数是再犯者〕戒烟天数是从0到他〔她〕退出戒烟或研究截止时间〔1年〕的天数。假定他们全部没有人中途退出研究。请答复以下问题：1〕试分析上述234人中再次吸烟的累加发病率分布情况〔如不同年龄段、不同性别等因素下的累加发病率分布情况〕。2〕你认为年龄、性别、每日抽烟支数及调整的CO浓度等因素会影响戒烟时间〔天数〕长短吗？如果影响请利用附录中的数据，分别给出戒烟时间与上述你认为有影响的因素之间的定量分析结果。3〕请利用附录中的数据建立适当的数学模型，讨论影响戒烟成功的主要因素有哪些，并对你的模型进展可靠性分析。4〕请根据你的模型，撰写一篇500字左右的短文，向有志于戒烟的人士提供戒烟对策和建议。二、模型的假设1、原烟民戒烟天数缺乏365天的数据都是可靠的；2、原烟民戒烟的可信度很低，可以说他们多数是再犯者；3、自愿者中全部没有人中途退出研究。三、符号说明符号符号说明年龄性别每日抽烟数调整的CO浓度戒烟天数权向量一致性指标随机一致性指标一致性比率最大特征根四、问题的分析4.1问题一的分析针对问题一，首先对缺失数据进展补充，然后求出在不同因素下，不同年龄段、不同性别、不同每日抽烟数、不同调整CO浓度条件下的累加发病率的分布情况，并作出图表，进展分析。4.2问题二的分析针对问题二，首先对其它各变量与戒烟天数的相关性分析，得到和戒烟天数相关的变量，然后根据题意，得到会对戒烟天数长短的变量。对相关变量的数据进展分组处理，求不同区间戒烟天数的平均值，把组距中间值作为那一组的数值，用spss和MATLAB对数据进展处理得到相关变量和戒烟天数的定量分析结果。4.3问题三的分析我们认为CO浓度和距抽最后一支烟的分钟数是来控制调整CO浓度的，这两个因素可以用调整CO浓度说明，因此在考虑戒烟成功主要因素时没有针对这两个因素讨论，则可以假定戒烟成功受年龄，性别，每日抽烟数以及调整CO浓度影响。就这四个因素建立层次分析，探究每个因素对戒烟成功的影响。五、模型的建立与求解5.1问题一的求解5.1.1问题一模型的建立观察附录可以发现，附录中所给数据缺失，为了更加全面的进展分析，首先我们使用ｓｐｓｓ中缺失值分析命令对于数据进展分析，结果如下。列表均值案例数目Min_lastLogCOadjCO224177.311397.48257.37列表协方差Min_lastLogCOadjCOMin_last91053.229LogCOadj-27.30741300.914CO-20993.63522198.70818998.861列表相关性Min_lastLogCOadjCOMin_last1LogCOadj.0001CO-.505.7921回归均值aMin_lastLogCOadjCO179.581399.76258.04a.将随机选中的案例的残差添加到各个估计。回归协方差aMin_lastLogCOadjCOMin_last94546.889LogCOadj-957.70840165.315CO-20909.16621649.54018918.215a.将随机选中的案例的残差添加到各个估计。回归相关性aMin_lastLogCOadjCOMin_last1LogCOadj-.0161CO-.494.7851a.将随机选中的案例的残差添加到各个估计。所给数据缺失率为0.0446，缺失率较小，可以利用spss进展补充，然后使用缺失值替换命令进展数据的补充，经过对已有数据的分析，决定使用线性差值法进展补充，补充后数据见附表。为了直观得了解所有234人的总体戒烟情况，作出如下散点图：由图5.1.1看出，被调查的234人多数经过很短时间后又再次抽烟，只有少数人戒烟天数到达365天，只占总人数的14.10％，本文假设在研究截止时间没有再抽烟的烟民戒烟成功，所以累加发病率的具体定义为戒烟天数小于365的烟民数量占研究样本总人数的比例，本文据此对在不同性别、年龄、每日抽烟数、调整后CO浓度的情况下对累加发病率进展比拟和分析。首先运用spss做出男性与女性吸烟者的戒烟天数频率分布图，并求出男性与女性的累加发病率，做出图像，图表如下。性别成功戒烟人数总人数累加发病率男171100.845女161240.871由图5.1.2看出男性与女性再次吸烟的累加发病率分布情况相差不大，而由表5.1.1可以发现男性累加发病率为0.845，略低于女性的0.871。年龄：通过对所给数据进展分析可以得到参与研究人员的年龄跨度为22—72，所以本文将所有参与研究人员的年龄分为6段，分别为20—29、30—39、40—49、50—59、60—69、70—79，经过统计可以得到不同年龄段的累加发病率的分布，如下图。年龄段成功戒烟人数总人数累加发病率20-294440.9130-3911640.8340-498570.8650-597470.8560-691170.9470-79150.80通过表5.1.3和图5.1.1可以看出30-59年龄段的累加发病率较低，70岁以上的被研究者累加发病率最低，60-69年龄段的被研究者的累加发病率最高。每日抽烟数：首先用Matlab做出每日抽烟数跟戒烟天数的散点图通过对数据的分析可发现每日抽烟数的跨度为2—90，所以将每日抽烟数分为7个档次，分别为1—10、11—20、31—40、41—50、51—60、60以上，求出其累加发病率，并做成图表，如下。每日抽烟数戒烟人数总人数累加发病率1—104170.76511—2010810.87721—3012720.83331—404370.89241—50180.87551—60071.00060以上220.000通过观察图表可得，每日抽烟数60以上的样本数过少具有偶然性，所以在分析时，可适当忽略，由图表可以看出，根本趋势为每日抽烟数越多则累加发病率越高。调整CO浓度：首先用E*cel做出调整CO浓度跟戒烟天数的散点图，如以下图。通过对数据的分析可以得到调整CO浓度的跨度为682—1951，所以将调整CO浓度分为7个区段，分别为600-799、800-999、1000-1199、1200-1399、1400-1599、1600-1799、1800-2000，求出其累加发病率,并做成统计图，图表如下。调整的CO浓度成功戒烟人数总人数累加发病率600-799110.000800-9992100.8001000-11992190.8951200-139910730.8631400-159915980.8471600-17993310.9031800-2000021.000通过图表可以发现累加发病率的趋势是随着调整后CO浓度的增加而提高。5.2问题二模型的建立与求解5.2.1问题二模型的建立由图5.1.2可以看出戒烟天数的长短与性别的关系不明显，由spss软件得到它们的相关系数为0.006，所以可以认为戒烟天数和性别无关。由spss中双变量相关系分析得到的每个变量之间的相关性数据〔见附表〕，得到变量Age、Gender、Cig_Day、CO、Min_last、LogCOadj和戒烟天数的显著性〔双侧〕值分别为0.458、0.925、0.686、0.002、0.019、0.037。表5.2.1相关变量与戒烟天数的相关性变量显著性〔双侧〕值相关性Age0.458不相关Gender0.925不相关Cig_Day0.686不相关CO0.0020.01水平显著相关Min_last0.0190.05水平显著相关LogCOadj0.0370.05水平显著相关根据题意变量CO和Min_last可得到LogCOadj。LogCOadj水平提供了一个自愿者们先前抽烟数量的客观标准，在一定程度上，LogCOadj和Cig_Day有一定的关系。由此，我们可以主观的认为每日抽烟数和调整的CO浓度能影响戒烟时间的长短。戒烟天数为365的虽然不能确定其以后会抽烟，但是，他们已经坚持了365天，我们要研究影响因素与戒烟时间的定量分析，就不能忽略365天的数据。我们分别对每日抽烟数和调整的CO浓度和戒烟天数的数据进展处理。由上表可得每日抽烟数的围为[2，90],我们对每日抽烟数以10为组距在[1，90]分组，处理后的数据为表5.2.2每日戒烟数和戒烟天数统计表每日抽烟数组距中间值戒烟总天数总人数戒烟天数均值1-105.5304317179.011-2015.553958166.621-3025.560047283.431-4035.526513771.641-5045.5689886.151-6055.5200728.661-7065.53651365.071-8075.5000.081-9085.53651365.0用E*cel可以得到以下图图5.2.1每日抽烟数和戒烟天数折线图由图5.2.1可以看出[61,90]的两个点均为戒烟天数为365天的，偏离了图形正常的趋势，这两个点只能说明个别情况，不能代表整体趋势，所以把这两个点删除后再处理数据。表5.2.3删除后每日戒烟数和戒烟天统计表每日抽烟数组距中间值戒烟天数均值1—105.5179.011—2015.566.621—3025.583.431—4035.571.641—5045.586.151—6055.528.6用spss中曲线估计进展拟合，分析后得到立方拟合效果最好。此时=0.948，F=12.038，Sig=0.078。用spss得到拟合的图形，如下：表5.2.4每日抽烟数和戒烟天数相关分析表模型汇总和参数估计值因变量:戒烟天数方程模型汇总参数估计值R方Fdf1df2Sig.常数b1b2b3三次.94812.03832.078278.633-22.711.760-.008自变量为每日抽烟数。由表5.2.4得到每日抽烟数和戒烟天数的回归分析方程式：〔5-1〕用Matlab求出参数的置信区间参数参数估计值参数置信区间常数278.633(115.2,442.1)b1-22.711(-46.27,0.8484)b20.760(-0.1359,1.656)b3-0.008(-0.01748,0.001903)R2=0.761314F=25.9115P=0.000653633可以看出每日抽烟数和戒烟天数在整体上呈负相关，调整的CO浓度与戒烟天数的定量分析图5.2.3调整的CO浓度值从小到大散点图由上图和调整的CO浓度频率表可得其在[682，1951]分布均匀。我们在[600,1999]按组距为200调整的CO浓度组距中间值戒烟总天数总人数平均值600-799699.53651365.0800-999899.58891088.91000-11991099.513001968.41200-13991299.567577392.61400-15991499.578159879.71600-17991699.515753150.81800-19991899.51125.5我们用spss曲线估计拟合得到二次的R方为0.771且为最大值。此时F值为5.044较小，Sig值为0.110较大，拟合的不好。模型汇总和参数估计值因变量:戒烟天数方程模型汇总参数估计值R方Fdf1df2Sig.常数b1b2二次.7715.04423.1101141.496-1.581.001自变量为调整的CO浓度。所以我们绝定取组距为100，在[681，1980]分13组进展处理。调整CO浓度组距中间值戒烟天数均值681—780730.5365781—880830.5124.67881—980930.533.25981—10801030.5104.631081—11801130.558.181181—12801230.5115.341281—13801330.581.941381—14801430.589.741481—15801530.573.281581—16801630.552.91681—17801730.52.131781—18801830.5351881—19801930.57此时R方为0.530，F值为5.646，Sig值为0.023，R方较小，我们考虑到可能是戒烟天数为365天数据的存在造成的影响比拟大，因为题目中假设原烟民戒烟的可信度是很低的，在234个数据中戒烟天数为365天的数据为33，其频率为0.141，所以戒烟天数为365天的只是一少局部人，不能代表整体，我们将[681,780]区间戒烟天数平均值为365的数据删除，再把偏离总趋势较大的点删除后进展处理。此时用spss曲线估计拟合得到三次的R方为0.785且为最大值。此时F值为12.799，Sig值为0.005较大，拟合的较好。由spss得到的参数估计值，可以得到调整的CO浓度和戒烟天数的回归分析方程式：〔5-2〕由Matlab得到参数的置信区间见下表：参数参数估计值参数置信区间常数117.306(-504.3,1693)b26.436×10-6(-0.001076,0.002752)b3-1.784×10-8(-6.782e-7,2.446e-7)R2=0.785F=12.799P=0.005可以看出调整的CO浓度和戒烟天数呈负相关。5.3问题三模型的建立与求解我们认为戒烟天数到达365天的为戒烟成功，抽取出所有戒烟到达365天的数据，根据成功戒烟的数据分析影响因素。根据一、二问结论，我们假定戒烟成功受年龄，性别，每日抽烟数以及调整CO浓度影响为消除不同变量的量纲的影响，且不改变变量的相关系数，首先需要对变量进展标准化处理，标准化处理可通过spss实现。我们忽略其他因素对戒烟成功的影响，认为年龄，性别，每日抽烟数以及调整CO浓度为影响戒烟成功的主成分，相应的主成分特征值和累计奉献利率如下表：解释的总方差成份初始特征值提取平方和载入旋转平方和载入合计方差的%累积%合计方差的%累积%合计方差的%累积%11.48637.15537.1551.48637.15537.1551.00125.01325.01321.00025.00262.1561.00025.00262.1561.00025.01050.0233.92923.22185.377.92923.22185.3771.00024.99675.0194.58514.623100.000.58514.623100.000.99924.981100.000提取方法：主成份分析。考虑到每一种成分对戒烟成功的影响不一样，为确定各指标对戒烟成功影响的权重，本文采用层次分析法进展影响程度评判。模型的层次构造如下：目标层A为“戒烟成功因素分析〞；准则层B包括年龄*1、性别*2、每天吸烟支数*3和调整的CO浓度*4；方案层C为“戒烟成功〞和“戒烟失败〞。对于层次分析法中的判断矩阵，根据每个成分的累计奉献率来确定各个指标之间相对重要程度，判断矩阵表如下：年龄性别每日抽烟数调整CO浓度年龄1.0001.5880.6250.929性别0.6301.0000.3940.585每日抽烟数1.6002.5411.0001.486调整CO浓度1.0761.7090.6731.000得到判断矩阵后，求其最大特征向量，将该特征向量归一化处理后即可得到个影响成分的权向量：为了能用成比照拟阵的对应于特征根的特征向量作为被比拟因素的权向量，我们需要对成比照拟阵进展一致性检验。检验过程中，我们需要借助Satty提出的一致性指标理论：其中，CI=0时，矩阵为一致阵；CI越大，矩阵的不一致程度越严重。为了衡量一致性指标CI的标准，我们还需要引入随机一致性指标RI，计算公式为：其中，RI数值可以查表得知〔见下表〕。满足CR<0.1时就能通过一致性检验。查阅相关资料，查找随机一致性指标RI的数值如下：n1234567891011RI000.580.891.121.241.321.411.451.491.51接下来借助Matlab运算易得，，，代入RI=0.89，得：。CR<0.1，说明成比照拟阵的不一致程度在允许围，判断矩阵具有满意一致性，可以作为评判的权向量。根据各因素权重向量可知戒烟成功主要因素有每天抽烟数，调整CO浓度和年龄。其比重分别为37.15%，25.00%，23.22%。5.4问题四的求解众所周知，抽烟可能诱发多种疾病，如尼古丁对心、脑血管的影响，还可致癌，对个体安康危害极大；同时被动吸烟者所吸入的有害物质浓度并不比吸烟者为低，受到的伤害甚至高于吸烟者，所以现如今戒烟未尝不是一个良好的选择。结合上面所建模型对问题一到问题三的求解，我们可以为戒烟人士做出如下的建议与指导。通过对问题一的探究发现，戒烟天数与性别显著性并不明显，所以对于戒烟人士来说性别并不是阻碍，无论男士或是女士只要心智坚决，是可以戒烟成功的。通过对年龄的分析得出，年龄越小者的累加发病率越高，且年龄对戒烟天数奉献的权重到达0.2322，可能与年轻人自制力不强有关，对于年青的人士来说戒烟可能会更困难些，所以更应注重体育锻炼和户外运动，坚决信念、早日戒烟成功。通过对问题二和三的分析不难发现，每日抽烟数少、调整CO浓度低的被研究人员戒烟成功率较高，可能与烟草的成瘾性有关，而每日抽烟数和调整的CO浓度对戒烟天数奉献的权重分别为0.3715和0.2499，所以建议戒烟人士戒烟时采取逐步减量法，并换抽CO含量低的烟草。总体来说，戒烟不是一朝一夕的事情，在戒烟期间要加强体育锻炼，多进展户外运动、呼吸新鲜空气，同时可辅助药疗和针灸，以到达更好的效果。戒烟期间最重要的因素是戒烟者本人，只有本人坚决信念才能最终戒烟成功。六、模型的评价及改良模型优点：1、分区段对数据对相关性不是很明显的数据进展在不同区段取平均值处理，可以是数据之间的关系更加清晰。2、运用层次分析法时建立的判断矩阵利用主成分分析中的奉献度来赋值，有效地防止了主观感受对两指标间影响程度进展赋值的人为因素。模型缺点：1、只考虑了不同因素和戒烟天数的影响，忽略了各因素之间的数据在处理时都具有相关的影响。2、认为戒烟天数到达365天的为戒烟成功，只考虑了各因素对成功戒烟的影响。模型改良：1、在对不同因素的数据进展处理时，可以以戒烟天数为定值，对其它因素的值取平均值，建立不同因素之间和戒烟天数有联系的方程，对方程进展相关性检验。2、可以设置一个戒烟成功程度，分析各因素对戒烟程度的影响，使分析更全面科学。七、参考文献[1]启源，金星，叶俊，数学模型〔第四版〕，，高等教育，2011年[2]邓维斌，唐兴艳，胡大权，spss统计分析教程，，电子工业，2012年[3]周品，新芬，MATLAB数学建模与仿真，，国防工业，2011年八、附录问题二(1)beta0=ones(4,1);[beta,R,J]=nlinfit(*,y,'myfun1',beta0)'myfun1'functionf=myfun1(beta0,*)*=[5.500015.500025.500035.500045.500055.5000]';*1=*;*2=*1.**1;*3=*1.**1.*1;a=beta0(1);b=beta0(2);c=beta0(3);d=beta0(4);f=a+b.**1+c.**2+d.**3(2)beta0=ones(4,1);[beta,R,J]=nlinfit(*,y,'myfun1',beta0)'myfun1'functionf=myfun1(beta0,*)*=[1.0e+003*0.83051.03051.13051.23051.33051.43051.53051.63051.83051.9305]';y=[124.67104.6358.18115.3481.9489.7473.2852.9357]';*1=*;*2=*1.**1;*3=*1.**1.*1;a=beta0(1);b=beta0(2);c=beta0(3);d=beta0(4);f=a+b.**1+c.**2+d.**3问题三A=[1.00001.58800.62500.92900.63001.00000.39400.58501.60002.54101.00001.48601.07601.70900.67301.0000][*,y]=eig(A)*=-0.44220.0824-0.4639i0.0824+0.4639i0.3616-0.2785-0.2746+0.1759i-0.2746-0.1759i-0.5297-0.7074-0.4012+0.2952i-0.4012-0.2952i0.7673-0.47590.65080.6508-0.0001y=4.000300000.0002+0.0002i00000.0002-0.0002i0000-0.0006得*为特征向量矩阵

y为特征值矩阵

%找到y中对应最大的特征值所在列m

w=*(:,m)/sum(*(:,m))

w就是权重！

w=*(:,1)/sum(*(:,1))处理过的数据Id年龄性别每日抽烟数CO浓度距离抽最后一支烟的分钟数调整的CO浓度戒烟天数118221141406011332285231402401013367141241202706014191416824114260701409332226120175321212253927120310489133412332712980945145536527727181809012623305271201401051162721322811539065158114166281203456515282130729135435101594433029116430601621151283012234050151216340301232706014191326231114953301137354132125580601751184321181959012967943214028090145445252321333401014873293321202752401541203283212545070164720443312031515145732643316501440156736523833128480551665133443311019011012986134341132501571447715135120270621420902533512020060128882663516035011494121207361132503313681422936140500101654124273361343907515882276361304106016001445371405401801796514637130450516052185371253555148630313371202605136736532137123851440179796334371204004515807345371208570013251965939135420901630365160391340205682365224391303201014603239391181909512882311391301904381507453253911519078017252249401112654131583213044113029045144036533941117210701316111014214041565160821274215027560142773294216016015116325343160450281315314943125310401466017043120240551364613644134330451496627444110852401031155451403006014643652154612020010125636524546123651020141330224746123255120143211251461403701151813178471232706014195834815041545159611295481402104012973300481232609014211308481201007099411864919170821232602845014070908513651751120520120174221775112850070169332915113017045120836514521201601051220471215214026099139230222521253609215647214531245980122742255319019051231365264531602402101463142675312043595164815354130220160139338541302557014004557541354205516073347541202906014507133551201406121486034155123265130145515055612525070139216102561202601651469471625612014012511753715814043085163736589591304207516201175591203006514682915601111609512143616113027060141962116112020091130836529462150440451621621963130215901339365626312030012015032235631262900141136528863120270301399130346411019085128216364120180551239742786412012510105243336414018560125467686616022060133042323661201501951249124671302705014123653327011329015142112372169071013563651047611514090115393062121512060106671143222447107518484210232181806612463236242604513801483143142422510060987167425291103010091410525215440971654113525230290751459157203252203001751538290262304506516441416126215401164129425241262404858016861246262332808014471312262109074013753355127220130225120636515927220180100126863189272670140883142792721515580119016124282828011014662139282172851001468106205282302907514593652312823540520156941242282915090118260244282521026128823726528218475201638430228240601035138731133728225240751377333292281459011682111229222375601561365179292203305014996335292143601201582499302202708514352109302201501511353651263022834050151236516430245480601668422630225150135121121320302403206014920763123043055161833153227401320139315320833225380215305530933211270401406365144342291708012306348342222651501468153035