RQNOJ521 sudoku

         这是noip2009 的第四题，也是一道让我心痛不已的题目……

         拿到题目发现这显然是一到搜索，这不难看出。可能有人用dancing links 来解这道题，我认为这确实是正解，但不得不承认dancing links是一种比较高级的数据结构了，并非那么轻易就掌握的。

         既然如此，单纯的但富有有技巧的搜索能不能解决这个问题呢？答案是肯定的！我今天的这个代码用到了位运算判重——这是可以过掉80%的测试数据的。如果注意搜索顺序，加以改进便在0.5s之内就解决了问题，远远小于2s的时间限制。

         当然，可以看出我的这份代码还有优化的空间。搜索顺序中我只考虑了行的因素，那同理，列与九宫格也可以一起考虑进来一起排序。另外，这个顺序是动态的，如果用堆来维护，效果会更好！而且，我还没有加上对于ans的最优化剪枝！可以说优化空间巨大！

         不过，抛开对题目的剖析，仅仅从应试的角度来看，下面的代码，不仅仅AC了这道题，也将编程复杂度降到了理想的水平。应该说在noip中，为这道压轴题敲出130行的满分代码不会占用其他题目过多时间，十分划算！

 

 

 

program sudoku; const all=511; map:array[1..9,1..9] of longint=((1,1,1,2,2,2,3,3,3), (1,1,1,2,2,2,3,3,3), (1,1,1,2,2,2,3,3,3), (4,4,4,5,5,5,6,6,6), (4,4,4,5,5,5,6,6,6), (4,4,4,5,5,5,6,6,6), (7,7,7,8,8,8,9,9,9), (7,7,7,8,8,8,9,9,9), (7,7,7,8,8,8,9,9,9)); scr:array[1..9,1..9] of longint=((6,6,6,6,6,6,6,6,6), (6,7,7,7,7,7,7,7,6), (6,7,8,8,8,8,8,7,6), (6,7,8,9,9,9,8,7,6), (6,7,8,9,10,9,8,7,6), (6,7,8,9,9,9,8,7,6), (6,7,8,8,8,8,8,7,6), (6,7,7,7,7,7,7,7,6), (6,6,6,6,6,6,6,6,6)); var data:array[1..9,1..9] of longint; exp,row,col,nine,line,sum:array[0..10] of longint; lob:array[0..1024] of longint; ans,tmpans:longint; procedure prework; var i:longint; begin ans:=-1; tmpans:=0; for i:=1 to 9 do begin exp[i]:=1 shl (i-1); lob[1 shl (i-1)]:=i; end; fillchar(row,sizeof(row),0); fillchar(nine,sizeof(nine),0); fillchar(col,sizeof(col),0); fillchar(line,sizeof(sum),0); end; procedure init; var i,j,k:longint; begin for i:=1 to 9 do for j:=1 to 9 do begin read(k); data[i,j]:=k; if k=0 then continue; inc(sum[i]); inc(row[i],exp[k]); inc(col[j],exp[k]); inc(nine[map[i,j]],exp[k]); inc(tmpans,k*scr[i,j]); end; for i:=1 to 9 do line[i]:=i; for i:=1 to 8 do for j:=i+1 to 9 do if sum[i]<sum[j] then begin k:=sum[i]; sum[i]:=sum[j]; sum[j]:=k; k:=line[i]; line[i]:=line[j]; line[j]:=k; end; end; procedure outit; begin writeln(ans); end; function lowbit(x:longint):longint; begin exit(x and -x); end; procedure dfs(ii,j,tmp:longint); var i,bit,k,p,t:longint; begin if ii=10 then begin if tmp>ans then ans:=tmp; exit; end; if j=10 then begin dfs(ii+1,1,tmp); exit; end; i:=line[ii]; if data[i,j]<>0 then begin dfs(ii,j+1,tmp); exit; end; t:=map[i,j]; bit:=all and not (row[i] or col[j] or nine[t]); while bit>0 do begin p:=lowbit(bit); dec(bit,p); k:=lob[p]; inc(row[i],p); inc(col[j],p); inc(nine[t],p); dfs(ii,j+1,tmp+k*scr[i,j]); dec(row[i],p); dec(col[j],p); dec(nine[t],p); end; end; begin assign(input,'sudoku.in'); reset(input); assign(output,'sudoku.out'); rewrite(output); prework; init; dfs(1,1,tmpans); outit; close(input); close(output); end.