两个程序比较，为什么第一个通不过呢？

bei_fang

程序1
model:
sets:
node /supply,res_1,res_2,simple_1,simple_2,simple_3,simple_4,simple_5,outlet/;
nl (node)/res_1,res_2/;
month /Jan, Feb, Mar,Apr, May, Jun, Jul, Aug, Sep, Oct, Nov, Dec/:sup1,sup2,demand_dom,demand_ind,demand_agr,divert_dom,divert_ind,divert_agr,out,int;

endsets
sets: net_flow(node,node)
/supply,res_1
supply,res_2
res_1,simple_1
res_2,simple_1
simple_1,simple_2
simple_2,simple_3
simple_3,simple_4
simple_4,simple_5
simple_5,outlet
outlet,supply/;
endsets
sets:
net_time(net_flow,month):waterflow;
endsets
data:
sup1=128 125 234 290 350 400 500 440 240 210 181 128;   !水源1的入流量;
sup2=39 39  52  121 168 144 105 78  49  44  45 39;      !水源2的入流量;
demand_dom=20 20 20 20 20 20 20 20 20 20 20 20;
demand_ind=40 40 40 40 40 40 40 40 40 40 40 40;
demand_agr=60 60 60 60 60 60 60 60 60 60 60 60;

enddata

@for(net_time(i,j,t)|i#eq#1 #and# j#eq#2:waterflow(i,j,t)=sup1(t));     
@for(net_time(i,j,t)|i#eq#1 #and# j#eq#3:waterflow(i,j,t)=sup2(t));     

max=@sum(month(t):divert_dom(t)+divert_ind(t)+divert_agr(t));

@for(month(t):divert_dom(t)>=demand_dom(t));
@for(month(t):divert_ind(t)>=demand_ind(t));
@for(month(t):divert_agr(t)>=demand_agr(t));
@for(month(t):int(t)=@sum(net_time(j,i,t)|i#eq#5:waterflow(j,i,t)));
@for(month(t):out(t)=divert_dom(t)+divert_ind(t)+divert_agr(t));
@for(month(t):int(t)>=out(t));


@for(month(t):
    @for(node(i):@sum(net_time(i,j,t):waterflow(i,j,t))=
            @if(i#lt#5,@sum(net_time(j,i,t):waterflow(j,i,t)),
            @if(i#eq#5,@sum(net_time(j,i,t):waterflow(j,i,t))-divert_dom(t)-divert_ind(t)-divert_agr(t),@sum(net_time(j,i,t):waterflow(j,i,t))
                  )
                )
                )
      );
end


程序2
model:
sets:
month /Jan, Feb, Mar,Apr, May, Jun, Jul, Aug, Sep, Oct, Nov, Dec/:sup1,sup2,demand_dom,demand_ind,demand_agr,divert_dom,divert_ind,divert_agr,out,int;
endsets

data:
sup1=128 125 234 290 350 400 500 440 240 210 181 128;   !水源1的入流量;
sup2=39 39  52  121 168 144 105 78  49  44  45 39;      !水源2的入流量;
demand_dom=20 20 20 20 20 20 20 20 20 20 20 20;
demand_ind=40 40 40 40 40 40 40 40 40 40 40 40;
demand_agr=60 60 60 60 60 60 60 60 60 60 60 60;
enddata

max=@sum(month(t):divert_dom(t)+divert_ind(t)+divert_agr(t));

@for(month(t):divert_dom(t)>=demand_dom(t));
@for(month(t):divert_ind(t)>=demand_ind(t));
@for(month(t):divert_agr(t)>=demand_agr(t));
@for(month(t):out(t)=divert_dom(t)+divert_ind(t)+divert_agr(t));
@for(month(t):int(t)=sup1(t)+sup2(t));
@for(month(t):int(t)>=out(t));
end

两个程序中，运行结果显示int（t）的值是相同的啊，为什么第一个通不过，而第二个却可以通过？

madio

可以看出来程序二是程序一去掉了一些约束的，程序一提示没有可行解，也就是这些被去掉的约束无法满足，你可以看看这些约束是什么意思，是不是和上面的约束相矛盾

bei_fang

程序1多了一个平衡条件，程序1中的int（t）与程序2中的int（t）的运行值是一样的，但不知道为什么程序1却没有可行解了？

magic2728

约束之间有相互矛盾的地方。