bzoj 2085: [Poi2010]Hamsters hash+倍增floyd

题意

给定n个长度总和不超过10W的字符串，求一个最短的母串，使所有字符串的出现次数之和=m 这n个字符串保证不互相包含

分析

一开始题目没有给出字符串总长，以为就是AC自动机上dp+矩阵乘法随便优化一下就好了。。。

既然不能对长度和字符集dp，那我们就可以考虑对字符串dp。
设f[i,j]表示j接到i后面最少增加多少个字符。
看到m辣么大，很容易想到倍增，那我们就可以用f[i,j,k]表示i后面接长度为2^k个字符串且最有一个为j时的最小长度。转移的话就像floyd那样即可。
一开始在预处理f[i,j,0]的时候，我用的是hash+二分，后来对拍发现原来这个东西是不满足二分性的，随手画一下就知道了。
其实f[i,j,0]是可以直接hash+枚举的，算了一下总复杂度是O(nL)，不会T。
注意倍增数组最好滚动，不然会T。

代码

#include<iostream>#include<cstdio>#include<cstdlib>#include<cstring>#include<algorithm>#include<cmath>using namespace std;typedef long long LL;const int N=205;const int L=100005;const int MOD1=233333333;const int MOD2=1000000007;const LL inf=1e17;int n,m,hash1[N][L],hash2[N][L],mi1[L],mi2[L],len[N];LL f[N][N],ans[N],tmp[N],g[N][N];char ch[L];int get_dis(int x,int y){    int lim=min(len[x],len[y]);    if (x==y) lim--;    for (int l=lim;l;l--)    {        int u1=(hash1[x][len[x]]-(LL)hash1[x][len[x]-l]*mi1[l]%MOD1+MOD1)%MOD1;        int v1=(hash2[x][len[x]]-(LL)hash2[x][len[x]-l]*mi2[l]%MOD2+MOD2)%MOD2;        int u2=hash1[y][l],v2=hash2[y][l];        if (u1==u2&&v1==v2) return len[y]-l;    }    return len[y];}int main(){    scanf("%d%d",&n,&m);    int mxl=0;    for (int i=1;i<=n;i++)    {        scanf("%s",ch+1);        len[i]=strlen(ch+1);        mxl=max(mxl,len[i]);        for (int j=1;j<=len[i];j++)        {            hash1[i][j]=((LL)hash1[i][j-1]*27%MOD1+ch[j]-'a'+1)%MOD1;            hash2[i][j]=((LL)hash2[i][j-1]*27%MOD2+ch[j]-'a'+1)%MOD2;        }    }    mi1[0]=mi2[0]=1;    for (int i=1;i<=mxl;i++) mi1[i]=(LL)mi1[i-1]*27%MOD1,mi2[i]=(LL)mi2[i-1]*27%MOD2;    for (int i=1;i<=n;i++)        for (int j=1;j<=n;j++)            f[i][j]=get_dis(i,j);    int T=log(m)/log(2)+0.1;    for (int i=1;i<=n;i++) ans[i]=len[i];    m--;    for (int l=0;l<=T;l++)    {        if ((m&(1<<l)))        {            for (int j=1;j<=n;j++) tmp[j]=inf;            for (int j=1;j<=n;j++)                for (int k=1;k<=n;k++)                    tmp[k]=min(tmp[k],ans[j]+f[j][k]);            for (int j=1;j<=n;j++) ans[j]=tmp[j];        }        for (int i=1;i<=n;i++)            for (int j=1;j<=n;j++)            {                g[i][j]=inf;                for (int k=1;k<=n;k++)                    g[i][j]=min(g[i][j],f[i][k]+f[k][j]);            }        for (int i=1;i<=n;i++)            for (int j=1;j<=n;j++)                f[i][j]=g[i][j];    }    LL mn=inf;    for (int i=1;i<=n;i++) mn=min(mn,ans[i]);    printf("%lld",mn);    return 0;}