【POJ 3415】Common Substrings
【链接】h在这里写链接
【题意】
相同的重复计数。
【题解】
然后,对于这一个字符串,求出它的Height和Sa数组。
然后,把height 数组【连续】大于等于k的,分成一组去考虑。
(没有大于等于k了的话,显然子串的长度就不符合要求)
如果我们在【这一组】里遇到了一个a串的后缀。
那么显然,我们会在这个组里再找若干个b串的后缀。
然后分别求出这个a串和这些b串的最长公共前缀temp,然后累加temp-k+1就好.
但是这样实现的话,时间复杂度为O(n^2)
有一个性质,就是任意两个后缀的最长公共前缀为它们之间(排名之间)的height的最小值。
根据这个,用单调队列来做一个优化。
对于连续的大于等于k的height值。
i顺序扫过去。
遇到的如果是a串,那么直接cnt累加height[i]-k+1,然后把它放到单调队列里面
去,之后,如果再遇到a串的话,就看看是不是这个a串的height比之前的a串的height
值更小,如果更小的话,就用新的,更小的height去取代旧的且比较大的height;
(因为遇到了更小的,就说明,后面的b串如果要和a串递增答案,用的temp值就
更小了,则cnt也会相应改变)
这个过程就对应了,把单调队列的队尾的大的height给去掉。
然后更新cnt;(cnt是下次遇到b串的时候累加答案的)
然后把这个height[i]直接加到队列尾巴.
因为遇到了一个新的a的后缀,所以后面的b如果要和这个a后缀做答案,
肯定要从这个height[i]开始算最小值,后面如果有最小值,则后面再用
上面的方法改就好。
(之前虽然有更小的height,但是他们在[l,r]这个区间的左边。
因此把它加到队列是没有问题的.
如果遇到的是b串的后缀,那么,也应该用这个b串的后缀,来尝试改变一下单调
队列里面的最小值,如果可改的话,那么他就取代了某些a串的height值了。还
是一样,弹出队列尾巴比它大的。把这个新的加进去(因为是b的后缀,所以这个
加入到队列尾巴以后,除非他取代了之前的a串的height,否则不计算个数),然
后更改cnt的值。因为是b串,所以把之前的a串的cnt值累加进最后的答案就好。
(所以这个cnt值,就相当于是预先算出来了前面的所有a的后缀,和下一个会遇
到的b的后缀的答案的和)
然后要反过来做一遍。
即把b放在前面,遇到a累加b的答案。
比如
AABBAA
(A和B代表了A和B的一个后缀)
如果只做第一种
3->5
4->5
3->6
4->6
这4种就会漏掉
【错的次数】
【反思】
给你n个数字ai,
每个数字属于A或属于B;
对于所有属于B的数字a[i]
求出j<i;
且a[j]属于A,
temp = ∑(min(a[j..i])-k+1)
对于所有的a[i]的temp值求和
【代码】
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std; int k; const int N = 2e5;
const int MAX_CHAR = 255;//每个数字的最大值。
char s[N + 10];//如果是数字,就写成int s[N+10]就好,从0开始存
int Sa[N + 10], T1[N + 10], T2[N + 10], C[N + 10];
int Height[N + 10], Rank[N + 10];
int cnt[N + 10], dl[N + 10], tail; void build_Sa(int n, int m) {
int i, *x = T1, *y = T2;
for (i = 0; i<m; i++) C[i] = 0;
for (i = 0; i<n; i++) C[x[i] = s[i]]++;
for (i = 1; i<m; i++) C[i] += C[i - 1];
for (i = n - 1; i >= 0; i--) Sa[--C[x[i]]] = i;
for (int k = 1; k <= n; k <<= 1)
{
int p = 0;
for (i = n - k; i<n; i++) y[p++] = i;
for (i = 0; i<n; i++) if (Sa[i] >= k) y[p++] = Sa[i] - k;
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];
swap(x, y);
p = 1; x[Sa[0]] = 0;
for (i = 1; i<n; i++)
x[Sa[i]] = y[Sa[i - 1]] == y[Sa[i]] && y[Sa[i - 1] + k] == y[Sa[i] + k] ? p - 1 : p++;
if (p >= n) break;
m = p;
}
} void getHeight(int n)
{
int i, j, k = 0;
for (i = 1; i <= n; i++) Rank[Sa[i]] = i;
for (i = 0; i<n; i++) {
if (k) k--;
j = Sa[Rank[i] - 1];
while (s[i + k] == s[j + k]) k++;
Height[Rank[i]] = k;
}
} int main()
{
//freopen("F:\\rush.txt", "r", stdin);
while (~scanf("%d", &k) && k)
{
scanf("%s", s);
int l1 = strlen(s);
s[l1] = 'z' + 1;//分隔符
scanf("%s", s + l1 + 1);
int n = strlen(s);
s[n] = 0;//一个字符就链接完成了 //然后开始求后缀数组
build_Sa(n + 1, MAX_CHAR);//调用n+1
getHeight(n); //后缀数组求完了
//开始做算法 tail = 0;
long long temp = 0,ans =0;
for (int i = 2; i <= n-1; i++)//求连续的大于等于k的组,a->b
if (Height[i]>=k){
int num = 0;
while (tail >= 1 && Height[i] < dl[tail])
{//比Height[i]小的都去掉,因为区间最值肯定变小了
num+=cnt[tail];//加上a的后缀个数
temp -= 1LL*cnt[tail] * (dl[tail] - k + 1);
temp += 1LL*cnt[tail] * (Height[i] - k + 1);//改变累加值
tail--;
}
tail++;//把height[i]放到队列尾巴
dl[tail] = Height[i];
if (Sa[i-1] < l1)//是一个a串
{
temp += (Height[i] - k + 1);//累加它对后面的b的影响
cnt[tail] = num + 1;//把它看成一个整体了,前面的被他
//"吃掉的"都和它一样了。
}
else
{
//是一个b串
cnt[tail] = num;//他本身没有影响
}
if (Sa[i] > l1)//height[i]的影响范围只限于...i
//height[i+1]的影响范围是..i+1
//所以如果i是b串的话,它的height值已经算进去了
//可以想一下i=2时的运行情况。
{
ans += temp;
//i和j的lcp
//为height[i+1]..height[j]里面的最小值
}
}
else
{
tail = 0, temp = 0;
} temp = 0, tail = 0;
for (int i = 2; i <= n; i++)//求连续的大于等于k的组,b->a
if (Height[i] >= k) {
int num = 0;
while (tail >= 1 && Height[i] < dl[tail])
{//比Height[i]小的都去掉,因为区间最值肯定变小了
num += cnt[tail];//加上b的后缀个数
temp -= 1LL*cnt[tail] * (dl[tail] - k + 1);
temp += 1LL*cnt[tail] * (Height[i] - k + 1);//改变累加值
tail--;
}
tail++;//把height[i]放到队列尾巴
dl[tail] = Height[i];
if (Sa[i - 1] > l1)//是一个b串
{
temp += (Height[i] - k + 1);//累加它对后面的a的影响
cnt[tail] = num + 1;//把它看成一个整体了,前面的被他
//"吃掉的"都和它一样了。
}
else
{
//是一个a串
cnt[tail] = num;//他本身没有影响
}
if (Sa[i] < l1)//height[i]的影响范围只限于...i
//height[i+1]的影响范围是..i+1
//所以如果i是a串的话,它的height值已经算进去了
//可以想一下i=2时的运行情况。
{
ans += temp;
//i和j的lcp
//为height[i+1]..height[j]里面的最小值
}
}
else
{
tail = 0, temp = 0;
}
printf("%lld\n", ans);
}
return 0;
}
【POJ 3415】Common Substrings的更多相关文章
- 【POJ 3415】Common Substrings 长度不小于k的公共子串的个数
长度不小于k的公共子串的个数,论文里有题解,卡了一上午,因为sum没开long long!!! 没开long long毁一生again--- 以后应该早看POJ里的Discuss啊QAQ #inclu ...
- 【POJ - 1458】Common Subsequence(动态规划)
Common Subsequence Descriptions: A subsequence of a given sequence is the given sequence with some e ...
- 【POJ 1201】 Intervals(差分约束系统)
[POJ 1201] Intervals(差分约束系统) 11 1716的升级版 把原本固定的边权改为不固定. Intervals Time Limit: 2000MS Memory Limit: ...
- bzoj 2295: 【POJ Challenge】我爱你啊
2295: [POJ Challenge]我爱你啊 Time Limit: 1 Sec Memory Limit: 128 MB Description ftiasch是个十分受女生欢迎的同学,所以 ...
- 【链表】BZOJ 2288: 【POJ Challenge】生日礼物
2288: [POJ Challenge]生日礼物 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 382 Solved: 111[Submit][S ...
- BZOJ2288: 【POJ Challenge】生日礼物
2288: [POJ Challenge]生日礼物 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 284 Solved: 82[Submit][St ...
- BZOJ2293: 【POJ Challenge】吉他英雄
2293: [POJ Challenge]吉他英雄 Time Limit: 1 Sec Memory Limit: 128 MBSubmit: 80 Solved: 59[Submit][Stat ...
- BZOJ2287: 【POJ Challenge】消失之物
2287: [POJ Challenge]消失之物 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 254 Solved: 140[Submit][S ...
- BZOJ2295: 【POJ Challenge】我爱你啊
2295: [POJ Challenge]我爱你啊 Time Limit: 1 Sec Memory Limit: 128 MBSubmit: 126 Solved: 90[Submit][Sta ...
随机推荐
- Android引入library失败的可能原因
eclipse环境,引入library一直失败 谷歌到的原因是 Windows下 目标工程必须和引用工程在同一磁盘盘符下 然后我将要引入的library项目移动和现在项目同一磁盘下,引入成功 至于为什 ...
- vue methods 方法中 方法 调用 另一个方法。
vue在同一个组件内: methods中的一个方法调用methods中的另外一个方法. 可以在调用的时候 this.$options.methods.test(); this.$options.met ...
- CodeForcesGym 100502H Clock Pictures
Clock Pictures Time Limit: 1000ms Memory Limit: 524288KB This problem will be judged on CodeForcesGy ...
- tomcat日志配置之一自带log
问题 tomcat每次启动时,自动在logs目录下生产以下日志文件,且每天都会生成对应日期的一个文件,造成日志文件众多: localhost.2012-07-05.txt catalina.2012- ...
- IOS - 获取UITextField的输入文本
当UITextField文本改变时, 依据内容更新数据, 通过写监听事件就可以. 加入监听: [timesField addTarget:self action:@selector(textField ...
- Android 上的 制表符(tab) —— 一个奇妙的字符 (二)
接到上回的说,主要是上回那个问题,我认为是android的bug,黎叔认为是cocos2dx的bug,叫我去提交bug.所以我又继续研究了下. 上回说到会调用java层的函数去创建一个image,然后 ...
- 终结者:借助pinyin4j相关jar包提取汉字的首字母
import net.sourceforge.pinyin4j.PinyinHelper; import net.sourceforge.pinyin4j.format.HanyuPinyinCase ...
- 对 hiren bootcd 15.2 中的 XP 系统作了汉化, 同时支持中文输入法。提供下载
对 hiren bootcd 15.2 中的 XP 系统作了汉化, 同时支持中文输入法.提供下载 对该PE 中的 XP 系统作了汉化, 由于一个 中文字库 就要 10M 多:加之原系统过于精简,对中文 ...
- Vue router的query对象里的值的问题
在使用 $router.push() 时,如果使用了query,则可以在跳转后从query中获取到对应的参数.如果传的是字符串自然没问题,但是如果传的其他类型的数据,在跳转之后是正常的,而跳转之后再刷 ...
- MFC窗口显隐
使用SetLayeredWindowAttributes可以方便的制作透明窗体,此函数在w2k以上才支持,而且如果希望直接使用的话,可能需要下载最新的SDK.不过此函数在w2k的user32.dll里 ...