2. 武汉大学全球健康研究中心 湖北武汉 430071;
3. 四川大学华西公共卫生学院 四川成都 610041;
4. 四川大学健康城市发展研究中心 四川成都 610041
2. Global Health Institute, Wuhan University, Wuhan Hubei 430071, China;
3. West China School of Public Health, Sichuan University, Chengdu Sichuan 610041, China;
4. Institute for Healthy Cities and West China Research Center for Rural Health Development, Sichuan University, Chengdu Sichuan 610041, China
癌症作为一种严重威胁人类生命健康的恶性疾病,是延长预期寿命的重要障碍,已成为全球性的重大公共卫生问题。中国作为世界上最大的发展中国家,随着经济社会快速发展和人口老龄化加剧,面临着更为紧迫和严峻的癌症防治形势。[1]近年来,中国在全球癌症新发病例和癌症死亡病例数据排行中均居于高位。[2]最新评估报告指出,2022年中国大约有超过480万的癌症新发病例和超过320万的癌症死亡病例[3],并且高昂的治疗费用给患者和医疗保障系统带来了极为沉重的经济负担和负面冲击[4]。
为破解这一难题,提高包括抗肿瘤药物等在内的一系列专利创新药及独家原研药的可负担性和可及性,中国政府在2015年启动了国家医保药品价格谈判(简称“国谈”)工作,以进一步确定纳入国家医保目录中的药品支付价格标准。[5]相关研究指出,经过价格谈判,进入国家医保目录的抗肿瘤药品(简称“国谈抗肿瘤药”)的价格明显下降,显著减轻了患者的经济负担。[4, 6-7]为进一步提升抗肿瘤等类别国谈药的可及性水平,2021年4月,国家医疗保障局、国家卫生健康委员会联合印发《关于建立完善国家医保谈判药品“双通道”管理机制的指导意见》,提出在定点医疗机构之外,增设定点零售药店作为国谈抗肿瘤药的第二购药渠道,执行与定点医疗机构相同的医保报销待遇,以更好的满足患者的合理用药需求。
不难看出,国谈政策先是有力地压缩了国谈抗肿瘤药的价格空间,同时借助有关的政策工具在继续强化医疗机构首要供应责任的基础上将社会零售药店作为重要供给主体。前者在较大程度上提升了患者的经济可得性,而后者则拓展了患者的空间可得性。作为国谈抗肿瘤药可得性的两种重要构面,经济可得性已经得到了学术界和政策界的大量关注;相比之下,其空间可得性(在本研究语境下,亦可称为“空间可达性”)却几乎无人问津,这也成为当前医疗卫生资源配置领域研究的明显缺憾。因此,本研究以武汉市为例,以区域内配备《国家基本医疗保险、工伤保险和生育保险药品目录(2021年)》(简称“医保目录”)所划分的国谈抗肿瘤药的配备机构(简称“配药机构”)作为研究对象,借助空间视角的科学量化评价方法,评估和测量研究区域内各需求单元的配药机构可达性水平和空间分布情况。
1 资料与方法 1.1 研究区域研究区域为武汉市,武汉市下辖共计13个行政区和6个功能区,本文的主要研究区域为其中心城区,包括江岸区、江汉区、硚口区、汉阳区、武昌区、青山区和洪山区。
1.2 数据来源 1.2.1 供给点数据本研究的供给点单位为配药机构,包含定点医疗机构和定点零售药店两类。武汉市中心城区的具体配药机构名单及其配备的抗肿瘤药物品种数量均来源于中国国家医保服务平台官网(https://fuwu.nhsa.gov.cn/nationalHallSt/#/home)。经整理,共获得42个定点医疗机构(含不同院区)和69个定点零售药店的相关信息。其中,配药机构经纬度地理坐标同时借助高德地图Web服务API (Application Programming Interface)(https://lbs.amap.com/)和高德地图坐标拾取器(https://lbs.amap.com/tools/picker)请求和解析数据,并对二者返回的数据加以比对,误差超过一定阈值时重新进行人工检索。度量配药机构的国谈抗肿瘤药供给能力的价格信息主要来自2021年医保目录中的支付标准数据和各省份的官方药械集中采购平台。
需求点的潜在需求水平由WorldPop网站(https://hub.worldpop.org/)提供的最新2020年100m分辨率开源人口分布栅格数据所表征。该数据是基于随机森林的dasymetric再分配方法所做出的人口分布预测,具有精度较高且较易获取的特点,适合本研究需求。我们根据第七次全国人口普查的武汉市分区数据对人口栅格数据进行修正,同时结合湖北省癌症防治中心于2022年公布的最新的全省肿瘤登记地区恶性肿瘤发病率数据进行估算,最后借助ArcGIS软件,创建1km×1km的渔网,合并小单元人口栅格,构建武汉市中心城区潜在肿瘤患病人群的空间分布预测数据集。此外,我们以需求点的几何质心即1km×1km的渔网中心点计算经纬度坐标。
1.2.3 OD成本数据根据供给点和需求点的经纬度坐标信息,我们通过高德地图路径规划API测算始发地和目的地(Origin-Destination, OD)两点之间的交通出行距离,作为对供给点和需求点通行成本的度量。考虑到潜在的需求群体为体质较为虚弱的肿瘤患者,且一般需求单元距离供给单元存在较长的通行距离,并且一般不会仅在最邻近的配药机构获取服务,因此我们仅采用最为常用的驾车路径规划API和公交路径规划API测算了两种交通方式的出行距离,并将患者的空间搜索阈值设置为所有中心城区供需点之间的最长路网极限距离,而将步行、骑行等交通出行方式排除在外。
1.2.4 其他数据本研究所使用的武汉市行政区划数据主要来自全国地理信息资源目录服务系统(https://www.webmap.cn/),路网数据来自OpenStreetMap的开源GIS矢量数据集(https://www.openstreetmap.org/),包括高速公路、国道、省道、县道以及城市快速路、主干路、次干路、支路等信息。
1.3 研究方法 1.3.1 基于高斯函数的两步移动搜索法两步移动搜索法(2-Step Floating Catchment Area, 2SFCA)可以同时考虑医药服务供给侧和需求侧两个方面的作用因素,以此计算搜索阈值范围内的供需比例,其可以较为全面和客观地测量医药服务可达性,具有较强的适用性和操作性。[8-11]该方法的主要计算思路可分为两步:
第一步,以配药机构j为中心,以人们前往配药机构的路网极限距离d0为搜索半径建立搜索域j,搜索供给阈值d0范围内的所有需求点k(本文的最小需求单元为1km×1km渔网),并汇总相应的域内人口数量Dk,以此计算j点的供给规模Sj和k点的需求规模Dk的供需比Rj:
(1) |
第二步,以需求点i为中心,搜索阈值d0范围内的所有供应点j,将搜索到的所有j点的供需比累加,即得到i点的国谈药配备医药机构的可达性评分Ai:
(2) |
至此,传统两步移动搜索法的展开思路和计算过程已经得到明确,但进一步来看,此种方法也存在一定的内生缺陷,即未能考虑距离衰减作用和不同配药机构服务能力差异对患者就医购药行为选择的影响。为弥补方法学上的第一个共性不足,有学者提出应该在原有算法上引入距离衰减“惩罚项”,即在形式上对2SFCA的搜索半径部分乘上距离衰减函数,其一般形式为[12-13]:
(3) |
式中,f(dij)为一般的距离衰减函数形式,当其处于搜索半径之内时为g(dij),超过搜索半径时则为常数0。本研究选择高斯函数作为2SFCA搜索半径内的距离衰减函数,其暗含的意义在于,配药机构的可达性衰减速率呈现出一个随距离的增加先边际递增而后边际递减过程,具体形式为:
(4) |
即在高斯衰减函数的基础上,式(2)可以改写为:
(5) |
此外,针对配药机构中的定点医疗机构的分析中,如果单一沿用前述改进后的2SFCA依旧还存在一定的适用性风险,即不同配药机构的服务供给能力存在差异,特别是在药物储存运输、使用指导和患者随访评估等方面高水平配药机构更具优势,对使用患者的吸纳能力较强。基于此,有必要再对2SFCA算法进行特定优化,即在算法中另外增加对应配药机构权重因子项Wj:
(6) |
式中,Pjl为第j个配药机构配备的第l种抗肿瘤药物的医保支付价格,如果有多个规格的价格信息,则对这些多规格价格取平均值;经过加权计算后的最终结果即为配药机构的可达性指数。
1.3.2 空间自相关分析本研究主要使用全局莫兰指数I(Global Morans'I)来测量研究区域内配药机构空间可达性的关联性和分布特征。我们对渔网的属性数据进行了标准化处理之后,使用如下全局Morans'I的计算公式开展分析:
(7) |
式中,Yi和Yj分别表示i和j需求单元的配药机构可达性指数,n为地区总数,Wij为空间邻接矩阵。如果地区i与地区j相邻,则W=1,否则为0,特定地区预期自身的邻近关系为0。Morans'I的值域为[-1, 1],一定的统计学显著性水平上的正(负)值表示邻近的两个地区之间具有空间正(负)相关性,绝对值越大表明相关性越强,即相邻地区具有更高的属性相似度;而绝对值越小则表明这种相关性越弱,即邻近地区之间的差异越大,整体呈离散分布状态。
1.3.3 K-means聚类分析K-means聚类分析是一种迭代求解算法,其主要目的是将具有相似特征的数据进行归类,以便揭示数据中的隐含结构。在本研究中,我们采用K-means聚类分析方法来探讨研究区域内配药机构的可达性及其与供需组合模式的关联。通过从供应、需求和交通等方面深入了解可达性差异的根本原因,能够有助于在各个空间单元中找出主要影响可达性的限制性因素。在进行聚类分析时,主要考虑了以下几个指标:需求单元的配药机构可达性指数、潜在需求人口数量以及以需求单元为起点到各个配药机构的平均OD成本。这些指标综合反映了研究区域内配药机构的可达性状况、现实供给和潜在需求情况,为进一步的分析提供有力支持。为确保分析的准确性和可比性,在进行聚类之前,对所有指标进行了标准化处理。
1.3.4 软件实现本研究采用ArcGIS 10.8和QGIS 3.28地理信息系统软件对地理空间数据加以处理、分析和可视化呈现。在K-means聚类分析时,选用Python 3.10编程语言,并借助于其中的scikit-learn库展开测算。
2 结果 2.1 整体可达性评价结果如图 2所示,武汉市中心城区中西部区域需求单元的可达性明显要高于东部区域,大体上呈现出由一环线向二环线、三环线圈层递减的分布模式。在驾车或者公交出行的情况下,武汉市中心城区各个需求单元的配药定点医疗机构的可达性评价结果的平均值分别为97.625(SD=10.040)和97.068(SD=11.040),配药定点零售药店的可达性评价结果的平均值分别为36.792(SD=4.972)和35.972(SD=5.665),具有较为明显的空间分异特征。我们借助自然间断点方法将评价结果划分为很好、较好、一般、较差、很差五等级。以驾车情景为例,共有671个需求单元的配药定点医疗机构的可达性处于较好水平之上,占总需求单元的65.2%;而有615个需求单元的配药定点零售药店达到同样的层次,占总需求单元的59.7%;从两种配药机构来看,可达性不佳的需求单元数量占比在17.8%~20.2%。公交出行情景下的两种配药机构的可达性评价结果总体与前者类似(表 1)。
根据各区可达性的描述性统计结果(表 2),江汉区、武昌区、硚口区、汉阳区和江岸区的配药机构可达性明显较好,而青山区和洪山区的可达性排名则靠后。具体而言,江汉区在四种评价场景中,辖区内所有需求单元的配药机构可达性指数均超过武汉市中心城区的总体平均水平,平均可达性均为最高,分别达到了107.074~107.510和41.956~41.962;并且,在同种评价场景的横向比较中的标准差最小,分别仅为1.576~1.681和0.362~0.395,表明其区域内的可达性均等化水平最好。相比之下,排名最后的洪山区,其辖区内需求单元的定点医疗机构和定点零售药店的平均可达性水平仅有92.885~93.885和33.498~34.636,对应的标准差则高于其他行政区,说明其辖区内部的可达性较低且存在明显的分化,均等化水平最差。
此部分沿用基于自然间断点的武汉市中心城区总体可达性五等级划分数值界限,对各区的可达性等级加以归类。从驾车情景下的测算结果等级分布情况来看(图 3),各区需求单元处于配药定点医疗机构较高可达性水平及以上的占比排名依次为江汉区(100.0%)、汉阳区(98.6%)、硚口区(98.5%)、武昌区(98.3%)、江岸区(89.4%)、青山区(72.9%)、洪山区(48.2%);而在配药定点零售药店方面,各区的具体排名依次为江汉区(100.0%)、硚口区(100.0%)、汉阳区(99.3%)、武昌区(96.7%)、江岸区(95.2%)、青山区(71.4%)、洪山区(38.7%)。从图 3亦可以看出,江汉区、汉阳区等内部等级分化程度低,而青山区和洪山区内部等级分化程度高,特别是洪山区的需求单元在各等级区间内均有分布,这印证了前述的推断。
表 3展示了4种情景下的各需求单元的可达性评价空间自相关检验结果。在这4种情景中,Moran's I均十分接近于1,表明存在高度明显的空间聚类模式,同时也发现,在驾车出行情景下的Moran's I和Z-score相比公共交通出行情景下的数值有所降低,表明前者需求单元的配药机构可达性分布会更为均匀,而后者需求单元的可达性分布可能会出现较大的差异。
表 4显示了K-means聚类结果,其对应的空间分布情况如图 4所示。根据各项指标的聚类中心数值所在的三分位数区间,可以将武汉市中心城区的配药定点医疗机构和配药定点零售药店的可达性及其供需状况划分为3种类型:高需求低成本高可达性区域、中需求中成本中可达性区域(或高需求中成本中可达性区域)、中需求高成本低可达性区域。
无论是使用何种出行方式,配药定点医疗机构和配药定点零售药店的高需求低成本高可达性区域基本上分布在武汉市中心城区二环线以内及其附近、东湖的西北部区域(青山区与洪山区交界处)、东南部小范围区域,需求单元个数分别为148~150个和131~140个;中需求高成本低可达性区域主要集中在武汉市中心城区的东部以及西南行政边界处,需求单元个数分别为306~309个和328~357个;各需求单元的配药定点医疗机构的中需求中成本中可达性区域与配药定点零售药店的高需求中成本中可达性区域面积较大并且二者的耦合性同样较高,需求单元个数分别达到了571~576个和542~562个。此外,经过比较发现驾车情景下前两种分布类型相比公交情景下的范围明显更大。
3 讨论 3.1 空间可达性模型在评估药物配置方面具有应用潜力目前,药物的经济可负担性往往成为研究者所关注的重点,其空间可达性却未能给予应有的关注。本研究以武汉市中心城区为研究区域,针对2021年医保目录中所涉及的国谈抗肿瘤药的配药机构的空间布局进行分析,以链接估算需求患者数量的1km×1km渔网表征对国谈抗肿瘤药的潜在需求水平,采用经过改进的2SFCA开展可达性评价,并通过使用空间自相关分析和K-means聚类分析进一步识别了具体的区域分布模式和类型归属,所获得的一系列发现能够为优化国谈抗肿瘤药的区域配置提供有益思考。从方法学上讲,空间可达性模型能够科学评估研究区域内特定医疗资源与目标人群之间的供需关系。该方法可以为药物经济学和健康地理学研究提供有益借鉴。
3.2 武汉市中心城区配药机构可达性总体评估情况分析抗肿瘤药物价格高昂一直是限制中国内地医疗机构的配备动力以及广大患者的使用意愿和数量的主要因素之一。[6, 14]自2015年起,原卫生部开始探索和组织药品价格谈判,以减轻患者沉重的治疗负担,并释放了抗肿瘤药物的市场需求。[4, 15]同时,伴随着2021年国谈药“双通道”管理政策的扩散和推广,各省份的定点医疗机构和定点零售药店可以同时面向患者提供国谈抗肿瘤药,使其空间可及性得到显著提升。整体而言,在武汉市中心城区,配药机构的可达性呈现出明显的圈层扩散形态,其中中西部可达性明显高于东部。原因在于中西部是武汉城市传统核心区域,人口密度和潜在需求水平明显高于其他区域,该处也集中了大量的优质医疗资源,诸多较高水平的定点医疗机构和资质良好的定点零售药店分布于此,加之完善的基础设施和健全的交通网络能够极大缩小患者的出行成本,使得该区域的可达性水平远高于东部等城市边缘区域。
将武汉市中心城区的整体可达性结果分解到各个具体的辖区后可以发现,不同辖区内部的可达性等级分布不尽相同。无论是采取何种出行方式,江汉区、汉阳区等的可达性水平较高且内部的可达性分化程度较低,而青山区和洪山区的可达性水平较低且内部存在明显的分层现象。洪山区东部较为靠近山区,地形起伏较大,人口密度较低,路网也比其他辖区更为稀疏,大片区域除仅有的1个配备国谈抗肿瘤药的定点医疗机构外,尚未有设置其他相应的定点医疗机构和定点零售药店,这使得具有潜在需求的人群不得不前往医疗资源集中的中西部区域寻求就医和购药服务,大大降低了该地区需求点的配药机构可达性。同时,研究还发现,驾车出行和公交出行情景下,配药定点医疗机构和配药定点零售药店的可达性评价结果具有相似的空间分布特征。然而,在驾车出行情景下的可达性相对较好,这可能与私家车出行具有更高的出行灵活性、更短的出行时间和更低的出行成本有关。此外,公共交通出行情景下的可达性分布可能受到公交线路和站点布局的影响,这使得部分需求单元的可达性出现了较大差异。[16]
3.3 配药机构可达性的空间分布模式及其供需关系组合本研究的空间自相关分析结果表明,武汉市中心城区配药定点医疗机构和配药定点零售药店的可达性具有显著正相关的空间分布模式,区域内高值集聚与分区内部分化现象并存。同时,K-means聚类分析将武汉市中心城区的配药定点医疗机构和配药定点零售药店的可达性及其供需状况之间的关系模式分别划分为三类。配药定点医疗机构和配药定点零售药店的高需求低成本高可达性区域范围较小,中需求高成本低可达性区域范围较大,而中需求中成本中可达性区域或高需求中成本中可达性区域的范围最大,这些分类结果与武汉市中心城区全域可达性评价结果具有较高的耦合度。根据所划分的关系模式,应着重关注中需求高成本低可达性区域,该类区域配药机构数量十分有限,人口密度低且交通条件差,各方面发展水平处于相对落后的状态,未来随着潜在需求人口的增加,该区域势必会面临更为严峻的配药机构可达性提升压力。
4 建议 4.1 完善和落实国谈药“双通道”管理机制持续推进国家医保谈判工作和完善国谈药“双通道”管理机制是提高抗肿瘤国谈药可获得性的基础性手段。一般而言,国谈抗肿瘤药的配备和使用对于医疗卫生机构的等级水平和服务能力有着更高的要求,需要有力压实定点医疗机构的谈判药品配备主体责任,鼓励其因时因势、适度适量地增加抗肿瘤药物的供给,特别是对于一些需求量较大、对医疗机构配备和使用无过多特殊要求的药品可以更多借助低等级或基层医疗机构延伸供给链条。同时,需要发挥定点零售药店分布广泛、市场化程度高、服务灵活的优势,遴选资质合规、管理规范、信誉良好的定点零售药店作为国谈抗肿瘤药供给的重要补充渠道,扩大定点零售药店的分布范围,按照便民、利民的原则推进配药机构布局的科学化和合理化,满足不同区域患者的就医和购药需求。
4.2 优化城市交通基础设施以提升通行效率城市交通基础设施的优化在增强患者对抗肿瘤药物的获取能力方面发挥着重要作用。一套高效、连贯且覆盖广泛的交通系统,能大大降低患者前往配药机构的时间和空间成本,从而提高抗肿瘤药物的空间可达性。需要进一步完善城市基础设施建设,形成多样化、互补性强的城市交通系统,特别是进一步优化公共交通,增强其延展性和连通性。尤其对于那些偏远、人口稀少的地区,提高公共交通的运行效率、覆盖范围和使用便利性,将有助于显著减少患者就医和购药的通行距离和时间成本。
4.3 提升全民健康水平以减少潜在药物需求政府部门和社会各界需要重视公众健康教育,强化医学知识和保健知识宣传,引导公众形成自我“健康第一责任人”的观念和健康规律的生活方式。同时,提升健康体检和筛查的覆盖率,有助于个体早期发现和处理潜在的健康问题。通过这些综合手段,控制和降低肿瘤的全人群患病率,进一步增强全民体质,从而减少对抗肿瘤药物的需求。
5 研究局限性一方面,本研究用于表征人口需求的数据来自开源的预测数据集,即使经过了两次修正,但与真实世界中的武汉市中心城区需求人口数量势必有所误差;另一方面,由于具体地划分各肿瘤病种及测量对应的人口分布情况在当前现有技术手段上存在高度的困难,所以本研究仅利用单一化的口径将配药机构所配备的所有病种的抗肿瘤药物作为其供给能力的测量,未能做到进一步的精细化处理。
作者声明本文无实际或潜在的利益冲突。
[1] |
Wei W Q, Zeng H M, Zheng R S, et al. Cancer registration in China and its role in cancer prevention and control[J]. Lancet Oncol, 2020, 21(7): E342-E349. DOI:10.1016/S1470-2045(20)30073-5 |
[2] |
Qiu H B, Cao S M, Xu R H. Cancer incidence, mortality, and burden in China: a time-trend analysis and comparison with the United States and United Kingdom based on the global epidemiological data released in 2020[J]. Cancer Commun, 2021, 41(10): 1037-1048. DOI:10.1002/cac2.12197 |
[3] |
Xia C F, Dong X S, Li H, et al. Cancer statistics in China and United States, 2022:profiles, trends, and determinants[J]. Chin Med J (Engl), 2022, 135(5): 584-590. DOI:10.1097/CM9.0000000000002108 |
[4] |
Zhang Y C, Wushouer H, Han S, et al. The impacts of government reimbursement negotiation on targeted anticancer medication price, volume and spending in China[J]. BMJ Glob Health, 2021, 6(7): 11. |
[5] |
Tang M, Song P, He J. Progress on drug pricing negotiations in China[J]. Bioscience Trends, 2019, 13(6): 464-68. DOI:10.5582/bst.2019.01339 |
[6] |
Zhang Y C, Wei Y X, Li H Q Y, et al. Prices and Clinical Benefit of National Price-Negotiated Anticancer Medicines in China[J]. Pharmacoeconomics, 2022, 40(7): 715-724. DOI:10.1007/s40273-022-01161-7 |
[7] |
刘辉, 鄢玉杭, 蒋倩. 国家谈判药品"双通道"管理模式及实践探析[J]. 中华医院管理杂志, 2022, 38(9): 651-656. |
[8] |
Luo W, Wang F H. Measures of spatial accessibility to health care in a GIS environment: synthesis and a case study in the Chicago region[J]. Environ Plan B-Plan Des, 2003, 30(6): 865-884. DOI:10.1068/b29120 |
[9] |
Guagliardo M F. Spatial accessibility of primary care: concepts, methods and challenges[J]. International journal of health geographics, 2004, 3(1): 1-13. DOI:10.1186/1476-072X-3-1 |
[10] |
Mohammadi A, Mollalo A, Bergquist R, et al. Measuring COVID-19 vaccination coverage: an enhanced age-adjusted two-step floating catchment area model[J]. Infect Dis Poverty, 2021, 10(1): 13. DOI:10.1186/s40249-021-00801-y |
[11] |
Bryant J, Delamater P L. Examination of spatial accessibility at micro-and macro-levels using the enhanced two-step floating catchment area (E2SFCA) method[J]. Ann GIS, 2019, 25(3): 219-229. DOI:10.1080/19475683.2019.1641553 |
[12] |
Wang F H. Measurement, Optimization, and Impact of Health Care Accessibility: A Methodological Review[J]. Ann Assoc Am Geogr, 2012, 102(5): 1104-1112. DOI:10.1080/00045608.2012.657146 |
[13] |
Hu S J, Song W, Li C G, et al. A multi-mode Gaussian-based two-step floating catchment area method for measuring accessibility of urban parks[J]. Cities, 2020, 105: 10. |
[14] |
Luo X, Du X, Huang L, et al. The price, efficacy, and safety of within-class targeted anticancer medicines between domestic and imported drugs in China: a comparative analysis[J]. The Lancet Regional Health-Western Pacific, 2023, 32. |
[15] |
Huang C, Ung C O L, Wushouer H, et al. Trends of Negotiated Targeted Anticancer Medicines Use in China: An Interrupted Time Series Analysis[J]. International Journal of Health Policy and Management (IJHPM), 2021, 47. |
[16] |
Stentzel U, Piegsa J, Fredrich D, et al. Accessibility of general practitioners and selected specialist physicians by car and by public transport in a rural region of Germany[J]. BMC Health Serv Res, 2016, 16: 10. |
(编辑 薛云)