关闭

[POJ 1743] Musical Theme (字符串哈希+二分)

标签: acmpoj字符串哈希二分
74人阅读 评论(0) 收藏 举报

链接

POJ 1743


题意

给出一个序列,元素范围1到88,求出最长连续子序列长度,满足该序列(或该序列的逆序)在序列中出现过至少2次并且不相交(没有公共部分)。


思路

又是子串+长度问题,可用哈希+二分解决。
把序列看做字符串,字符集是1到88,对字符串正序逆序求哈希,二分长度len,从左向右枚举长度为len的子串。这时候需要对正序的子串和相同位置的逆序子串同时枚举。枚举后判断是否出现过且不相交,没出现过的话就存储该子串的哈希值到哈希表内,出现过则不重复存储,只记录最左侧的位置即可。


代码

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
typedef unsigned long long ulint;
const ulint seed = 30007uLL;
//const ulint mod = 1e9 + 1107uLL;
#define maxn 20020
ulint H[maxn], xp[maxn];
int s[maxn], N;

void initHash()
{
    H[0] = s[0];
    for(int i = 1; i < N; i++)
        H[i] = H[i - 1]*seed + s[i];
}

ulint askHash(int l, int r)
{
    if(l == 0) return H[r];
    return H[r] - H[l - 1]*xp[r - l + 1];
}


#define mod 100003
ulint h[mod];
int bg[mod], nx[mod], pos[mod];

bool check(int len)
{
    memset(bg, 0, sizeof(bg));

    ulint ht;
    int e = 0;

    for(int i = 0, l, r; i + len - 1 < N; i++)
    {
        l = i, r = i + len - 1;
        ht = askHash(l, r);

        for(int p = bg[ht % mod]; p; p = nx[p])
            if(h[p] == ht && i - pos[p] >= len)
                return true;

        h[++e] = ht;
        nx[e] = bg[ht % mod];
        bg[ht % mod] = e;
        pos[e] = i;
    }

    return false;
}

int main()
{
    //freopen("1743.txt", "r", stdin);

    xp[0] = 1;
    for(int i = 1; i < maxn; i++)
    {
        xp[i] = xp[i-1] * seed;
    }

    while(scanf("%d", &N) && N)
    {
        for(int i = 0; i < N; i++)
        {
            scanf("%d", &s[i]);
        }
        N--;
        for(int i = 0; i < N; i++)
        {
            s[i] = 89 + s[i+1] - s[i];
        }

        initHash();

        int ans = 0, l = 1, r = N / 2, m;
        while(l <= r)
        {
            m = (l + r) >> 1;

            if(check(m))
            {
                l = m + 1;
                ans = m;
            }
            else r = m - 1;
        }
        cout << ((ans >= 4)? ans+1: 0) << endl;
    }
    return 0;
}
0
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:7389次
    • 积分:667
    • 等级:
    • 排名:千里之外
    • 原创:63篇
    • 转载:0篇
    • 译文:0篇
    • 评论:0条