回溯算法

http://www.cnblogs.com/hustcat/archive/2008/04/09/1144645.html

http://www.cnblogs.com/hustcat/archive/2008/04/13/1151379.html

引言
      寻找问题的解的一种可靠的方法是首先列出所有候选解,然后依次检查每一个,在检查完所有或部分候选解后,即可找到所需要的解。理论上,当候选解数量有限并且通过检查所有或部分候选解能够得到所需解时,上述方法是可行的。不过,在实际应用中,很少使用这种方法,因为候选解的数量通常都非常大(比如指数级,甚至是大数阶乘),即便采用最快的计算机也只能解决规模很小的问题。对候选解进行系统检查的方法有多种,其中回溯和分枝定界法是比较常用的两种方法。按照这两种方法对候选解进行系统检查通常会使问题的求解时间大大减少(无论对于最坏情形还是对于一般情形)。事实上,这些方法可以使我们避免对很大的候选解集合进行检查,同时能够保证算法运行结束时可以找到所需要的解。因此,这些方法通常能够用来求解规模很大的问题。

算法思想
      回溯(backtracking)是一种系统地搜索问题解答的方法。为了实现回溯,首先需要为问题定义一个解空间(solution space),这个空间必须至少包含问题的一个解(可能是最优的)。

     下一步是组织解空间以便它能被容易地搜索。典型的组织方法是图(迷宫问题)或树(N皇后问题)。
      一旦定义了解空间的组织方法,这个空间即可按深度优先的方法从开始节点进行搜索。

回溯方法的步骤如下:
     1) 定义一个解空间,它包含问题的解。
     2) 用适于搜索的方式组织该空间。
     3) 用深度优先法搜索该空间,利用限界函数避免移动到不可能产生解的子空间。
回溯算法的一个有趣的特性是在搜索执行的同时产生解空间。在搜索期间的任何时刻,仅保留从开始节点到当前节点的路径。因此,回溯算法的空间需求为O(从开始节点起最长路径的长度)。这个特性非常重要,因为解空间的大小通常是最长路径长度的指数或阶乘。所以如果要存储全部解空间的话,再多的空间也不够用。

算法应用
回溯算法的求解过程实质上是一个先序遍历一棵"状态树"的过程,只是这棵树不是遍历前预先建立的,而是隐含在遍历过程中<<数据结构>>(严蔚敏).
(1) 幂集问题(组合问题) (参见《数据结构》(严蔚敏))
     求含N个元素的集合的幂集。
     如对于集合A={1,2,3},则A的幂集为
     p(A)={{1,2,3},{1,2},{1,3},{1},{2,3},{2},{3},Φ}
幂集的每个元素是一个集合,它或是空集,或含集合A中的一个元素,或含A中的两个元素,或者等于集合A。反之,集合A中的每一个元素,它只有两种状态:属于幂集的元素集,或不属于幂集元素集。则求幂集P(A)的元素的过程可看成是依次对集合A中元素进行“取”或“舍”的过程,并且可以用一棵状态树来表示。求幂集元素的过程即为先序遍历这棵状态树的过程。

下面给出程序:

#include  < stdio.h >
#include 
< malloc.h >
#define  ERROR 0
#define  OK 1

typedef 
int  ElemType;
typedef 
struct  LNode
{
    ElemType data;
    
struct LNode *next;
}
 LNode, * LinkList;

// 初始化
LinkList ListInit()
{

    LNode 
*base=(LinkList)malloc(sizeof(LNode));
    
base->data=0;
    
base->next=NULL;
    
return base;
}

// 插入一个元素
int  ListInsert(LinkList L, int  i,ElemType e)
{
    LNode 
*p,*s;
    
int j=0;
    p
=(LNode *)L;
    
while(p&&j<i-1)
    
{
        p
=p->next;
        
++j;
    }

    
if(!p||j>i-1)
        
return ERROR;
    s
=(LNode *)malloc(sizeof(LNode));
    s
->data=e;
    s
->next=p->next;
    p
->next=s;
    
return OK;
}

// 删除一个结点
int  ListDelete(LinkList  & L, int  i,ElemType  & e)
{
    LinkList p
=L,q;
    
int j=0;
    
while(p->next&&j<i-1)
    
{
        p
=p->next;
        
++j;
    }

    
if(!(p->next)||j>i-1)
        
return ERROR;
    q
=p->next;
    p
->next=q->next;
    e
=q->data;
    free(q);
}

