[poi2010]Hamsters

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

思路：首先可以处理出g[a][b]表示b接在a后面的长度(即重复部分之前都不算)

那么就可以转化成求长度m-1的最短路。

假设dis[i][a][b]为经过i步a走到b的最短路，那么dis[i][a][b] = min(dis[i/2][a][c] + dis[i-i/2][c][b]) (1<=c<=n)

那么是不是就可以利用快速幂的思想了。。

前面求g[a][b]可利用hash判断是否一样

code：

 1 /*
 2  * Author:  Yzcstc
 3  * Created Time:  2014/11/12 20:54:48
 4  * File Name: hamsters.cpp
 5  */
 6 #include<cstdio>
 7 #include<iostream>
 8 #include<cstring>
 9 #include<cstdlib>
10 #include<cmath>
11 #include<algorithm>
12 #include<string>
13 #include<map>
14 #include<set>
15 #include<vector>
16 #include<queue>
17 #include<stack>
18 #include<ctime>
19 #define M0(x)  memset(x, 0, sizeof(x))
20 #define repf(i, a, b) for (int i = (a); i <= (b); ++i)
21 #define Inf  0x7fffffff
22 #define M 10007
23 using namespace std;
24 typedef long long ll;
25 const int maxn = 204;
26 int n, m;
27 char str[120000];
28 vector<int> hs[210];
29 int sz[maxn], pw[120000];
30 ll A[maxn][maxn], B[maxn][maxn];
31 ll s[maxn], tmp[maxn];
32
33 inline int Hash(const int &p, const int& l, const int& r){
34        return hs[p][r] - hs[p][l-1] * pw[r-l+1];
35 }
36
37 void init(){
38      int len, v;
39      repf(i, 1, n){
40            scanf("%s", str + 1);
41            sz[i] = len = strlen(str + 1);
42            hs[i].push_back(v = 0);
43            repf(j, 1, len)
44                 v = v * M + str[j], hs[i].push_back(v);
45      }
46 }
47
48 void modify(ll *s, ll A[][maxn]){
49      memset(tmp, 0x3f, sizeof(tmp));
50 //     cout << tmp[0] << endl;
51      for (int i = 1; i <= n; ++i)
52             for (int j = 1; j <= n; ++j)
53                    tmp[j] = min(tmp[j], s[i] + A[i][j]);
54      repf(i, 1, n) s[i] = tmp[i];
55 }
56
57 void modify(ll A[][maxn]){
58      memset(B, 0x3f, sizeof(B));
59      repf(k, 1, n) repf(i, 1, n) repf(j, 1, n)
60            B[i][j] = min(B[i][j], A[i][k] + A[k][j]);
61      repf(i, 1, n) repf(j, 1, n) A[i][j] = B[i][j];
62 }
63
64 void quick(int n){
65      for (;n>0; n>>=1){
66            if (n & 1) modify(s, A);
67            modify(A);
68      }
69 }
70
71 void solve(){
72      repf(i, 1, n) repf(j, 1, n){
73            A[i][j] = sz[j];
74            for (int k = max(2, sz[i]-sz[j]+2); k <= sz[i]; ++k)
75                 if (Hash(i, k, sz[i]) == Hash(j, 1, sz[i]-k+1)){
76                       A[i][j] = sz[j] - (sz[i]-k+1);  break;
77                 }
78      }
79      repf(i, 1, n) s[i] = sz[i];
80      quick(--m);
81      ll ans = 1LL<<62;
82 //     cout << ans << endl;
83      repf(i, 1, n) ans = min(ans, s[i]);
84      cout << ans << endl;
85 }
86
87 int main(){
88     pw[0] = 1;
89     repf(i, 1, 100000) pw[i] = pw[i-1] * M;
90     while (scanf("%d%d", &n, &m) != EOF){
91         init();
92         solve();
93     }
94     return 0;
95 }