bzoj2085 [Poi2010]Hamsters 哈希+倍增

Description

---

Tz养了一群仓鼠，他们都有英文小写的名字，现在Tz想用一个字母序列来表示他们的名字，只要他们的名字是字母序列中的一个子串就算，出现多次可以重复计算。现在Tz想好了要出现多少个名字，请你求出最短的字母序列的长度是多少。

1<=n<=200
1<=m<=10^9

Solution

---

一开始完全没看到m的范围于是傻了，被文化课吊打后还要在这里翻船

可以用f[i,j,k]表示以i开始j结束出现了2^k个名字的最短序列，推的时候倍增就行了
初始条件f[i,j,0]可以用哈希+二分，实际上暴力扫是不会T的
为了保证空间f可以滚动

Code

---

#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <math.h>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define drp(i,st,ed) for (int i=st;i>=ed;--i)
#define fill(x,t) memset(x,t,sizeof(x))
#define copy(x,t) memcpy(x,t,sizeof(x))

typedef long long LL;
const LL INF=100000000000000000LL;
const int MOD1=1000000007;
const int MOD2=1000000009;
const int N=205;
const int L=100005;

LL pow1[L],pow2[L],tmp[N];
LL f[N][N],g[N][N],ans[N];

int len[N],hash1[N][L],hash2[N][L];

char str[L];

LL get_dis(int x,int y) {
    int lim=std:: min(len[x],len[y])-(x==y);
    drp(i,lim,1) {
        LL u1=(hash1[x][len[x]]-(LL)hash1[x][len[x]-i]*pow1[i]%MOD1+MOD1)%MOD1;
        LL u2=hash1[y][i];
        if (u1!=u2) continue;
        LL v1=(hash2[x][len[x]]-(LL)hash2[x][len[x]-i]*pow2[i]%MOD2+MOD2)%MOD2;
        LL v2=hash2[y][i];
        if (v1==v2) return len[y]-i;
    }
    return len[y];
}

int main(void) {
    int n,m,max=0; scanf("%d%d",&n,&m); m-=1;
    rep(i,1,n) {
        scanf("%s",str+1);
        len[i]=strlen(str+1);
        max=std:: max(len[i],max);
        rep(j,1,len[i]) {
            hash1[i][j]=(hash1[i][j-1]*27LL%MOD1+str[j]-'a'+1)%MOD1;
            hash2[i][j]=(hash2[i][j-1]*27LL%MOD2+str[j]-'a'+1)%MOD2;
        }
    }
    pow1[0]=1; rep(i,1,max) pow1[i]=pow1[i-1]*27LL%MOD1;
    pow2[0]=1; rep(i,1,max) pow2[i]=pow2[i-1]*27LL%MOD2;
    rep(i,1,n) ans[i]=len[i];
    rep(i,1,n) rep(j,1,n) {
        f[i][j]=get_dis(i,j);
    }
    int T=log2(m)+0.1;
    rep(pp,0,T) {
        LL now=1<<pp;
        if (now&m) {
            rep(i,1,n) tmp[i]=INF;
            rep(i,1,n) {
                rep(j,1,n) if (f[i][j]+ans[i]<tmp[j]) {
                    tmp[j]=f[i][j]+ans[i];
                }
            }
            rep(i,1,n) ans[i]=tmp[i];
        }
        rep(i,1,n) rep(j,1,n) {
            g[i][j]=INF;
            rep(k,1,n) {
                if (f[i][k]+f[k][j]<g[i][j]) {
                    g[i][j]=f[i][k]+f[k][j];
                }
            }
        }
        rep(i,1,n) rep(j,1,n) f[i][j]=g[i][j];
    }
    LL prt=INF;
    rep(i,1,n) prt=std:: min(prt,ans[i]);
    printf("%lld\n", prt);
    return 0;
}