codevs 1174 靶形数独

去题面的传送门
QAQ这道题已经调了很久了，一直是TLE
第一次交，codevs 35分，洛谷 70分（好吧看到这样子我还以为洛谷的数据水来着）
后来加了一个优化，codevs 95分，洛谷 80分（这才知道其实是codevs数据水）
优化是这样的：
一开始读入之后，把所有没有填的空位，它的所在行、列、九宫格里的已经填了的数字统计一下，放到结构体里，然后按照数字数量由大到小排序（所在区域已经填了的数字越多，越优先搜索），dfs时按照结构体顺序搜索，改掉原来一行一行地搜索
但是这样肯定不行，因为我还是T了（╮(╯▽╰)╭）
想想这样做有什么缺点：
①只统计了数字个数，而不是已经填了的数字的种类
②随着搜索的进行，下一个搜索的最优目标可能会改变，这时，初始时设定的已经排列好的顺序就不是最优的了，还是会有很多冗杂的操作

所以，最后的优化就是，在搜索过程中，用函数计算下一个要搜索的点（其实是启发式搜索的启发函数）。计算的标准就是，改点所在的行、列、九宫格内，还没有没使用的数字的数量，越少越好

（虽然一开始觉得这样很暴力，但是可以让搜索树变得超级小）

一开始T掉的点跑了大约6秒（手动测），优化之后变成717ms，时间限制是4s

ps：启发式搜索还是要选好启发函数

代码：

#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>using namespace std;int ans,cnt,sum;int ju[10][10],X[10][10],Y[10][10],S[10][10];bool vis[10][10],used[10][10];bool flg=0;struct hh{    int x,y;};int muti[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}};hh getit(){    hh ret;    int cm=15,sum=0;    for(int i=1;i<=9;++i)      for(int j=1;j<=9;++j)        if(!ju[i][j])        {            sum=0;            int wh=(i-1)/3*3+(j-1)/3+1;            for(int k=1;k<=9;++k)              if(!X[i][k]&&!Y[j][k]&&!S[wh][k])                sum++;            if(sum<cm)            {                cm=sum;                ret.x=i;                ret.y=j;            }        }    return ret;}void dfs(int x,int y,int k,int tot){    if(x>=1&&y<=9&&x<=9&&y>=1&&!vis[x][y]&&!used[x][y])    {        int wh=(x-1)/3*3+(y-1)/3+1;        for(int i=1;i<=9;++i)          if(!X[x][i]&&!Y[y][i]&&!S[wh][i])          {            X[x][i]=Y[y][i]=S[wh][i]=vis[x][y]=1;            ju[x][y]=i;             if(k==cnt) ans=max(ans,tot+i*muti[x][y]);            else            {                hh xx=getit();                dfs(xx.x,xx.y,k+1,tot+ju[x][y]*muti[x][y]);            }            ju[x][y]=0;            X[x][i]=Y[y][i]=S[wh][i]=vis[x][y]=0;          }    }}int main(){    for(int i=1;i<=9;++i)      for(int j=1;j<=9;++j)      {         scanf("%d",&ju[i][j]);         if(!ju[i][j])          {            cnt++;            continue;         }         sum+=ju[i][j]*muti[i][j];         used[i][j]=vis[i][j]=true;         int wh=(i-1)/3*3+(j-1)/3+1;         X[i][ju[i][j]]=S[wh][ju[i][j]]=Y[j][ju[i][j]]=1;      }    hh xx=getit();    dfs(xx.x,xx.y,1,sum);     if(ans==0) ans--;     printf("%d",ans);    return 0;}