找回密码
 立即注册→加入我们

QQ登录

只需一步,快速开始

搜索
热搜: 下载 VB C 实现 编写
查看: 2642|回复: 0

qq空间老贴转载-骑士周游问题我的做法

[复制链接]
发表于 2014-2-19 19:20:25 | 显示全部楼层 |阅读模式

欢迎访问技术宅的结界,请注册或者登录吧。

您需要 登录 才可以下载或查看,没有账号?立即注册→加入我们

×
自己写的,也做了优化但效果不大,求出所有结果是世界级的难题,貌似没什么好办法
SIN<8时很快 >=8很慢,无结果,

改进法在后面
  1. [/align][align=left]#include <conio.h>
  2. #include <iostream>
  3. #include <time.h>[/align][align=left]using namespace std;[/align][align=left]#define SIN 8
  4. #define MAX SIN*SIN
  5. int sum=0;[/align][align=left]
  6. struct POINT
  7. {
  8. int x;
  9. int y;
  10. };[/align]
  11. [align=left]class HorsePuzzle
  12. {
  13. POINT horse[MAX];
  14. bool  permit[SIN+4][SIN+4];//防越界
  15. public:
  16. HorsePuzzle();
  17. void printOut();
  18. bool IsValid(int n);
  19. void placeHorse(int n,int x,int y);
  20. };[/align][align=left]HorsePuzzle::HorsePuzzle()
  21. {
  22. int i,j;
  23. for(i=0;i<MAX;i++)
  24. {
  25.   horse[i].x=0;
  26.   horse[i].y=0;
  27.   *((bool*)permit+i)=true;
  28. }
  29. for(i=2;i<SIN+2;i++)
  30. {
  31.   for(j=0;j<SIN+2;j++)
  32.   {
  33.    permit[i][j]=true;//中间位置允许放置
  34.   }
  35. }
  36. //0/1/SIN+2/SIN+3行列不允许放置
  37. for(i=0;i<SIN+4;i++)
  38. {
  39.   permit[0][i]=false;
  40.   permit[1][i]=false;
  41.   permit[SIN+2][i]=false;
  42.   permit[SIN+3][i]=false;
  43.   permit[i][0]=false;
  44.   permit[i][1]=false;
  45.   permit[i][SIN+2]=false;
  46.   permit[i][SIN+3]=false;
  47. }
  48. }[/align][align=left]void HorsePuzzle::printOut()
  49. {
  50. int data[SIN][SIN],i,j;
  51. cout<<"解法"<<++sum<<':'<<endl;
  52. for(i=0;i<MAX;i++)
  53. {
  54.   data[horse[i].x][horse[i].y]=i;
  55. }
  56. for(i=0;i<SIN;i++)
  57. {
  58.   for(j=0;j<SIN;j++)
  59.   {
  60.    cout<<data[i][j]<<'\t';
  61.   }
  62.   cout<<endl;
  63. }
  64. cout<<endl;
  65. }[/align][align=left]void HorsePuzzle::placeHorse(int n,int x,int y)
  66. {
  67. horse[n].x=x;
  68. horse[n].y=y;
  69. if(permit[horse[n].x+2][horse[n].y+2])
  70. {
  71.   permit[x+2][y+2]=false;
  72.   n++;
  73.         if(n == MAX)
  74.   {
  75.    printOut();
  76.    permit[x+2][y+2]=true;
  77.    return;
  78.   }[/align][align=left]  placeHorse(n,x-2,y-1);
  79.   placeHorse(n,x-1,y-2);
  80.   placeHorse(n,x+1,y-2);
  81.   placeHorse(n,x+2,y-1);
  82.   placeHorse(n,x-2,y+1);
  83.   placeHorse(n,x-1,y+2);
  84.   placeHorse(n,x+1,y+2);
  85.   placeHorse(n,x+2,y+1);[/align][align=left]  permit[x+2][y+2]=true;
  86. }
  87. }[/align][align=left]void main()
  88. {
  89. time_t begin=time(NULL);
  90. HorsePuzzle horse;
  91. horse.placeHorse(0,0,0);
  92. cout<<"共"<<sum<<"组解,时间:"<<time(NULL)-begin<<endl;
  93. }[/align][align=left]
