靶型数独.docVIP

type arr1=array[1..80,1..9] of boolean; 　　 arr2=array[1..80] of integer; 　　 rec=record 　　　　 r,len:integer; 　　　　 h:arr1; 　　　　 l:arr2; 　　 end; var h,ht:arr1; 　　l,lt:arr2; 　　h1,h2,h3:array[1..9,1..9] of boolean; 　　h4:array[1..80] of boolean; 　　a,b,fen:array[1..9,1..9] of integer; 　　v:array[1..80,1..2] of integer; 　　n,t,max,c:integer; 　　zh:array[0..80] of rec; function sma(x,y:integer):integer; begin if xy then sma:=x else sma:=y; end; procedure put; var p,i,j:integer; begin p:=c; for i:=1 to n do inc(p,a[v[i,1],v[i,2]]*fen[v[i,1],v[i,2]]); if pmax then max:=p; end; procedure ready; var i,j,x,y,min,p:integer; begin for i:=1 to 9 do for j:=1 to 9 do begin read(a[i,j]); 　　　　 b[i,j]:=(i-1) div 3*3+(j-1) div 3+1; 　　　　 fen[i,j]:=sma(sma(9-i+1,i),sma(9-j+1,j))+5; if a[i,j]=0 then begin 　 inc(n); v[n,1]:=i; v[n,2]:=j; 　　 end　else begin 　　　　 h1[i,a[i,j]]:=true; 　　　　 h2[j,a[i,j]]:=true; 　　　　 h3[b[i,j],a[i,j]]:=true; 　　　　 inc(c,a[i,j]*fen[i,j]); 　　　　 end; end; min:=9; p:=0; for i:=1 to n do begin x:=v[i,1]; y:=v[i,2]; for j:=1 to 9 do if (not h1[x,j]) and (not h2[y,j]) and (not h3[b[x,y],j]) then 　　 begin 　　 inc(l[i]); h[i,j]:=true; 　　 end; if l[i]min then begin min:=l[i]; p:=i; end; end; t:=1; zh[t].r:=p; zh[0].h:=h; zh[0].l:=l; h4[p]:=true; end; procedure zhao; var min,p,i:integer; begin min:=9; p:=0; for i:=1 to n do if not h4[i] then if lt[i]min then begin min:=lt[i]; p:=i; end; inc(t); zh[t].r:=p; zh[t].len:=0; h4[p]:=true; end; procedure expend(z:integer); var x,y,i:integer; begin x:=v[zh[t].r,1]; y:=v[zh[t].r,2]; ht:=zh[t-1].h; lt:=zh[t-1].l; for i:=1 to n do if not h4[i] and ht[i,z] then 　　 if (v[i,1]=x) or (v[i,2]=y) or (b[v[i,1],v[i,2]]=b[x,y]) then 　　 begin 　　 dec(lt[i]); ht[i,z]:=false; 　　 if lt[i]=0 then exit; 　　 end; a[x,y]:=z; if t=n then begin put; exit; end; zh[t].h:=ht; zh[t].l:=lt; zhao; end; procedure work; begin while t0 do begin inc(zh[t].len); if zh[t].len9 then 　　 beg