lingo 求修改,没有可行解,

1个回答

  • sets:

    x05dy/1..21/;

    x05nl/1..7/:qz;

    x05dw/1..7/:d;

    x05link(dy,nl):C;

    links(nl,dw):Z;

    x05DD(dy,nl,dw):X;

    x05link2(dy,dw):y;

    endsets

    data:

    qz=0.350396277663283 0.237472604788531 0.158965769190678 0.105558337909803 0.069645448792050 0.046163178982660 0.031798382672994;

    C=86 90 71 80 78 88 93

    82 87 91 70 85 91 90

    80 75 88 93 86 78 78

    90 88 95 79 77 82 90

    87 83 90 86 65 84 91

    90 91 73 67 77 89 81

    88 87 85 86 90 94 85

    91 88 70 65 61 78 84

    73 89 91 87 88 78 80

    87 88 90 91 83 84 89

    80 87 72 90 86 84 79

    87 90 80 82 79 90 80

    95 78 88 90 76 78 63

    93 92 77 82 84 90 88

    78 85 84 89 73 88 83

    88 90 79 83 87 76 80

    95 78 90 91 78 83 85

    88 77 89 85 83 90 78

    92 89 70 75 78 92 88

    89 93 75 80 87 89 89

    95 88 76 80 74 86 78;

    enddata

    max=@sum(nl(j):

    x05x05qz(j)*@sum(dy(i):

    x05x05x05x05x05C(i,j)*@sum(dw(k):

    x05x05x05x05x05x05x05X(i,j,k))));

    @for(DD:@bin(X));

    @for(dw(k):@sum(link(i,j):X(i,j,k))=7);

    @for(dy(i):@sum(dw(m):y(i,m))