NYOJ 492 King (状态压缩)

举报
Linux猿 发表于 2021/08/05 00:12:23 2021/08/05
【摘要】 做题感悟:做完这题发现状态压缩有许多需要优化的地方。 解题思路:状态压缩                开始自己用的一般的思路,就和炮兵阵地,郑厂长等题类似的方法做的,开始超时,然后把数组开到了最小的极限就险过。然后看了别人的代码感觉需要优化(注意)的地方太多了。   &nb...


做题感悟:做完这题发现状态压缩有许多需要优化的地方。

解题思路:状态压缩

               开始自己用的一般的思路,就和炮兵阵地,郑厂长等题类似的方法做的,开始超时,然后把数组开到了最小的极限就险过。然后看了别人的代码感觉需要优化(注意)的地方太多了。

               首先我们这题可以预处理出来上下两行对应的合法状态,这样我们就确定下来上下两行对应的状态了,这是第一步的优化 。因为当前行只与上一行有关(这里只考虑当前行对上一行的影响就可以),我们可以把 dp 第一维开成 2 的就可以了,这是第二步的优化,虽然这步作用有点小。最后 ,我们还可以抛弃一些不合法的输入 ,假如是 n * n 的棋盘,一行最多可以防止(n+1)/ 2 个,然后 n 行最多可以放置 ( n + 1 ) / 2 * ( n + 1 ) / 2 个。 为什么是 ( n + 1 ) / 2 行放置呢 ? 因为每个格子攻击相邻的八个格子,不要用n / 2 ,因为如果 奇数行的话就少算了一行。

               这题也可以用DFS递推和铺方格差不多,这样也相当于预处理出上下两行对应的状态。

代码:


  
  1. #include<iostream>
  2. #include<sstream>
  3. #include<map>
  4. #include<cmath>
  5. #include<fstream>
  6. #include<queue>
  7. #include<vector>
  8. #include<sstream>
  9. #include<cstring>
  10. #include<cstdio>
  11. #include<stack>
  12. #include<bitset>
  13. #include<ctime>
  14. #include<string>
  15. #include<cctype>
  16. #include<iomanip>
  17. #include<algorithm>
  18. using namespace std ;
  19. #define INT long long int
  20. #define L(x) (x * 2)
  21. #define R(x) (x * 2 + 1)
  22. const int INF = 0x3f3f3f3f ;
  23. const double esp = 0.0000000001 ;
  24. const double PI = acos(-1.0) ;
  25. const int mod = 1e9 + 7 ;
  26. const int MY = 1400 + 5 ;
  27. const int MX = 150 ;
  28. int n ,m ,num ,top ;
  29. INT dp[2][MX][101] ;
  30. int key[MX] ,h[MX] ;
  31. struct node
  32. {
  33. int x ,y ,c ;
  34. }P[MY] ;
  35. void init() // 预处理合法状态
  36. {
  37. for(int S = 0 ;S < (1<<n) ; ++S)
  38. if(!(S&(S>>1)) && !(S&(S<<1)))
  39. {
  40. int nx = 0 ;
  41. for(int i = 0 ;i < n ; ++i)
  42. if(S&(1<<i))
  43. nx++ ;
  44. key[num] = S ;
  45. h[num++] = nx ;
  46. }
  47. }
  48. int main()
  49. {
  50. while(~scanf("%d%d" ,&n ,&m))
  51. {
  52. if(((n+1)/2)*((n+1)/2) < m) // 抛弃不合法的状态
  53. {
  54. puts("0") ;
  55. continue ;
  56. }
  57. num = 0 ; top = 0 ;
  58. init() ; // 预处理合法状态
  59. for(int i = 0 ;i < num ; ++i) // 预处理上下两行对应的合法状态
  60. for(int j = 0 ;j < num ; ++j)
  61. if(!(key[i]&(key[j]>>1)) && !(key[i]&(key[j]<<1)) && !(key[i]&key[j]))
  62. {
  63. P[top].x = i ;
  64. P[top].y = j ;
  65. P[top++].c = h[j] ;
  66. }
  67. memset(dp[0] ,0 ,sizeof(dp[0])) ;
  68. dp[0][0][0] = 1 ;
  69. for(int i = 1 ;i <= n ; ++i)
  70. {
  71. memset(dp[i&1] ,0 ,sizeof(dp[i&1])) ;
  72. for(int j = 0 ;j < top ; ++j) // 上一行的合法状态
  73. for(int t = 0 ;t <= m ; ++t) // 安放的个数
  74. if(t + P[j].c <= m) // 个数限制
  75. dp[i&1][P[j].y][t+P[j].c] += dp[(i+1)&1][P[j].x][t] ;
  76. }
  77. INT ans = 0 ;
  78. for(int i = 0 ;i < num ; ++i)
  79. ans += dp[n&1][i][m] ;
  80. printf("%lld\n" ,ans) ;
  81. }
  82. return 0 ;
  83. }


文章来源: blog.csdn.net,作者:Linux猿,版权归原作者所有,如需转载,请联系作者。

原文链接:blog.csdn.net/nyist_zxp/article/details/40142521

【版权声明】本文为华为云社区用户转载文章,如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱: cloudbbs@huaweicloud.com
  • 点赞
  • 收藏
  • 关注作者

评论(0

0/1000
抱歉,系统识别当前为高风险访问,暂不支持该操作

全部回复

上滑加载中

设置昵称

在此一键设置昵称,即可参与社区互动!

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。