搜索
ACM程序设计杭州电子科技大学刘春英acm@hdu.edu.cn2020/1/12上周，你了吗？2020/1/13每周一星（8）：gaojie2020/1/14第九讲一招制敌之搜索题2020/1/15根据“信息学初学者之家”网站的统计，Ural（俄罗斯的Ural州立大学的简称，有名的UralOnlineProblemSet就是该校的系统）的题目类型大概呈如下的分布：搜索动态规划贪心构造图论约10%约15%约5%约5%约10%计算几何纯数学题数据结构其它约5%约20%约5%约25%统计信息：2020/1/16——摘自《ACM竞赛之新人向导》“算法中最基本和常用的是搜索，这里要说的是，有些初学者在学习这些搜索基本算法是不太注意剪枝，这是十分不可取的，因为所有搜索的题目给你的测试用例都不会有很大的规模，你往往察觉不出程序运行的时间问题，但是真正的测试数据一定能过滤出那些没有剪枝的算法。实际上参赛选手基本上都会使用常用的搜索算法，题目的区分度往往就是建立在诸如剪枝之类的优化上了。”引言2020/1/17什么是搜索算法呢？搜索算法是利用计算机的高性能来有目的地穷举一个问题的部分或所有的可能情况，从而求出问题的解的一种方法。搜索过程实际上是根据初始条件和扩展规则构造一棵解答树并寻找符合目标状态的节点的过程。2020/1/18预热一下：二分查找2345681220324565748695100headmidtail2020/1/19查找示意图：A[1]~A[15]A[1]~A[7]A[9]~A[15]A[1]~A[3]A[5]~A[7]A[1]~A[1]A[3]~A[3]……2020/1/110思考：1、在一百万个元素里查找某个元素大约需要比较多少次？2、时间复杂度：O(logN)2020/1/111举例分析从简单的字符串搜索讲起2020/1/112SampleInput23ABCDBCDFFBRCD2roseorchidSampleOutput22HDOJ_1238Substrings2020/1/113题目分析：这是一道入门级别的搜索题，基本思想比较简单，但是如果用最朴素的算法，可能会超时如何降低算法的复杂度呢？下面的算法如何：先将字符串按长度从短到长排序，枚举最短的字符串的子串，判断是否都是别的字符串的子串，求出最大长度即可。2020/1/114说明：本题除了可以练习基本搜索算法，也是练习字符串处理的好题目，题中用到的相关知识点有：求反串求子串字符串查找求字符串长度强烈推荐！！2020/1/115再来一道数值型搜索题2020/1/116SampleInput5129999999999916805161970112002411000SampleOutput22313313237343433753HDOJ_12392020/1/117获取有用信息a.给定整数m,a,b(4m=100000and1=a=b=1000)b.需要找到两个数(不妨设为p,q)满足以下条件：p,q均为质数；p*q=m;a/b=p/q=1;c.输出所有满足以上条件的p,q中乘积最大的一对p,q2020/1/118算法分析1.典型的搜索从所有可能的p,q中寻找满足条件的一对2.p,q的要求p,q均为质数,且p=q=100000;3.按上述思想流程应为：a.从1—100000中搜出质数b.两层循环，试遍所有的组合(p,q可能相等）c.每种组合去判断是否符合条件，如是，将p*q与当前最大值比较，判断，保存2020/1/119面临的问题：超时！从1—100000的质数运算约为1e+8,而这只是准备工作。因此，如不加以分析简化此题无法在规定时间内出解2020/1/120深入分析：p,q的范围其实可在2—50000(why?)然而，这是最小的范围吗？考虑大于10000的某个质数，不妨设为Q，另一个质数为P，则：如果P10，P/Q0.001如果P10，P*Q100000而考虑到a,b的取值范围(1=a=b=1000)可知min(a/b)=0.001同时，要求：p*q=m=100000所以无论如何质数都不能超过100002020/1/121搜索时的技巧：搜索顺序很重要。建议从大往小搜（num:质数的个数）for(i=num-1;i=0;i--)for(j=i;j=num-1;j++)……注意剪枝：If(a[j]m||a[j]*a[i]m||((double)a[i]/a[j])s)……2020/1/122典型搜索——迷宫搜索2020/1/123SampleInput445S.X...X...XD....345S.X...X....D000SampleOutputNOYESHDOJ_1010TempteroftheBone2020/1/124要点分析：典型的迷宫搜索，熟练掌握该题将具有里程碑式的意义！每个block只能走一次要求恰好某个给定的时间到达出口2020/1/125剪枝条件？如果可走的block的总数小于时间，将会产生什么情况？还想到什么剪枝？2020/1/126奇偶性剪枝可以把map看成这样：010101101010010101101010010101从为0的格子走一步，必然走向为1的格子从为1的格子走一步，必然走向为0的格子即：0-1或1-0必然是奇数步0-0走1-1必然是偶数步结论：所以当遇到从0走向0但是要求时间是奇数的，或者，从1走向0但是要求时间是偶数的都可以直接判断不可达！2020/1/127这个题目没问题了吧？2020/1/128思考：求某给定时间以内能否找到出口找到出口的最短时间条件变为可以停留2020/1/129四、深度优先搜索基本思想：从初始状态S开始，利用规则生成搜索树下一层任一个结点，检查是否出现目标状态G，若未出现，以此状态利用规则生成再下一层任一个结点，再检查是否为目标节点G，若未出现，继续以上操作过程，一直进行到叶节点（即不能再生成新状态节点），当它仍不是目标状态G时，回溯到上一层结果，取另一可能扩展搜索的分支。生成新状态节点。若仍不是目标状态，就按该分支一直扩展到叶节点，若仍不是目标，采用相同的回溯办法回退到上层节点，扩展可能的分支生成新状态，…，一直进行下去，直到找到目标状态G为止。2020/1/130DFS算法（1）把起始节点S线放到OPEN表中。（2）如果OPEN是空表，则失败退出，否则继续。（3）从OPEN表中取最前面的节点node移到CLOSED表中。（4）若node节点是叶结点（若没有后继节点），则转向（2）。（5）扩展node的后继节点，产生全部后继节点，并把他们放在OPEN表的前面。各后继结点指针指向node节点。（6）若后继节点中某一个是目标节点，则找到一个解，成功退出。否则转向（2）循环。2020/1/131三、广度优先搜索基本思想：从初始状态S开始，利用规则，生成所有可能的状态。构成树的下一层节点，检查是否出现目标状态G，若未出现，就对该层所有状态节点，分别顺序利用规则。生成再下一层的所有状态节点，对这一层的所有状态节点检查是否出现G，若未出现，继续按上面思想生成再下一层的所有状态节点，这样一层一层往下展开。直到出现目标状态为止。2020/1/132BFS算法：（1）把起始节点S线放到OPEN表中（2）如果OPEN是空表，则失败退出，否则继续。（3）在OPEN表中取最前面的节点node移到CLOSED表中。（4）扩展node节点。若没有后继（即叶节点），则转向（2）循环。（5）把node的所有后继节点放在OPEN表的末端。各后继结点指针指向node节点。（6）若后继节点中某一个是目标节点，则找到一个解，成功退出。否则转向（2）循环。2020/1/133小结：广度和深度优先搜索有一个很大的缺陷,就是他们都是在一个给定的状态空间中穷举。这在状态空间不大的情况下是很合适的算法，可是当状态空间十分大，且不预测的情况下就不可取了。他的效率实在太低，甚至不可完成。所以，在这里再次强调“剪枝”！2020/1/134附录：推荐搜索题：1010、1240、1241、12421072、1253、1312、1372(以上题目类似于1010)1238、1239、1015、10161401、1515、15482020/1/135附:参考源码（HDOJ_1010）附录：hdoj_1010月下版#includeiostream.h#includestring.h#includestdlib.hcharmap[9][9];intn,m,t,di,dj;boolescape;intdir[4][2]={{0,-1},{0,1},{1,0},{-1,0}};voiddfs(intsi,intsj,intcnt){inti,temp;if(sin||sjm||si=0||sj=0)return;if(cnt==t&&si==di&&sj==dj)escape=1;if(escape)return;temp=(t-cnt)-abs(si-di)-abs(sj-dj);if(temp0||temp&1)return;for(i=0;i4;i++){if(map[si+dir[i][0]][sj+dir[i][1]]!='X'){map[si+dir[i][0]][sj+dir[i][1]]='X';dfs(si+dir[i][0],sj+dir[i][1],cnt+1);map[si+dir[i][0]][sj+dir[i][1]]='.';}}return;}intmain(){inti,j,si,sj;while(cinnmt){if(n==0&&m==0&&t==0)break;intwall=0;for(i=1;i=n;i++)for(j=1;j=m;j++){cinmap[i][j];if(map[i][j]=='S'){si=i;sj=j;}elseif(map[i][j]=='D'){di=i;dj=j;}elseif(map[i][j]=='X')wall++;}if(n*m-wall=t){coutNOendl;continue;}escape=0;map[si][sj]='X';dfs(si,sj,0);if(escape)coutYESendl;elsecoutNOendl;}return0;}2020/1/136下一讲：？？？2020/1/137WelcometoHDOJThankYou~