294. 計算重複（算法競賽進階指南，倍增優化 DP）

一.題目鏈接：

計算重複

二.題目大意：

題目簡單易懂

三.分析：

題目要求最大的 m，使得 conn(s2, n2 * m) 是 conn(s1, n1) 的子序列.

等價於求最大的 h，使得 conn(s2, h) 是 conn(s1, n1) 的子序列

易知，m = h / n2.

第一想法是二分答案，可是 check 不好寫，因此可以通過 預處理 + 倍增 求解.

由於在倍增枚舉 k 時，如果能湊出來，那麼需要知道應從當前點 s1[p] 跳多遠

因此可以預處理出 f[i][j]: 從 s1[i] 開始，最少走多少，能包含 2^j 個 s2.

說的有點亂，詳見代碼.

四.代碼實現：

/**
狀態表示 f[i][j]: 從 s1[i] 開始，最少要走多少步，才能包含 2^j 個 s2
狀態計算 f[i][j] = f[i][j - 1] + f[i + f[i][j - 1]][j - 1], j > 0
**/

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;

const int M = (int)1e2;

int n1, n2;
char s1[M + 5], s2[M + 5];
ll f[M + 5][27];

ll work()
{
    int len1 = strlen(s1);
    int len2 = strlen(s2);
    for(int i = 0; i < len1; ++i)
    {
        f[i][0] = 0;
        int p = i, c = 0;
        for(int j = 0; j < len2; ++j)
        {
            while(s1[p] != s2[j])
            {
                p = (p + 1) % len1;
                ++c;
                if(c == len1)
                {
                    return 0;
                }
            }
            f[i][0] += c + 1;
            c = 0;
            p = (p + 1) % len1;
        }
    }
    for(int j = 1; j <= 26; ++j)
    {
        for(int i = 0; i < len1; ++i)
        {
            f[i][j] = f[i][j - 1] + f[(i + f[i][j - 1]) % len1][j - 1];
        }
    }
    ll p = 0, m = 0;
    for(int k = 26; k >= 0; --k)
    {
        if(p + f[p % len1][k] <= n1 * len1)
        {
            m |= (1<<k);
            p += f[p % len1][k];
        }
    }
    return m / n2;
}

int main()
{
//    freopen("input.txt", "r", stdin);
    while(~scanf("%s %d %s %d", s2, &n2, s1, &n1))
    {
        printf("%lld\n", work());
    }
    return 0;
}