// 长度
int  ListLength(LinkList L)
{
    LinkList p
=L;
    
int j=0;
    
if(!L)
        
return ERROR;
    
while(p->next)
    
{
        p
=p->next;
        
++j;
    }

    
return j;
}

// 查找一个元素
int  GetElem(LinkList L, int  i,ElemType  & e)
{
    LNode 
*p=L;
    
int j=0;
    
while(p->next&&j<i)
    
{
        p
=p->next;
        
++j;
    }

    
if(!p||j>i)
        
return ERROR;
    e
=p->data;
    
return OK;
}

// 输出链表元素
void  Display(LinkList L)
{
    LNode 
*p=L;
    
if(!(p->next))
    
{
        printf(
"NULL,");
        
return;
    }

    
else
        p
=p->next;
    
while(p)
    
{

        printf(
"%d,",p->data);
        p
=p->next;    
    }

}

// 求幂集
void  PowerSet( int  i,LinkList A,LinkList  & B)
{
    
int k=0;
    ElemType e
=0;
    
if(i>ListLength(A))
    
{
        Display(B);
        printf(
"\n");
    }

    
else
    
{
        GetElem(A,i,e);
        k
=ListLength(B);
        ListInsert(B,k
+1,e);
        PowerSet(i
+1,A,B);

        ListDelete(B,k
+1,e);
        PowerSet(i
+1,A,B);
    }

}


int  main()
{
    LinkList list
=ListInit(); //初始化
    LinkList list2=ListInit();//初始化

    ListInsert(list,
1,1);//插入元素
    ListInsert(list,2,2);
    ListInsert(list,
3,3);

    Display(list);
//输出元素
    printf("\npower set is:\n");
    PowerSet(
1,list,list2);//求幂集
}

 

(2)迷宫问题(参见《数据结构》(严蔚敏))
计算机解迷宫时,通常用的是"试探和回溯"的方法,即从入口出发,顺某一方向向前探索,若能走通,则继续往前走;否则沿原路退回,换一个方向再继续探索,直至所有可能的通路都探索到为止,如果所有可能的通路都试探过,还是不能走到终点,那就说明该迷宫不存在从起点到终点的通道。

  1.从入口进入迷宫之后,不管在迷宫的哪一个位置上,都是先往东走,如果走得通就继续往东走,如果在某个位置上往东走不通的话,就依次试探往南、往西和往北方向,从一个走得通的方向继续往前直到出口为止;

  2.如果在某个位置上四个方向都走不通的话,就退回到前一个位置,换一个方向再试,如果这个位置已经没有方向可试了就再退一步,如果所有已经走过的位置的四个方向都试探过了,一直退到起始点都没有走通,那就说明这个迷宫根本不通;
   
   3.所谓"走不通"不单是指遇到"墙挡路",还有"已经走过的路不能重复走第二次",它包括"曾经走过而没有走通的路"。显然为了保证在任何位置上都能沿原路退回,需要用一个"后进先出"的结构即栈来保存从入口到当前位置的路径。并且在走出出口之后,栈中保存的正是一条从入口到出口的路径。

由此,求迷宫中一条路径的算法的基本思想是:
若当前位置"可通",则纳入"当前路径",并继续朝"下一位置"探索;若当前位置"不可通",则应顺着"来的方向"退回到"前一通道块",然后朝着除"来向"之外的其他方向继续探索;若该通道块的四周四个方块均"不可通",则应从"当前路径"上删除该通道块。

设定当前位置的初值为入口位置; 
  do{
    若当前位置可通, 
    则{
     将当前位置插入栈顶;       // 纳入路径 
     若该位置是出口位置,则算法结束; 
      // 此时栈中存放的是一条从入口位置到出口位置的路径
     否则切换当前位置的东邻方块为新的当前位置; 
     }
    否则
    {
    若栈不空且栈顶位置尚有其他方向未被探索, 
    则设定新的当前位置为: 沿顺时针方向旋转找到的栈顶位置的下一相邻块;
    若栈不空但栈顶位置的四周均不可通, 
    则{ 删去栈顶位置;         // 从路径中删去该通道块
      若栈不空,则重新测试新的栈顶位置, 
      直至找到一个可通的相邻块或出栈至栈空; 
     } 
   } 
} while (栈不空);

程序如下:

 

#include  < stdio.h >

#define  WALL   0   //
#define  CORRIDOR 1  // 通道
#define  PATH  9  // 为路径上的一块
#define  TRIED 2  //

#define  ROW_NUM    7  // 迷宫数组行数
#define  COL_NUM   13  // 列数