复制代码
今天看了《c/c++常用算法》这本书,10.8节也是骑士周游问题,我把他的代码改成通用的,数组维数可以变化,代码如下:
  1. [/align]#include <stdio.h>
  2. #include <stdlib.h>
  3. #define MYSIZE 7
  4. typedef struct
  5. {
  6. int x;
  7. int y;
  8. }Coordinate;
  9. int chessboard[MYSIZE][MYSIZE];
  10. int curstep;
  11. Coordinate fangxiang[8]={{-2,1},{-1,2},{1,2},{2,1},{2,-1},{1,-2},{-1,-2},{-2,-1}};
  12. void Move(Coordinate curpos)
  13. {
  14. Coordinate next;
  15. int i,j;
  16. if(curpos.x < 0 || curpos.x > MYSIZE-1 || curpos.y < 0 || curpos.y > MYSIZE-1)
  17. return;
  18. if(chessboard[curpos.x][curpos.y])
  19. return;
  20. chessboard[curpos.x][curpos.y]=curstep;
  21. curstep++;
  22. if(curstep > MYSIZE*MYSIZE)
  23. {
  24. for(i=0;i<MYSIZE;i++)
  25. {
  26. for(j=0;j<MYSIZE;j++)
  27. {
  28. printf("%5d",chessboard[i][j]);
  29. }
  30. printf("\n");
  31. }
  32. exit(0);
  33. }
  34. else
  35. {
  36. for(i=0;i<8;i++)
  37. {
  38. next.x=curpos.x+fangxiang[i].x;
  39. next.y=curpos.y+fangxiang[i].y;
  40. if(next.x < 0 || next.x > MYSIZE-1 || next.y < 0 || next.y > MYSIZE-1)
  41. {

  42. }
  43. else
  44. {
  45. Move(next);
  46. }
  47. }
  48. }
  49. chessboard[curpos.x][curpos.y]=0;
  50. curstep--;
  51. }

  52. void main()
  53. {
  54. int i,j;
  55. Coordinate start;
  56. scanf("%d%d",&start.x,&start.y);
  57. if(start.x < 1 || start.y < 1 || start.x > MYSIZE || start.y > MYSIZE)
  58. {
  59. printf("error");
  60. exit(0);
  61. }
  62. for(i=0;i<MYSIZE;i++)
  63. {
  64. for(j=0;j<MYSIZE;j++)
  65. {
  66. chessboard[i][j]=0;
  67. }
  68. }
  69. start.x--;
  70. start.y--;
  71. curstep=1;
  72. Move(start);
  73. }
  74. [align=left]
复制代码

