scut 125. 笔芯回文

https://scut.online/p/125

看数据量，这题可能是O（n^2）的dp

也可能是区间dp，但是区间dp一般复杂度是O(n^3)，虽然也可以优化，但是比赛的时候那么多人“秒”了，应该不会是那么麻烦的。

套路：设dp[i]表示前i个字符中能拿到的最大贡献。dp[len]就是答案。

如果[L, R]这段区间能组成回文串，那么就有两种决策，删除或则不删除。

删除：dp[R] = dp[L - 1] + a[R - L + 1]

不：  dp[R] = dp[R];

取个max就行。

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <assert.h>
#define IOS ios::sync_with_stdio(false)
using namespace std;
#define inf (0x3f3f3f3f)
typedef long long int LL;

#include <iostream>
#include <sstream>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <string>
#include <bitset>
const int maxn = 1e5 + ;
int a[maxn];
LL dp[maxn];
char str[maxn], sub[maxn];
int p[maxn];
int manacher(char str[], int lenstr) {
    str[] = '*'; //表示一个不可能的值
    //目标要插入lenstr+1个'#'，所以长度变成2*lenstr+1
    for (int i = lenstr; i >= ; i--) { //str[lenstr+1]是'\0'
        //i=lenstr时，i+i+2那个值要赋为'\0';
        //总长度只能是lenstr+lenstr+2,所以i从lenstr开始枚举
        str[i + i + ] = str[i + ];
        str[i + i + ] = '#';
    }
    int id = , maxlen = ; //现在开始在str[2]了
    for (int i = ; i <=  * lenstr + ; i++) { //2*lenstr+1是'#'没用
        if (p[id] + id > i) { //没取等号的，只能去到p[id]+id-1
            //p[id]+id是越界的，减去i即为区间长度
            //p[id]+id-i,这个是所有可能中的最大值了
            p[i] = min(p[id] + id - i, p[ * id - i]);
        } else p[i] = ; //记得修改的是p[i]
        while (str[i + p[i]] == str[i - p[i]]) ++p[i];
        if (p[id] + id < p[i] + i) id = i;
        maxlen = max(maxlen, p[i]);
    }
    return maxlen - ;
}
bool isok(int be, int en) {
    int len = en - be + ;
    be <<= ;
    en <<= ;
    int mid = (be + en) >> ;
    return p[mid] -  >= len;
}
void work() {
    int n;
    cin >> n;
    for (int i = ; i <= n; ++i) cin >> a[i];
    memset(dp, , sizeof dp);
    cin >> str + ;
    int lenstr = strlen(str + );
    strcpy(sub + , str + );
    manacher(str, lenstr);
    for (int i = ; i <= lenstr; ++i) {
        for (int j = ; j <= i; ++j) {
            if (i - j +  > n) continue;
            if (isok(j, i)) {
                dp[i] = max(dp[i], dp[j - ] + a[i - j + ]);
            }
        }
    }
    cout << dp[lenstr] << endl;
}

int main() {
#ifdef local
    freopen("data.txt", "r", stdin);
//    freopen("data.txt", "w", stdout);
#endif
    int t;
    cin >> t;
    while (t--) work();
    return ;
}