数读9*9 个人训练赛/ G/F ——POJ - 3074 (DLK)

35 篇文章 0 订阅
17 篇文章 0 订阅


本来以为跟之前的做的一个4*4一样的一个题目,用普通的dfs去解决,但是开始dfs会输出各个答案,在加入flag之后使他停止在输出一组上,但是还是超时

超时代码

#include <iostream>
#include<cstring>
using namespace std;
int vis[15][15];
char tu[15][15];
int ju(int r,int l)
{
    for(int i=0;i<9;i++)
    {
        if(tu[i][l]==tu[r][l]&&i!=r)
            return 0;
    }
    for(int i=0;i<9;i++)
    {
        if(tu[r][i]==tu[r][l]&&i!=l)
            return 0;
    }
    int mr=r;
    int ml=l;
    while(mr%3!=0)
        mr--;
    while(ml%3!=0)
        ml--;
    for(int i=mr;i<=mr+2;i++)
    {
        for(int j=ml;j<ml+2;j++)
        {
            if(tu[i][j]==tu[r][l]&&i!=r&&j!=l)
                return 0;
        }
    }
    return 1;
}
int flag=0;
void dfs(int x)
{
    if(x==9*9)
    {    if(flag)
        return ;
        for(int i=0;i<9;i++)
        {
            for(int j=0;j<9;j++)
                cout<<tu[i][j];

        }
        flag=1;
        cout<<endl;
        return ;
    }
    else
    {
        int r=x/9;
        int l=x%9;
        if(tu[r][l]=='.')
        {
            if(flag)
              return;
            for(int k=1;k<=9;k++)
            {   if(flag)
                   return;
                tu[r][l]='0'+k;
                if(ju(r,l))
                {
                    dfs(x+1);
                }
            }
            tu[r][l]='.';
        }
        else
        {
           dfs(x+1);
        }
    }
    if(flag)
      return;
}
int main()
{  int flagg=0;
   while(1)
    {    flagg=0;
      memset(tu,0,sizeof(tu));
        for(int i=0;i<9;i++)
      {
          for(int j=0;j<9;j++)
        {

             cin>>tu[i][j];
             if(tu[i][j]=='d')
                {
                flagg=1;break;
                }
        }
        if(flagg)
            break;
      }
       if(flagg)
            break;
        flag=0;
      dfs(0);
    }
   // cout << "Hello world!" << endl;
    return 0;
}
搜了一下题解之后发现是一种dance links 的一种算法

点击打开链接

这个超好的 博客将DLX的!点击打开链接

这里引用一个题解


#include<iostream>  
    #include<cstring>  
    #include<string>  
    #include<cstdio>  
    #include<algorithm>  
    #include<vector>  
    #include<algorithm>  
      
    using namespace std;  
    //   列:(行+列+块)*9种可能+9*9个格子  
    //   行: 9*9*9  表示第i行第j列填k  
    const int MAXN=(9+9+9)*9+9*9+9*9*9*9*9*4+10;  
    #define INF 0xFFFFFF  
    int size;  
    int head,sz;  
    int U[MAXN],D[MAXN],L[MAXN],R[MAXN];  
    int H[MAXN],ROW[MAXN],C[MAXN],S[MAXN],O[MAXN];  
      
    void remove(int c)  
    {  
        L[R[c]]=L[c];  
        R[L[c]]=R[c];  
        for(int i=D[c];i!=c;i=D[i])  
        {  
            for(int j=R[i];j!=i;j=R[j])  
            {  
                U[D[j]]=U[j];  
                D[U[j]]=D[j];  
                --S[C[j]];  
            }  
         }  
    }  
      
    void resume(int c)  
    {  
        for(int i=U[c];i!=c;i=U[i])  
        {  
            for(int j=L[i];j!=i;j=L[j])  
            {  
                ++S[C[j]];  
                U[D[j]]=j;  
                D[U[j]]=j;  
              }  
         }  
         L[R[c]]=c;  
         R[L[c]]=c;  
    }  
      
    bool dfs(int k)  
    {  
        if(R[head]==head)  
        {  
            sort(O,O+9*9);  
            int p=0;  
            for(int i=0;i<9;i++)  
            {  
                for(int j=0;j<9;j++)  
                {  
                    int num=O[p++];  
                    //cout<<num<<endl;  
                    num=num-(i*9+j)*9;  
                    printf("%d",num);  
                }  
            }  
            printf("\n");  
            return  true;  
        }  
        int s=INF,c;  
        for (int t=R[head];t!=head;t=R[t])  
        {  
            if (S[t]<s)  
            {  
                s=S[t];  
                c=t;  
            }  
         }  
         remove(c);  
         for(int i=D[c];i!=c;i=D[i])  
         {  
              O[k]=ROW[i];  
              for(int j=R[i];j!=i;j=R[j])  
                  remove(C[j]);  
              if(dfs(k+1))  
                   return  true;  
              for(int j=L[i];j!=i;j=L[j])  
                   resume(C[j]);  
         }  
         resume(c);  
         return  false;  
    }  
      
    void initDL(int n)  
    {  
        head=0;  
        for(int i=0;i<=n;i++)  
        {  
            U[i]=i;D[i]=i;  
            L[i]=i-1;R[i]=i+1;  
            S[i]=0;  
        }  
        R[n]=0;L[0]=n;S[0]=INF+1;  
        sz=n+1;  
        memset(H,0,sizeof(H));  
    }  
      
    void insert(int i, int j)  
    {  
        if(H[i])  
        {  
            L[sz]=L[H[i]];  
            R[sz]=H[i];  
            L[R[sz]]=sz;  
            R[L[sz]]=sz;  
        }  
        else  
        {  
            L[sz]=sz;  
            R[sz]=sz;  
            H[i]=sz;  
        }  
        U[sz]=U[j];  
        D[sz]=j;  
        U[D[sz]]=sz;  
        D[U[sz]]=sz;  
        C[sz]=j;  
        ROW[sz]=i;  
        ++S[j];  
        ++sz;  
    }  
      
    char str[200];  
      
    void build()  
    {  
        int p=0;  
        initDL(9*9*4);  
        for(int i=0;i<9;i++)  
            for(int j=1;j<=9;j++,p++)  
            {  
                int base=(i*9+j-1)*9;  
                if(str[p]=='.')  
                {  
                    for(int k=1;k<=9;k++)  
                    {  
                        int r;  
                        r=base+k;  
                        //第i行有数字k  
                        insert(r,i*9+k);  
                        //第j列有数字k  
                        insert(r,9*9+(j-1)*9+k);  
                        //第k块有数字k  
                        int block=(j-1)/3*3+i/3;  
                        insert(r,9*9*2+block*9+k);  
                        //第i行j列有一个数字(限制一个格子只填一个数)  
                        insert(r,9*9*3+i*9+j);  
                    }  
                }  
                else  
                {  
                    int k=str[p]-'0';  
                    int r=base+k;  
                    //第i行有数字k  
                    insert(r,i*9+k);  
                    //第j列有数字k  
                    insert(r,9*9+(j-1)*9+k);  
                    //第k块有数字k  
                    int block=(j-1)/3*3+i/3;  
                    insert(r,9*9*2+block*9+k);  
                    //第i行j列有一个数字(限制一个格子只填一个数)  
                    insert(r,9*9*3+i*9+j);  
                }  
            }  
    }  
      
    int main()  
    {  
        size=9; //9*9数独  
        while(~scanf("%s",str))  
        {  
            if(strcmp(str,"end")==0)  
                break;  
            build();  
            dfs(0);  
        }  
        return 0;  
    }  



  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
