0047算法笔记——【随机化算法】拉斯维加斯(Las Vegas)算法和n后问题

      1、拉斯维加斯(Las Vegas)算法

     拉斯维加斯算法不会得到不正确的解。一旦用拉斯维加斯算法找到一个解，这个解就一定是正确解。但有时用拉斯维加斯算法找不到解。与蒙特卡罗算法类似，拉斯维加斯算法找到正确解的概率随着它所用的计算时间的增加而提高。对于所求解问题的任一实例，用同一拉斯维加斯算法反复对该实例求解足够多次，可使求解失败的概率任意小。拉斯维加斯算法的一个显著特征是它所作的随机性决策有可能导致算法找不到所需的解。

[cpp] view plaincopyprint?

1. void obstinate(Object x, Object y)  
2. {// 反复调用拉斯维加斯算法LV(x,y)，直到找到问题的一个解y   
3.     bool success= false;  
4.     while (!success) success=lv(x,y);  
5. }  

void obstinate(Object x, Object y){// 反复调用拉斯维加斯算法LV(x,y)，直到找到问题的一个解y    bool success= false;    while (!success) success=lv(x,y);}

      设p(x)是对输入x调用拉斯维加斯算法获得问题的一个解的概率。一个正确的拉斯维加斯算法应该对所有输入x均有p(x)>0。设t(x)是算法obstinate找到具体实例x的一个解所需的平均时间 ,s(x)和e(x)分别是算法对于具体实例x求解成功或求解失败所需的平均时间，则有。解此方程得：

     2、n后问题

     问题描速：在n×n格的棋盘上放置彼此不受攻击的n个皇后。按照国际象棋的规则，皇后可以攻击与之处在同一行或同一列或同一斜线上的棋子。n后问题等价于在n×n格的棋盘上放置n个皇后，任何2个皇后不放在同一行或同一列或同一斜线上。

     1)纯拉斯维加斯随机算法求解思路

      对于n后问题的任何一个解而言，每一个皇后在棋盘上的位置无任何规律，不具有系统性，而更象是随机放置的。在棋盘上相继的各行中随机地放置皇后，并注意使新放置的皇后与已放置的皇后互不攻击，直至n个皇后均已相容地放置好，或已没有下一个皇后的可放置位置时为止。

      具体实现代码如下：

     1、RandomNumber.h

[cpp] view plaincopyprint?

1. #include"time.h"   
2. //随机数类   
3. const unsigned long maxshort = 65536L;  
4. const unsigned long multiplier = 1194211693L;  
5. const unsigned long adder = 12345L;  
6.   
7. class RandomNumber  
8. {  
9.     private:  
10.         //当前种子   
11.         unsigned long randSeed;  
12.     public:  
13.         RandomNumber(unsigned long s = 0);//构造函数，默认值0表示由系统自动产生种子  
14.         unsigned short Random(unsigned long n);//产生0：n-1之间的随机整数  
15.         double fRandom(void);//产生[0,1)之间的随机实数  
16. };  
17.   
18. RandomNumber::RandomNumber(unsigned long s)//产生种子  
19. {  
20.     if(s == 0)  
21.     {  
22.         randSeed = time(0);//用系统时间产生种子   
23.     }  
24.     else  
25.     {  
26.         randSeed = s;//由用户提供种子   
27.     }  
28. }  
29.   
30. unsigned short RandomNumber::Random(unsigned long n)//产生0：n-1之间的随机整数  
31. {  
32.     randSeed = multiplier * randSeed + adder;//线性同余式  
33.     return (unsigned short)((randSeed>>16)%n);  
34. }  
35.   
36. double RandomNumber::fRandom(void)//产生[0,1)之间的随机实数  
37. {  
38.     return Random(maxshort)/double(maxshort);  
39. }  

