搜索 素数方阵

Posted 2020-09-02 Lazer

一个比较麻烦的数独

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


 

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





 
由于求胜心切，小城找到了善于编程的你，让你帮他求出，对于给定的靶形数独，能够得到的最高分数。

暴力并不难，并且可以获得一个不错的90分（不累赘的话）；

但是100分的话貌似需要些位运算的知识

因为数据水的原因.......改变一下搜索顺序就A了....(那天去学学位运算....)

代码如下

#include<iostream>
#include<cstring>
#include<string>
#include<algorithm>
using namespace std;
int h[10][10];
int l[10][10];
int kaka[10][10]={
    0,0,0,0,0,0,0,0,0,0,
    0,1,1,1,4,4,4,7,7,7,
    0,1,1,1,4,4,4,7,7,7,
    0,1,1,1,4,4,4,7,7,7,
    0,2,2,2,5,5,5,8,8,8,
    0,2,2,2,5,5,5,8,8,8,
    0,2,2,2,5,5,5,8,8,8,
    0,3,3,3,6,6,6,9,9,9,
    0,3,3,3,6,6,6,9,9,9,
    0,3,3,3,6,6,6,9,9,9,
};
int grade[10][10]={
    0,0,0,0,0,0,0,0,0,0,
    0,6,6,6,6,6,6,6,6,6,
    0,6,7,7,7,7,7,7,7,6,
    0,6,7,8,8,8,8,8,7,6,
    0,6,7,8,9,9,9,8,7,6,
    0,6,7,8,9,10,9,8,7,6,
    0,6,7,8,9,9,9,8,7,6,
    0,6,7,8,8,8,8,8,7,6,
    0,6,7,7,7,7,7,7,7,6,
    0,6,6,6,6,6,6,6,6,6,
};
struct zz
{
    int x,y;
}q[1000000];
int xiao[10][10];
int topx=0;
int f[1000];
int top[1000];
int toop[10][10];
int num[10][10];
int ff[20][20];
int ans[1000][100];
int chick(int x,int y,int z)
{
    int k=kaka[x][y];
    if(ff[k][z]==1)return 0;
    for(int i=1;i<=9;i++)
    if(num[i][y]==z||num[x][i]==z)return 0;
    return 1;
}
void chuli()
{
    for(int i=1;i<=topx;i++)
    {
        int x=q[i].x,y=q[i].y;
        for(int u=1;u<=9;u++)
        {
            if(chick(x,y,u)==1)
            ans[i][++top[i]]=u;
        }
    }
}
bool my(zz a,zz b)
{
    return grade[a.x][a.y]<grade[b.x][b.y];
}
int maxx=-1;
int dfs(int i,int man)
{
    if(i==topx+1){maxx=max(maxx,man);return 0;}
    int x=q[i].x,y=q[i].y;
    int k=kaka[x][y];
    for(int z=1;z<=top[i];z++)
    {
        int u=ans[i][z];
        if(ff[k][u]==0&&h[x][u]==0&&l[y][u]==0)
        {
            h[x][u]=1;
            l[y][u]=1;
            num[x][y]=u;
            ff[k][u]=1;
            dfs(i+1,man+num[x][y]*grade[x][y]);
            h[x][u]=0;
            l[y][u]=0;
            num[x][y]=0;
            ff[k][u]=0;
        }
    }
    return 0;
}
int main()
{
    int man=0;
    for(int i=1;i<=9;i++)
        for(int u=1;u<=9;u++)
        {
            cin>>num[i][u];
            if(num[i][u]==0)
                q[++topx].x=i,q[topx].y=u;
            if(num[i][u]!=0)
            {
                int k=num[i][u];
                int w=kaka[i][u];
                ff[w][k]=1;
                h[i][k]=1;
                l[u][k]=1;
                
                
            }
            man+=num[i][u]*grade[i][u];
        }
    //sort(q+1,q+topx+1,my);
        chuli();
    dfs(1,man);
    cout<<maxx<<endl;
    return 0;
}