Matlab 0-1规划求解八皇后问题

这是我以前写的一个程序,发布在一些Matlab的论坛,现在搬过来吧

clear;clc;close all
N=8;
c=ones(N);
% 行求和
blkele=num2cell(c,2);
A1=blkdiag(blkele{:});
% 列求和
A2=repmat(eye(N),1,N);
Aeq=[A1;A2];
beq=ones(2*N,1);
% 斜线情况判断
M=N-2;
A=zeros(4*M+2,N*N);
d{1}=arrayfun(@(i)find(diag(diag(c,i),i)),-M:M,'UniformOutput',0);
d{2}=arrayfun(@(i)find(fliplr(diag(diag(c,i),i))),-M:M,'UniformOutput',0);
d0=cat(1,d{:});
linearIndex=arrayfun(@(x)sub2ind(size(A),x+zeros(1,length(d0{x})),...
    d0{x}'),1:numel(d0),'UniformOutput',0);
A(cell2mat(linearIndex))=1;
b=ones(4*M+2,1);
x=bintprog([],A,b,Aeq,beq);
Res=reshape(x,[],N);


 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值