#define  TRUE 1
#define  FALSE 0
#define  MAXSIZE 50
typedef 
struct  
{
    
int row;
    
int col;
}
PosType;

typedef 
struct  
{
    
int ord;      //通道块在路径上的"序号"
    PosType seat; //通道块在迷宫中的坐标
    int di;       //当前通道块的方向
}
SElemType;
typedef 
struct  
{
    SElemType S[MAXSIZE];
    
int top;
}
MazeType;
// 迷宫
int  grid[ROW_NUM][COL_NUM] = {{1110111001111},
                            
{1011101111101},
                            
{1000001010101},
                            
{1000111010111},
                            
{1111100001000},
                            
{0000100000000},
                            
{0000111111111}}
;
// 当前位置是否可以通过
bool  Valid(PosType pos)
{
    
if(pos.row>=0&&pos.row<=ROW_NUM&&pos.col>=0&&pos.col<=COL_NUM&&grid[pos.row][pos.col]==CORRIDOR)
        
return TRUE;
    
else
        
return FALSE;
}

void  FootPrint(PosType pos) // 留下足迹
{
    grid[pos.row][pos.col]
=PATH;
}

void  Undo(PosType pos)  // 留下不能通过的标识
{
    grid[pos.row][pos.col]
=TRIED;
}

// 当前位置的下一个位置
PosType NextPos(PosType cur, int  di)
{
    PosType next;
    
switch(di)
    
{
    
case 0//
        next.row=cur.row;
        next.col
=cur.col+1;
        
break;
    
case 1//
        next.row=cur.row+1;
        next.col
=cur.col;
        
break;
    
case 2:  //西
        next.row=cur.row;
        next.col
=cur.col-1;
        
break;
    
case 3:  //
        next.row=cur.row-1;
        next.col
=cur.col;
        
break;
    }

    
return next;
}

// 是否到达终点
bool  Done(PosType cur,PosType end)
{
    
if(cur.row==end.row&&cur.col==end.col)
        
return TRUE;
    
else
        
return FALSE;
}

// 寻找迷宫路径
bool  MazePath(MazeType  & path,PosType start,PosType end)
{
    SElemType e;
    path.top
=-1;
    
int step=1;
    PosType curpos
=start;
    
do
    
{
        
if(Valid(curpos))
        
{
            FootPrint(curpos);
            e.ord
=step;
            e.di
=0;
            e.seat
=curpos;
            path.S[
++path.top]=e;
            
if(Done(curpos,end))
                
return TRUE;
            curpos
=NextPos(curpos,0);
            step
++;
        }

        
else
        
{
            
if(path.top>-1)//棧不空
            {
                e
=path.S[path.top--];
                
while(e.di==3&&path.top>-1)
                
{
                    Undo(e.seat);
                    e
=path.S[path.top--];
                }

                
if(e.di<3)
                
{
                    e.di
++;
                    path.S[
++path.top]=e;
                    curpos
=NextPos(e.seat,e.di);
                }

            }
//if
        }
//else
    }
while(path.top>-1);
    
return FALSE;
}

// 输出路径
void  PrintPath(MazeType path)
{
    
int i=0;
    
while(i<=path.top)
    
{
        printf(
"第%d步:(%d,%d)\n",path.S[i].ord,path.S[i].seat.row,path.S[i].seat.col);
        i
++;
    }

}

// 输出路径
void  PrintPath2()
{
    
for(int i=0;i<ROW_NUM;i++)
        
for(int j=0;j<COL_NUM;j++)
        
if(grid[i][j]==PATH)
            printf(
"(%d,%d)\n",i,j);
}

int  main()
{
    MazeType path;
    PosType start
={0,0},end={6,12};
    
if(MazePath(path,start,end))
        PrintPath(path);
    
else
        printf(
"not reachable!\n");

    PrintPath2();
}

(3)N皇后问题:
在一个N*N的棋盘上放置N个皇后,且使得每两个之间不能互相攻击,也就是使得每两个不在同一行,同一列和同一斜角线上。
对于N=1,问题的解很简单,而且我们很容易看出对于N=2和N=3来说,这个问题是无解的。所让我们考虑4皇后问题并用回溯法对它求解。因为每个皇后都必须分别占据—行,我们需要做的不过是为图1棋盘上的每个皇后分配一列。

 
     我们从空棋盘开始,然后把皇后1放到它所在行的第一个可能位置上,也就是第一行第—列。对于皇后2,在经过第一列和第二列的失败尝试之后,我们把它放在第一个可能的位置,就是格子〔23),位于第二行第二列的格子。这被证明是一个死胡同,因为皇后:将没有位置可放。所以,该算法进行回溯,把皇后2放在下一个可能位置(24)上。然后皇后3就可以放在(32),这被证明是另一个死胡同。该算法然后就回溯到底,把皇后1移到(12)。接着皇后2(24),皇后3(31),而皇后4(43),这就是该问题的一个解。图2给出了这个查找的状态空间树。

 程序如下:

 

