[NOIP2009]靶形数独

 1 #include <algorithm>
 2 #include <cstdio>
 3 #include <cctype>
 4 #include <cmath>
 5 #include <vector>
 6 #include <queue>
 7 
 8 #define maxn ((int)1e5 + 4)
 9 #if defined DEBUG
10 FILE *in = fopen("test", "r");
11 #define out stdout
12 #else
13 FILE *in = fopen("sudoku.in", "r");
14 FILE *out = fopen("sudoku.out", "w");
15 #endif
16 using namespace std;
17 inline void getint(int &x){
18     char c = fgetc(in);
19     while(!isdigit(c)) c = fgetc(in);
20     x = c - ‘0‘;
21     while(isdigit(c = fgetc(in)))x = x * 10 - ‘0‘ + c;
22 }
23 
24 /*===========================================*/
25 #define lowbit(x) (x & -x)
26 #define gs(x, y) (x / 3 * 3 + y / 3)
27 #define f(a) (bitcnt(RNG &~(row[a.x] | col[a.y] | squ[gs(a.x,a.y)])))
28 const int RNG = 1022/*状态的范围*/,val[9][9] = {
29     {6,6,6,6,6,6,6,6,6},
30     {6,7,7,7,7,7,7,7,6},
31     {6,7,8,8,8,8,8,7,6},
32     {6,7,8,9,9,9,8,7,6},
33     {6,7,8,9,10,9,8,7,6},
34     {6,7,8,9,9,9,8,7,6},
35     {6,7,8,8,8,8,8,7,6},
36     {6,7,7,7,7,7,7,7,6},
37     {6,6,6,6,6,6,6,6,6}};
38 struct P{
39     int x, y;
40     P(){}
41 }blank[80];
42 int bcnt = 0, ord[80];
43 int row[9] = {0}, col[9] = {0}, squ[9] = {0};//状态
44 int ans = 0;
45 bool known = 0;
46 inline int bitcnt(int x){
47     int a = 0;
48     while(x){
49         if(x & 1)++a;
50         x >>= 1;
51     }
52     return a;
53 }
54 inline void init(){
55     int i, j, k, t;
56     for(i = 0;i < 9;++i)for(j = 0;j < 9;++j){
57         getint(k);
58         if(k){
59             t = 1 << k;
60             row[i] ^= t, col[j] ^= t, squ[gs(i,j)] ^= t;
61             ans += k * val[i][j];
62         }
63         else{
64             ord[bcnt] = bcnt;
65             blank[bcnt].x = i;
66             blank[bcnt++].y = j;
67         }
68     }
69 }
70 inline void dfs(int cur, int now){
71     if(cur == bcnt){
72         if(now > ans)ans = now,known = 1;
73         return;
74     }
75     int mincnt = 10, min_it, i, t;
76     for(i = cur;i < bcnt;++i)
77         if((t = f(blank[ord[i]])) < mincnt)
78             mincnt = t, min_it = i;
79     P x = blank[ord[min_it]];
80     if(cur!=min_it)ord[cur] ^= ord[min_it] ^= ord[cur] ^= ord[min_it];
81     int bits = RNG &~(row[x.x] | col[x.y] | squ[gs(x.x,x.y)]);
82     while(bits){
83         t = lowbit(bits); bits ^= t;
84         i = bitcnt(t - 1);
85         row[x.x] ^= t, col[x.y] ^= t, squ[gs(x.x,x.y)] ^= t;
86         dfs(cur + 1, now + i * val[x.x][x.y]);
87         row[x.x] ^= t, col[x.y] ^= t, squ[gs(x.x,x.y)] ^= t;
88     }
89 }
90 int main(){
91     init();
92     dfs(0, ans);
93     if(known)
94         fprintf(out, "%d\n", ans);
95     else fprintf(out, "-1\n");
96     return 0;
97 }