#include"time.h"//随机数类const unsigned long maxshort = 65536L;const unsigned long multiplier = 1194211693L;const unsigned long adder = 12345L;class RandomNumber{private://当前种子unsigned long randSeed;public:RandomNumber(unsigned long s = 0);//构造函数，默认值0表示由系统自动产生种子unsigned short Random(unsigned long n);//产生0：n-1之间的随机整数double fRandom(void);//产生[0,1)之间的随机实数};RandomNumber::RandomNumber(unsigned long s)//产生种子{if(s == 0){randSeed = time(0);//用系统时间产生种子}else{randSeed = s;//由用户提供种子}}unsigned short RandomNumber::Random(unsigned long n)//产生0：n-1之间的随机整数{randSeed = multiplier * randSeed + adder;//线性同余式return (unsigned short)((randSeed>>16)%n);}double RandomNumber::fRandom(void)//产生[0,1)之间的随机实数{return Random(maxshort)/double(maxshort);}

     2、7d4d1.cpp

[cpp] view plaincopyprint?

1. //随机化算法 拉斯维加斯算法 n后问题   
2. #include "stdafx.h"   
3. #include "RandomNumber.h"  
4. #include <cmath>   
5. #include <iostream>   
6. using namespace std;  
7.   
8. class Queen  
9. {  
10.     friend void nQueen(int);  
11.     private:  
12.         bool Place(int k);      //测试皇后k置于第x[k]列的合法性  
13.         bool QueensLv(void);    //随机放置n个皇后拉斯维加斯算法  
14.         int n;                  //皇后个数  
15.         int *x,*y;              //解向量  
16. };  
17.   
18. //测试皇后k置于第x[k]列的合法性   
19. bool Queen::Place(int k)  
20. {  
21.     for(int j=1; j<k; j++)  
22.     {  
23.         if((abs(k-j)==abs(x[j]-x[k]))||(x[j]==x[k]))  
24.         {  
25.             return false;  
26.         }  
27.     }  
28.     return true;  
29. }  
30.   
31. //随机放置n个皇后的拉斯维加斯算法   
32. bool Queen::QueensLv(void)  
33. {  
34.     RandomNumber rnd;       //随机数产生器   
35.     int k = 1;              //下一个皇后的编号  
36.     int count = 1;          //在一列中，可以放置皇后的个数  
37.   
38.     while((k<=n)&&(count>0))  
39.     {  
40.         count = 0;  
41.         for(int i=1; i<=n; i++)  
42.         {  
43.             x[k] = i;//位置   
44.             if(Place(k))  
45.             {  
46.                 y[count++] = i;  
47.             }  
48.         }  
49.         if(count>0)  
50.         {  
51.             x[k++] = y[rnd.Random(count)];      //随机位置  
52.         }  
53.     }  
54.   
55.     return (count>0);        //count>0 表示放置成功  
56. }  
57.   
58. //解n后问题的拉斯维加斯算法   
59. void nQueen(int n)  
60. {  
61.     Queen X;  
62.     X.n = n;  
63.   
64.     int *p = new int[n+1];  
65.     for(int i=0; i<=n; i++)  
66.     {  
67.         p[i] = 0;  
68.     }  
69.     X.x = p;  
70.     X.y = new int[n+1];  
71.   
72.     //反复调用随机放置n个皇后的拉斯维加斯算法，直到放置成功   
73.     while(!X.QueensLv());  
74.   
75.     for(int i=1; i<=n; i++)  
76.     {  
77.         cout<<p[i]<<" ";  
78.     }  
79.     cout<<endl;  
80.     delete []p;  
81. }  
82.   
83. int main()  
84. {  
85.     int n=8;    
86.     cout<<n<<"皇后问题的解为："<<endl;    
87.     nQueen(n);    
88.     return 0;    
89. }  