#include  < stdio.h >
#include 
< math.h >
#define  N  4
int  col[N + 1 ];
// 输出结果
void  Output()
{
    
for(int i=1;i<=N;i++)
    
{
        printf(
"(%d,%d)\n",i,col[i]);
    }

    printf(
"\n");
}

// 求解函数
void  Queen( int  i, int  n)
{
    
if(i>n)
        Output();
    
else
    
{
        
for(int j=1;j<=n;++j)
        
{
            
int k=1;
            col[i]
=j;
            
while(k<i)
            
{
                
if((col[k]-col[i])*(fabs(col[k]-col[i])-fabs(k-i))!=0)
                
{
                    k
++;
                    
if(k==i)
                        Queen(i
+1,n);
                }

                
else
                
{
                    
break;
                }

            }

        }

    }

}

int  main()
{
    printf(
"the answer is:\n");
    
for(int i=1;i<=N;i++)
    
{
        col[
1]=i; //设置第一行
        Queen(2,N);
    }

}

结果如下:

 

引言
      在<<算法系列---回溯算法>>一节,讨论回溯算法及其应用,回溯能够在可以接受的时间内解决某些规模的组合问题,这节再讨论它的一个非常有意思的应用---跳马问题(骑士周游问题)。

问题
跳马问题也称为骑士周游问题,是算法设计中的经典问题。其一般的问题描述是:

     考虑国际象棋棋盘上某个位置的一只马,它是否可能只走63步,正好走过除起点外的其他63个位置各一次?如果有一种这样的走法,则称所走的这条路线为一条马的周游路线。试设计一个算法找出这样一条马的周游路线。

此题实际上是一个Hamilton回路问题,和迷宫问题很相似,可以用回溯算法来解决.
考虑到马在每一个位置,最多有8种跳法,如下图所示:
 

               
   

K7

 

K0

     
 

K6

     

K1

   
     

K

       
 

K5

     

K2

   
   

K4

 

K3

     
               
               

可以使用N皇后问题的算法模板。
算法如下:

#include  < stdio.h >
#include 
< stdlib.h >
/**/
// 棋盘行数
const   int  N  =   8 ;

