求解数独Fortran程序
!*******************************************************************
!编写：飞行棋 编写时间：2013年儿童节
!程序说明：
!主程序采用回溯法，先按顺序将每个未知数的位置坐标存入x、y数组中，
!共有step个未知数。每个位置从1开始试验，若满足条件则t=t+1，试验
!下一个未知数；若1~9均不能满足条件，则回溯到上一个未知数t=t-1，
!将试验值+1，再进行循环。当循环到最后一个未知数时，若满足条件则输
!出解；若回溯到第一个未知数时，填入1~9均不能满足条件，则数独无解。
!子程序用来判断试验值是否满足条件。
!m,n=未知数在校宫格中的位置，step=未知数个数，t=试验值在未知数中
!的次序,A=存放数独值的数组，x,y=存放未知数坐标的数组
!********************************************************************
program sudoku
implicit none
integer i,j,m,n,k,step,t
integer A(1:9,1:9),x(81),y(81)
logical flag
t=0
open(8,file='soku.txt')
do i=1,9 !读入数独初值，未知数赋0
do j=1,9
read(8,'(i2)') A(i,j)
if(A(i,j)==0)then
t=t+1
step=t !计算未知数个数
x(t)=i
y(t)=j
endif
enddo
enddo
close(8)
print '(2x,9i2)',((A(i,j),j=1,9),i=1,9)
t=1
print*,step
do while(t<=step)
i=x(t)
j=y(t)
if(A(i,j)==9)then !若回溯到的上一个未知数值为9，则继续往前回溯
t=t-2
A(i,j)=0
goto 200
endif
do k=A(i,j)+1,9
A(i,j)=k
call judge(i,j,k,flag,A) !判断是否满足条件,若满足，退出循环，试验下一个未知数
if(flag.eqv..true.)then
exit
endif
if(k==9.and.(flag.eqv..false.))then !若1~9均不能满足条件，回溯上一个未知数
t=t-2
A(i,j)=0
endif
enddo
200 t=t+1
if(t==step+1)then
print*,'本数独解为：'
exit
endif
if(t<=0.and.(flag.eqv..false.))then
print*,'本数独无解,再检查一下题目吧！'
exit
endif
enddo
print '(2x,9i2)',((A(i,j),j=1,9),i=1,9)
10 format(1x,3I3)
End
子程序：
subroutine judge(vi,vj,k,flag,A)
implicit none
integer i,j,vi,vj,k,m,n
logical flag
integer A(1:9,1:9)
flag=.true.
m=mod(vi,3)+3*int(1-mod(vi,3)/3.)
n=mod(vj,3)+3*int(1-mod(vj,3)/3.)
do i=1,9 !判断行条件
if(A(i,vj)==k.and.i/=vi)then
flag=.false.
endif
enddo
do j=1,9 !判断列条件
if(A(vi,j)==k.and.j/=vj)then
flag=.false.
endif
enddo
do i=vi-m+1,vi-m+3 !判断小宫格条件
do j=vj-n+1,vj-n+3
if((i/=vi.or.j/=vj).and.A(i,j)==k)then
flag=.false.
end if
enddo
enddo
end
输入数据格式示例：
0
0
0
0
0
0
8
0
0
4
0
0
2
0
8
0
5
1
0
8
3
9
0
0
0
0
7
0
4
0
5
0
0
0
8
2
0
0
5
0
0
0
4
0
0
8
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
每行一个数据，共82行。
输出示例：
1 2 6 3 5 7 8 4 9
4 7 9 2 6 8 3 5 1
5 8 3 9 1 4 2 6 7
3 4 1 5 7 6 9 8 2
2 6 5 8 9 1 4 7 3
8 9 7 4 2 3 5 1 6
6 3 2 7 4 5 1 9 8
7 5 8 1 3 9 6 2 4
9 1 4 6 8 2 7 3 5
本文档为【数独求解程序】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。
该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。
[版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
浙公网安备 33021202002483