经验首页 前端设计 程序设计 Java相关 移动开发 数据库/运维 软件/图像 大数据/云计算 其他经验
当前位置:技术经验 » 程序设计 » C++ » 查看文章
BZOJ4598: [Sdoi2016]模式字符串(点分治 hash)
来源:cnblogs  作者:自为风月马前卒  时间:2018/12/7 9:33:28  对本文有异议

题意

题目链接

Sol

直接考虑点分治+hash匹配

\(up[i]\)表示\(dep \% M = i\)的从下往上恰好与前\(i\)位匹配的个数

\(down\)表示\(dep \% M = i\)的从上往下恰好与后\(i\)位匹配的个数

暴力转移即可

复杂度:\(O(nlog^2n)??\)

代码写起来有一车边界

  1. #include<bits/stdc++.h>
  2. #define ull unsigned long long
  3. #define LL long long
  4. #define int long long
  5. #define siz(v) ((int)v.size())
  6. using namespace std;
  7. const int MAXN = 1e6 + 10, INF = 1e10 + 10;
  8. inline int read() {
  9. char c = getchar(); int x = 0, f = 1;
  10. while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
  11. while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
  12. return x * f;
  13. }
  14. int N, M, Root, siz[MAXN], mx[MAXN], Siz, dep[MAXN], up[MAXN], down[MAXN], su[MAXN], sd[MAXN];
  15. LL ans;
  16. bool det[MAXN];
  17. char a[MAXN], b[MAXN];
  18. vector<int> v[MAXN];
  19. ull hs[MAXN], hp[MAXN], po[MAXN], base = 1331;
  20. map<ull, bool> mp;
  21. void FindRoot(int x, int fa) {
  22. siz[x] = 1; mx[x] = 0;
  23. for(int i = 0; i < siz(v[x]); i++) {
  24. int to = v[x][i];
  25. if(to == fa || det[to]) continue;
  26. FindRoot(to, x);
  27. siz[x] += siz[to]; mx[x] = max(mx[x], siz[to]);
  28. }
  29. mx[x] = max(mx[x], Siz - siz[x]);
  30. if(mx[x] < mx[Root]) Root = x;
  31. }
  32. int dfs(int x, int fa, ull now) {
  33. siz[x] = 1;
  34. dep[x] = dep[fa] + 1;
  35. now = now * base + a[x];
  36. if(hp[dep[x]] == now) up[(dep[x] - 1) % M + 1]++, ans += sd[M - (dep[x] - 1) % M];
  37. if(hs[dep[x]] == now) down[(dep[x] - 1) % M + 1]++, ans += su[M - (dep[x] - 1) % M];
  38. // printf("%d %d\n", x, ans);
  39. int td =1;
  40. for(int i = 0; i < siz(v[x]); i++) {
  41. int to = v[x][i];
  42. if(to == fa || det[to]) continue;
  43. td = max(td, dfs(to, x, now) + 1);
  44. siz[x] += siz[to];
  45. }
  46. return td;
  47. }
  48. void work(int x) {
  49. int tk = 0, tmp = 0;
  50. det[x] = 1; dep[x] = 1; su[1] = sd[1] = 1;//tag;
  51. for(int i = 0; i < siz(v[x]); i++) {
  52. int to = v[x][i];
  53. if(det[to]) continue;
  54. tk = min(M, dfs(to, x, a[x]) + 1), tmp = max(tmp, tk);
  55. for(int j = 1; j <= tk; j++) su[j] += up[j], sd[j] += down[j], up[j] = down[j] = 0;
  56. }
  57. for(int i = 1; i <= tmp; i++) su[i] = sd[i] = 0;
  58. for(int i = 0; i < siz(v[x]); i++) {
  59. int to = v[x][i];
  60. if(to == x || det[to]) continue;
  61. Siz = siz[to]; Root = 0; FindRoot(to, x);
  62. work(Root);
  63. }
  64. }
  65. void init() {
  66. for(int i = 1; i <= N; i++) v[i].clear();
  67. memset(det, 0, sizeof(det));
  68. memset(siz, 0, sizeof(siz));
  69. memset(mx, 0, sizeof(mx));
  70. ans = 0;
  71. }
  72. void solve() {
  73. N = read(); M = read();
  74. init();
  75. scanf("%s", a + 1);
  76. for(int i = 1; i <= N - 1; i++) {
  77. int x = read(), y = read();
  78. v[x].push_back(y); v[y].push_back(x);
  79. }
  80. for(int i = 1; i <= N; i++) reverse(v[i].begin(), v[i].end());
  81. scanf("%s", b + 1); po[0] = 1;
  82. for(int i = 1; i <= N; i++) {
  83. hp[i] = hp[i - 1] + b[(i - 1) % M + 1] * po[i - 1];
  84. hs[i] = hs[i - 1] + b[M - (i - 1) % M] * po[i - 1];
  85. po[i] = base * po[i - 1];
  86. }
  87. Siz = N; mx[0] = INF; Root = 0; FindRoot(1, 0);
  88. work(1);
  89. printf("%d\n", ans);
  90. }
  91. signed main() {
  92. freopen("a.in", "r", stdin);
  93. for(int T = read(); T; T--, solve());
  94. return 0;
  95. }
 友情链接:直通硅谷  点职佳  北美留学生论坛

本站QQ群:前端 618073944 | Java 606181507 | Python 626812652 | C/C++ 612253063 | 微信 634508462 | 苹果 692586424 | C#/.net 182808419 | PHP 305140648 | 运维 608723728

W3xue 的所有内容仅供测试,对任何法律问题及风险不承担任何责任。通过使用本站内容随之而来的风险与本站无关。
关于我们  |  意见建议  |  捐助我们  |  报错有奖  |  广告合作、友情链接(目前9元/月)请联系QQ:27243702 沸活量
皖ICP备17017327号-2 皖公网安备34020702000426号