结果发现他的程序居然在MYSIZE=8时很快找到了答案。。。百思不得其解,遂看了一下他的程序。
乍一看,觉得写得很一般,效率应该没我的高,验证了一下MYSIZE=7的时候都卡,证明了我的想法
,因为我用的判断直接是判断数组,不需要函数,于是我找原因,
开始的时候我发现他传的参数类型是coordinate,在底层实现的时候,其实传的参数是2个int了。考虑到压栈可能耗费时间,
所以我把我的程序中的placeHorse递归也从3个参数改成2个参数,结果发现还是没结果,我又找,发现书中代码有fangxiang这个数组,
用来存储8个方向,在函数中以for循环遍历,我想,我直接手写的,你用for循环永远不可能比我快啊,,不解,肯定还有猫腻,
果然不出我所料,后来我发现 数组为{-2,1},{-1,2},{1,2},{2,1},{2,-1},{1,-2},{-1,-2},{-2,-1},这8个方向为顺时针排列,再看我的,明白了,
原来出结果的速度与这个顺序也是有关的,打个比方,起始点选0,0的话,如果前几个搜索方向为横向和纵向数轴的相反方向,则由于方向问题导致需要遍历到根部,而如果为数轴正方向的话,开始的搜索方向总是正确的,循环是从开始遍历的,所以开始错的越少,效率越高,从根部循环越容易找到。。。 最后换用他的数组,很快得到结果。
不过话说到这里,要想遍历所有解,还是遍历不出来的,或者说很慢很慢,世界难题没那么好解。
还有,给定初始位置,其实可以选最合适的初始方向的。。。不过我没继续研究了到此为止,最后最优代码如下:
  1. [/align][align=left]#include <conio.h>[/align]#include <iostream>
  2. #include <time.h>
  3. using namespace std;

  4. #define SIN 8
  5. #define MAX SIN*SIN
  6. int sum=0;

  7. struct POINT
  8. {
  9. int x;
  10. int y;
  11. };

  12. class HorsePuzzle
  13. {
  14. POINT horse[MAX];
  15. bool  permit[SIN+4][SIN+4];//防越界
  16. public:
  17. HorsePuzzle();
  18. void printOut();
  19. void placeHorse(/*int n,*/int x,int y);
  20. };

  21. HorsePuzzle::HorsePuzzle()
  22. {
  23. int i,j;
  24. for(i=0;i<MAX;i++)
  25. {
  26. horse[i].x=0;
  27. horse[i].y=0;
  28. *((bool*)permit+i)=true;
  29. }
  30. for(i=2;i<SIN+2;i++)
  31. {
  32. for(j=0;j<SIN+2;j++)
  33. {
  34. permit[i][j]=true;//中间位置允许放置
  35. }
  36. }
  37. //0/1/SIN+2/SIN+3行列不允许放置
  38. for(i=0;i<SIN+4;i++)
  39. {
  40. permit[0][i]=false;
  41. permit[1][i]=false;
  42. permit[SIN+2][i]=false;
  43. permit[SIN+3][i]=false;
  44. permit[i][0]=false;
  45. permit[i][1]=false;
  46. permit[i][SIN+2]=false;
  47. permit[i][SIN+3]=false;
  48. }
  49. }

  50. void HorsePuzzle::printOut()
  51. {
  52. int data[SIN][SIN],i,j;
  53. cout<<"解法"<<++sum<<':'<<endl;
  54. for(i=0;i<MAX;i++)
  55. {
  56. data[horse[i].x][horse[i].y]=i;
  57. }
  58. for(i=0;i<SIN;i++)
  59. {
  60. for(j=0;j<SIN;j++)
  61. {
  62. cout<<data[i][j]<<'\t';
  63. }
  64. cout<<endl;
  65. }
  66. cout<<endl;
  67. }

  68. int n=0;

  69. void HorsePuzzle::placeHorse(int x,int y)
  70. {
  71. if(permit[x+2][y+2])
  72. {
  73. horse[n].x=x;
  74. horse[n].y=y;
  75. n++;
  76. permit[x+2][y+2]=false;

  77.         if(n != MAX)
  78. {
  79. placeHorse(x-2,y+1);
  80. placeHorse(x-1,y+2);
  81. placeHorse(x+1,y+2);
  82. placeHorse(x+2,y+1);
  83. placeHorse(x+2,y-1);
  84. placeHorse(x+1,y-2);
  85. placeHorse(x-1,y-2);
  86. placeHorse(x-2,y-1);
  87. }
  88. else
  89. {
  90. printOut();
  91. }

  92. permit[x+2][y+2]=true;
  93. n--;
  94. }
  95. }

  96. void main()
  97. {
  98. time_t begin=time(NULL);
  99. HorsePuzzle horse;
  100. horse.placeHorse(0,0);
  101. cout<<"共"<<sum<<"组解,时间:"<<time(NULL)-begin<<endl;
  102. }

  103. [align=left]
复制代码
如果你还能提出效率高出很多的建议,并可以实现求所有解的更快方法,欢迎和我讨论!!!


回复

使用道具 举报

本版积分规则

QQ|Archiver|小黑屋|技术宅的结界 ( 滇ICP备16008837号 )|网站地图

GMT+8, 2024-11-21 21:51 , Processed in 0.034024 second(s), 25 queries , Gzip On.

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表