找回密码
 立即注册

QQ登录

只需一步,快速开始

手机号码,快捷登录

手机号码,快捷登录

Flux电机有限元分析教程西莫团购入口 | 奖励入口当当网购物入口 | 奖励入口欢迎使用!西莫论坛App开放下载Motor-CAD电机多物理域设计教材购买入口 | 奖励入口
★新会员论坛须知★《西莫电机技术》第39期发售火热进行中Flux电机电磁阀有限元分析教程团购入口 | 奖励入口论坛微信公众平台欢迎入驻
西莫电机及相关产品供需交流群开放邀请★ 论坛VIP会员申请 ★Motor-CAD.MANATEE电磁热振动噪声教程 | 奖励入口西莫团队欢迎您的加盟!
宣传推广合作请联系QQ:25941174西莫电机论坛微信群正式开放Flux变压器与电抗器有限元分析团购入口 | 奖励入口西莫电机论坛技术版区QQ群汇总
查看: 663|回复: 0

[原创] 自己编写的八皇后问题解

[复制链接]

该用户从未签到

发表于 2016-3-15 22:28 | 显示全部楼层 |阅读模式 来自: 中国湖南湘潭

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?立即注册

×
本帖最后由 hhbhhy 于 2016-3-15 22:29 编辑

八皇后问题描述:在国际象棋盘中,8*8格,放置8个皇后,使得每行,每列,每对角斜线上都只有一个皇后。问一共有多少种可能的排布方法。

程序如下,高手指点!
%N个皇后排座问题
N=8;
%方案个数
M=0;

%以一为不可使用位置。
for i1=1:N
       %初始化
       huanghou=zeros(N,N);
       %设定第一位皇后位置
       huanghou(i1,1)=1;
       %列可以不管,设置本行后面元素全为1,
       huanghou(i1,:)=1;
       %对角排除不可放的位置
       %斜向上
       if i1>1
           for g1=1:i1-1;
               huanghou(i1-g1,1+g1)=1;
           end
       end
       %斜向下
       for g1=1:N-i1
           huanghou(i1+g1,1+g1)=1;
       end

       %第二步,与第一步类似。
       for i2=1:N
          %承接第一步的结果
          huanghou2=huanghou;
          %如果本列相加为维数N,说明没有位置可放皇后,则这种方案放弃,直接返回第一步。
          if sum(huanghou2(:,2))==N
              break;
          end
          %如果(i2,2)位置上为1,不可放皇后,则跳到下一个i2进行尝试。
          if huanghou2(i2,2)==1
              continue;
          end
          %找到第二个皇后的位置。下面与第一步相同,不再解释。
          huanghou2(i2,2)=2;
          huanghou2(i2,3:N)=1;
           if i2>1
               for g1=1:min(N-2,i2-1);
                  huanghou2(i2-g1,2+g1)=1;
               end
           end
          for g1=1:min(N-i2,N-2)
           huanghou2(i2+g1,2+g1)=1;
          end

         

         %第三步         
         for i3=1:N
          huanghou3=huanghou2;
          if sum(huanghou3(:,3))==N
              break;
          end
          if huanghou3(i3,3)==1
              continue;
          end
          huanghou3(i3,3)=3;
          huanghou3(i3,4:N)=1;
           if i3>1
               for g1=1:min(N-3,i3-1);
                  huanghou3(i3-g1,3+g1)=1;
               end
           end
          for g1=1:min(N-i3,N-3)
           huanghou3(i3+g1,3+g1)=1;
          end

         %第四步         
         for i4=1:N
          huanghou4=huanghou3;
          if sum(huanghou4(:,4))==N
              break;
          end
          if huanghou4(i4,4)==1
              continue;
          end
          huanghou4(i4,4)=4;
          huanghou4(i4,5:N)=1;
           if i4>1
               for g1=1:min(N-4,i4-1);
                  huanghou4(i4-g1,4+g1)=1;
               end
           end
          for g1=1:min(N-i4,N-4)
           huanghou4(i4+g1,4+g1)=1;
          end

         %第五步         
         for i5=1:N
          huanghou5=huanghou4;
          if sum(huanghou5(:,5))==N
              break;
          end
          if huanghou5(i5,5)==1
              continue;
          end
          huanghou5(i5,5)=5;
          huanghou5(i5,6:N)=1;
           if i5>1
               for g1=1:min(N-5,i5-1);
                  huanghou5(i5-g1,5+g1)=1;
               end
           end
          for g1=1:min(N-i5,N-5)
           huanghou5(i5+g1,5+g1)=1;
          end


         %第六步         
         for i6=1:N
          huanghou6=huanghou5;
          if sum(huanghou6(:,6))==N
              break;
          end
          if huanghou6(i6,6)==1
              continue;
          end
          huanghou6(i6,6)=6;
          huanghou6(i6,7:N)=1;
           if i6>1
               for g1=1:min(N-6,i6-1);
                  huanghou6(i6-g1,6+g1)=1;
               end
           end
          for g1=1:min(N-i6,N-6)
           huanghou6(i6+g1,6+g1)=1;
          end

         %第七步         
         for i7=1:N
          huanghou7=huanghou6;
          if sum(huanghou7(:,7))==N
              break;
          end
          if huanghou7(i7,7)==1
              continue;
          end
          huanghou7(i7,7)=7;
          huanghou7(i7,8:N)=1;
           if i7>1
               for g1=1:min(N-7,i7-1);
                  huanghou7(i7-g1,7+g1)=1;
               end
           end
          for g1=1:min(N-i7,N-7)
           huanghou7(i7+g1,7+g1)=1;
          end
         

         %最后一步         
         for i8=1:N
          huanghou8=huanghou7;
          if sum(huanghou8(:,8))==N
              break;
          end
          if huanghou8(i8,8)==1
              continue;
          end
          %找到最后一个皇后的位置
          huanghou8(i8,8)=8;
          M=M+1;
          %输出:
          disp(['第' num2str(M) '个方案'])
          huanghou8

%最后一步结束
end


%第七步结束
end


%第六步结束
end



%第五步结束
end

%第四步结束

end



%第三步结束
end
      
         
  %第二步结束        
         
end



  %第一步结束  
end


西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

西莫电机论坛微信公众平台欢迎您的关注!

QQ|Archiver|手机版|小黑屋|西莫电机圈 ( 浙ICP备10025899号-3|浙公网安备:33028202000436号 )

GMT+8, 2024-12-23 02:23 , Processed in 0.040673 second(s), 22 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表