//随机化算法 拉斯维加斯算法 n后问题#include "stdafx.h"#include "RandomNumber.h"#include <cmath>#include <iostream>using namespace std;class Queen{friend void nQueen(int);private:bool Place(int k);//测试皇后k置于第x[k]列的合法性bool QueensLv(void);//随机放置n个皇后拉斯维加斯算法int n;//皇后个数int *x,*y;//解向量};//测试皇后k置于第x[k]列的合法性bool Queen::Place(int k){for(int j=1; j<k; j++){if((abs(k-j)==abs(x[j]-x[k]))||(x[j]==x[k])){return false;}}return true;}//随机放置n个皇后的拉斯维加斯算法bool Queen::QueensLv(void){RandomNumber rnd;//随机数产生器int k = 1;//下一个皇后的编号int count = 1;//在一列中，可以放置皇后的个数while((k<=n)&&(count>0)){count = 0;for(int i=1; i<=n; i++){x[k] = i;//位置if(Place(k)){y[count++] = i;}}if(count>0){x[k++] = y[rnd.Random(count)];//随机位置}}return (count>0);//count>0 表示放置成功}//解n后问题的拉斯维加斯算法void nQueen(int n){Queen X;X.n = n;int *p = new int[n+1];for(int i=0; i<=n; i++){p[i] = 0;}X.x = p;X.y = new int[n+1];//反复调用随机放置n个皇后的拉斯维加斯算法，直到放置成功while(!X.QueensLv());for(int i=1; i<=n; i++){cout<<p[i]<<" ";}cout<<endl;delete []p;}int main(){int n=8;      cout<<n<<"皇后问题的解为："<<endl;      nQueen(n);      return 0;  }

     程序运行结果如下：

     2)与回溯法结合的拉斯维加斯随机算法求解思路

     如果将上述随机放置策略与回溯法相结合，可能会获得更好的效果。可以先在棋盘的若干行中随机地放置皇后，然后在后继行中用回溯法继续放置，直至找到一个解或宣告失败。随机放置的皇后越多，后继回溯搜索所需的时间就越少，但失败的概率也就越大。 

　 算法具体代码如下：

    1、RandomNumber.h如上

    2、7d4d1-2.cpp

[cpp] view plaincopyprint?

1. //随机化算法 拉斯维加斯算法 n后问题   
2. #include "stdafx.h"   
3. #include "RandomNumber.h"  
4. #include <cmath>   
5. #include <iostream>   
6. using namespace std;  
7.   
8. class Queen  
9. {  
10.     friend bool nQueen(int);  
11.     private:  
12.         bool Place(int k);                  //测试皇后k置于第x[k]的合法性  
13.         bool Backtrack(int t);              //解n后问题的回溯法  
14.         bool QueensLV(int stopVegas);       //随机放置n个皇后拉斯维加斯算法  
15.         int n,*x,*y;  
16. };  
17.   
18. //测试皇后k置于第x[k]列的合法性   
19. bool Queen::Place(int k)  
20. {  
21.     for(int j=1; j<k; j++)  
22.     {  
23.         if((abs(k-j)==abs(x[j]-x[k]))||(x[j]==x[k]))  
24.         {  
25.             return false;  
26.         }  
27.     }  
28.     return true;  
29. }  
30.   
31. //解n后问题的回溯法   
32. bool Queen::Backtrack(int t)  
33. {  
34.     if(t>n)  
35.     {  
36.         for(int i=1; i<=n; i++)  
37.         {  
38.             y[i] = x[i];//问题的解   
39.         }  
40.         return true;  
41.     }  
42.     else  
43.     {  
44.         for(int i=1; i<=n; i++)  
45.         {  
46.             x[t] = i;  
47.             if(Place(t)&&Backtrack(t+1))  
48.             {  
49.                 return true;  
50.             }  
51.         }  
52.     }  
53.     return false;  
54. }  
55.   
56.   
57. //随机放置n个皇后拉斯维加斯算法   
58. bool Queen::QueensLV(int stopVegas)  
59. {  
60.     RandomNumber rnd;       //随机数产生器   
61.     int k = 1;  
62.     int count = 1;  
63.   
64.     //1<=stopVegas<=n   
65.     while((k<=stopVegas)&&(count>0))  
66.     {  
67.         count = 0;  
68.         for(int i=1; i<=n; i++)  
69.         {  
70.             x[k] = i;  
71.             if(Place(k))  
72.             {  
73.                 y[count++] = i;  
74.             }  
75.         }  
76.   
77.         if(count>0)  
78.         {  
79.             x[k++] = y[rnd.Random(count)];//随机位置  
80.         }  
81.     }  
82.     return (count>0);        //count>0表示放置成功  
83. }  
84.   
85. //与回溯法相结合的解n后问题的拉斯维加斯算法   
86. bool nQueen(int n)  
87. {  
88.     Queen X;  
89.       
90.     //初始化X   
91.     X.n = n;  
92.   
93.     int *p = new int[n+1];  
94.     int *q = new int[n+1];  
95.   
96.     for(int i=0; i<=n; i++)  
97.     {  
98.         p[i] = 0;  
99.         q[i] = 0;  
100.     }  
101.   
102.     X.y = p;  
103.     X.x = q;  
104.   
105.     int stop = 3;  
106.     if(n>15)  
107.     {  
108.         stop = n-15;  
109.     }  
110.   
111.     bool found = false;  
112.     while(!X.QueensLV(stop));  
113.   
114.     //算法的回溯搜索部分   
115.     if(X.Backtrack(stop+1))  
116.     {  
117.         for(int i=1; i<=n; i++)  
118.         {  
119.             cout<<p[i]<<" ";  
120.         }  
121.         found = true;  
122.         cout<<endl;  
123.     }  
124.       
125.     delete []p;  
126.     delete []q;  
127.     return found;  
128. }  
129.   
130. int main()  
131. {  
132.     int n=8;    
133.     cout<<n<<"皇后问题的解为："<<endl;    
134.     while(!nQueen(n));    
135.     return 0;    
136. }  