这是一道比较经典的计数问目描述如下: 给定一个 $n \times n$ 的网格图,其中一些格子被标记为障碍。一个连通块是指一些被标记为障碍的格子的集合,满足这些格子在网格图中连通。一个格子是连通的当且仅当它与另一个被标记为障碍的格子在网格图中有公共边。 现在,你需要计算在这个网格图中,有多少个不同的连通块,满足这个连通块的大小(即包含的格子数)恰好为 $k$。 这是一道比较经典的计数问,一般可以通过计算生成函数的方法来解决。具体来说,我们可以定义一个生成函数 $F(x)$,其中 $[x^k]F(x)$ 表示大小为 $k$ 的连通块的个数。那么,我们可以考虑如何计算这个生成函数。 对于一个大小为 $k$ 的连通块,我们可以考虑它的形状。具体来说,我们可以考虑以该连通块的最左边、最上边的格子为起点,从上到下、从左到右遍历该连通块,把每个格子在该连通块中的相对位置记录下来。由于该连通块的大小为 $k$,因此这些相对位置一定是 $(x,y) \in [0,n-1]^2$ 中的 $k$ 个不同点。 现在,我们需要考虑如何计算这些点对应的连通块是否合法。具体来说,我们可以考虑从左到右、从上到下依次处理这些点,对于每个点 $(x,y)$,我们需要考虑它是否能够与左边的点和上边的点连通。具体来说,如果 $(x-1,y)$ 和 $(x,y)$ 都在该连通块中且它们在网格图中有公共边,那么它们就是连通的;同样,如果 $(x,y-1)$ 和 $(x,y)$ 都在该连通块中且它们在网格图中有公共边,那么它们也是连通的。如果 $(x,y)$ 与左边和上边的点都不连通,那么说明这个点不属于该连通块。 考虑到每个点最多只有两个方向需要检查,因此时间复杂度为 $O(n^2 k)$。不过,我们可以使用类似于矩阵乘法的思想,将这个过程优化到 $O(k^3)$ 的时间复杂度。 具体来说,我们可以设 $f_{i,j,k}$ 表示状态 $(i,j)$ 所代表的点在连通块中,且连通块的大小为 $k$ 的方案数。显然,对于一个合法的 $(i,j,k)$,我们可以考虑 $(i-1,j,k-1)$ 和 $(i,j-1,k-1)$ 这两个状态,然后把点 $(i,j)$ 加入到它们所代表的连通块中。因此,我们可以设计一个 $O(k^3)$ 的 DP 状态转移,计算 $f_{i,j,k}$。 具体来说,我们可以考虑枚举连通块所包含的最右边和最下边的格子的坐标 $(x,y)$,然后计算 $f_{x,y,k}$。对于一个合法的 $(x,y,k)$,我们可以考虑将 $(x,y)$ 所代表的点加入到 $(x-1,y,k-1)$ 和 $(x,y-1,k-1)$ 所代表的连通块中。不过,这里需要注意一个细节:如果 $(x-1,y)$ 和 $(x,y)$ 在网格图中没有相邻边,那么它们不能算作连通的。因此,我们需要特判这个情况。 最终,$f_{n,n,k}$ 就是大小为 $k$ 的连通块的个数,时间复杂度为 $O(n^2 k + k^3)$。 参考代码:

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值