【NOIP模拟题】【DFS】【位运算】【舞蹈链】2016.11.12第三题题解

C

C
noip2009原题
【问题描述】
小城和小华都是热爱数学的好学生，最近，他们不约而同地迷上了数独游戏，好胜的他们想用数独来一比高低。但普通的数独对他们来说都过于简单了，于是他们向Z博士请教，Z博士拿出了他最近发明的“靶形数独”，作为这两个孩子比试的题目。
靶形数独的方格同普通数独一样，在9格宽×9格高的大九宫格中有9个3格宽×3格高的小九宫格（用粗黑色线隔开的）。在这个大九宫格中，有一些数字是已知的，根据这些数字，利用逻辑推理，在其他的空格上填入1到9的数字。每个数字在每个小九宫格内不能重复出现，每个数字在每行、每列也不能重复出现。但靶形数独有一点和普通数独不同，即每一个方格都有一个分值，而且如同一个靶子一样，离中心越近则分值越高。（如图）

上图具体的分值分布是：最里面一格（黄色区域）为10分，黄色区域外面的一圈（红色区域）每个格子为9分，再外面一圈（蓝色区域）每个格子为8分，蓝色区域外面一圈（棕色区域）每个格子为7分，最外面一圈（白色区域）每个格子为6分，如上图所示。比赛的要求是：每个人必须完成一个给定的数独（每个给定数独有可能有不同的填法），而且要争取更高的总分数。而这个总分数即每个方格上的分值和完成这个数独时填在相应格上的数字的乘积的总和。如图，在以下这个已经填完数字的靶形数独游戏中，总分为2829。游戏规定，将以总分数的高低决出胜负。

由于求胜心切，小城找到了善于编程的你，让你帮他求出，对于给定的靶形数独，能够得到的最高分数。
【输入】
输入文件名为c.in。
一共9行，每行9个整数（每个数都在0—9的范围内），表示一个尚未填满的数独方格，未填满的空格用“0”表示。每两个数字之间用一个空格隔开。
【输出】
输出文件c.out共1行。
输出可以得到的靶形数独的最高分数。如果这个数独无解，则输出整数-1。
【输入输出样例1】
c.in
c.out
7 0 0 9 0 0 0 0 1
1 0 0 0 0 5 9 0 0
0 0 0 2 0 0 0 8 0
0 0 5 0 2 0 0 0 3
0 0 0 0 0 0 6 4 8
4 1 3 0 0 0 0 0 0
0 0 7 0 0 2 0 9 0
2 0 1 0 6 0 8 0 4
0 8 0 5 0 4 0 1 2
2829
【输入输出样例2】
c.in
c.out
0 0 0 7 0 2 4 5 3
9 0 0 0 0 8 0 0 0
7 4 0 0 0 5 0 1 0
1 9 5 0 8 0 0 0 0
0 7 0 0 0 0 0 2 5
0 3 0 5 7 9 1 0 8
0 0 0 6 0 1 0 0 0
0 6 0 9 0 0 0 0 1
0 0 0 0 0 0 0 0 6
2852
【数据范围】
40%的数据，数独中非0数的个数不少于30。
80%的数据，数独中非0数的个数不少于26。
100%的数据，数独中非0数的个数不少于24。
明显的DFS，关键是怎么优化。
不妨用位运算，当然我们也会想到直接用舞蹈链搞。
附位运算代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<set>
#include<queue>
#include<algorithm>
#include<vector>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<stack>
#define rez(i,x,y) for(int i=x;i>=y;i--)
#define res(i,x,y) for(int i=x;i<=y;i++)
#define INF 2100000000
#define ll long long
#define clr(x)  memset(x,0,sizeof(x))
#define NAME "c"
#define con(i) (1<<(i))
#define id(i,j) ((i)/3*3+(j)/3)
#define order(i,j) (3*((i)-(i)/3*3)+(j)-(j)/3*3)
using namespace std;
int mp[11][11],row[9],lie[9],line[9],ma[9];
int f[512],ans,node[9],cnt[9],fin;
int sore[9][9]={{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}};
void grade(){
    int t=0;
    for(int i=0;i<9;i++){
        for(int j=0;j<9;j++){
            t+=mp[i][j]*sore[i][j];
        }
    }
    if(ans<t)ans=t;
}
void search(int ti){
    int i,j,p;
    int pos,k;
    if(ti==9){
        grade();return;
    }
    i=node[ti];
    if(cnt[i]==0){
        search(ti+1);return;
    }
    cnt[i]--;
    p=511^row[i];p=p&-p;
    row[i]|=p;j=f[p];
    pos=511^(lie[i]|line[j]|ma[id(i, j)]);
    while(pos>0){
        k=pos&-pos;pos^=k;
        lie[i]|=k;line[j]|=k;
        ma[id(i,j)]|=k;mp[i][j]=f[k]+1;
        search(ti);
        lie[i]^=k;line[j]^=k;
        ma[id(i,j)]^=k;
    }
    cnt[i]++;
    row[i]^=p;
}
int main(){ 
    freopen(NAME".in","r",stdin);
    freopen(NAME".out","w",stdout);
    int t=0;
    for(int i=1,j=0;i<=511;i<<=1,j++)f[i]=j;
    for(int i=0;i<9;i++){
        for(int j=0;j<9;j++){
            cin>>mp[i][j];
            if(mp[i][j]!=0){
                row[i]|=con(j);
                t=con(mp[i][j]-1);
                if(((lie[i]&t)!=0)||((line[j]&t)!= 0)||((ma[id(i,j)]&t)!=0)){
                    cout<<-1<<endl;return 0;
                }
                lie[i]|=t;line[j]|=t;
                ma[id(i,j)]|=t;
            }else cnt[i]++;
        }
    }
    for(int i=0;i<9;i++)node[i]=i;
    for(int i=0;i<9;i++){
        for(int j=i+1;j<9;j++){
            if(cnt[node[i]]>cnt[node[j]]){
                node[i]^=node[j];
                node[j]^=node[i];
                node[i]^=node[j];
            }
        }
    }
    int tot=0;
    while(cnt[node[tot]]==0)tot++;
    search(tot);
    if(ans==0){
        cout<<-1<<endl;return 0;
    }   
    cout<<ans<<endl;
    return 0;
}