int  step[N  *  N]  =   {-1} // 保存每一步做出的选择
int  chess[N][N]  =   {0} // 棋盘
// 下一个方向
int  Jump[ 8 ][ 2 =   {{-21}{-12}{12}{21}{2-1}{1-2}{-1-2}{-2-1}} ;

int  p  =   0 ; // 对解的个数计数

// 判断是否合法
int  canJump( int  x,  int  y)
{
    
if (x >= 0 && x < N && y >= 0 && y < N && chess[x][y] == 0)
        
return 1;
    
return 0;
}

// 输出结果
void  OutPutChess()
{
    
int i;
    
for (i = 1; i <= N * N - 1++i)
        
{
            printf(
"%d ", step[i]);
        }

    printf(
"\n");
    
for(i=0;i<N;i++)
    
{
        
for(int j=0;j<N;j++)
            printf(
"%3d ",chess[i][j]);
        printf(
"\n");
    }

}


// 回溯算法
void  BackTrace( int  t,  int  x,  int  y)
{
    
if (t >= N * N) 
    
//if(t>=37)
    {
        p
++;
        OutPutChess();
//输出结果
        exit(1);   //求得一个解时,退出程序
        
//return;  //求得所有解
    }

    
else
    
{
        
for (int i = 0; i < 8++i)
        
{
            
if (canJump(x + Jump[i][0], y + Jump[i][1]))
            
{   //求下一个结点
                x += Jump[i][0];
                y 
+= Jump[i][1];

                printf(
"(%2d,%2d)",x,y);//打印当结点
                chess[x][y] = t + 1;
                step[t] 
= i;
                BackTrace(t 
+ 1, x, y);//递归调用
                
//回溯
                chess[x][y] = 0;
                x 
-= Jump[i][0];
                y 
-= Jump[i][1];
            }

        }

        
    }

}

/**/
int  main()
{
    
int x = 0;
    
int y = 0;
    chess[x][y] 
= 1;
    BackTrace(
1, x, y);
    printf(
"All Results Number = %d ", p);
}

该算法最坏的时间复杂度为O(8 N*N )。这是一个相当大的数字,不能接受,但实际情况好得多。并且在37步的时候,开
发生回溯,可通过改BackTrace中的第一个if中的参数发现 据说,总的回溯次数达300多百万次.
我的机子运行20分钟也没运行出结果.我们可以考虑,当N=8时,为2192,即使对于2100=1.3*1030,对于一台每秒1万亿(1012)次操作的计算机,也需要4*1010才能完成,超过了45亿年---地球的估计年龄.
但是,该算法可以适当改进,考虑到:
即向前看两步,当每准备跳一步时,设准备跳到(x, y)点,计算(x, y)这一点可能往几个方向跳(即向前看两步),将这个数目设为(x, y)点的权值,将所  有可能的(x, y)按权值排序,从最小的开始,循环遍历所有可能的(x, y),回溯求出结果。算法可以求出所有可能的马跳棋盘路径,算出一个可行 的结果很快,但在要算出所有可能结果时,仍然很慢,因为最坏时间复杂度本质上并没有改变,仍为O(8^(N * N)),但实际情况很好,在瞬间即可得到一个解,当然,要求得所有解,也需要很长的时间.
下面是实现这一思想的代码:

#include  < stdio.h >
#include 
< stdlib.h >
/**/
// 棋盘行数
const   int  N  =   8 ;

int  step[N  *  N]  =   {-1} // 保存每一步做出的选择
int  chess[N][N]  =   {0} // 棋盘

int  Jump[ 8 ][ 2 =   {{-21}{-12}{12}{21}{2-1}{1-2}{-1-2}{-2-1}} ;

int  p  =   0 ; // 对解的个数计数

// 判断是否合法
int  canJump( int  x,  int  y)
{
    
if (x >= 0 && x < N && y >= 0 && y < N && chess[x][y] == 0)
        
return 1;
    
return 0;
}

// 求权值
int  weightStep( int  x,  int  y)
{
    
int count = 0;
    
for (int i = 0; i < 8++i)
    
{
        
if (canJump(x + Jump[i][0], y + Jump[i][1]))
            count
++;
    }

    
return count;
}

// 权值排序
void  inssort( int  a[],  int  b[],  int  n)
{
    
if (n <= 0return;
    
for (int i = 0; i < n; ++i)
    
{
        
for (int j = i; j > 0--j)
        
{
            
if (a[j] < a[j - 1])
            
{
                
int temp = a[j - 1];
                a[j 
- 1= a[j];
                a[j] 
= temp;
                
                temp 
= b[j - 1];
                b[j 
- 1= b[j];
                b[j] 
= temp;
            }

        }

    }

}

// 输出结果
void  OutPutChess()
{
    
int i;
    
for (i = 1; i <= N * N - 1++i)
        
{
            printf(
"%d ", step[i]);
        }

    printf(
"\n");
    
for(i=0;i<N;i++)
    
{
        
for(int j=0;j<N;j++)
            printf(
"%3d ",chess[i][j]);
        printf(
"\n");
    }

}

// 回溯算法
void  BackTrace( int  t,  int  x,  int  y)
{
    
if (t >= N * N)
    
{
        p
++;
        OutPutChess();
//输出结果
        exit(1);   //求得一个解时,退出程序
        
//return;  //求得所有解
    }

    
else
    
{
        
int i;
        
int count[8], possibleSteps[8];
        
int k = 0;
        
for (i = 0; i < 8++i)
        
{
            
if (canJump(x + Jump[i][0], y + Jump[i][1]))
            
{
                count[k] 
= weightStep(x + Jump[i][0], y + Jump[i][1]); //求权值
                possibleSteps[k++= i;   //保存下一个结点的序号
            }

        }

        
        inssort(count, possibleSteps, k);
//排序

        
for (i = 0; i < k; ++i)
        
{
            
int d = possibleSteps[i];
            
//跳向下一个结点
            x += Jump[d][0];
            y 
+= Jump[d][1];

            chess[x][y] 
= t + 1;
            step[t] 
= d;
            BackTrace(t 
+ 1, x, y);  //递归调用
            
//回溯
            chess[x][y] = 0;
            x 
-= Jump[d][0];
            y 
-= Jump[d][1];
            
        }

    }

}


int  main()
{
    
    
int x = 0;
    
int y = 0;
    chess[x][y] 
= 1;
    BackTrace(
1, x, y);
    printf(
"All Results Number = %d ", p);
}


如果如下:

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值