NOIP 2009靶形数独解题报告.docVIP

靶形数独解题报告 题目来源:NOIP2009 TG 4 题目概述最大。 P.S 这份题解只有50%原创，代码都不是我的…… 主要考察知识点 解题思路’的做法。30分做法的话就是有一个弊端，每次检查每行每列以及每个九宫格有没填过k都要进行一次扫描（一般27次），这样本来状态量就大，在这里耗费了一个状态大量时间就有劣势。我们可以用类似状态压缩的思路来解决更新状态的难处。 设有数组hang[1..9] of longint记录9个longint，表示1~9行的填数情况。 假设我现在要在第一行填一个数字8(不用具体到第几行第几列，我们关心的是行而不是一格)。 那么就 hang[1]:=hang[1]+(1 shl (8-1)); 相当于在这个数的二进制形式第八位+1。 40分做法实际上是回溯法做的，那么消去这位的方法也和上面类似 hang[1]:=hang[1]-(1 shl (8-1)); 实际操作中为了减少代码量用个子过程。 procedure change(i,j,p:longint); //在行i摆数字j begin hang[i]:=hang[i]+p*(1 shl (j-1)); end; 举例：行5放8: change(5,8,1) 行5消去8: change(5,8,-1) 这样就实现了O(1)的判断和状态生成。拿到40分。 【45分做法】第二次交，在上面的基础上加个卡时（玩RP啊），我设的是timelimit=2700000 【100分做法，参程做法】至今没懂也没用过。使用M67神的八皇后位运算解题思想。 这回位运算就不限于我那么2的用来检测了…… 这个位运算方法还能一次性找出所有能填的数（45分做法里面还有一个for i:= 1 - 9来枚举检查）。注意一下C/C++的位运算 分别为pas的shr, shl。 : and |: or ^: xor /blog/archives/263 【100分做法，草根做法】第五次交……（直接被PIA），由于数独越里面权值越大，我们可以不按部就班的一个一个搜，而是按照从外圈向内圈（或从内圈向外圈）的顺序进行搜索。这样能保证最早的方案权值尽可能的大。如果有两个格子的选择范围一样大, 那么优先搜索尽量靠里的格子/blog/static/164280893201057102350461/ P.S：其实它很短。见参程2； 好了就这些了…… 至于lqs牛的90分解法也可以问下。 参考程序#includestdio.h#includemath.h int h[10]={},hs[10]={},zs[10]={},xj[5][5]={},hq[10]={};int ans=-1,st[10],a[10][10];void make(){　　int sum=0,i,j;　　fo r (i=1;i5;i++)　　 ?{fo r (j=i;j11-i;j++)　　　　sum+=(a[i][j]+a[10-i][j])*(5+i);　　 ? fo r (j=i+1;j10-i;j++)　　　　sum+=(a[j][i]+a[j][10-i])*(5+i);　　　　}　　sum+=a[5][5]*10;　　　　if (sumans) ans=sum;}void dfs(int k){　　 if (k==10) make();　　 else　　 {　　 ?int x,y,j,pos,p,i=st[k];　　 ?x=511-h[i];　　 ?y=x-x;　　 ?h[i]|=y;　　 ?j=(int)log2(y)+1;　　 ?pos=511-(hs[i]|zs[j]|xj[(i-1)/3][(j-1)/3]);　　 ?while (pos0)　　 ?{p=pos-pos;　　 ? pos-=p;　　 ? a[i][j]=(int)log2(p)+1;　　 ? hs[i]|=p;　　 ? zs[j]|=p;　　 ? xj[(i-1)/3][(j-1)/3]|=p;　　 ? if (x==y) dfs(k+1);　　 ? else dfs(k);　　 ? hs[i]-=p;　　 ? zs[j]-=p;　　 ? xj[(i-1)/3][(j-1)/3]-=p;　　 ? };　　 ? h[i]-=y;　　 ? };}int main(){　　int i,j,p0;　　fo r (i=1;i10;i++)　　 ?fo r (j=1;j10;j++)　　　　{scanf(%d,a[i][j]);　　　　 if (a[i][j]0)　　　　 ? {h[i]|=1(j-1