//随机化算法 拉斯维加斯算法 n后问题#include "stdafx.h"#include "RandomNumber.h"#include <cmath>#include <iostream>using namespace std;class Queen{friend bool nQueen(int);private:bool Place(int k);//测试皇后k置于第x[k]的合法性bool Backtrack(int t);//解n后问题的回溯法bool QueensLV(int stopVegas);//随机放置n个皇后拉斯维加斯算法int n,*x,*y;};//测试皇后k置于第x[k]列的合法性bool Queen::Place(int k){for(int j=1; j<k; j++){if((abs(k-j)==abs(x[j]-x[k]))||(x[j]==x[k])){return false;}}return true;}//解n后问题的回溯法bool Queen::Backtrack(int t){if(t>n){for(int i=1; i<=n; i++){y[i] = x[i];//问题的解}return true;}else{for(int i=1; i<=n; i++){x[t] = i;if(Place(t)&&Backtrack(t+1)){return true;}}}return false;}//随机放置n个皇后拉斯维加斯算法bool Queen::QueensLV(int stopVegas){RandomNumber rnd;//随机数产生器int k = 1;int count = 1;//1<=stopVegas<=nwhile((k<=stopVegas)&&(count>0)){count = 0;for(int i=1; i<=n; i++){x[k] = i;if(Place(k)){y[count++] = i;}}if(count>0){x[k++] = y[rnd.Random(count)];//随机位置}}return (count>0);//count>0表示放置成功}//与回溯法相结合的解n后问题的拉斯维加斯算法bool nQueen(int n){Queen X;//初始化XX.n = n;int *p = new int[n+1];int *q = new int[n+1];for(int i=0; i<=n; i++){p[i] = 0;q[i] = 0;}X.y = p;X.x = q;int stop = 3;if(n>15){stop = n-15;}bool found = false;while(!X.QueensLV(stop));//算法的回溯搜索部分if(X.Backtrack(stop+1)){for(int i=1; i<=n; i++){cout<<p[i]<<" ";}found = true;cout<<endl;}delete []p;delete []q;return found;}int main(){int n=8;      cout<<n<<"皇后问题的解为："<<endl;      while(!nQueen(n));      return 0;  }

     程序运行结果如图：