【POJ 3415】Common Substrings 长度不小于k的公共子串的个数
2016-04-22 11:46
465 查看
长度不小于k的公共子串的个数,论文里有题解,卡了一上午,因为sum没开long long!!!
没开long long毁一生again~~~
以后应该早看POJ里的Discuss啊QAQ
没开long long,毁我青春,耗我钱财,颓我精神==
没开long long毁一生again~~~
以后应该早看POJ里的Discuss啊QAQ
#include<cstdio> #include<cstring> #include<iostream> #include<algorithm> using namespace std; const int N = 200003; int t1 , t2 , c ; void st(int *x, int *y, int *sa, int n, int m) { int i; for(i = 0; i < m; ++i) c[i] = 0; for(i = 0; i < n; ++i) ++c[x[y[i]]]; for(i = 1; i < m; ++i) c[i] += c[i - 1]; for(i = n - 1; i >= 0; --i) sa[--c[x[y[i]]]] = y[i]; } void mkhz(int *a, int *sa, int n, int m) { int *t, *x = t1, *y = t2, i, j, p; for(i = 0; i < n; ++i) x[i] = a[i], y[i] = i; st(x, y, sa, n, m); for(p = 1, j = 1; p < n; j <<= 1, m = p) { for(p = 0, i = n - j; i < n; ++i) y[p++] = i; for(i = 0; i < n; ++i) if (sa[i] >= j) y[p++] = sa[i] - j; st(x, y, sa, n, m); for(t = x, x = y, y = t, p = 1, x[sa[0]] = 0, i = 1; i < n; ++i) x[sa[i]] = y[sa[i]] == y[sa[i - 1]] && y[sa[i] + j] == y[sa[i - 1] + j] ? p - 1 : p++; } } void mkh(int *r, int *sa, int *rank, int *h, int n) { int i, j, k = 0; for(i = 1; i <= n; ++i) rank[sa[i]] = i; for(i = 1; i <= n; h[rank[i++]] = k) for(k ? --k : 0, j = sa[rank[i] - 1]; r[i + k] == r[j + k]; ++k); } char s ; int a , sa , rank , h , n, k, sta , top, tmp, mark ; long long sum ; int main() { while (scanf("%d", &k), k) { scanf("%s", s + 1); tmp = strlen(s + 1); s[++tmp] = 1; scanf("%s", s + tmp + 1); n = strlen(s + 1); for(int i = 1; i <= n; ++i) a[i] = s[i]; mkhz(a, sa, n + 1, 130); mkh(a, sa, rank, h, n); long long ans = 0; for(int i = 1; i <= n; ++i) { if (h[i] < k) { top = sum[0] = sum[1] = 0; } else { for(int j = top; sta[j] > h[i] - k + 1 && j; --j) { sum[mark[j]] += h[i] - k + 1 - sta[j]; sta[j] = h[i] - k + 1; } sta[++top] = h[i] - k + 1; if (sa[i - 1] < tmp) mark[top] = 0; if (sa[i - 1] > tmp) mark[top] = 1; sum[mark[top]] += h[i] - k + 1; if (sa[i] < tmp) ans += sum[1]; if (sa[i] > tmp) ans += sum[0]; } } printf("%lld\n", ans); } return 0; }
没开long long,毁我青春,耗我钱财,颓我精神==
相关文章推荐
- WEB系统架构
- python urllib、urlparse、urllib2、cookielib
- hibernate---联合主键关联
- 自定义初学5——自定义View显示图片
- 23种设计模式分类+SOLID设计原则+从设计模式角度看MVC框架
- 小米引爆用户参与感有哪些值得借签的法则?
- C++第4次实验项目2-学生类
- android layout_gravity和gravity的区别
- 《疯狂Java讲义(第3版)》.(李刚)——运算符
- 从APP产品原型七步搞定APP设计
- Hibernate单表操作
- java 连接 linux redis 失
- IsNumeric 函数
- php 多页面间共享变量
- JVM——GC
- 关于测试工具以及前端性能测试的一些思考
- Linux CentOS 7 & Apache Tomcat 7 安装与配置
- Ubuntu 12.04下安装OpenCV 2.4.2 和 ffmpeg-0.11.1
- 知道这20个正则表达式,能让你少写1,000行代码
- 数位DP