肾脏供者交换的周期检测（onete02）
/****************************************************************//* *//*S A S S A M P L E L I B R A R Y*//* *//*姓名：onete02*//*标题：肾脏供者交换的周期检测（onete02）*//*产品：或*//*系统：所有*//*关键字：或*//*程序：OPTNET、OPTMODEL、PRINT*//*数据：*//* *//*支持：更新：*//*参考：*//*MISC：OPTNET文档中的示例2*//* *//****************************************************************//*用弧概率p在n个节点上创建随机图和均匀（0,1）重量*/%设n＝100；%设p=0.02；数据链接设置；do从=0到&n-1；do to=0 to&n-1；如果从等式到，则继续；否则，如果拉努尼（1）<&p则会；重量=拉努尼（2）；产出；结束；结束；结束；运行；/*生成具有2<=长度<=最大长度的所有循环*/%设max_length=10；程序optnetloglevel=中等graph_direction=定向data_links=LinkSetIn；周期最小长度=2最大长度=最大长度（&M）out=循环模式=所有循环；运行；%放置&_ROPTNET_；%放置&_ROPTNET_CYCLE_；/*将周期转换为每个弧的一个观测值*/数据周期0（保持=c i j）；设置周期；保留最后一个；c=循环；i=最后；j=节点；最后=j；如果命令ne为1，则输出；运行；/*求解集合装箱问题寻找最大权节点不交并短定向循环*/过程操作模型；/*声明索引集和参数，并读取数据*/设置<num，num>ARCS；num权重{ARCS}；将数据LinkSetIn读入ARCS=[from to]weight；设置<num，num，num>TRIPLES；将数据Cycles0读入TRIPLES=[c i j]；集合CYCLES=TRIPLES}c中的集合{c，i，j>；在CYCLES}中设置ARCS_c{c=在TRIPLES}中的{<（c），i，j>集；在CYCLES}中设置NODES_c{c=ARCS_c[c]}{i，j}中的并集{；设置NODES=联合{c in CYCLES}NODES_c[c]；num cycle_weight｛c in CYCLES｝=ARCS_c[c]}中的sum｛＜i，j＞weight〔i，j〕；/*如果使用循环c，则UseCycle[c]=1，否则为0*/var UseCycle{CYCLES}二进制；/*宣布目标*/最大总重量=CYCLES}中的sum{c cycle_weight[c]*UseCycle[c]；/*每个节点最多出现一个周期*/con node_packing{i在NODES}中：CYCLES中的sum{c：NODES_c[c]}中的i UseCycle[c]<=1；/*调用解算器*/用milp求解；/*输出最优解*/从创建数据解决方案[c] ={CYCLES中的c：UseCycle[c].sol>0.5}cycle_weight；退出；%放置&_OROPTMODEL_；proc print data=解决方案noobs标签；运行；