poj3167&hdu4749 KMP+BIT

        一道比较好玩的KMP,他在匹配的时候不是按该位字符是否想到,而是按当前字符的名次来判断是否匹配,转化一下就是对i,j,若i之前中比s[i]小的和j之前中比s[j]小的数的个数相等并且i之前小于等于s[i]的数和j之前小于等于s[j]的个数相等,那么就判定s[i]和s[j]匹配。以此方法构造状态转移图后,拿母串跑一遍即可。维护的过程可以用树状数组实现。hdu4749那个和poj3167基本一样,最后求出所有答案后加一个贪心就过了....

            贴一个poj3167的代码

           

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <string>
#include <cstring>
#include <queue>
#include <stack>
#include <map>
using namespace std;
typedef long long ll;
const int maxn=200000+1000;
int a[maxn],b[maxn];
int b1[maxn],b2[maxn];
int f[maxn];
bool ok[maxn];
int in[maxn];
int n,m,t,tt,k,s;
struct BIT
{
    int tree[50];
    void init()
    {
        memset(tree,0,sizeof tree);
    }
    int lowbit(int x)
    {
        return x&(-x);
    }
    void modify(int x,int c=1)
    {
        for (int i=x; i<=s; i+=lowbit(i))
        tree[i]+=c;
    }
    int query(int x)
    {
        int res=0;
        for (int i=x; i>0; i-=lowbit(i))
        res+=tree[i];
        return res;
    }
}bit;
void getFail()
{
    memset(in,0,sizeof in);
    memset(bit.tree,0,sizeof bit.tree);
    memset(f,0,sizeof f);
    f[0]=f[1]=0;
    for (int i=1; i<k; i++)
    {
        int j=f[i];
        bit.modify(b[i]); in[i]++;

        while (j && (bit.query(b[i]-1)!=b1[j] || bit.query(b[i])!=b2[j]))
        {
            for (int kk=i-j; kk<i-f[j]; kk++)
            bit.modify(b[kk],-1),in[kk]--;
            j=f[j];
        }
        if (bit.query(b[i]-1)==b1[j] && bit.query(b[i])==b2[j]) f[i+1]=++j;
        else f[i+1]=0;
    }
}
int kmp()
{
    memset(in,0,sizeof in);
    memset(bit.tree,0,sizeof bit.tree);
    int j=0;
    int res=0;
    for (int i=0; i<n; i++)
    {
        bit.modify(a[i],1); in[i]++;
        while(j && (bit.query(a[i]-1)!=b1[j] || bit.query(a[i])!=b2[j]))
        {
            for (int kk=i-j; kk<i-f[j]; kk++)
            bit.modify(a[kk],-1),in[kk]--;
            j=f[j];
        }
        if (bit.query(a[i]-1)==b1[j] && bit.query(a[i])==b2[j]) ++j;
        if (j==k) res++,ok[i]=true;
    }
    return res;
}
int main()
{
    memset(b1,0,sizeof b1);
    memset(b2,0,sizeof b2);
    memset(ok,false,sizeof ok);
    scanf("%d%d%d",&n,&k,&s);
    bit.init();
    for (int i=0; i<n; i++)
    scanf("%d",&a[i]);
    for (int i=0; i<k; i++)
    {
        scanf("%d",&b[i]);
        bit.modify(b[i]);
        b1[i]=bit.query(b[i]-1);
        b2[i]=bit.query(b[i]);
    }
    getFail();
//    for (int i=0; i<=k; i++)
//    printf("%d ",f[i]);
//    puts("");
    printf("%d\n",kmp());
    for (int i=0; i<n; i++)
    {
        if (ok[i]) printf("%d\n",i-k+2);
    }


    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值