数独问题求解三:数独终局生成

数独问题求解三:数独终局生成

1、主要过程

由上一篇文章可知,本题中解决数独终局的问题的关键是对1e6个数独进行生成
在解决的过程中,我设置了一个全局变量field数组来存储数独
由于最多生成1e6个数独,每个数独由81个数字组成,每一行的最后一位是换行符,因此设置了数组大小为1e8.
数组中数独的数字用字符的形式保存

在这里插入图片描述

以下为主要代码段
这段代码为对数独的4到6行进行全排列,共有6种排列方式,因此循环6次,每次进入循环时都要判断是否已经生成了足够的数独。

for (int i = 0; i < 6; i++)//4到6行的六种排列方式
 {
  if (flag >= n)
  {
   break;
  }
  if (i)
  {
   //交换4到6行
   next_permutation(trans + 3, trans + 6);     //全排列函数,输出所有比当前排列大的排列,顺序是从小到大。
   trans[6] = 2, trans[7] = 5, trans[8] = 8;
  }

这段与上段思路相同,为对7到9行进行排列

for (int j = 0; j < 6; j++)//7到9行的六种排列方式
  {
   if (flag >= n)
    break;
   if (j)
    //交换7到9行
    next_permutation(trans + 6, trans + 9);   //全排列函数

先手动输入第一行,根据题设,我生成的数独第一个数字应该是2,然后对第一行进行全排列,共有8!=40320种方式

char row[10] = "213456789";
   //先生成第一种情况
   for (int k = 0; k < 40320; k++)//第一行的8!种排列方式
   {
    if (flag >= n)
     break;
    if (k)
     next_permutation(row + 1, row + 9);//第一个数字不能换

将每一种排列方式进行变换,利用trans数组(trans[9] = { 0,3,6,1,4,7,5,2,8 })进行,并将数独存放到数组中,flag++。

for (int i = 0; i < 9; i++)
    {
     for (int j = 0; j < 9; j++)
     {
      field[m + flag] = row[(j + trans[i]) % 9];
      m += 1;
     }
     field[(m - 1) + flag] = '\n';
     //m += 1;
    }
    field[m + flag] = '\n';
    flag++;

2、结果

为了知道程序运行的时间,我加入了clock()函数

printf("\n用时:%.2lfs", (double)clock() / CLOCKS_PER_SEC);//调用clock()函数得出运行时间

在这里插入图片描述
代码运行后,可以在约2s的时间内生成1e6个数独终局
通过对比,程序运行效率良好

3、一些问题

在编码的过程中,我遇到了诸多问题,虽然都很蠢但是还是在此列举出来,避免踩坑
1)在声明field数组时没有考虑到数组过大造成的栈溢出,后改为申请全局变量
2)在使用Visual Studio2017编码时,遇到了无法查找或打开 PDB 文件的报错,该问题可通过对编译器参数修改解决
3)在解决输出格式时,忽略了程序结束后的换行问题,导致输出field数组时,总是只输出第一个数组

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
一个简单的解数独的小程序 //***求数独的解,参数mod为0或1,time0为搜索开始时的时间,mod=0时仅检查Data1中数独是否有解,有解则抛出1,mod=1时求出所有解并输出*** { int i,j,im=-1,jm,min=10; int mark[10]; for(i=0;i<9;i++) { for(j=0;j<9;j++) { if(Data1[i][j]) //如果该位置有数据则跳过 { continue; } int c=Uncertainty(i,j,mark); //如果该位置为空则先求不确定度 if(c==0) //如果不确定度为0则表示该数独无解 { return; } if(c<min) //得到不确定度最小的位置(第im行 第jm列 不确定度为min) { im=i; jm=j; min=c; } } } if(im==-1) //所有位置都已经确定,数独已经解出,按要求输出解 { if(mod==1) //显示所有解 { if(IsSolved()==true) { if(Solutions++<MAXANSNUM) { cout<<"第 "<<Solutions<<" 个 解:"<<endl; Display(1); } if((time(NULL)-time0)>TIMEOUT) { throw(Solutions); } } return; } else //只给出一个解 { throw(1); //跳出所有递归调用,返回1 } } Uncertainty(im,jm,mark); //从不确定度最小的位置开始解 for(i=1;i<=9;i++) { if(mark[i]==0) { Data1[im][jm]=i; //对不确定度最小的位置尝试可能的赋值 Search(mod,time0); //递归调用 } } Data1[im][jm]=0; } void Csudoku::Set(int n) //***随机生成数独,参数n表示数独中待填元素个数*** { srand((unsigned)time(NULL)); int i,j,k; do { for(i=0;i<9;i++) //随机给每行的某一个位置赋值 { for(j=0;j<9;j++) { Data1[i][j]=0; } j=rand()%9; Data1[i][j]=i+1; } } while(!Solve(0)); //按照随机赋的值给出一个解 for(k=0;k<n;) //从中随机去掉n个数据 { i=rand()%81; j=i%9; i=i/9; if(Data1[i][j]>0) { Data1[i][j]=0; k++; } } for(i=0;i<9;i++) //将生成数独存入Data0数组 { for(j=0;j<9;j++) { Data0[i][j]=Data1[i][j]; } } }
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值